mersenneforum.org  

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

Reply
 
Thread Tools
Old 2019-12-08, 01:03   #23
a1call
 
a1call's Avatar
 
"Rashid Naimi"
Oct 2015
Remote to Here/There

3×5×137 Posts
Default

Well, I do have a detailed algorithm which may or may not speed thins up, but in short Pari-code does the sieving and if the candidates don't have small factors it writes the input file for the PFGW and launches it using the system command. I even have counter outputs for multi-threading, but at this stage the PFGW is too fast to use it. It takes longer to read and write to harddisk than finish the primality test,
a1call is offline   Reply With Quote
Old 2019-12-08, 01:18   #24
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

3,739 Posts
Default

Code:
273249*6*(2^4423-1)-1  is 3-PRP! (0.0079s+0.0031s)                                                                 
273249*6*(2^4423-1)+1 is 3-PRP! (0.0079s+0.0029s)
Found in 10 minutes with ./pfgw64 -N -f mersenne_twin where mersenne_twin contains:

Code:
ABC2 $a*6*(2^4423-1)-1 & $a*6*(2^4423-1)+1
a: from 1 to 100000000
paulunderwood is offline   Reply With Quote
Old 2019-12-08, 01:42   #25
a1call
 
a1call's Avatar
 
"Rashid Naimi"
Oct 2015
Remote to Here/There

3×5×137 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
Code:
273249*6*(2^4423-1)-1  is 3-PRP! (0.0079s+0.0031s)                                                                 
273249*6*(2^4423-1)+1 is 3-PRP! (0.0079s+0.0029s)
Found in 10 minutes with ./pfgw64 -N -f mersenne_twin where mersenne_twin contains:

Code:
ABC2 $a*6*(2^4423-1)-1 & $a*6*(2^4423-1)+1
a: from 1 to 100000000

I am still at 196998 since last night and that's for M4253.

There are ranges where sieving for small factors on PARI would be faster than PFGW, I am pretty sure.

Last fiddled with by a1call on 2019-12-08 at 01:49
a1call is offline   Reply With Quote
Old 2019-12-08, 02:12   #26
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

36·13 Posts
Talking

Quote:
Originally Posted by a1call View Post

I am still at 196998 since last night and that's for M4253.

There are ranges where sieving for small factors on PARI would be faster than PFGW, I am pretty sure.
1. Just write a simple sieve using libgmp. ...or use Polysieve? it will be way faster.

2. Any number less than 10,000 digits won't be worth the paper it is written on. Why not just start with M44497 (...and that's only for a warm up) ?
Batalov is offline   Reply With Quote
Old 2019-12-08, 02:55   #27
a1call
 
a1call's Avatar
 
"Rashid Naimi"
Oct 2015
Remote to Here/There

205510 Posts
Default

Quote:
Originally Posted by Batalov View Post
1. Just write a simple sieve using libgmp. ...or use Polysieve? it will be way faster.

2. Any number less than 10,000 digits won't be worth the paper it is written on. Why not just start with M44497 (...and that's only for a warm up) ?
I am aware that any twins less than 50k dd won't make it to top 10 listing, but there are 3 points. I remember I read an article somewhere that Tibetans believed that if they finish counting to some number the world would come to an end. I can't find that article but this is interesting:

Quote:
if you ask a Tibetan to do finger-counting, they won't bend down one finger at a time, rather they will use their thumb to count the phalanges of the finger (the three bones that make up every digit).
https://theculturetrip.com/asia/chin...etan-language/

* There's is a value to a complete set even if it includes single digit primes.
* I am doing this as a hobby
* There is already indications towards the facts that there is a bias for small k values which is an indication for infinitude of twin primes.

Plus I have 3 cores working for the 80k Mersenne for days now and on last check they only found a single non-twin prime so far.

Thank you for the sieving tips. I will look into them.

Last fiddled with by a1call on 2019-12-08 at 02:58
a1call is offline   Reply With Quote
Old 2019-12-08, 03:18   #28
a1call
 
a1call's Avatar
 
"Rashid Naimi"
Oct 2015
Remote to Here/There

3×5×137 Posts
Default

407635.6.M4253 +/- 1 are a pair of minimal Twin-Primes for M4253 with 1287 dd each.

Found Using PFGW and PaulUnderwood.


ETA:
Looks like someone fed these to factorDB more than a year ago:

http://factordb.com/index.php?query=...4253-1%29%2B1+
http://factordb.com/index.php?query=...5E4253-1%29-1+

