mp_sqrmod
mp_err mp_sqrmod(const mp_int *a, const mp_int *m, mp_int *c)
#define mp_sqrmod(a, m, c) mp_mulmod(a, a, m, c)