mp_div(0, &sr, sr, g->r);
mpbarrett_create(&b, g->r);
ss = mp_mul(ss, sr, c->u); ss = mpbarrett_reduce(&b, ss, ss);
ss = mp_add(ss, ss, mr); mp_div(0, &ss, ss, g->r);
mp_div(0, &sr, sr, g->r);
mpbarrett_create(&b, g->r);
ss = mp_mul(ss, sr, c->u); ss = mpbarrett_reduce(&b, ss, ss);
ss = mp_add(ss, ss, mr); mp_div(0, &ss, ss, g->r);