Last fiddled with by a1call on 2019-12-08 at 03:22
a1call is offline   Reply With Quote
Old 2019-12-08, 03:25   #29
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

E9B16 Posts
Default

Quote:
Originally Posted by a1call View Post
407635.6.M4253 +/- 1 are a pair of minimal Twin-Primes for M4253 with 1287 dd each.

Found Using PFGW and PaulUnderwood.


ETA:
Looks like someone fed these to factorDB more than a year ago:

http://factordb.com/index.php?query=...4253-1%29%2B1+
http://factordb.com/index.php?query=...5E4253-1%29-1+
Quote:
Primality proving
Proven by certificate
Size 202,645 bytes
Status Verified
Uploaded by Edwin Hall
Date July 4, 2012, 4:39 am
A BLS (Brillhart-Lehmer-Selfridge) proof is quicker than ECPP. This can be done by putting the prime factors in a helper file:

Code:
cat > mersenne_twin_4423.helper
2
3
3
3
97
113
2^4423-1
Code:
./pfgw64 -N -tp -q"273249*6*(2^4423-1)-1" -h"mersenne_twin_4423.helper"
PFGW Version 4.0.1.64BIT.20191203.x86_Dev [GWNUM 29.8]

Primality testing 273249*6*(2^4423-1)-1 [N+1, Brillhart-Lehmer-Selfridge]                                    
Reading factors from helper file mersenne_twin_4423.helper
Running N+1 test using discriminant 3, base 3+sqrt(3)                                         
273249*6*(2^4423-1)-1 is prime! (0.1872s+0.0158s)
Code:
./pfgw64 -N -t -q"273249*6*(2^4423-1)+1" -h"mersenne_twin_4423.helper"
PFGW Version 4.0.1.64BIT.20191203.x86_Dev [GWNUM 29.8]

Primality testing 273249*6*(2^4423-1)+1 [N-1, Brillhart-Lehmer-Selfridge]                                    
Reading factors from helper file mersenne_twin_4423.helper
Running N-1 test using base 2                                                                 
273249*6*(2^4423-1)+1 is prime! (0.0355s+0.0001s)

Last fiddled with by paulunderwood on 2019-12-08 at 03:43
paulunderwood is offline   Reply With Quote
Old 2019-12-08, 06:10   #30
a1call
 
a1call's Avatar
 
"Rashid Naimi"
Oct 2015
Remote to Here/There

3·5·137 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
Quote:

Primality proving
Proven by certificate
Size 202,645 bytes
Status Verified
Uploaded by Edwin Hall
Date July 4, 2012, 4:39 am

I have no idea where that quote is from (except that it seems to be from the "The Prime Database"), but it certainly does not seem to be from any source that is indexed by Google.
a1call is offline   Reply With Quote
Old 2019-12-08, 06:26   #31
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

1110100110112 Posts
Default

Quote:
Originally Posted by a1call View Post
I have no idea where that quote is from (except that it seems to be from the "The Prime Database"), but it certainly does not seem to be from any source that is indexed by Google.
It comes from factorDB for that number.
paulunderwood is offline   Reply With Quote
Old 2019-12-09, 01:31   #32
a1call
 
a1call's Avatar
 
"Rashid Naimi"
Oct 2015
Remote to Here/There

3×5×137 Posts
Default

FTR:

For M9689 I rewrote my Pari-GP code so that it does sieving in Pari-GP but without multitthreading (writing out a counter where multiple instances can read a common counter ), and it is much faster to sieve in Pari and feeding it to PFGW, than running ABC2 in PFGW alone.

FTR:
I still don't know how Paul got the July 2012 information.
in this link:
http://factordb.com/index.php?query=...4253-1%29%2B1+

If i click on More information I get:
Quote:
Others:
Create time
Before November 4, 2018, 12:20 am

Last fiddled with by a1call on 2019-12-09 at 01:32
a1call is offline   Reply With Quote
Old 2019-12-09, 01:37   #33
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

72338 Posts
Default

Quote:
Originally Posted by a1call View Post
I still don't know how Paul got the July 2012 information.
in this link:
http://factordb.com/index.php?query=...4253-1%29%2B1+

If i click on More information I get:
It seems to have vanished.

Last fiddled with by paulunderwood on 2019-12-09 at 01:40
paulunderwood is offline   Reply With Quote
Reply

Thread Tools


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


Fri Jul 16 17:09:33 UTC 2021 up 49 days, 14:56, 1 user, load averages: 1.95, 1.63, 1.54

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.