[Gmp-commit] /var/hg/gmp: 3 new changesets

[mercurial at gmplib.org]

details:   /var/hg/gmp/rev/6b1cd8580be6
changeset: 18567:6b1cd8580be6
user:      Marco Bodrato <bodrato at mail.dm.unipi.it>
date:      Tue Jul 21 18:42:58 2026 +0200
description:
Whitespace correction

details:   /var/hg/gmp/rev/ae5d56fd7e47
changeset: 18568:ae5d56fd7e47
user:      Marco Bodrato <bodrato at mail.dm.unipi.it>
date:      Tue Jul 21 18:43:31 2026 +0200
description:
-typo

details:   /var/hg/gmp/rev/ed17cc6e34b9
changeset: 18569:ed17cc6e34b9
user:      Marco Bodrato <bodrato at mail.dm.unipi.it>
date:      Tue Jul 21 18:44:16 2026 +0200
description:
mpn/generic/bsqrtinv.c: Adapt assert to usage, the first limb is overwritten anyway

diffstat:

 mpn/generic/bsqrt.c    |  2 +-
 mpn/generic/bsqrtinv.c |  2 +-
 mpz/perfsqrt.c         |  2 +-
 3 files changed, 3 insertions(+), 3 deletions(-)

diffs (36 lines):

diff -r d2888762e87d -r ed17cc6e34b9 mpn/generic/bsqrt.c
--- a/mpn/generic/bsqrt.c	Mon Jul 20 14:30:32 2026 +0200
+++ b/mpn/generic/bsqrt.c	Tue Jul 21 18:44:16 2026 +0200
@@ -59,7 +59,7 @@
     }
   else
     {
-      mp_size_t nn= 1 + nb / GMP_NUMB_BITS;
+      mp_size_t nn = 1 + nb / GMP_NUMB_BITS;
       if (nn <= 4) /* THINK: should we tune this? */
 	{
 	  mp_ptr sp = tp + nn;
diff -r d2888762e87d -r ed17cc6e34b9 mpn/generic/bsqrtinv.c
--- a/mpn/generic/bsqrtinv.c	Mon Jul 20 14:30:32 2026 +0200
+++ b/mpn/generic/bsqrtinv.c	Tue Jul 21 18:44:16 2026 +0200
@@ -88,7 +88,7 @@
   ASSERT (bnb > 0);
 
 #ifndef BSQRTINV_RP_NOT_ZEROED
-  ASSERT (mpn_zero_p (rp, 1 + bnb / GMP_NUMB_BITS));
+  ASSERT ((bnb <= GMP_NUMB_BITS) || mpn_zero_p (rp, 1 + bnb / GMP_NUMB_BITS));
 #endif
   if (UNLIKELY (bnb == 1))
     {
diff -r d2888762e87d -r ed17cc6e34b9 mpz/perfsqrt.c
--- a/mpz/perfsqrt.c	Mon Jul 20 14:30:32 2026 +0200
+++ b/mpz/perfsqrt.c	Tue Jul 21 18:44:16 2026 +0200
@@ -52,7 +52,7 @@
   if (! mpn_probab_perfect_square_p (PTR (op), u_size))
     return 0;
 
-  /* This is the precise size of sqrt(op) because leading limp is non-zero. */
+  /* This is the precise size of sqrt(op) because leading limb is non-zero. */
   mp_size_t rop_size = (SIZ (op) + 1) / 2;
 
   /* mpn_sqrtrem doesn't allow sp == np */