mersenneforum.org  

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

Reply
 
Thread Tools
Old 2015-02-03, 12:23   #1
TheCount
 
TheCount's Avatar
 
Sep 2013
Perth, Au.

2·72 Posts
Thumbs up fpcsieve version 0.1.1 released

I've taken Geoffrey Reynolds fpsieve from https://sites.google.com/site/geoffr...ograms/testing and ported it to VS2012 and added Compositorial sieving support to create fpcsieve.
Windows version 0.1.1 binaries are available at:
https://github.com/RogerKarpin/fpcsieve/tree/master/bin
GitHub is new to me but looks really good.

The only changes from fpsieve to include Compositorial sieving are additions to app.c and app.h.
If someone could compile the original code base with these changes and produce linux binaries that would be appreciated.

This is a CPU sieve that has multithreaded support.
Won't port well to GPU as needs to refer to list of Composites in memory (just like Primorial does to list of Primes).
See here for more info:
http://www.primegrid.com/forum_thread.php?id=5952
I've started sieving n<1M. Currently up to 3.8G.
PM me if your interested in Compositorial prime searching.
This is one of the baseless prime searches like Factorial and Primorial.
Will need to reach n>80k to get in the top 5,000 list. Currently only searched to n=16k which is ridiculous for this day and age:
https://primes.utm.edu/primes/search...&Comment=orial
This is my first new thread post on mersenneforum. Hopefully many more to come.
TheCount is offline   Reply With Quote
Old 2015-02-06, 22:34   #2
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

2×5×23×41 Posts
Default

96835!/96835# + 1 is prime
Batalov is offline   Reply With Quote
Old 2015-02-07, 06:11   #3
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

2×11×167 Posts
Default

Quote:
Originally Posted by Batalov View Post
Congrats

Are there any gaps in this table?

Code:
-----  -------------------------------- ------- ----- ---- --------------
 rank           description              digits  who year comment
-----  -------------------------------- ------- ----- ---- --------------
 2209  96835!/96827#+1                   398821  p381 2015 Compositorial (**)
 ----  41400!/41400#-1                   155301  p381 2015 Compositorial
 ----  33684!/33684#-1                   123321  p381 2015 Compositorial
40882  20943!/20939#-1                    72401   p16 2000 Compositorial
42827  18160!/18149#-1                    61635   p16 2000 Compositorial
44240  17258!/17257#+1                    58207   p16 2000 Compositorial
47617  15250!/15241#-1                    50623   p16 2000 Compositorial
48663  13917!/13913#+1                    45631   p16 2000 Compositorial
49214  13724!/13723#-1                    44919   p16 2000 Compositorial
51963  10977!/10973#-1                    34876   p16 2000 Compositorial
56655  8711!/8707#-1                      26816    p7 1999 Compositorial
59086  6851!/6841#+1                      20374   p34 2000 Compositorial
60838  6187!/6173#+1                      18135    p7 1999 Compositorial
68706  3692!/3691#+1                       9996    p7 1999 Compositorial
-----  -------------------------------- ------- ----- ---- --------------

Last fiddled with by Batalov on 2015-02-10 at 07:10
paulunderwood is offline   Reply With Quote
Old 2015-02-07, 07:16   #4
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

100100110101102 Posts
Default

There must be!
...'cause I started from 90,200, really, and went up, only until one on either side... well, maybe I'll hunt for the minus side, as well.

Found in gaps:
Feb 9 18:39:34 PST: 33684!/33684#-1 is prime (123321 digits)
Feb 9 23:11:33 PST: 41400!/41400#-1 is prime (155301 digits)

Last fiddled with by Batalov on 2015-02-10 at 07:12 Reason: added two primes
Batalov is offline   Reply With Quote
Old 2015-02-07, 14:04   #5
TheCount
 
TheCount's Avatar
 
Sep 2013
Perth, Au.

2·72 Posts
Default

Quote:
Originally Posted by paulunderwood View Post
Are there any gaps in this table?
Was previously tested to 14000 on both the Plus and Minus sides according to OEIS.
https://oeis.org/A049420

Before implementing the sieving program I double checked this initial range and continued on to:
- Compositorial Plus side: 21706
- Compositorial Minus side: 20621

I've now completed sieving n<1M up to 10G.
969528 terms remain (total of both sides).
TheCount is offline   Reply With Quote
Old 2015-02-07, 17:46   #6
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

2·5·23·41 Posts
Default

I've briefly searched google for "Daniel Heuer compositorial" and I think some traces (and possibly DH's limits) can be found in the archives of the primeform yahoo group. E.g. this (about C(20493)-1 which you will soon double-check) ...and this appears to be all of his messages (API of yahoo groups changed so many times over years; the new API is so inconvenient).

