What ( if tracked ) is the error rate for Trial Factoring
Is an error rate for Trial Factoring tracked ?
I realize that a false negative ( small factor not found/reported, when one really exists) will eventually be covered even if it takes doing a LL. What about a false positive ( small factor found/reported, when none exist) ? Is it even possible for this to happen ? Or is there some sort of very quick verification ( just do the TF calc using the factor found on the mersenne ) ? If a factor is found and no quick verification will that be the end of work on that mersenne ? Possibly missing a mersenne prime. 
Re: What ( if tracked ) is the error rate for Trial Factorin
[quote]What about a false positive ( small factor found/reported, when none exist) ?
Is it even possible for this to happen ? [/quote] There are no false positives in the factors database. This can easily be verified, I just did so today on the Sep 1 data. If you have a Linux box, you can do this verification yourself. Basically, you just need to verify that 2[sup]p[/sup] mod f == 1 using a multipleprecision arithmetic library like [url=http://www.swox.com/gmp/]GNU MP[/url]. Here p is the exponent and f is the factor. [font=courier]unzip factors.zip [/font] # creates the file factors.cmp [font=courier]decomp f 1 100000000 [/font] # creates the file [font=courier]factors[/font] (about 58 MB) [font=courier]./factorverify factors[/font] # where the [url=http://www.mersenneforum.org/showthread.php?s=&postid=12917#post12917]program factorverify.c[/url] is attached in a later post in this thread. # Expected output = none (no lines with false factors) On a 2.8 GHz P4 this takes less than 20 seconds for all 2.6 million factors in the GIMPS database. Edit: previously instead of using the factorverify.c program, this post used: [font=courier] cat factors  sed e 's/^\(.*\),\(.*\)$/(2^\11)%\2/'  xargs ./pexpr  grep n v '^0$'[/font] # Where pexpr is a small program from the [font=courier]demos[/font] directory of the GMP source code. # Previous line outputs linenumber of any falsepositives in the "factors" file. # Expected output = none. This was clumsier and slower. 
There is backend verification so you cannot get false positives. False negatives I'm sure have happened but are probably relatively rare.

Another question for those worried about missing Mersenne primes might be: have any exponents fallen through the cracks? That is, untested, untrialfactored, and unassigned.
This can be verified by combining all unique exponents from:  BAD (probably redundant, all exps in BAD should already be in LUCAS_V.TXT too)  LUCAS_V.TXT  HRF3.TXT  factors (the file produced by decomp f 1 100000000 on FACTORS.CMP)  nofactor (the file produced by decomp n 1 100000000 on NOFACTOR.CMP)  status.txt (the uptothehour file of all current assignments of Primenet)  cleared.txt (the uptothehour file of all current cleared exponents, probably entirely redundant unless you insist on catching any exponents cleared in the few days since the last database update of BAD, LUCAS_V.TXT etc. which occurs more or less weekly)  a 39line file of known Mersenne prime exponents and comparing this combined list of exponents to a file of all prime numbers less than 100,000,000. It is very easy to generate the latter file by the way... a program like sieve2310.c by John Moyer does this in a few seconds on a 2.8 GHz P4. The comparison can be done with standard Linux utilities like cut, sort, uniq, comm, sed. It turns out that there are no "missing" exponents less than 79.3 M. If we omit "nofactors" from the set of files above, we can discover if there are any exponents which have been trialfactored, but not LL tested or currently assigned by Primenet. It turns out there are indeed some, but this is easily explained as manual tests or people using nonPrimenet programs like Glucas or Mlucas on ranges that George has reserved for them. Here is the complete list of "trialfactored but notLLtested or currently assigned by Primenet" exponents less than M39. Nothing alarming here. On the other hand, if there was a "missing" exponent somewhere well behind the current leading edge of doublechecks, it would be time for conspiracy theories. :) 12495941 12808997 12809701 12809717 12822037 12822581 12822619 12822713 12823441 12824093 12827161 12827189 12827329 12827539 12827561 12827743 12827897 12827989 12828031 12829651 13110371 There are 20 "trialfactored but notLLtested or currently assigned by Primenet" exponents in the 12M range, 14 in the 13M range, 85 in the 14M range, etc. 
These are likely exponents released for manual testing by George. If they are in the nofactors database, then they have not slipped through the cracks.

1 Attachment(s)
Attached is a C program that can be used to verify factors of Mersenne exponents. It's intended for use under Unix (and uses the [url=http://www.swox.com/gmp/]GMP multipleprecision library[/url] that comes with most Linux distributions).
On a 2.8GHz Pentium 4, it verifies all 2.6 million factors in the GIMPS database in less than 20 seconds. 
[QUOTE]
Here is the complete list of "trialfactored but notLLtested or currently assigned by Primenet" exponents less than M39. Nothing alarming here. On the other hand, if there was a "missing" exponent somewhere well behind the current leading edge of doublechecks, it would be time for conspiracy theories. :) 12495941 . . . [/QUOTE] Actually, 12495941 has already been tested (by user Amperio) and (hopefully) successfully DCed by me. I grabbed that exponent a while ago, immediately after it had expired, but some days later it was removed from my worktodo.ini and was no longer assigned to me on the Primenet status.txt report. I suppose that the original owner finished the test, and, because I had not yet reported any work, Primenet unassigned it. As I had already started working, I kept on crunching and then sent the result by email to George so he could credit it as an early DC. Hopefully, in the next release of lucas_v.txt we will see it... The hrf3.txt file from the 21st of September has the exponent. It had expired a couple of days before. 
All times are UTC. The time now is 06:14. 
Powered by vBulletin® Version 3.8.11
Copyright ©2000  2021, Jelsoft Enterprises Ltd.