[genius] Fri Jan 12 23:27:23 2018 Jiri (George) Lebl <jirka 5z com>



commit 4737ce8f41ed2fba6f47b76fbc7ec93e4f86e8ae
Author: Jiri (George) Lebl <jiri lebl gmail com>
Date:   Fri Jan 12 23:27:31 2018 -0800

    Fri Jan 12 23:27:23 2018  Jiri (George) Lebl <jirka 5z com>
    
        * lib/number_theory/primes.gel: update for 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 ae216d3..9fd357d 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,7 @@
+Fri Jan 12 23:27:23 2018  Jiri (George) Lebl <jirka 5z com>
+
+       * lib/number_theory/primes.gel: update for newest data from GIMPS
+
 Thu Oct 12 17:52:41 2017  Jiri (George) Lebl <jirka 5z com>
 
        * examples/plane-curves.gel: Simple brute force algorithm for
diff --git a/lib/number_theory/primes.gel b/lib/number_theory/primes.gel
index b948ead..ae9d20a 100644
--- a/lib/number_theory/primes.gel
+++ b/lib/number_theory/primes.gel
@@ -107,7 +107,8 @@ MersennePrimeExponents = [2
 42643801
 43112609
 57885161
-74207281];
+74207281
+77232917];
 
 # See http://www.mersenne.org
 SetHelp("IsMersennePrimeExponent","number_theory","Test if Mp is a Mersenne prime using a table");
@@ -122,10 +123,10 @@ function IsMersennePrimeExponent(p) = (
                return true;
 
        # http://www.mersenne.org / GIMPS doublechecked everything up
-       # to 38,201,927 (on Dec 20, 2016)
-       if p < 38201927 then
+       # to 42 120 083 (on Jan 12, 2018)
+       if p < 42120083 then
                return false;
 
-       error("IsMersennePrimeExponent: Number too large (known values up to: " + 38201927 + ")");
+       error("IsMersennePrimeExponent: Number too large (known values up to: " + 42120083 + ")");
        bailout
 );


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