From bdd13ff8c07fa2a535c6e66b4df7a1832e5974ff Mon Sep 17 00:00:00 2001 From: Paul Eggert Date: Fri, 4 Jul 2025 11:52:36 -0700 Subject: [PATCH] factor: port back to mini-gmp mini-gmp lacs mpn_tdiv_qr, so supply an emulation of it when using mini-gmp. * src/factor.c (copy_mpn_from_mpz, mpn_tdiv_qr) [!mpn_tdiv_qr]: New functions. --- src/factor.c | 24 ++++++++++++++++++++++++ 1 file changed, 24 insertions(+) diff --git a/src/factor.c b/src/factor.c index 4049297298..766de24925 100644 --- a/src/factor.c +++ b/src/factor.c @@ -615,6 +615,30 @@ mpz_va_init (void (*mpz_single_init)(mpz_t), ...) } #endif +#ifndef mpn_tdiv_qr + +static void +copy_mpn_from_mpz (mp_limb_t *p, mp_size_t n, mpz_t z) +{ + mp_size_t zsize = mpz_size (z); + mpn_copyi (p, mpz_limbs_read (z), zsize); + mpn_zero (p + zsize, n - zsize); +} + +static void +mpn_tdiv_qr (mp_limb_t *qp, mp_limb_t *rp, MAYBE_UNUSED mp_size_t qxn, + mp_limb_t const *np, mp_size_t nn, + mp_limb_t const *dp, mp_size_t dn) +{ + mpz_t q, r, n, d; + mpz_inits (q, r, nullptr); + mpz_tdiv_qr (q, r, mpz_roinit_n (n, np, nn), mpz_roinit_n (d, dp, dn)); + copy_mpn_from_mpz (qp, nn - dn + 1, q); + copy_mpn_from_mpz (rp, dn, r); + mpz_clears (q, r, nullptr); +} +#endif + static struct mp_factors mp_factor (mpz_t); static struct mp_factors -- 2.47.3