mp_sqrmod
return mp_sqrmod(a, &meth->irr, r);
mp_err mp_sqrmod(const mp_int *a, const mp_int *m, mp_int *c);