| /netbsd/src/external/lgpl3/gmp/dist/mpn/generic/ |
| D | divrem.c | 36 mpn_divrem (mp_ptr qp, mp_size_t qxn, in mpn_divrem() argument
40 ASSERT (qxn >= 0); in mpn_divrem()
45 ASSERT (! MPN_OVERLAP_P (qp, nn-dn+qxn, np, nn) || qp==np+dn+qxn); in mpn_divrem()
46 ASSERT (! MPN_OVERLAP_P (qp, nn-dn+qxn, dp, dn)); in mpn_divrem()
58 q2p = TMP_ALLOC_LIMBS (nn + qxn); in mpn_divrem()
60 np[0] = mpn_divrem_1 (q2p, qxn, np, nn, dp[0]); in mpn_divrem()
61 qn = nn + qxn - 1; in mpn_divrem()
70 return mpn_divrem_2 (qp, qxn, np, nn, dp); in mpn_divrem()
80 if (UNLIKELY (qxn != 0)) in mpn_divrem()
83 TMP_ALLOC_LIMBS_2 (n2p, nn + qxn, in mpn_divrem()
[all …]
|
| D | divrem_2.c | 55 mpn_divrem_2 (mp_ptr qp, mp_size_t qxn, in mpn_divrem_2() argument
65 ASSERT (qxn >= 0); in mpn_divrem_2()
67 ASSERT (! MPN_OVERLAP_P (qp, nn-2+qxn, np, nn) || qp >= np+2); in mpn_divrem_2()
92 qp += qxn; in mpn_divrem_2()
103 if (UNLIKELY (qxn != 0)) in mpn_divrem_2()
105 qp -= qxn; in mpn_divrem_2()
106 for (i = qxn - 1; i >= 0; i--) in mpn_divrem_2()
|
| D | divrem_1.c | 89 mpn_divrem_1 (mp_ptr qp, mp_size_t qxn, in mpn_divrem_1() argument
97 ASSERT (qxn >= 0); in mpn_divrem_1()
101 ASSERT (MPN_SAME_OR_SEPARATE_P (qp+qxn, up, un)); in mpn_divrem_1()
103 n = un + qxn; in mpn_divrem_1()
136 for (i = qxn - 1; i >= 0; i--) in mpn_divrem_1()
157 for (i = qxn - 1; i >= 0; i--) in mpn_divrem_1()
216 for (i = qxn - 1; i >= 0; i--) in mpn_divrem_1()
245 for (i = qxn - 1; i >= 0; i--) in mpn_divrem_1()
|
| D | tdiv_qr.c | 47 mpn_tdiv_qr (mp_ptr qp, mp_ptr rp, mp_size_t qxn, in mpn_tdiv_qr() argument
50 ASSERT_ALWAYS (qxn == 0); in mpn_tdiv_qr()
55 ASSERT (! MPN_OVERLAP_P (qp, nn - dn + 1 + qxn, np, nn)); in mpn_tdiv_qr()
56 ASSERT (! MPN_OVERLAP_P (qp, nn - dn + 1 + qxn, dp, dn)); in mpn_tdiv_qr()
|
| /netbsd/src/external/lgpl3/mpfr/dist/src/ |
| D | mpfr-mini-gmp.c | 180 mpn_divrem_1 (mp_limb_t *qp, mp_size_t qxn, mp_limb_t *np, mp_size_t nn, in mpn_divrem_1() argument
191 if (qxn == 0) in mpn_divrem_1()
198 mpz_init2 (n, (nn + qxn) * GMP_NUMB_BITS); in mpn_divrem_1()
199 mpn_copyi (n->_mp_d + qxn, np, nn); in mpn_divrem_1()
200 mpn_zero (n->_mp_d, qxn); in mpn_divrem_1()
201 n->_mp_size = nn + qxn; in mpn_divrem_1()
206 if (q->_mp_size < nn + qxn) in mpn_divrem_1()
207 mpn_zero (qp + q->_mp_size, nn + qxn - q->_mp_size); in mpn_divrem_1()
211 if (qxn != 0) in mpn_divrem_1()
249 mpn_tdiv_qr (mp_limb_t *qp, mp_limb_t *rp, mp_size_t qxn, in mpn_tdiv_qr() argument
[all …]
|
| /netbsd/src/external/lgpl3/gmp/dist/tests/mpn/ |
| D | t-divrem_1.c | 34 mp_size_t qxn; in check_data() member
73 qsize = data[i].nsize + data[i].qxn; in check_data()
76 got_r = mpn_divrem_1 (got_q, data[i].qxn, data[i].n, data[i].nsize, in check_data()
86 printf (" qxn=%ld\n", (long) data[i].qxn); in check_data()
98 got_r = mpn_preinv_divrem_1 (got_q, data[i].qxn, in check_data()
|
| /netbsd/src/external/lgpl3/gmp/dist/mpn/ia64/ |
| D | divrem_1.asm | 61 C qxn = r33
95 shladd r32 = r15, 3, r32 C r32 = rp + n + qxn
98 adds r32 = -8, r32 C r32 = rp + n + qxn - 1
155 shladd r14 = r15, 3, r32 C r14 = rp + n + qxn
158 adds r32 = -8, r14 C r32 = rp + n + qxn - 1
|
| /netbsd/src/external/lgpl3/gmp/dist/tests/ |
| D | refmpn.c | 1654 refmpn_divrem_2 (mp_ptr qp, mp_size_t qxn, in refmpn_divrem_2() argument
1661 tp = refmpn_malloc_limbs (nn + qxn); in refmpn_divrem_2()
1662 refmpn_zero (tp, qxn); in refmpn_divrem_2()
1663 refmpn_copyi (tp + qxn, np, nn); in refmpn_divrem_2()
1664 qh = refmpn_sb_div_qr (qp, tp, nn + qxn, dp, 2); in refmpn_divrem_2()
2340 refmpn_tdiv_qr (mp_ptr qp, mp_ptr rp, mp_size_t qxn, in refmpn_tdiv_qr() argument
2344 ASSERT (qxn == 0); in refmpn_tdiv_qr()
|
| /netbsd/src/external/lgpl3/gmp/dist/doc/ |
| D | gmp.texi | 5359 @deftypefun void mpn_tdiv_qr (mp_limb_t *@var{qp}, mp_limb_t *@var{rp}, mp_size_t @var{qxn}, const …
5367 @var{qxn} operand must be zero.
5370 @deftypefun mp_limb_t mpn_divrem (mp_limb_t *@var{r1p}, mp_size_t @var{qxn}, mp_limb_t *@var{rs2p},…
5379 In addition to an integer quotient, @var{qxn} fraction limbs are developed, and
5380 stored after the integral limbs. For most usages, @var{qxn} will be zero.
5390 The area at @var{r1p} needs to be @var{rs2n} @minus{} @var{s3n} + @var{qxn}
5394 @deftypefn Function mp_limb_t mpn_divrem_1 (mp_limb_t *@var{r1p}, mp_size_t @var{qxn}, @w{mp_limb_t…
5399 The integer quotient is written to @{@var{r1p}+@var{qxn}, @var{s2n}@} and in
5400 addition @var{qxn} fraction limbs are developed and written to @{@var{r1p},
5401 @var{qxn}@}. Either or both @var{s2n} and @var{qxn} can be zero. For most
[all …]
|
| /netbsd/src/external/lgpl3/gmp/dist/ |
| D | ChangeLog | 14487 * mpn/generic/divrem.c: Rate qxn test as UNLIKELY.
31992 * mpn/generic/divrem.c: qextra_limbs => qxn.
|