[genius] Wed Apr 20 15:46:47 2022 Jiri (George) Lebl <jirka 5z com>



commit a21eaed555ba0edf2d27c47b85899975f6ae917f
Author: Jiri (George) Lebl <jiri lebl gmail com>
Date:   Wed Apr 20 15:46:50 2022 -0500

    Wed Apr 20 15:46:47 2022  Jiri (George) Lebl <jirka 5z com>
    
            * lib/number_theory/primes.gel: Update to newest data from GIMPS

 ChangeLog                    | 4 ++++
 lib/number_theory/primes.gel | 9 +++++----
 2 files changed, 9 insertions(+), 4 deletions(-)
---
diff --git a/ChangeLog b/ChangeLog
index 2061099d..b2388a90 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,7 @@
+Wed Apr 20 15:46:47 2022  Jiri (George) Lebl <jirka 5z com>
+
+       * lib/number_theory/primes.gel: Update to newest data from GIMPS
+
 Wed Oct 27 15:47:14 2021  Jiri (George) Lebl <jirka 5z com>
 
        * Release 1.0.27
diff --git a/lib/number_theory/primes.gel b/lib/number_theory/primes.gel
index e436b695..a7293378 100644
--- a/lib/number_theory/primes.gel
+++ b/lib/number_theory/primes.gel
@@ -108,7 +108,8 @@ MersennePrimeExponents = [2
 43112609
 57885161
 74207281
-77232917];
+77232917
+82589933];
 
 # See http://www.mersenne.org
 SetHelp("IsMersennePrimeExponent","number_theory","Test if Mp is a Mersenne prime using a table");
@@ -123,10 +124,10 @@ function IsMersennePrimeExponent(p) = (
                return true;
 
        # http://www.mersenne.org / GIMPS doublechecked everything up
-       # to 42 866 599 (on Apr 7, 2018)
-       if p < 42866599 then
+       # to 60 141 677 (on Apr 20, 2022)
+       if p < 60141677 then
                return false;
 
-       error("IsMersennePrimeExponent: Number too large (known values up to: " + 42866599 + ")");
+       error("IsMersennePrimeExponent: Number too large (known values up to: " + 60141677 + ")");
        bailout
 );


[Date Prev][Date Next]   [Thread Prev][Thread Next]   [Thread Index] [Date Index] [Author Index]