Home | History | Annotate | Download | only in generic

Lines Matching refs:qn

45       q = -n * d^{-1} mod 2^{qn * GMP_NUMB_BITS},
47       r = (n + q * d) * 2^{-qn * GMP_NUMB_BITS}
98   mp_size_t qn;
111   qn = nn - dn;
113   if (qn > dn)
115       /* Reduce qn mod dn without division, optimizing small operations.  */
117 	qn -= dn;
118       while (qn > dn);
121       if (BELOW_THRESHOLD (qn, DC_BDIV_QR_THRESHOLD))
122 	cy = mpn_sbpi1_bdiv_qr (qp, np, 2 * qn, dp, qn, dinv);
124 	cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, qn, dinv, tp);
127       if (qn != dn)
129 	  if (qn > dn - qn)
130 	    mpn_mul (tp, qp, qn, dp + qn, dn - qn);
132 	    mpn_mul (tp, dp + qn, dn - qn, qp, qn);
133 	  mpn_incr_u (tp + qn, cy);
135 	  rr = mpn_add (np + qn, np + qn, nn - qn, tp, dn);
139       np += qn;
140       qp += qn;
142       qn = nn - dn - qn;
145 	  rr += mpn_add_1 (np + dn, np + dn, qn, cy);
149 	  qn -= dn;
151       while (qn > 0);
156   if (BELOW_THRESHOLD (qn, DC_BDIV_QR_THRESHOLD))
157     cy = mpn_sbpi1_bdiv_qr (qp, np, 2 * qn, dp, qn, dinv);
159     cy = mpn_dcpi1_bdiv_qr_n (qp, np, dp, qn, dinv, tp);
162   if (qn != dn)
164       if (qn > dn - qn)
165 	mpn_mul (tp, qp, qn, dp + qn, dn - qn);
167 	mpn_mul (tp, dp + qn, dn - qn, qp, qn);
168       mpn_incr_u (tp + qn, cy);
170       rr = mpn_add (np + qn, np + qn, nn - qn, tp, dn);