mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Riesel Prime Search (https://www.mersenneforum.org/forumdisplay.php?f=59)
-   -   Choose your own K and work on finding a top-5000 prime! (https://www.mersenneforum.org/showthread.php?t=4963)

kar_bon 2020-11-27 08:20

@diep
I've created the [url='https://www.rieselprime.de/ziki/Riesel_prime_32767']Riesel k=32767[/url] in the Wiki and reserved it for you (including tested to n=50k).
I've also updated the k=89 reservation. If you got any further results for k=32767 please fill them in.

diep 2020-11-27 10:02

Kar, c'est bon!

You can be sure i will report it if there is a new one above n=1569 :)

Please note for the few who might not know: testspeed of LLR with larger k's is a lot slower than small k's. k = 89 and k=69 here tests roughly factor 2 faster than k=32767, under the same conditions.

At the CPU's with many cores here what matters most is that LLR FFT-size stays within combined L3 cachesize : #cores * L3_cache_a_core

I do not know details yet the FFT-size seems dependant upon which k-factor you use.

(so if you know the FFT size for that specific k-factor and n-size in units multiply it with 16 bytes - because the size is in doubles and you need 1 double which has t he information and 1 to write new information to,
and a double is 64 bits so that's 8 bytes and you need 2 of that makes 16 bytes. The blabla from manufacturers that you can add up the L2 and L1 cachesize to it seems utter blabla from the manufacturers to me by the way yet that's not relevant here - just L3 cachesize matters there in that comparision - except if you do not have a L3 cache then L2 cache size matters. So to say it technical the algorithm doesn't seem to work in-situ regrettably)

Happy5214 2020-11-27 23:32

[QUOTE=diep;564551]I do not know details yet the FFT-size seems dependant upon which k-factor you use.[/QUOTE]

According to that LLRTools tool I linked, to go from the max [I]n[/I] for [I]k[/I]=1 of a given FFT length (which I'll call [$]n_m[/$]) to the max [I]n[/I] for any [I]k[/I] for that same FFT length (which I'll call [$]n_k[/$]), the formula is [TEX]n_k = n_m - (\log_2 k + \frac{1}{2}(\log_2 k)(fftlen))[/TEX].

Happy5214 2020-11-28 07:33

1 Attachment(s)
Here is the promised script. It's a little slow, and I haven't tested it on Windows (I don't currently have access to a Windows box), but it's completely automated once started. It uses a binary search (which it prints out) to find the minimum [I]n[/I] for every FFT length below a given limit (provided on the command line; you can give the exact number or a "K" or "M" suffixed one like LLR prints out; the actual length doesn't have to be present on that platform). It dynamically finds the relevant FFT lengths for that platform through the search. Create a "temp" subfolder as a working directory for it to dump all of the temp files (it creates a ton of them). It uses the CLI version of LLR (assuming its name is "llr" in the current directory). Let me know if anyone finds this useful.

diep 2020-11-30 20:45

[QUOTE=Happy5214;564652]Here is the promised script. It's a little slow, and I haven't tested it on Windows (I don't currently have access to a Windows box), but it's completely automated once started. It uses a binary search (which it prints out) to find the minimum [I]n[/I] for every FFT length below a given limit (provided on the command line; you can give the exact number or a "K" or "M" suffixed one like LLR prints out; the actual length doesn't have to be present on that platform). It dynamically finds the relevant FFT lengths for that platform through the search. Create a "temp" subfolder as a working directory for it to dump all of the temp files (it creates a ton of them). It uses the CLI version of LLR (assuming its name is "llr" in the current directory). Let me know if anyone finds this useful.[/QUOTE]

How to use the script?
What do i type on the command line?

something like this?

python happy.py 32767 5800000

Nescro 2020-12-01 07:04

k=2805
 
Releasing k=8847 and reserving k=2805.

Happy5214 2020-12-01 08:23

1 Attachment(s)
[QUOTE=diep;564875]How to use the script?
What do i type on the command line?

something like this?

python happy.py 32767 5800000[/QUOTE]
To be honest, I use Linux, so I'd just set the executable flag and the shebang (the first comment line with the Unix-style interpreter path) would do the rest. You appear to be using Windows, though, given that invocation. It uses Python 3, so I think you'd call "python3" instead of "python" (which calls Python 2). At least that's how it works on Ubuntu.

You have the arguments flipped. The [I]max FFT length[/I] ([I]not[/I] the max [I]n[/I]) is the first one, since it's the only required one. The [I]k[/I] you choose for the optional second argument shouldn't change the result significantly unless it starts zero-padding the FFTs (high 6-digits), and I chose [I]k[/I]=100005 to avoid zero-padding while still shortening the timing tests (since the [I]n[/I]'s will be lower). [I]n[/I]=5.8M is at the 480k FFT length for your [I]k[/I] on my Core 2 Quad, and the Xeon won't be too much different, if at all, so run it to that or 512k. Something like "python3 happy.py 480K".

I've attached another script that reads in the generated maxlen.txt and converts the max [I]n[/I]'s (which are adjusted in the file to match Mersenne values) for any chosen [I]k[/I]. They should closely match the test results from the verbose printouts of the first script. It's called using "python3 fftlen.py [I]k[/I]".

Happy5214 2020-12-02 09:05

I posted those two scripts (renamed to generate_maxlen.py and adjusted_fftlen.py) and a few others to a new GitHub repo: [url]https://github.com/happy5214/rps-scripts[/url]. I'll post updates there. You'll need the fftlen.py library module for both scripts. I plan on adding a new argument interface to generate_maxlen.py to allow you to specify a max [I]n[/I] instead of a max FFT length, and also to specify verbosity via a flag instead of by editing the code.

Happy5214 2020-12-08 08:10

I posted an update to the generate_maxlen.py script which revamps the CLI interface. It now uses the Python argparse standard library module, which makes it a lot nicer. I've also added support for specifying a max [I]n[/I] in lieu of a max FFT length and moving the verbose flag to the CLI. It now has a built-in help message generated by the module, so check that for the details.

kar_bon 2020-12-12 10:01

Perhaps an admin should make a new thread for such script-related posts, better finding later.

Happy5214 2021-01-04 07:20

The (near-)Woodall [I]k[/I]'s listed in [url]https://www.mersenneforum.org/showpost.php?p=550539&postcount=363[/url] have been completed to [I]n[/I]=425k. The following primes were found:[list][*]1993191*2^405251-1[*]582833*2^417560-1[*]197673*2^421871-1[/list]

Happy5214 2021-01-14 05:11

1 Attachment(s)
I've completed the remaining RPS 9th and 10th Drive [I]k[/I]'s with missing ranges from [I]n[/I]=325k to 350k. 14 primes were found, which are attached. I brought an old laptop into the fold, and I have a new Kubuntu Focus on the way, so this project (to [I]n[/I]=400k) may be done by this summer.

bur 2021-03-01 09:52

[QUOTE=Happy5214;568280]The (near-)Woodall [I]k[/I]'s listed in [url]https://www.mersenneforum.org/showpost.php?p=550539&postcount=363[/url] have been completed to [I]n[/I]=425k. The following primes were found:[/QUOTE]What are (near-)Woodall k's? I know about near-Woodall primes where k is n-1, but that doesn't seem to be the case here.

kar_bon 2021-03-01 10:11

(General) Near Woodalls are primes of type (n-1)*b^n-1 or (n+1)*b^n-1, examples [url='https://www.rieselprime.de/ziki/Gen_Woodall_prime_M_2_table']here[/url] or [url='https://www.rieselprime.de/ziki/Gen_Woodall_prime_P_2_table']here[/url].
Most of the are not spotted immediately, because the normal form in for example the Top5000 page is displayed different:
Example: 1060*2^1059-1 is normalized written as 265*2^1061-1

Happy5214 2021-03-02 03:41

To be clear, these are [I]k[/I]'s where near-Woodall and Woodall primes were found for larger [I]n[/I] values, and I'm filling in the untested missing ranges. I'm not actually finding new primes of those forms, which are searched for largely by Steven Harvey and PrimeGrid, respectively.

manrus 2021-03-29 09:30

Reserving k=3317

Happy5214 2021-04-03 16:44

As of March 27, the (near-)Woodall [I]k[/I]'s listed in [url]https://www.mersenneforum.org/showpost.php?p=550539&postcount=363[/url] have been completed to [I]n[/I]=450k. The following primes were found:[list][*]1268979*2^426051-1[*]1195203*2^444882-1[*]197673*2^447700-1[/list]

gd_barnes 2021-04-06 03:52

[QUOTE=Happy5214;540140]I'll also tack on k=14549535 for 125k < n < 300k as a double-check.[/QUOTE]


Happy,

From a year ago...was this done? I'd like to go ahead and double-check k=14549535 for all n<=300K and then test it up to n=500K. Is that OK?

Gary

Happy5214 2021-04-06 19:02

[QUOTE=gd_barnes;575291]Happy,

From a year ago...was this done? I'd like to go ahead and double-check k=14549535 for all n<=300K and then test it up to n=500K. Is that OK?

Gary[/QUOTE]
Funny you mention that. I had started work on that [I]k[/I] in late February, but I had suspended work on my local PRPNet servers for a bit in March to work on the GIMPS PRP-CF-DC backlog. I have since resumed work on both clients. That [I]k[/I] is currently at [I]n[/I]=288k, and I should get to [I]n[/I]=300k by the end of this week. You can take over from there.

gd_barnes 2021-04-06 19:24

[QUOTE=Happy5214;575349]Funny you mention that. I had started work on that [I]k[/I] in late February, but I had suspended work on my local PRPNet servers for a bit in March to work on the GIMPS PRP-CF-DC backlog. I have since resumed work on both clients. That [I]k[/I] is currently at [I]n[/I]=288k, and I should get to [I]n[/I]=300k by the end of this week. You can take over from there.[/QUOTE]


OK I'll start testing within a couple of days for n=300k-500k. Thanks.

gd_barnes 2021-04-09 04:51

k=14549535 is complete for n=300K-500K, 5 primes were found shown below, the k is released. Happy is still working on a double-check for n=125K-300K.

Primes:
14549535*2^376607-1
14549535*2^390393-1
14549535*2^430198-1
14549535*2^446171-1
14549535*2^453217-1

Happy5214 2021-04-10 05:53

The aforementioned double-check of [I]k[/I]=14549535 for [I]n[/I]=125k-300k is complete. In addition to the known primes, the following new primes were found:[list][*]14549535*2^149070-1[*]14549535*2^160953-1[*]14549535*2^161852-1[*]14549535*2^180548-1[*]14549535*2^199042-1[*]14549535*2^204177-1[/list]

gd_barnes 2021-05-10 20:48

2 Attachment(s)
k=8847:

Double-check complete to n=500K. No problems found.

Hole filled for n=500K-700K. In addition to confirming the 3 previously found primes, one additional prime was found:
8847*2^514298-1

Results attached.

gd_barnes 2021-05-10 20:57

2 Attachment(s)
k=32767:

This one is currently reserved by Diep for n=5M-30M. Since it is so low weight I thought it would be interesting to see if there were any more primes up to n=1M.

Double-check complete to n=50K. 4 previously found primes confirmed. No problems found.

First-time check complete for n=50K-1M. No primes found.

This is a very low-weight k! :-)

Results attached.

Happy5214 2021-05-13 02:33

As of May 9, the (near-)Woodall [I]k[/I]'s listed in [url]https://www.mersenneforum.org/showpost.php?p=550539&postcount=363[/url] have been completed to [I]n[/I]=475k. The following primes were found:[list][*]667071*2^450745-1[*]938237*2^455888-1[*]667071*2^465478-1[*]582833*2^467596-1[/list]

gd_barnes 2021-05-19 19:15

4 Attachment(s)
k=10005 has a large prime gap so I double-checked it for n=0-1.25M.

34 primes confirmed; no problems found.

Results attached.

Happy5214 2021-06-02 04:42

1 Attachment(s)
I've completed the remaining RPS 9th and 10th Drive [I]k[/I]'s with missing ranges from [I]n[/I]=350k to 375k. 21 primes were found, which are attached. Of the laptops mentioned in my last post on this set of [I]k[/I]'s, the Kubuntu Focus has been dedicated to GIMPS and is not working on this project, but the older laptop is actively working on this. The expected finish date is currently October, though I may dedicate some additional cycles to it once it gets closer to the finish line.

gd_barnes 2021-06-23 07:04

5 Attachment(s)
Back about 13-14 years ago Beyond (SB2) did a lot of searching on many low-weight k's. Although he gave details about what he was working on his details were frequently incomplete and he did not complete a lot of what he was sieving. One of the main issues is that we never knew where he started his searches. Sometimes they started at n=1, other times at n=300K, many times we never knew. I found these 4 main status posts of his:

post link / date
[URL]https://www.mersenneforum.org/showpost.php?p=100670[/URL] 2007-03-13
[URL]https://www.mersenneforum.org/showpost.php?p=101240[/URL] 2007-03-18
[URL]https://www.mersenneforum.org/showpost.php?p=110378[/URL] 2007-07-14
[URL]https://www.mersenneforum.org/showpost.php?p=153292[/URL] 2008-12-14

I took all of the k's found in those posts and removed some k's that were already checked by others since then. This left 98 k's all of which are shown below.

All 98 k's have been double-checked by me for n=1 to 1M. 5 missing primes were found for ranges that he said were done and 13 new primes were found for ranges that he sieved but never searched or failed to report.

Missing primes:
9613*2^734907-1
124679*2^12-1
239857*2^161901-1
334147*2^103669-1
36231101*2^34482-1

New primes:
273613*2^397407-1
313979*2^243184-1
412843*2^107943-1
412843*2^142983-1
412843*2^377223-1
434729*2^941408-1
451049*2^281984-1
464909*2^136620-1
469933*2^330667-1
473279*2^103840-1
700477*2^492061-1
783073*2^402635-1
1034503*2^174283-1

All k's searched with all primes listed below.
[B][COLOR=Red]Red = missing[/COLOR][/B]
[COLOR=Red][COLOR=Black][B]Bold = new[/B][/COLOR][/COLOR]
[code]
k / n-primes
2039: 24, 64, 319920, 841312
2171: 14, 202, 1594, 6754, 101974
2279: 900, 6572, 23772, 35252
2591: 134, 326, 830, 2174, 17270
2741: 50, 138, 1610, 1850, 3954, 10202, 23178, 84858, 468594, 737586
2759: 388, 405820, 569812
2797: 21, 6093, 7389, 9861
2983: 11, 671, 767, 2891, 771455
3343: 483, 1299
3389: 92, 848, 2684, 23864, 209996, 504584, 749828
3539: 12, 192, 4236, 869868
3817: 40381, 161197, 218917
3821: 70, 15814, 439006
3823: 7, 287, 335, 599, 1015, 2407, 384223
4327: 13, 37, 133, 1189
4597: 17, 37, 161, 221, 721, 112009, 799781
4813: 515, 11379, 13155
5077: 117, 333, 490629, 881829
5099: 184, 1216, 9936, 19072, 27648
5323: 27, 83, 8651, 107643
6119: 4544, 678080
7331: 34, 87154, 704842
9179: 4120, 10408, 41680
9613: 1479, 34239, 464439, 636159, [B][COLOR=Red]734907[/COLOR][/B]
9913: 1187, 8387
10949: 1492, 1516
11519: 164444
12791: 2818, 29398, 308026
17077: 677
19919: 888, 2688, 8460, 170820
21547: 1, 121, 2233, 16001, 310625, 936473
21673: 7, 775, 129175
22183: 727, 2471, 773447
24067: 1753, 2069, 2881, 700073
24671: 1290, 2250, 11698, 108906, 382066
25229: 238652
30727: 1121, 3353, 15929
43541: 507098, 850394
48973: 43, 115, 259, 32119
65279: 1640, 4952, 66632, 625704
70079: 1428
106303: 1939
114613: 211, 1687, 461551
117037: 1029
124679: [B][COLOR=red]12[/COLOR][/B], 948, 2460, 10308, 497244, 769908
132217: 15281
171607: 21117, 73245
179743: 99731
192509: 52
236377: 19693, 111973, 178285, 824773
239857: 33, 7233, [B][COLOR=red]161901[/COLOR][/B]
248047: 11613
249127: 5189, 31301
273613: 39, 87, 159, [B]397407[/B]
278713: 17159, 57719
279731: 1782
284579: 56, 4736
294907: 545, 203921, 221409
306251: 22, 29038, 37462
313979: 4, 184, [B]243184[/B]
320107: 41, 257, 8249
323563: 243, 9099
332159: 4416
334147: 37, 7981, [B][COLOR=red]103669[/COLOR][/B]
334331: 78, 294, 7998
362737: 153
370421: 201442
384029: 24
403993: 379
412843: 15, 16047, [B]107943, 142983, 377223[/B]
434729: 8, 20, 80, [B]941408[/B]
441907: 837, 42861
451049: 7748, 22784, [B]281984[/B]
464353: 191, 527, 1919, 371279
464909: 73740, [B]136620[/B]
465187: 3949, 16441
469933: 51, [B]330667[/B]
473279: 160, 92032, [B]103840[/B]
492787: 553, 1865, 6961, 58849
516569: 20
648433: 75
700477: 973, [B]492061[/B]
780427:
783073: 11, [B]402635[/B]
808477: 11673
842711: 226, 3034, 3754
844559:
893887: 257
1005449: 204, 41100
1034503: 43, 5971, [B]174283[/B]
1210421: 28674, 655938
10013593: 3847
10108837: 33657
10247561: 730
10284899: 868, 2596
10346561: 32714
36231101: [B][COLOR=red]34482[/COLOR][/B], 469002
37616317:
[/code] Karsten, this cleans up a former mess of k's, many of which need to be entered into the Wiki but couldn't because their accuracy was in question. I will slowly work on getting them all into the Wiki going in descending order from the largest k. Larger k's have less history and so are easier to enter. I will be referencing Beyond's work, my own current and former work, the old low-weight stats page done by Joss, and other misc. work that has popped up by other people or projects that I find from searches on the k's in this forum including work done by NPLB for k=2000-3000.

Attached are primes and residues from the effort.

Gary

Happy5214 2021-06-28 02:03

The (near-)Woodall [I]k[/I]'s listed in [url]https://www.mersenneforum.org/showpost.php?p=550539&postcount=363[/url] have been completed to [I]n[/I]=500k. The following primes were found:[list][*]197673*2^475399-1[*]1467763*2^480087-1[*]665127*2^483394-1[*]1268979*2^485869-1[*]1268979*2^487073-1[*]665127*2^488224-1[*]938237*2^491576-1[*]1183953*2^496687-1[/list]
I want to reiterate my reservation of [I]k[/I]=454483, which will resume work from this point.


All times are UTC. The time now is 15:48.

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