From: Paul Eggert <eggert@cs.ucla.edu>
Date: Fri, 4 Jul 2025 18:52:36 +0000 (-0700)
Subject: factor: port back to mini-gmp
X-Git-Tag: v9.8~214
X-Git-Url: http://git.ipfire.org/?a=commitdiff_plain;h=bdd13ff8c07fa2a535c6e66b4df7a1832e5974ff;p=thirdparty%2Fcoreutils.git

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.
---

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