mersenneforum.org  

Go Back   mersenneforum.org > Prime Search Projects > And now for something completely different

Reply
 
Thread Tools
Old 2015-12-12, 00:13   #12
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

100101000001012 Posts
Default

ok (for independence sake, constructed it with gp from scratch)
...running a 2-PRP test now...
...done
Generic modular reduction using generic reduction FMA3 FFT length 96K, Pass1=384, Pass2=256 on A 971486-bit number
2357111317...619511619537619543619561619 is 2-PRP! (3880.7502s+4.8723s)
...running a 11-PRP test now...
...done
Generic modular reduction using generic reduction FMA3 FFT length 96K, Pass1=384, Pass2=256 on A 971486-bit number
2357111317...619511619537619543619561619 is 11-PRP! (3552.4437s+2.4134s)
Attached Files
File Type: zip CE294k.zip (113.2 KB, 128 views)

Last fiddled with by Batalov on 2015-12-12 at 02:42
Batalov is offline   Reply With Quote
Old 2015-12-12, 01:04   #13
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

3,739 Posts
Default

Please attach the decimal representation of the PRP
paulunderwood is offline   Reply With Quote
Old 2015-12-12, 02:02   #14
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

11·577 Posts
Default

Quote:
Originally Posted by Batalov View Post
ok (for independence sake, constructed it with gp from scratch)
...running a 2-PRP test now...
...done
2357111317...19511619537619543619561619 is 2-PRP! (3880.7502s+4.8723s)
Thank you Serge!
rogue is online now   Reply With Quote
Old 2015-12-12, 03:03   #15
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

373910 Posts
Default

Thanks for the attachment, Serge. After about 50 mins using 80% of a 3.5GHz core, my GWNUM based program (using a=0, since Jacobi(-1,n)==-1) reports "Likely prime!" meaning (x+2)^(n+1)==5 (mod n, x^2+1).

Last fiddled with by paulunderwood on 2015-12-12 at 03:04
paulunderwood is offline   Reply With Quote
Old 2015-12-12, 13:36   #16
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

11×577 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
Thanks for the attachment, Serge. After about 50 mins using 80% of a 3.5GHz core, my GWNUM based program (using a=0, since Jacobi(-1,n)==-1) reports "Likely prime!" meaning (x+2)^(n+1)==5 (mod n, x^2+1).
Thanks to you as well Paul!
rogue is online now   Reply With Quote
Old 2015-12-17, 11:50   #17
ericw
 
ericw's Avatar
 
Oct 2015

238 Posts
Default

Quote:
Originally Posted by rogue View Post
If I've written my code correctly then CE(292447) is PRP. Even if my code isn't correct, it is a fairly large PRP. Could someone please independently verify?
Verified also using Mathematica's PrimeQ (returns True -> probable prime)

Congratulations!
ericw is offline   Reply With Quote
Old 2015-12-17, 12:07   #18
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

11000110010112 Posts
Default

Quote:
Originally Posted by ericw View Post
Verified also using Mathematica's PrimeQ (returns True -> probable prime)

Congratulations!
How long did that take?
rogue is online now   Reply With Quote
Old 2015-12-17, 12:54   #19
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

72·131 Posts
Default

I am confused about the decimal representation, it seems to stop half-way through a number

619511
619537
619543
619561
619

I had thought you were looking at the number as a series of concatenated primes, rather than as a digit-string (made up of concatenated primes) whose prefixes you were testing for primality.
fivemack is offline   Reply With Quote
Old 2015-12-17, 14:31   #20
R. Gerbicz
 
R. Gerbicz's Avatar
 
"Robert Gerbicz"
Oct 2005
Hungary

22×7×53 Posts
Default

Quote:
Originally Posted by rogue View Post
How long did that take?
Couldn't be that too slow as Mathematica is using gmp: http://www.wolfram.com/LGPL/GMP/ .
R. Gerbicz is offline   Reply With Quote
Old 2015-12-17, 14:32   #21
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

11·577 Posts
Default

Quote:
Originally Posted by fivemack View Post
I am confused about the decimal representation, it seems to stop half-way through a number

619511
619537
619543
619561
619

I had thought you were looking at the number as a series of concatenated primes, rather than as a digit-string (made up of concatenated primes) whose prefixes you were testing for primality.
No. You are thinking about the Smarandache-Wellin search.
rogue is online now   Reply With Quote
Old 2015-12-21, 13:56   #22
ericw
 
ericw's Avatar
 
Oct 2015

1910 Posts
Default

Quote:
Originally Posted by rogue View Post
How long did that take?
32 CPU-hours on moderately old hardware. (PrimeQ uses a multiple Rabin-Miller test in bases 2 and 3 combined with a Lucas pseudoprime test.)
ericw is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Erdos numbers xilman Lounge 83 2020-12-19 15:03
Erdos-Turan Conjecture. wildrabbitt Math 4 2015-09-16 02:40
Mills' Primes, Mills' Constant Andrew Computer Science & Computational Number Theory 1 2013-02-17 22:22
Constant n Search kar_bon Riesel Prime Data Collecting (k*2^n-1) 5 2009-06-22 23:00
Erdos-Kac Theorem flouran Math 84 2009-05-20 21:37

All times are UTC. The time now is 17:11.


Fri Jul 16 17:11:40 UTC 2021 up 49 days, 14:58, 1 user, load averages: 2.95, 1.96, 1.67

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, Jelsoft Enterprises Ltd.

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.