mersenneforum.org  

Go Back   mersenneforum.org > Prime Search Projects > Riesel Prime Search

Reply
 
Thread Tools
Old 2014-09-26, 03:00   #1
Citrix
 
Citrix's Avatar
 
Jun 2003

22·397 Posts
Default Calcprimes program

I have been using the calcprimes.jar program I found on the forum to calculate the odds of finding primes. Overall I have found it very accurate. Are there any sequences that anyone knows of that beats the odds predicted by the program?

Do all the Riesel k<300 follow the odds predicted by the program? Looking at the statistics of the various drives the program looks accurate.
Citrix is offline   Reply With Quote
Old 2014-09-27, 14:07   #2
pepi37
 
pepi37's Avatar
 
Dec 2011
After milion nines:)

32·52·7 Posts
Default

Hi!
I just search for that program and in "my range" show 1.79 expected primes.
I should hope that will be 2 primes :))
Thanks for program, and if it is accurate as you say, that is good news for me :)

Quick sieve test ( on known set of primes)
REPORTED PREDICTED
2*10^n-1 28 23.84
3*10^n-1 19 20.45
5*10^n-1 25 24.74
6*10^n-1 31 27.63
8*10^n-1 23 23.42
9*10^n-1 20 18.05


so it looks pretty accurate to me :)

Last fiddled with by pepi37 on 2014-09-27 at 14:38 Reason: add more info
pepi37 is offline   Reply With Quote
Old 2014-09-27, 16:25   #3
VBCurtis
 
VBCurtis's Avatar
 
"Curtis"
Feb 2005
Riverside, CA

2×3×11×83 Posts
Default

Citrix-
A sequence that "beats the odds" is a matter of statistical distribution- there are surely a few that substantially beat the estimate so far, but we have no reason to think they will continue to do so. Like rolling dice, past performance is not connected to future results.

If you are asking if any sequence exceeds the odds and we have a reason to think it will continue to do so, that 'reason' would be publication-worthy in most cases. So, no.

Pepi-
The theory that program uses is well-known, and any inaccuracies are due to the irregular nature of the distribution of primes rather than any error/inaccuracy in the program.
VBCurtis is online now   Reply With Quote
Old 2014-09-27, 17:04   #4
Batalov
 
Batalov's Avatar
 
"Serge"
Mar 2008
Phi(4,2^7658614+1)/2

100110110101012 Posts
Default

Quote:
Originally Posted by Citrix View Post
Do all the Riesel k<300 follow the odds predicted by the program? Looking at the statistics of the various drives the program looks accurate.
I would expect those forms that have a partial algebraic factorization to need a correction coefficient (f < 1, not in the direction of 'beating the odds'). Consider 27*2^n-1.
Batalov is offline   Reply With Quote
Old 2014-09-27, 21:40   #5
Citrix
 
Citrix's Avatar
 
Jun 2003

22·397 Posts
Default

Quote:
Originally Posted by VBCurtis View Post

If you are asking if any sequence exceeds the odds and we have a reason to think it will continue to do so, that 'reason' would be publication-worthy in most cases. So, no.
Yes that is what I was wondering. I believe (though do not have proof) that sequences that have special factors like mersenne numbers, generalized fermat, GM, GQ etc are slightly more likely than regular numbers to produce primes.

The likelihood increases in the case of mersenne numbers/GM/GQ as the size of the p increases by a factor of ln(p).

It would remain constant for a Generalized fermat series depending on the exponent. The increased likelihood for generalized fermat numbers would be ln(2^x) (where 2^x is the exponent).

Is there any way of confirming this empirically. I do not have access to the mersenne factor database/generalized fermat search database to test this.

Can anyone disprove this or give a counterexample to this?
Citrix is offline   Reply With Quote
Old 2014-09-27, 21:52   #6
Mini-Geek
Account Deleted
 
Mini-Geek's Avatar
 
"Tim Sorbera"
Aug 2006
San Antonio, TX USA

2×3×23×31 Posts
Default

Quote:
Originally Posted by Batalov View Post
I would expect those forms that have a partial algebraic factorization to need a correction coefficient (f < 1, not in the direction of 'beating the odds'). Consider 27*2^n-1.
If you've eliminated those candidates from the sieve file already, I'd expect that the program is still accurate.
Mini-Geek is online now   Reply With Quote
Old 2014-09-27, 22:47   #7
VBCurtis
 
VBCurtis's Avatar
 
"Curtis"
Feb 2005
Riverside, CA

2·3·11·83 Posts
Default

Quote:
Originally Posted by Citrix View Post
The likelihood increases in the case of mersenne numbers/GM/GQ as the size of the p increases by a factor of ln(p).

Is there any way of confirming this empirically. I do not have access to the mersenne factor database/generalized fermat search database to test this.

Can anyone disprove this or give a counterexample to this?
For Mersennes, the number of expected primes in some future range has been listed on mersenne.org for many years. The incidence of primes is too small to "prove" or "disprove" any estimate method, but I believe the rate of mersenne discovery is in rough agreement with this standard method, while your factor of ln(p) increase would estimate vastly more primes than have been discovered. At the size of exponent p currently in testing (50M+), ln(p) is over 17! That would be a massive deviation from the usual estimates.
VBCurtis is online now   Reply With Quote
Old 2014-09-28, 01:20   #8
Citrix
 
Citrix's Avatar
 
Jun 2003

22×397 Posts
Default