Daniel was active in many curious projects. I stumble over his decade old contributions in many areas.

Maybe Paul has Daniel's contacts? Paul "was always there", as they say.
Batalov is offline   Reply With Quote
Old 2015-02-07, 17:58   #7
Antonio
 
Antonio's Avatar
 
"Antonio Key"
Sep 2011
UK

32·59 Posts
Default

There is a bug in the on screen status display. On my machine (Windows 7 64bit i5 3570k @ 4.4GHz) it displays: -

p=309370883, 18.19K p/sec, 3.94x0MHz CPU, 0.0% done

(This is with the 64bit version)

Question - wouldn't it sieve a lot faster if it didn't print the factors found file - fpfactors.txt?

Last fiddled with by Antonio on 2015-02-07 at 17:59 Reason: spelling correction
Antonio is offline   Reply With Quote
Old 2015-02-07, 19:18   #8
paulunderwood
 
paulunderwood's Avatar
 
Sep 2002
Database er0rr

2·11·167 Posts
Default

Quote:
Originally Posted by Batalov View Post

Daniel was active in many curious projects. I stumble over his decade old contributions in many areas.

Maybe Paul has Daniel's contacts? Paul "was always there", as they say.
http://primes.utm.edu/bios/page.php?id=223 has his email address, but Daniel may have moved on.

Last fiddled with by paulunderwood on 2015-02-07 at 19:21
paulunderwood is offline   Reply With Quote
Old 2015-02-07, 23:36   #9
rogue
 
rogue's Avatar
 
"Mark"
Apr 2003
Between here and the

142438 Posts
Default

I see no reason why this can't be built upon fpsievecl, the OpenCL implementation of fpsieve. The reason is that this form only requires the exclusion of certain numbers from the multiplication of the factorial. It is possible that global memory access might cause slowdowns though.
rogue is offline   Reply With Quote
Old 2015-02-08, 08:03   #10
TheCount
 
TheCount's Avatar
 
Sep 2013
Perth, Au.

2×72 Posts
Default

Quote:
Originally Posted by Antonio View Post
There is a bug in the on screen status display. On my machine (Windows 7 64bit i5 3570k @ 4.4GHz) it displays: -

p=309370883, 18.19K p/sec, 3.94x0MHz CPU, 0.0% done

(This is with the 64bit version)
I already posted this issue on GitHub:
https://github.com/RogerKarpin/fpcsieve/issues/2
"Speed of CPU in periodic report not working for 64bit. Need separate asm file to read CPU speed as VS2012 does not compile C file with assembler inline."
Its a low priority issue.

Quote:
Question - wouldn't it sieve a lot faster if it didn't print the factors found file - fpfactors.txt?
You can use -c option to set the checkpoint period, default is every 5 minutes. Set to 0 if you don't want it to periodically checkpoint.
I wouldn't think writing a checkpoint file would affect performance much.
TheCount is offline   Reply With Quote
Old 2015-02-08, 11:36   #11
TheCount
 
TheCount's Avatar
 
Sep 2013
Perth, Au.

1428 Posts
Default

Quote:
Originally Posted by rogue View Post
I see no reason why this can't be built upon fpsievecl, the OpenCL implementation of fpsieve. The reason is that this form only requires the exclusion of certain numbers from the multiplication of the factorial. It is possible that global memory access might cause slowdowns though.
Worth a try. I'd love to get a GPU sieve for Compositorial. I could have my GPU sieving while PRP testing on my CPU.

fpsievecl? I didn't know you'd extended it to Primorial. If you hadn't yet I would extend to Primorial while adding Compositorial. We can see if it's worthwhile after its implemented.

Is v1.0.8 the latest version? If so I already have it downloaded and know how to fix a bug:
http://www.primegrid.com/forum_threa...rap=true#80600

If Compositorial implemented would the new program be called fpcsievecl? I'd add a new GitHub repo.
TheCount is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
LLR Version 3.8.20 released Jean Penné Software 30 2018-08-13 20:00
LLR Version 3.8.15 released Jean Penné Software 28 2015-08-04 04:51
LLR Version 3.8.11 released Jean Penné Software 37 2014-01-29 16:32
LLR Version 3.8.9 released Jean Penné Software 37 2013-10-31 08:45
llr 3.8.2 released as dev-version opyrt Prime Sierpinski Project 11 2010-11-18 18:24

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

Thu May 13 13:16:43 UTC 2021 up 35 days, 7:57, 2 users, load averages: 3.08, 2.96, 2.80

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.