| /netbsd-src/external/lgpl3/gmp/dist/mpn/generic/ |
| H A D | div_q.c | 107 gmp_pi1_t dinv; in mpn_div_q() local 154 invert_pi1 (dinv, new_dp[dn - 1], new_dp[dn - 2]); in mpn_div_q() 155 qh = mpn_sbpi1_div_q (qp, new_np, new_nn, new_dp, dn, dinv.inv32); in mpn_div_q() 162 invert_pi1 (dinv, new_dp[dn - 1], new_dp[dn - 2]); in mpn_div_q() 163 qh = mpn_dcpi1_div_q (qp, new_np, new_nn, new_dp, dn, &dinv); in mpn_div_q() 188 invert_pi1 (dinv, dh, dp[dn - 2]); in mpn_div_q() 189 qh = mpn_sbpi1_div_q (qp, new_np, nn, dp, dn, dinv.inv32); in mpn_div_q() 196 invert_pi1 (dinv, dh, dp[dn - 2]); in mpn_div_q() 197 qh = mpn_dcpi1_div_q (qp, new_np, nn, dp, dn, &dinv); in mpn_div_q() 241 invert_pi1 (dinv, new_dp[qn], new_dp[qn - 1]); in mpn_div_q() [all …]
|
| H A D | dcpi1_bdiv_q.c | 58 mp_limb_t dinv, mp_ptr tp) in mpn_dcpi1_bdiv_q_n() argument 68 cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, lo, dinv, tp); in mpn_dcpi1_bdiv_q_n() 82 mpn_sbpi1_bdiv_q (qp, np, n, dp, n, dinv); in mpn_dcpi1_bdiv_q_n() 95 mp_limb_t dinv) in mpn_dcpi1_bdiv_q() argument 121 cy = mpn_sbpi1_bdiv_qr (qp, np, 2 * qn, dp, qn, dinv); in mpn_dcpi1_bdiv_q() 123 cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, qn, dinv, tp); in mpn_dcpi1_bdiv_q() 144 cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, dn, dinv, tp); in mpn_dcpi1_bdiv_q() 149 mpn_dcpi1_bdiv_q_n (qp, np, dp, dn, dinv, tp); in mpn_dcpi1_bdiv_q() 154 mpn_sbpi1_bdiv_q (qp, np, qn, dp, qn, dinv); in mpn_dcpi1_bdiv_q() 156 mpn_dcpi1_bdiv_q_n (qp, np, dp, qn, dinv, tp); in mpn_dcpi1_bdiv_q()
|
| H A D | dcpi1_bdiv_qr.c | 62 mp_limb_t dinv, mp_ptr tp) in mpn_dcpi1_bdiv_qr_n() argument 72 cy = mpn_sbpi1_bdiv_qr (qp, np, 2 * lo, dp, lo, dinv); in mpn_dcpi1_bdiv_qr_n() 74 cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, lo, dinv, tp); in mpn_dcpi1_bdiv_qr_n() 82 cy = mpn_sbpi1_bdiv_qr (qp + lo, np + lo, 2 * hi, dp, hi, dinv); in mpn_dcpi1_bdiv_qr_n() 84 cy = mpn_dcpi1_bdiv_qr_n (qp + lo, np + lo, dp, hi, dinv, tp); in mpn_dcpi1_bdiv_qr_n() 96 mp_srcptr dp, mp_size_t dn, mp_limb_t dinv) in mpn_dcpi1_bdiv_qr() argument 122 cy = mpn_sbpi1_bdiv_qr (qp, np, 2 * qn, dp, qn, dinv); in mpn_dcpi1_bdiv_qr() 124 cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, qn, dinv, tp); in mpn_dcpi1_bdiv_qr() 146 cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, dn, dinv, tp); in mpn_dcpi1_bdiv_qr() 157 cy = mpn_sbpi1_bdiv_qr (qp, np, 2 * qn, dp, qn, dinv); in mpn_dcpi1_bdiv_qr() [all …]
|
| H A D | dcpi1_divappr_q.c | 44 gmp_pi1_t *dinv, mp_ptr tp) in mpn_dcpi1_divappr_q_n() argument 53 qh = mpn_sbpi1_div_qr (qp + lo, np + 2 * lo, 2 * hi, dp + lo, hi, dinv->inv32); in mpn_dcpi1_divappr_q_n() 55 qh = mpn_dcpi1_div_qr_n (qp + lo, np + 2 * lo, dp + lo, hi, dinv, tp); in mpn_dcpi1_divappr_q_n() 70 ql = mpn_sbpi1_divappr_q (qp, np + hi, 2 * lo, dp + hi, lo, dinv->inv32); in mpn_dcpi1_divappr_q_n() 72 ql = mpn_dcpi1_divappr_q_n (qp, np + hi, dp + hi, lo, dinv, tp); in mpn_dcpi1_divappr_q_n() 86 mp_srcptr dp, mp_size_t dn, gmp_pi1_t *dinv) in mpn_dcpi1_divappr_q() argument 145 udiv_qr_3by2 (q, n1, n0, n2, n1, n0, d1, d0, dinv->inv32); in mpn_dcpi1_divappr_q() 177 qh = mpn_sbpi1_div_qr (qp, np - qn, 2 * qn, dp - qn, qn, dinv->inv32); in mpn_dcpi1_divappr_q() 179 qh = mpn_dcpi1_div_qr_n (qp, np - qn, dp - qn, qn, dinv, tp); in mpn_dcpi1_divappr_q() 204 mpn_dcpi1_div_qr_n (qp, np - dn, dp - dn, dn, dinv, tp); in mpn_dcpi1_divappr_q() [all …]
|
| H A D | div_qr_1.c | 63 mp_limb_t dinv, q; in mpn_div_qr_1() local 82 invert_limb (dinv, d); in mpn_div_qr_1() 83 return mpn_div_qr_1n_pi1 (qp, up, n, uh, d, dinv); in mpn_div_qr_1() 88 mp_limb_t dinv, ul; in mpn_div_qr_1() local 120 invert_limb (dinv, d); in mpn_div_qr_1() 122 udiv_qrnnd_preinv (*qh, uh, uh, ul, d, dinv); in mpn_div_qr_1() 123 return mpn_div_qr_1n_pi1 (qp, qp, n, uh, d, dinv) >> cnt; in mpn_div_qr_1()
|
| H A D | divrem_1.c | 147 mp_limb_t dinv; in mpn_divrem_1() local 148 invert_limb (dinv, d); in mpn_divrem_1() 153 udiv_qrnnd_preinv (*qp, r, r, n0, d, dinv); in mpn_divrem_1() 159 udiv_qrnnd_preinv (*qp, r, r, CNST_LIMB(0), d, dinv); in mpn_divrem_1() 226 mp_limb_t dinv, nshift; in mpn_divrem_1() local 227 invert_limb (dinv, d); in mpn_divrem_1() 236 udiv_qrnnd_preinv (*qp, r, r, nshift, d, dinv); in mpn_divrem_1() 241 udiv_qrnnd_preinv (*qp, r, r, n1 << cnt, d, dinv); in mpn_divrem_1() 247 udiv_qrnnd_preinv (*qp, r, r, CNST_LIMB(0), d, dinv); in mpn_divrem_1()
|
| H A D | dcpi1_div_qr.c | 44 gmp_pi1_t *dinv, mp_ptr tp) in mpn_dcpi1_div_qr_n() argument 53 qh = mpn_sbpi1_div_qr (qp + lo, np + 2 * lo, 2 * hi, dp + lo, hi, dinv->inv32); in mpn_dcpi1_div_qr_n() 55 qh = mpn_dcpi1_div_qr_n (qp + lo, np + 2 * lo, dp + lo, hi, dinv, tp); in mpn_dcpi1_div_qr_n() 70 ql = mpn_sbpi1_div_qr (qp, np + hi, 2 * lo, dp + hi, lo, dinv->inv32); in mpn_dcpi1_div_qr_n() 72 ql = mpn_dcpi1_div_qr_n (qp, np + hi, dp + hi, lo, dinv, tp); in mpn_dcpi1_div_qr_n() 93 gmp_pi1_t *dinv) in mpn_dcpi1_div_qr() argument 151 udiv_qr_3by2 (q, n1, n0, n2, n1, n0, d1, d0, dinv->inv32); in mpn_dcpi1_div_qr() 184 qh = mpn_sbpi1_div_qr (qp, np - qn, 2 * qn, dp - qn, qn, dinv->inv32); in mpn_dcpi1_div_qr() 186 qh = mpn_dcpi1_div_qr_n (qp, np - qn, dp - qn, qn, dinv, tp); in mpn_dcpi1_div_qr() 212 mpn_dcpi1_div_qr_n (qp, np - dn, dp - dn, dn, dinv, tp); in mpn_dcpi1_div_qr() [all …]
|
| H A D | sbpi1_bdiv_q.c | 50 mp_limb_t dinv) in mpn_sbpi1_bdiv_q() argument 58 ASSERT (-(dp[0] * dinv) == 1); in mpn_sbpi1_bdiv_q() 66 q = dinv * up[0]; in mpn_sbpi1_bdiv_q() 78 q = dinv * up[0]; in mpn_sbpi1_bdiv_q() 87 mp_limb_t q = dinv * up[0]; in mpn_sbpi1_bdiv_q() 95 *qp = dinv * up[0]; in mpn_sbpi1_bdiv_q()
|
| H A D | pre_divrem_1.c | 63 mp_limb_t dinv, int shift) in mpn_preinv_divrem_1() argument 80 ASSERT (dinv == want_dinv); in mpn_preinv_divrem_1() 102 udiv_qrnnd_preinv (*qp, r, r, n0, d, dinv); in mpn_preinv_divrem_1() 127 d, dinv); in mpn_preinv_divrem_1() 131 udiv_qrnnd_preinv (*qp, r, r, n1 << shift, d, dinv); in mpn_preinv_divrem_1() 138 udiv_qrnnd_preinv (*qp, r, r, CNST_LIMB(0), d, dinv); in mpn_preinv_divrem_1()
|
| H A D | div_qr_1n_pi1.c | 165 mp_limb_t d, mp_limb_t dinv) in mpn_div_qr_1n_pi1() argument 177 udiv_qrnnd_preinv (q, uh, uh, ul, d, dinv); in mpn_div_qr_1n_pi1() 189 mp_limb_t d, mp_limb_t dinv) in mpn_div_qr_1n_pi1() argument 204 udiv_qrnnd_preinv (qp[0], u1, u1, up[0], d, dinv); in mpn_div_qr_1n_pi1() 209 B2 = -d*dinv; in mpn_div_qr_1n_pi1() 211 umul_ppmm (q1, q0, dinv, u1); in mpn_div_qr_1n_pi1() 241 umul_ppmm (p1, t, u1, dinv); in mpn_div_qr_1n_pi1() 242 add_ssaaaa (q2, q1, -u2, u2 & dinv, CNST_LIMB(0), u1); in mpn_div_qr_1n_pi1() 266 udiv_qrnnd_preinv (t, u0, u1, u0, d, dinv); in mpn_div_qr_1n_pi1()
|
| H A D | sbpi1_bdiv_r.c | 55 mp_srcptr dp, mp_size_t dn, mp_limb_t dinv) in mpn_sbpi1_bdiv_r() argument 63 ASSERT (-(dp[0] * dinv) == 1); in mpn_sbpi1_bdiv_r() 67 mp_limb_t q = dinv * up[0]; in mpn_sbpi1_bdiv_r()
|
| H A D | sbpi1_bdiv_qr.c | 56 mp_srcptr dp, mp_size_t dn, mp_limb_t dinv) in mpn_sbpi1_bdiv_qr() argument 64 ASSERT (-(dp[0] * dinv) == 1); in mpn_sbpi1_bdiv_qr() 69 mp_limb_t q = dinv * up[0]; in mpn_sbpi1_bdiv_qr()
|
| H A D | tdiv_qr.c | 109 gmp_pi1_t dinv; in mpn_tdiv_qr() local 141 invert_pi1 (dinv, d2p[dn - 1], d2p[dn - 2]); in mpn_tdiv_qr() 143 mpn_sbpi1_div_qr (qp, n2p, nn, d2p, dn, dinv.inv32); in mpn_tdiv_qr() 148 mpn_dcpi1_div_qr (qp, n2p, nn, d2p, dn, &dinv); in mpn_tdiv_qr() 271 invert_pi1 (dinv, d2p[qn - 1], d2p[qn - 2]); in mpn_tdiv_qr() 273 mpn_sbpi1_div_qr (qp, n2p, 2 * qn, d2p, qn, dinv.inv32); in mpn_tdiv_qr() 275 mpn_dcpi1_div_qr (qp, n2p, 2 * qn, d2p, qn, &dinv); in mpn_tdiv_qr()
|
| H A D | sbpi1_divappr_q.c | 47 mp_limb_t dinv) in mpn_sbpi1_divappr_q() argument 96 udiv_qr_3by2 (q, n1, n0, n1, np[1], np[0], d1, d0, dinv); in mpn_sbpi1_divappr_q() 142 udiv_qr_3by2 (q, n1, n0, n1, np[1], np[0], d1, d0, dinv); in mpn_sbpi1_divappr_q() 186 udiv_qr_3by2 (q, n1, n0, n1, np[1], np[0], d1, d0, dinv); in mpn_sbpi1_divappr_q()
|
| H A D | div_qr_2.c | 279 gmp_pi1_t dinv; in mpn_div_qr_2() local 294 gmp_pi1_t dinv; in mpn_div_qr_2() local 295 invert_pi1 (dinv, d1, d0); in mpn_div_qr_2() 296 return mpn_div_qr_2n_pi1 (qp, rp, np, nn, d1, d0, dinv.inv32); in mpn_div_qr_2() 311 invert_pi1 (dinv, d1, d0); in mpn_div_qr_2() 312 return mpn_div_qr_2u_pi1 (qp, rp, np, nn, d1, d0, shift, dinv.inv32); in mpn_div_qr_2()
|
| H A D | sec_pi1_div.c | 81 mp_limb_t dinv, in FNAME() argument 124 umul_ppmm (q1h, dummy, nh, dinv); in FNAME() 132 umul_ppmm (q0h, dummy, nh, dinv); in FNAME()
|
| H A D | pre_mod_1.c | 43 mpn_preinv_mod_1 (mp_srcptr up, mp_size_t un, mp_limb_t d, mp_limb_t dinv) in mpn_preinv_mod_1() argument 58 udiv_rnnd_preinv (r, r, n0, d, dinv); in mpn_preinv_mod_1()
|
| H A D | sbpi1_div_q.c | 46 mp_limb_t dinv) in mpn_sbpi1_div_q() argument 99 udiv_qr_3by2 (q, n1, n0, n1, np[1], np[0], d1, d0, dinv); in mpn_sbpi1_div_q() 145 udiv_qr_3by2 (q, n1, n0, n1, np[1], np[0], d1, d0, dinv); in mpn_sbpi1_div_q() 189 udiv_qr_3by2 (q, n1, n0, n1, np[1], np[0], d1, d0, dinv); in mpn_sbpi1_div_q()
|
| /netbsd-src/external/lgpl3/gmp/dist/mpz/ |
| H A D | powm_ui.c | 59 mod (mp_ptr np, mp_size_t nn, mp_srcptr dp, mp_size_t dn, gmp_pi1_t *dinv, mp_ptr tp) in mod() argument 69 mpn_div_qr_2n_pi1 (qp, np, np, nn, dp[1], dp[0], dinv->inv32); in mod() 74 mpn_sbpi1_div_qr (qp, np, nn, dp, dn, dinv->inv32); in mod() 81 mpn_dcpi1_div_qr (qp, np, nn, dp, dn, dinv); in mod() 107 reduce (mp_ptr tp, mp_srcptr ap, mp_size_t an, mp_srcptr mp, mp_size_t mn, gmp_pi1_t *dinv) in reduce() argument 115 mod (rp, an, mp, mn, dinv, scratch); in reduce() 131 gmp_pi1_t dinv; in mpz_powm_ui() local 167 invert_pi1 (dinv, mp[mn - 1], m2); in mpz_powm_ui() 176 reduce (new_bp, bp, bn, mp, mn, &dinv); in mpz_powm_ui() 215 mod (tp, tn, mp, mn, &dinv, scratch); in mpz_powm_ui() [all …]
|
| /netbsd-src/external/lgpl3/gmp/dist/tests/mpn/ |
| H A D | t-divrem_1.c | 67 mp_limb_t dinv, got_r, got_q[numberof(data[0].want_q)]; in check_data() local 97 dinv = refmpn_invert_limb (data[i].d << shift); in check_data() 100 data[i].d, dinv, shift); in check_data() 106 mp_limb_trace (" dinv", dinv); in check_data()
|
| H A D | t-bdiv.c | 138 mp_limb_t dinv; in main() local 224 binvert_limb (dinv, dp[0]); in main() 245 rh = mpn_sbpi1_bdiv_qr (qp, rp, nn, dp, dn, -dinv); in main() 253 mpn_sbpi1_bdiv_q (qp, rp, nn - dn, dp, MIN(dn,nn-dn), -dinv); in main() 260 mpn_sbpi1_bdiv_r (rp, nn, dp, dn, -dinv); in main() 271 rh = mpn_dcpi1_bdiv_qr (qp, rp, nn, dp, dn, -dinv); in main() 282 mpn_dcpi1_bdiv_q (qp, rp, nn - dn, dp, MIN(dn,nn-dn), -dinv); in main()
|
| H A D | t-div.c | 134 gmp_pi1_t dinv; in main() local 244 invert_pi1 (dinv, dnp[dn - 1], dnp[dn - 2]); in main() 265 qp[nn - dn] = mpn_sbpi1_div_qr (qp, rp, nn, dnp, dn, dinv.inv32); in main() 275 qp[nn - dn] = mpn_sbpi1_divappr_q (qp, rp, nn, dnp, dn, dinv.inv32); in main() 285 qp[nn - dn] = mpn_sbpi1_div_q (qp, rp, nn, dnp, dn, dinv.inv32); in main() 347 qp[nn - dn] = mpn_dcpi1_div_qr (qp, rp, nn, dnp, dn, &dinv); in main() 359 qp[nn - dn] = mpn_dcpi1_divappr_q (qp, rp, nn, dnp, dn, &dinv); in main() 371 qp[nn - dn] = mpn_dcpi1_div_q (qp, rp, nn, dnp, dn, &dinv); in main()
|
| /netbsd-src/external/lgpl3/gmp/dist/mpn/sparc64/ |
| H A D | divrem_1.c | 180 mp_limb_t d, n1, n0, q, r, dinv, nshift, norm_rmask; in mpn_divrem_1() local 217 invert_limb (dinv, d); in mpn_divrem_1() 227 udiv_qrnnd_preinv (q, r, r, nshift, d, dinv); in mpn_divrem_1() 232 udiv_qrnnd_preinv (q, r, r, nshift, d, dinv); in mpn_divrem_1() 237 udiv_qrnnd_preinv (q, r, r, CNST_LIMB(0), d, dinv); in mpn_divrem_1()
|
| H A D | mod_1.c | 144 mp_limb_t n1, n0, r, dinv, dummy_q, nshift, norm_rmask; in mpn_mod_1_anynorm() local 170 invert_limb (dinv, d_limb); in mpn_mod_1_anynorm() 176 udiv_qrnnd_preinv (dummy_q, r, r, nshift, d_limb, dinv); in mpn_mod_1_anynorm() 182 udiv_qrnnd_preinv (dummy_q, r, r, nshift, d_limb, dinv); in mpn_mod_1_anynorm()
|
| /netbsd-src/external/lgpl3/gmp/dist/mpn/x86_64/zen/ |
| H A D | sbpi1_bdiv_r.asm | 71 push dinv C keep dinv on stack 78 imul dinv, %rdx C first quotient limb 104 mov %r8, %rdx C dinv 143 mov %r8, %rdx C dinv 182 mov (%rsp), %rdx C dinv 234 mov %r8, %rdx C dinv 272 mov (%rsp), %rdx C dinv 333 mov %r8, %rdx C dinv 373 mov (%rsp), %rdx C dinv 405 mov %r8, %rdx C dinv [all …]
|