1062 powm (uintmax_t b, uintmax_t e, uintmax_t n, uintmax_t ni, uintmax_t one)  in powm()  argument
1071       b = mulredc (b, b, n, ni);  in powm()
1075         y = mulredc (y, b, n, ni);  in powm()
1084        uintmax_t ni, const uintmax_t *one)  in powm2()  argument
1102           r0 = mulredc2 (r1m, r1, r0, b1, b0, n1, n0, ni);  in powm2()
1105       b0 = mulredc2 (r1m, b1, b0, b1, b0, n1, n0, ni);  in powm2()
1112           r0 = mulredc2 (r1m, r1, r0, b1, b0, n1, n0, ni);  in powm2()
1115       b0 = mulredc2 (r1m, b1, b0, b1, b0, n1, n0, ni);  in powm2()
1124 millerrabin (uintmax_t n, uintmax_t ni, uintmax_t b, uintmax_t q,  in millerrabin()  argument
1127   uintmax_t y = powm (b, q, n, ni, one);  in millerrabin()
1136       y = mulredc (y, y, n, ni);  in millerrabin()
1147 millerrabin2 (const uintmax_t *np, uintmax_t ni, const uintmax_t *bp,  in millerrabin2()  argument
1152   y0 = powm2 (&r1m, bp, qp, np, ni, one);  in millerrabin2()
1165       y0 = mulredc2 (&r1m, y1, y0, y1, y0, np[1], np[0], ni);  in millerrabin2()
1203   uintmax_t a_prim, one, ni;  in prime_p()  local
1219   binv (ni, n);                 /* ni <- 1/n mod B */  in prime_p()
1224   if (!millerrabin (n, ni, a_prim, q, k, one))  in prime_p()
1243                 = powm (a_prim, (n - 1) / factors.p[i], n, ni, one) != one;  in prime_p()
1272       if (!millerrabin (n, ni, a_prim, q, k, one))  in prime_p()
1286   uintmax_t ni;  in prime2_p()  local
1310   binv (ni, n0);  in prime2_p()
1318   if (!millerrabin2 (na, ni, a_prim, q, k, one))  in prime2_p()
1343               y[0] = powm2 (&y[1], a_prim, e, na, ni, one);  in prime2_p()
1355               y[0] = powm2 (&y[1], a_prim, e, na, ni, one);  in prime2_p()
1371       if (!millerrabin2 (na, ni, a_prim, q, k, one))  in prime2_p()
1464   uintmax_t x, z, y, P, t, ni, g;  in factor_using_pollard_rho()  local
1477       binv (ni, n);             /* FIXME: when could we use old 'ni' value?  */  in factor_using_pollard_rho()
1483               x = mulredc (x, x, n, ni);  in factor_using_pollard_rho()
1487               P = mulredc (P, t, n, ni);  in factor_using_pollard_rho()
1503               x = mulredc (x, x, n, ni);  in factor_using_pollard_rho()
1512           y = mulredc (y, y, n, ni);  in factor_using_pollard_rho()
1550   uintmax_t x1, x0, z1, z0, y1, y0, P1, P0, t1, t0, ni, g1, g0, r1m;  in factor_using_pollard_rho2()  local
1562       binv (ni, n0);  in factor_using_pollard_rho2()
1568               x0 = mulredc2 (&r1m, x1, x0, x1, x0, n1, n0, ni);  in factor_using_pollard_rho2()
1573               P0 = mulredc2 (&r1m, P1, P0, t1, t0, n1, n0, ni);  in factor_using_pollard_rho2()
1591               x0 = mulredc2 (&r1m, x1, x0, x1, x0, n1, n0, ni);  in factor_using_pollard_rho2()
1601           y0 = mulredc2 (&r1m, y1, y0, y1, y0, n1, n0, ni);  in factor_using_pollard_rho2()