mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Sierpinski/Riesel Base 5 (https://www.mersenneforum.org/forumdisplay.php?f=54)
-   -   Could something be wrong with Linux sr5sieve? (https://www.mersenneforum.org/showthread.php?t=7321)

jasong 2007-03-10 18:54

Could something be wrong with Linux sr5sieve?
 
1315 factors found 5700-5750, expected 1440.

That's a pretty extreme difference, in my opinion. What do you guys think?

masser 2007-03-10 19:00

No, that's not too bad. I wouldn't be too worried.

geoff 2007-03-10 23:53

[QUOTE=jasong;100466]1315 factors found 5700-5750, expected 1440.

That's a pretty extreme difference, in my opinion. What do you guys think?[/QUOTE]

Don't take the expected number of factors too seriously, the formula is not very accurate. It is calculated as N*(1-log(pmin)/log(pmax)) where N is the number of remaining terms, and pmin <= p <= pmax the range of factors tested.

Does anyone know a better estimate?

geoff 2007-03-12 03:13

I you are concerned that factors are being missed, there is a small test sieve sr5check.tar.gz at [url]http://www.geocities.com/g_w_reynolds/sr5sieve/tests[/url]. There are instructions in the enclosed README, but briefly if you download sr5check.tar.gz and run `sr5sieve -i sr5check.txt -p 100e6 -P 150e6' it should find 10533 new factors, and the factors.txt file should match the checked-factors.txt file (which has been verified with NewPGen 2.82).

Another test in included in the sr5sieve source (see INSTALL): running `make realclean' then `make ARCH=k8 check' will test the modular arithmetic functions against GMP (for the -k8 build in this case). You need to have libgmp installed for this test.

Of course there are many corner cases that are not covered, but if both of these tests pass then that should give a reasonable assurance that most factors are being found in normal use.


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

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