Home
last modified time | relevance | path

Searched refs:qn (Results 1 – 25 of 103) sorted by relevance

12345

/netbsd-src/external/lgpl3/gmp/dist/mpn/generic/
H A Ddcpi1_bdiv_qr.c98 mp_size_t qn; in mpn_dcpi1_bdiv_qr() local
111 qn = nn - dn; in mpn_dcpi1_bdiv_qr()
113 if (qn > dn) in mpn_dcpi1_bdiv_qr()
117 qn -= dn; in mpn_dcpi1_bdiv_qr()
118 while (qn > dn); in mpn_dcpi1_bdiv_qr()
121 if (BELOW_THRESHOLD (qn, DC_BDIV_QR_THRESHOLD)) in mpn_dcpi1_bdiv_qr()
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()
127 if (qn != dn) in mpn_dcpi1_bdiv_qr()
129 if (qn > dn - qn) in mpn_dcpi1_bdiv_qr()
[all …]
H A Ddcpi1_divappr_q.c88 mp_size_t qn; in mpn_dcpi1_divappr_q() local
99 qn = nn - dn; in mpn_dcpi1_divappr_q()
100 qp += qn; in mpn_dcpi1_divappr_q()
104 if (qn >= dn) in mpn_dcpi1_divappr_q()
106 qn++; /* pretend we'll need an extra limb */ in mpn_dcpi1_divappr_q()
109 qn -= dn; in mpn_dcpi1_divappr_q()
110 while (qn > dn); in mpn_dcpi1_divappr_q()
112 qp -= qn; /* point at low limb of next quotient block */ in mpn_dcpi1_divappr_q()
113 np -= qn; /* point in the middle of partial remainder */ in mpn_dcpi1_divappr_q()
118 if (qn == 1) in mpn_dcpi1_divappr_q()
[all …]
H A Ddcpi1_div_qr.c95 mp_size_t qn; in mpn_dcpi1_div_qr() local
108 qn = nn - dn; in mpn_dcpi1_div_qr()
109 qp += qn; in mpn_dcpi1_div_qr()
113 if (qn > dn) in mpn_dcpi1_div_qr()
117 qn -= dn; in mpn_dcpi1_div_qr()
118 while (qn > dn); in mpn_dcpi1_div_qr()
120 qp -= qn; /* point at low limb of next quotient block */ in mpn_dcpi1_div_qr()
121 np -= qn; /* point in the middle of partial remainder */ in mpn_dcpi1_div_qr()
124 if (qn == 1) in mpn_dcpi1_div_qr()
181 if (qn == 2) in mpn_dcpi1_div_qr()
[all …]
H A Ddcpi1_bdiv_q.c97 mp_size_t qn; in mpn_dcpi1_bdiv_q() local
110 qn = nn; in mpn_dcpi1_bdiv_q()
112 if (qn > dn) in mpn_dcpi1_bdiv_q()
116 qn -= dn; in mpn_dcpi1_bdiv_q()
117 while (qn > dn); in mpn_dcpi1_bdiv_q()
120 if (BELOW_THRESHOLD (qn, DC_BDIV_QR_THRESHOLD)) in mpn_dcpi1_bdiv_q()
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()
125 if (qn != dn) in mpn_dcpi1_bdiv_q()
127 if (qn > dn - qn) in mpn_dcpi1_bdiv_q()
[all …]
H A Dmu_div_q.c72 mp_size_t qn; in mpn_mu_div_q() local
78 qn = nn - dn; in mpn_mu_div_q()
80 tp = TMP_BALLOC_LIMBS (qn + 1); in mpn_mu_div_q()
82 if (qn >= dn) /* nn >= 2*dn + 1 */ in mpn_mu_div_q()
103 for (i = 0; i < qn + 1; i++) in mpn_mu_div_q()
111 MPN_COPY (qp, tp + 1, qn); in mpn_mu_div_q()
119 mpn_mul (pp, tp + 1, qn, dp, dn); in mpn_mu_div_q()
121 cy = (qh != 0) ? mpn_add_n (pp + qn, pp + qn, dp, dn) : 0; in mpn_mu_div_q()
124 qh -= mpn_sub_1 (qp, tp + 1, qn, 1); in mpn_mu_div_q()
126 MPN_COPY (qp, tp + 1, qn); in mpn_mu_div_q()
[all …]
H A Dmu_bdiv_qr.c74 mp_size_t qn; in mpn_mu_bdiv_qr_old() local
79 qn = nn - dn; in mpn_mu_bdiv_qr_old()
82 ASSERT (qn >= 2); in mpn_mu_bdiv_qr_old()
84 if (qn > dn) in mpn_mu_bdiv_qr_old()
96 b = (qn - 1) / dn + 1; /* ceil(qn/dn), number of blocks */ in mpn_mu_bdiv_qr_old()
97 in = (qn - 1) / b + 1; /* ceil(qn/b) = ceil(qn / ceil(qn/dn)) */ in mpn_mu_bdiv_qr_old()
113 while (qn > in) in mpn_mu_bdiv_qr_old()
132 qn -= in; in mpn_mu_bdiv_qr_old()
150 mpn_mullo_n (qp, rp, ip, qn); in mpn_mu_bdiv_qr_old()
152 if (BELOW_THRESHOLD (qn, MUL_TO_MULMOD_BNM1_FOR_2NXN_THRESHOLD)) in mpn_mu_bdiv_qr_old()
[all …]
H A Dtdiv_qr.c201 mp_size_t qn; in mpn_tdiv_qr() local
209 qn = nn - dn; in mpn_tdiv_qr()
210 qp[qn] = 0; /* zero high quotient limb */ in mpn_tdiv_qr()
211 qn += adjust; /* qn cannot become bigger */ in mpn_tdiv_qr()
213 if (qn == 0) in mpn_tdiv_qr()
220 in = dn - qn; /* (at least partially) ignored # of limbs in ops */ in mpn_tdiv_qr()
229 d2p = TMP_ALLOC_LIMBS (qn); in mpn_tdiv_qr()
230 mpn_lshift (d2p, dp + in, qn, cnt); in mpn_tdiv_qr()
233 n2p = TMP_ALLOC_LIMBS (2 * qn + 1); in mpn_tdiv_qr()
234 cy = mpn_lshift (n2p, np + nn - 2 * qn, 2 * qn, cnt); in mpn_tdiv_qr()
[all …]
H A Ddiv_q.c106 mp_size_t new_nn, qn; in mpn_div_q() local
128 qn = nn - dn + 1; /* Quotient size, high limb might be zero */ in mpn_div_q()
130 if (qn + FUDGE >= dn) in mpn_div_q()
172 qp[qn - 1] = qh; in mpn_div_q()
212 tp = TMP_ALLOC_LIMBS (qn + 1); in mpn_div_q()
215 new_nn = 2 * qn + 1; in mpn_div_q()
231 new_dp = TMP_ALLOC_LIMBS (qn + 1); in mpn_div_q()
232 mpn_lshift (new_dp, dp + dn - (qn + 1), qn + 1, cnt); in mpn_div_q()
233 new_dp[0] |= dp[dn - (qn + 1) - 1] >> (GMP_NUMB_BITS - cnt); in mpn_div_q()
235 if (qn + 1 == 2) in mpn_div_q()
[all …]
H A Dmu_div_qr.c105 mp_size_t qn; in mpn_mu_div_qr() local
108 qn = nn - dn; in mpn_mu_div_qr()
109 if (qn + MU_DIV_QR_SKEW_THRESHOLD < dn) in mpn_mu_div_qr()
121 qh = mpn_mu_div_qr2 (qp, rp + nn - (2 * qn + 1), in mpn_mu_div_qr()
122 np + nn - (2 * qn + 1), 2 * qn + 1, in mpn_mu_div_qr()
123 dp + dn - (qn + 1), qn + 1, in mpn_mu_div_qr()
127 if (dn - (qn + 1) > qn) in mpn_mu_div_qr()
128 mpn_mul (scratch, dp, dn - (qn + 1), qp, qn); /* prod is dn-1 limbs */ in mpn_mu_div_qr()
130 mpn_mul (scratch, qp, qn, dp, dn - (qn + 1)); /* prod is dn-1 limbs */ in mpn_mu_div_qr()
133 cy = mpn_add_n (scratch + qn, scratch + qn, dp, dn - (qn + 1)); in mpn_mu_div_qr()
[all …]
H A Dmu_divappr_q.c85 mp_size_t qn, in; in mpn_mu_divappr_q() local
91 qn = nn - dn; in mpn_mu_divappr_q()
94 if (qn + 1 < dn) in mpn_mu_divappr_q()
96 np += dn - (qn + 1); in mpn_mu_divappr_q()
97 nn -= dn - (qn + 1); in mpn_mu_divappr_q()
98 dp += dn - (qn + 1); in mpn_mu_divappr_q()
99 dn = qn + 1; in mpn_mu_divappr_q()
103 in = mpn_mu_divappr_q_choose_in (qn, dn, 0); in mpn_mu_divappr_q()
171 mp_size_t qn; in mpn_preinv_mu_divappr_q() local
180 qn = nn - dn; in mpn_preinv_mu_divappr_q()
[all …]
H A Dmu_bdiv_q.c74 mp_size_t qn; in mpn_mu_bdiv_q_old() local
79 qn = nn; in mpn_mu_bdiv_q_old()
82 ASSERT (qn >= 2); in mpn_mu_bdiv_q_old()
84 if (qn > dn) in mpn_mu_bdiv_q_old()
97 b = (qn - 1) / dn + 1; /* ceil(qn/dn), number of blocks */ in mpn_mu_bdiv_q_old()
98 in = (qn - 1) / b + 1; /* ceil(qn/b) = ceil(qn / ceil(qn/dn)) */ in mpn_mu_bdiv_q_old()
115 qn -= in; in mpn_mu_bdiv_q_old()
117 while (qn > in) in mpn_mu_bdiv_q_old()
148 qn -= in; in mpn_mu_bdiv_q_old()
180 mpn_sub_nc (rp + dn - in, np, tp + dn, qn - (dn - in), cy); in mpn_mu_bdiv_q_old()
[all …]
H A Dremove.c90 mp_size_t pn, nn, qn, i; in mpn_remove() local
108 qn = un; in mpn_remove()
111 while (qn >= pn) in mpn_remove()
113 qp[qn] = 0; in mpn_remove()
114 mpn_bdiv_qr_wrap (qp2, tp, qp, qn + 1, pp, pn); in mpn_remove()
122 qn = qn - pn; in mpn_remove()
123 mpn_neg (qp, qp, qn+1); in mpn_remove()
125 qn += qp[qn] != 0; in mpn_remove()
135 if (nn > qn) in mpn_remove()
139 np = TMP_ALLOC_LIMBS (qn + LOG); /* powers of V */ in mpn_remove()
[all …]
H A Dsbpi1_div_q.c49 mp_size_t qn, i; in mpn_sbpi1_div_q() local
66 qn = nn - dn; in mpn_sbpi1_div_q()
67 if (qn + 1 < dn) in mpn_sbpi1_div_q()
69 dp += dn - (qn + 1); in mpn_sbpi1_div_q()
70 dn = qn + 1; in mpn_sbpi1_div_q()
77 qp += qn; in mpn_sbpi1_div_q()
88 for (i = qn - (dn + 2); i >= 0; i--) in mpn_sbpi1_div_q()
230 if (qn + 1 < dn) in mpn_sbpi1_div_q()
232 dp += dn - (qn + 1); in mpn_sbpi1_div_q()
233 dn = qn + 1; in mpn_sbpi1_div_q()
[all …]
H A Ddivexact.c50 mp_size_t qn; in mpn_divexact() local
75 qn = nn + 1 - dn; in mpn_divexact()
82 ss = (dn > qn) ? qn + 1 : dn; in mpn_divexact()
90 wp = TMP_ALLOC_LIMBS (qn + 1); in mpn_divexact()
91 mpn_rshift (wp, np, qn + 1, shift); in mpn_divexact()
95 if (dn > qn) in mpn_divexact()
96 dn = qn; in mpn_divexact()
98 tp = TMP_ALLOC_LIMBS (mpn_bdiv_q_itch (qn, dn)); in mpn_divexact()
99 mpn_bdiv_q (qp, np, qn, dp, dn, tp); in mpn_divexact()
103 mpn_neg (qp, qp, qn); in mpn_divexact()
[all …]
H A Drootrem.c211 mp_size_t qn, rn, sn, wn, nl, bn; in mpn_rootrem_internal() local
342 qn = wn + sn; in mpn_rootrem_internal()
343 qn -= qp[qn - 1] == 0; in mpn_rootrem_internal()
347 if (qn > rn || (qn == rn && (perf_pow=mpn_cmp (qp, rp, rn)) > 0)) in mpn_rootrem_internal()
357 ASSERT_ALWAYS (rn >= qn); in mpn_rootrem_internal()
379 mpn_sub (rp, rp, rn, qp, qn); in mpn_rootrem_internal()
448 qn = rn - wn; /* expected quotient size */ in mpn_rootrem_internal()
449 if (qn <= bn) { /* Divide only if result is not too big. */ in mpn_rootrem_internal()
451 qn += qp[qn] != 0; in mpn_rootrem_internal()
461 if (qn > bn || (qn == bn && (b % GMP_NUMB_BITS != 0) && in mpn_rootrem_internal()
[all …]
H A Ddcpi1_div_q.c47 mp_size_t qn; in mpn_dcpi1_div_q() local
60 qn = nn - dn; in mpn_dcpi1_div_q()
61 wp = TMP_ALLOC_LIMBS (qn + 1); in mpn_dcpi1_div_q()
69 if (qn > dn) in mpn_dcpi1_div_q()
70 mpn_mul (tp, wp + 1, qn, dp, dn); in mpn_dcpi1_div_q()
72 mpn_mul (tp, dp, dn, wp + 1, qn); in mpn_dcpi1_div_q()
74 cy = (qh != 0) ? mpn_add_n (tp + qn, tp + qn, dp, dn) : 0; in mpn_dcpi1_div_q()
77 qh -= mpn_sub_1 (qp, wp + 1, qn, 1); in mpn_dcpi1_div_q()
79 MPN_COPY (qp, wp + 1, qn); in mpn_dcpi1_div_q()
82 MPN_COPY (qp, wp + 1, qn); in mpn_dcpi1_div_q()
H A Ddivrem.c54 mp_size_t qn; in mpn_divrem() local
61 qn = nn + qxn - 1; in mpn_divrem()
62 MPN_COPY (qp, q2p, qn); in mpn_divrem()
63 ret = q2p[qn]; in mpn_divrem()
76 mp_size_t qn; in mpn_divrem() local
88 qn = nn - dn + qxn; in mpn_divrem()
89 MPN_COPY (qp, q2p, qn); in mpn_divrem()
90 qhl = q2p[qn]; in mpn_divrem()
96 qn = nn - dn; in mpn_divrem()
97 MPN_COPY (qp, q2p, qn); in mpn_divrem()
[all …]
/netbsd-src/external/cddl/dtracetoolkit/dist/
H A Ddexplorer149 if [[ `$dtrace -b1k -qn 'BEGIN { trace(pid); exit(0); }'` == "" ]]; then
244 $dtrace -qn "$header"'
254 $dtrace -qn "$header"'
275 $dtrace -qn "$header"'
285 $dtrace -qn "$header"'
299 $dtrace -qn "$header"'
309 $dtrace -qn "$header"'
324 $dtrace -qn "$header"'
333 $dtrace -qn "$header"'
344 $dtrace -qn "$header"'
[all …]
/netbsd-src/external/cddl/dtracetoolkit/dist/Bin/
H A Ddexplorer149 if [[ `$dtrace -b1k -qn 'BEGIN { trace(pid); exit(0); }'` == "" ]]; then
244 $dtrace -qn "$header"'
254 $dtrace -qn "$header"'
275 $dtrace -qn "$header"'
285 $dtrace -qn "$header"'
299 $dtrace -qn "$header"'
309 $dtrace -qn "$header"'
324 $dtrace -qn "$header"'
333 $dtrace -qn "$header"'
344 $dtrace -qn "$header"'
[all …]
/netbsd-src/external/lgpl3/gmp/dist/mpz/
H A Dcdiv_q_ui.c40 mp_size_t ns, nn, qn; in mpz_cdiv_q_ui() local
67 qn = 1; /* a white lie, fixed below */ in mpz_cdiv_q_ui()
75 qn = nn - 2 + 1; in mpz_cdiv_q_ui()
84 qn -= qp[qn - 1] == 0; qn -= qn != 0 && qp[qn - 1] == 0; in mpz_cdiv_q_ui()
97 qn = nn - (qp[nn - 1] == 0); in mpz_cdiv_q_ui()
100 SIZ(quot) = ns >= 0 ? qn : -qn; in mpz_cdiv_q_ui()
H A Dfdiv_q_ui.c38 mp_size_t ns, nn, qn; in mpz_fdiv_q_ui() local
65 qn = 1; /* a white lie, fixed below */ in mpz_fdiv_q_ui()
73 qn = nn - 2 + 1; in mpz_fdiv_q_ui()
82 qn -= qp[qn - 1] == 0; qn -= qn != 0 && qp[qn - 1] == 0; in mpz_fdiv_q_ui()
95 qn = nn - (qp[nn - 1] == 0); in mpz_fdiv_q_ui()
98 SIZ(quot) = ns >= 0 ? qn : -qn; in mpz_fdiv_q_ui()
H A Dtdiv_q_ui.c38 mp_size_t ns, nn, qn; in mpz_tdiv_q_ui() local
72 qn = nn - 2 + 1; qn -= qp[qn - 1] == 0; qn -= qn != 0 && qp[qn - 1] == 0; in mpz_tdiv_q_ui()
78 qn = nn - (qp[nn - 1] == 0); in mpz_tdiv_q_ui()
81 SIZ(quot) = ns >= 0 ? qn : -qn; in mpz_tdiv_q_ui()
H A Dcdiv_qr_ui.c40 mp_size_t ns, nn, qn; in mpz_cdiv_qr_ui() local
71 qn = 1; /* a white lie, fixed below */ in mpz_cdiv_qr_ui()
81 qn = nn - 2 + 1; in mpz_cdiv_qr_ui()
92 qn -= qp[qn - 1] == 0; qn -= qn != 0 && qp[qn - 1] == 0; in mpz_cdiv_qr_ui()
113 qn = nn - (qp[nn - 1] == 0); in mpz_cdiv_qr_ui()
116 SIZ(quot) = ns >= 0 ? qn : -qn; in mpz_cdiv_qr_ui()
H A Dfdiv_qr_ui.c38 mp_size_t ns, nn, qn; in mpz_fdiv_qr_ui() local
69 qn = 1; /* a white lie, fixed below */ in mpz_fdiv_qr_ui()
79 qn = nn - 2 + 1; in mpz_fdiv_qr_ui()
90 qn -= qp[qn - 1] == 0; qn -= qn != 0 && qp[qn - 1] == 0; in mpz_fdiv_qr_ui()
111 qn = nn - (qp[nn - 1] == 0); in mpz_fdiv_qr_ui()
114 SIZ(quot) = ns >= 0 ? qn : -qn; in mpz_fdiv_qr_ui()
H A Dtdiv_qr_ui.c39 mp_size_t ns, nn, qn; in mpz_tdiv_qr_ui() local
80 qn = nn - 2 + 1; qn -= qp[qn - 1] == 0; qn -= qn != 0 && qp[qn - 1] == 0; in mpz_tdiv_qr_ui()
95 qn = nn - (qp[nn - 1] == 0); in mpz_tdiv_qr_ui()
98 SIZ(quot) = ns >= 0 ? qn : -qn; in mpz_tdiv_qr_ui()

12345