Lines Matching defs:n2p
71 mp_ptr n2p;
83 n2p = TMP_ALLOC_LIMBS (nn + 1);
84 cy = mpn_lshift (n2p, np, nn, cnt);
85 n2p[nn] = cy;
86 qhl = mpn_divrem_2 (qp, 0L, n2p, nn + (cy != 0), d2p);
89 rp[0] = (n2p[0] >> cnt)
90 | ((n2p[1] << (GMP_NUMB_BITS - cnt)) & GMP_NUMB_MASK);
91 rp[1] = (n2p[1] >> cnt);
95 n2p = TMP_ALLOC_LIMBS (nn);
96 MPN_COPY (n2p, np, nn);
97 qhl = mpn_divrem_2 (qp, 0L, n2p, nn, dp);
99 rp[0] = n2p[0];
100 rp[1] = n2p[1];
115 mp_ptr n2p, d2p;
126 n2p = TMP_ALLOC_LIMBS (nn + 1);
127 cy = mpn_lshift (n2p, np, nn, cnt);
128 n2p[nn] = cy;
135 n2p = TMP_ALLOC_LIMBS (nn + 1);
136 MPN_COPY (n2p, np, nn);
137 n2p[nn] = 0;
143 mpn_sbpi1_div_qr (qp, n2p, nn, d2p, dn, dinv.inv32);
148 mpn_dcpi1_div_qr (qp, n2p, nn, d2p, dn, &dinv);
153 mpn_mu_div_qr (qp, rp, n2p, nn, d2p, dn, scratch);
154 n2p = rp;
158 mpn_rshift (rp, n2p, dn, cnt);
160 MPN_COPY (rp, n2p, dn);
202 mp_ptr n2p, d2p;
233 n2p = TMP_ALLOC_LIMBS (2 * qn + 1);
234 cy = mpn_lshift (n2p, np + nn - 2 * qn, 2 * qn, cnt);
237 n2p[2 * qn] = cy;
238 n2p++;
242 n2p[0] |= np[nn - 2 * qn - 1] >> (GMP_NUMB_BITS - cnt);
250 n2p = TMP_ALLOC_LIMBS (2 * qn + 1);
251 MPN_COPY (n2p, np + nn - 2 * qn, 2 * qn);
254 n2p[2 * qn] = 0;
255 n2p++;
263 udiv_qrnnd (q0, r0, n2p[1], n2p[0] << GMP_NAIL_BITS, d2p[0] << GMP_NAIL_BITS);
264 n2p[0] = r0 >> GMP_NAIL_BITS;
268 mpn_divrem_2 (qp, 0L, n2p, 4L, d2p); /* FIXME: obsolete function */
273 mpn_sbpi1_div_qr (qp, n2p, 2 * qn, d2p, qn, dinv.inv32);
275 mpn_dcpi1_div_qr (qp, n2p, 2 * qn, d2p, qn, &dinv);
283 mpn_mu_div_qr (qp, r2p, n2p, 2 * qn, d2p, qn, scratch);
284 MPN_COPY (n2p, r2p, qn);
312 if (n2p[qn - 1] < h)
317 cy = mpn_add_n (n2p, n2p, d2p, qn);
321 n2p[qn] = cy;
333 cy1 = mpn_lshift (n2p, n2p, rn, GMP_NUMB_BITS - cnt);
334 n2p[0] |= np[in - 1] & (GMP_NUMB_MASK >> cnt);
337 cy2 = mpn_submul_1 (n2p, qp, qn, dp[in - 1] & (GMP_NUMB_MASK >> cnt));
340 ASSERT_ALWAYS (n2p[qn] >= cy2);
341 n2p[qn] -= cy2;
345 n2p[qn] = cy1 - cy2; /* & GMP_NUMB_MASK; */
360 MPN_COPY (rp, n2p, rn);
369 cy = mpn_sub (n2p, n2p, rn, tp + in, qn);
370 MPN_COPY (rp + in, n2p, dn - in);