Quote:
Originally Posted by VBCurtis View Post
For Mersennes, the number of expected primes in some future range has been listed on mersenne.org for many years. The incidence of primes is too small to "prove" or "disprove" any estimate method, but I believe the rate of mersenne discovery is in rough agreement with this standard method, while your factor of ln(p) increase would estimate vastly more primes than have been discovered. At the size of exponent p currently in testing (50M+), ln(p) is over 17! That would be a massive deviation from the usual estimates.
Let me clarify:-
Odds of a regular number being prime= 1/ln(N)
Odds of a mersenne being prime (where only prime exponents are used)= [1+ln(p)]/ln(2^p-1)

For a 50M number this would be 18/50,000,000 ~ 1 in 3 million prime exponents or close to 1 prime in a 50 million n range --> which seems about right.

I agree that there is not enough data to prove or disprove this.

Odds for a GFN number being prime=[1+ln(2^x)]/ln(b^2^x)
Citrix is offline   Reply With Quote
Old 2014-09-28, 07:49   #9
VBCurtis
 
VBCurtis's Avatar
 
"Curtis"
Feb 2005
Riverside, CA

2·3·11·83 Posts
Default

I think your observation is an obfuscated rephrasing of the structure of the factors of mersenne numbers, while not altering the odds any mersenne trial factored to n bits is prime using the formula applied by calcprimes.

Since mersenne factors are bigger than 2p, you can apply the "odds of prime" as if trial factoring has been done to 2p, which would increase the probability produced by the formula. However, it doesn't mean mersennes are any easier to find per primality test; it merely alters the prefactoring effort/effectiveness.
VBCurtis is online now   Reply With Quote
Old 2015-06-26, 17:03   #10
pepi37
 
pepi37's Avatar
 
Dec 2011
After milion nines:)

32·52·7 Posts
Default

Does same size of sieve(in MB) ( same base, but different K) suggested same number of primes in same range?
Why I asking this: I made test sieves until 1000000
for K 8 sieve size is 1.906 MB, and for K 96 is 1.907 MB
When using calcprimes. jar , I got predicted number of 27 primes for K=8 and only 21 primes for K=96
pepi37 is offline   Reply With Quote
Old 2015-06-26, 18:50   #11
Mini-Geek
Account Deleted
 
Mini-Geek's Avatar
 
"Tim Sorbera"
Aug 2006
San Antonio, TX USA

2×3×23×31 Posts
Default

Quote:
Originally Posted by pepi37 View Post
Does same size of sieve(in MB) ( same base, but different K) suggested same number of primes in same range?
Why I asking this: I made test sieves until 1000000
for K 8 sieve size is 1.906 MB, and for K 96 is 1.907 MB
When using calcprimes. jar , I got predicted number of 27 primes for K=8 and only 21 primes for K=96
Are you referring to the size of the sieved file (Newpgen or ABC format) containing the remaining candidates? I would only expect this if the K=96 file took up significantly more bytes per candidate.
Mini-Geek is online now   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
GPU LLR program jasong GPU Computing 19 2011-08-23 03:32
So you think you can program rogue Lounge 5 2009-10-02 15:02
Program for GPU tribal Information & Answers 5 2009-03-19 20:54
which program? drakkar67 Prime Sierpinski Project 14 2005-11-29 06:25
program P-1 for K*2^n-1 jocelynl 15k Search 19 2004-01-11 17:24

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


Fri Sep 30 16:28:05 UTC 2022 up 43 days, 13:56, 0 users, load averages: 1.37, 1.54, 1.55

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

≠ ± ∓ ÷ × · − √ ‰ ⊗ ⊕ ⊖ ⊘ ⊙ ≤ ≥ ≦ ≧ ≨ ≩ ≺ ≻ ≼ ≽ ⊏ ⊐ ⊑ ⊒ ² ³ °
∠ ∟ ° ≅ ~ ‖ ⟂ ⫛
≡ ≜ ≈ ∝ ∞ ≪ ≫ ⌊⌋ ⌈⌉ ∘ ∏ ∐ ∑ ∧ ∨ ∩ ∪ ⨀ ⊕ ⊗ 𝖕 𝖖 𝖗 ⊲ ⊳
∅ ∖ ∁ ↦ ↣ ∩ ∪ ⊆ ⊂ ⊄ ⊊ ⊇ ⊃ ⊅ ⊋ ⊖ ∈ ∉ ∋ ∌ ℕ ℤ ℚ ℝ ℂ ℵ ℶ ℷ ℸ 𝓟
¬ ∨ ∧ ⊕ → ← ⇒ ⇐ ⇔ ∀ ∃ ∄ ∴ ∵ ⊤ ⊥ ⊢ ⊨ ⫤ ⊣ … ⋯ ⋮ ⋰ ⋱
∫ ∬ ∭ ∮ ∯ ∰ ∇ ∆ δ ∂ ℱ ℒ ℓ
𝛢𝛼 𝛣𝛽 𝛤𝛾 𝛥𝛿 𝛦𝜀𝜖 𝛧𝜁 𝛨𝜂 𝛩𝜃𝜗 𝛪𝜄 𝛫𝜅 𝛬𝜆 𝛭𝜇 𝛮𝜈 𝛯𝜉 𝛰𝜊 𝛱𝜋 𝛲𝜌 𝛴𝜎𝜍 𝛵𝜏 𝛶𝜐 𝛷𝜙𝜑 𝛸𝜒 𝛹𝜓 𝛺𝜔