mersenneforum.org  

Go Back   mersenneforum.org > Great Internet Mersenne Prime Search > Software

Reply
 
Thread Tools
Old 2017-07-12, 21:27   #1
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

2×2,969 Posts
Default Prime constellations?

What's the best software for searching for prime constellations? I have constellations of the form kn + {k0, k1, ..., kj} for constants k, k0, ..., kj, where k0, ..., kj are all small (in my case, k = 23# and 1 = k0 < ... < kj < 100).

I seem to remember PFGW had a feature to do something like this, but maybe there is something more specialized? Any suggestions?
CRGreathouse is offline   Reply With Quote
Old 2017-07-13, 00:40   #2
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

1000010112 Posts
Default

APSieve can do that or write a own sieve. APSieve is not the fastest. What is the exact term ?

Last fiddled with by Cybertronic on 2017-07-13 at 00:44
Cybertronic is offline   Reply With Quote
Old 2017-07-13, 01:45   #3
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

134628 Posts
Default

Incidentally what I meant above was that kj - k0 is small, not that (say) kj is small.

Quote:
Originally Posted by Cybertronic View Post
APSieve can do that or write a own sieve. APSieve is not the fastest. What is the exact term ?
Well, I actually have ~200 constellations. For example,

223092870n + 217153 + {0, 4, 6, 10, 16, 18, 24, 28, 30, 34, 40, 46, 48, 54, 58, 60, 66, 70, 76, 84, 88, 94, 96}
CRGreathouse is offline   Reply With Quote
Old 2017-07-13, 05:55   #4
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

3·89 Posts
Default

"223092870n + 217153 + {0, 4, 6, 10, 16, 18, 24, 28, 30, 34, 40, 46, 48, 54, 58, 60, 66, 70, 76, 84, 88, 94, 96}"

Is it an AP project ?
I can send to you a link for downloading APSieve. Maybe it is useful.

Last fiddled with by Cybertronic on 2017-07-13 at 06:19
Cybertronic is offline   Reply With Quote
Old 2017-07-13, 10:51   #5
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

3×89 Posts
Default

Okay, looks like a tuplet project, but a 23 -tuplet have d=94, not 96 ?!

23 94 : 0 4 6 10 12 16 24 30 34 40 42 46 52 54 60 66 70 72 76 82 84 90 94
23 94 : 0 4 10 12 18 22 24 28 34 40 42 48 52 54 60 64 70 78 82 84 88 90 94

( ktpatt.txt at T.Forbes site )

I wrote for my 15 tuplet a sieve.
To find an example is not easy. Maybe n up to 10^30 one case.
I wrote a test sieve for a 18-tuplet. No chance. All primes up to 10M into the sieve, maybe 10 n's left for PRPing in an 100 billion -n-block.
BTW ,are you Mr. Stocker ?

Last fiddled with by Cybertronic on 2017-07-13 at 11:34
Cybertronic is offline   Reply With Quote
Old 2017-07-13, 16:58   #6
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

2×2,969 Posts
Default

Quote:
Originally Posted by Cybertronic View Post
Is it an AP project ?
It's just something I'm doing on my own.

Quote:
Originally Posted by Cybertronic View Post
I can send to you a link for downloading APSieve. Maybe it is useful.
OK, great!

Quote:
Originally Posted by Cybertronic View Post
Okay, looks like a tuplet project, but a 23 -tuplet have d=94, not 96 ?!

23 94 : 0 4 6 10 12 16 24 30 34 40 42 46 52 54 60 66 70 72 76 82 84 90 94
23 94 : 0 4 10 12 18 22 24 28 34 40 42 48 52 54 60 64 70 78 82 84 88 90 94
Yes, it's a little different from those. I have other constraints that give me a slightly larger final term (hence my collection is not a prime 23-tuplet by the standard definition).

Quote:
Originally Posted by Cybertronic View Post
To find an example is not easy. Maybe n up to 10^30 one case.
I only need one case between all 200, so 0.01 or 0.02 solutions on average would be OK. But without specialized sieving software I can't even get that high (it took 8 hours for my naive GP script to search all forms up to 3e16).

Quote:
Originally Posted by Cybertronic View Post
BTW ,are you Mr. Stocker ?
No, I'm some other guy.
CRGreathouse is offline   Reply With Quote
Old 2017-07-13, 17:21   #7
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

3·89 Posts
Default

Here is the download link

https://www.sendspace.com/file/idqrzc

Read readme.txt , there are useful instructions.

Problems, ask me. There are some bugs in APSieve....

Best
Cybertronic is offline   Reply With Quote
Old 2017-07-13, 18:16   #8
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

2·2,969 Posts
Default

Well, I'm running it now with
Code:
apsieve -v -npb -ex23# -b23,217153,217157,217159,217163,217169,217171,217177,217181,217183,217187,217193,217199,217201,217207,217211,217213,217219,217223,217229,217237,217241,217247,217249 -k100000000 -K4000000000
but it's a pity that it's limited to 4 billion.
CRGreathouse is offline   Reply With Quote
Old 2017-07-13, 18:35   #9
CRGreathouse
 
CRGreathouse's Avatar
 
Aug 2006

2×2,969 Posts
Default

Oh wait, this is better because it PRPs as well:
Code:
apsieve -v 134473145 4000000000 23# 23 217153 217157 217159 217163 217169 217171 217177 217181 217183 217187 217193 217199 217201 217207 217211 217213 217219 217223 217229 217237 217241 217247 217249 1000000000

Last fiddled with by CRGreathouse on 2017-07-13 at 18:56
CRGreathouse is offline   Reply With Quote
Old 2017-07-13, 19:50   #10
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

10B16 Posts
Default

Can you find a command - line for this ?

for n=1000000 to n=10000000 , 1009# with 2 b's 16061 , 16063 and stop maxP=100000000 without PRP-Test -> so left n write APSieve in output.txt automatic

Never I found a solution for this with P<2^31 . Program ever crashed. Better, reload an output.txt file , set maxP <2^31

...and yes, 4 billion is maximum but 40 x 100 Million in order is faster via bat file

Last fiddled with by Cybertronic on 2017-07-13 at 20:42
Cybertronic is offline   Reply With Quote
Old 2017-07-14, 09:45   #11
Cybertronic
 
Cybertronic's Avatar
 
Jan 2007
DEUTSCHLAND !

3×89 Posts
Default

Found a way:

apsieve -npb -k1000000000 -K5000000000 -ex1200# -b4,16061,16063,16067,16069 -mP100000000

Stop at 100M without PRP -> write the output.txt file

PRPing with PFGW is twice faster totay, so its better to start pfgwxx -f0 output.txt after sieving

But reload an output.txt and set maxP < 2^31 is not possible.
Cybertronic is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
(M48) NEW MERSENNE PRIME! LARGEST PRIME NUMBER DISCOVERED! dabaichi News 571 2020-10-26 11:02
Prime Constellations 2 MattcAnderson MattcAnderson 25 2018-07-31 23:11
Prime Constellations MattcAnderson MattcAnderson 119 2018-03-14 20:22
disk died, prime work lost forever? where to put prime? on SSD or HDD? emily PrimeNet 3 2013-03-01 05:49
The 40th known Mersenne prime, 220996011-1 is not PRIME! illman-q Miscellaneous Math 33 2004-09-19 05:02

All times are UTC. The time now is 13:38.

Thu Nov 26 13:38:39 UTC 2020 up 77 days, 10:49, 3 users, load averages: 1.35, 1.50, 1.52

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2020, 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.