[genius] * lib/number_theory/primes.gel: update data from GIMPS (new prime)



commit 95b164231fdcfd7fe8401201bb94a10bdbd47480
Author: Jiri (George) Lebl <jiri lebl gmail com>
Date:   Wed Jan 20 09:37:42 2016 -0600

        * lib/number_theory/primes.gel: update data from GIMPS (new prime)

 ChangeLog                    |    4 ++++
 lib/number_theory/primes.gel |    9 +++++----
 2 files changed, 9 insertions(+), 4 deletions(-)
---
diff --git a/ChangeLog b/ChangeLog
index ad728b0..6528d19 100644
--- a/ChangeLog
+++ b/ChangeLog
@@ -1,3 +1,7 @@
+Wed Jan 20 09:37:38 2016  Jiri (George) Lebl <jirka 5z com>
+
+       * lib/number_theory/primes.gel: update data from GIMPS (new prime)
+
 Wed Jan 06 13:49:05 2016  Jiri (George) Lebl <jirka 5z com>
 
        * Release 1.0.21
diff --git a/lib/number_theory/primes.gel b/lib/number_theory/primes.gel
index e067db1..b3db114 100644
--- a/lib/number_theory/primes.gel
+++ b/lib/number_theory/primes.gel
@@ -106,7 +106,8 @@ MersennePrimeExponents = [2
 37156667
 42643801
 43112609
-57885161];
+57885161
+74207281];
 
 # See http://www.mersenne.org
 SetHelp("IsMersennePrimeExponent","number_theory","Test if Mp is a Mersenne prime using a table");
@@ -121,10 +122,10 @@ function IsMersennePrimeExponent(p) = (
                return true;
 
        # http://www.mersenne.org / GIMPS doublechecked everything up
-       # to 33,560,269 (on Feb 26, 2015)
-       if p < 33560269 then
+       # to 34,969,871 (on Jan 20, 2016)
+       if p < 34969871 then
                return false;
 
-       error("IsMersennePrimeExponent: Number too large (known values up to: " + 33560269 + ")");
+       error("IsMersennePrimeExponent: Number too large (known values up to: " + 34969871 + ")");
        bailout
 );


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