20210730, 00:25  #661  
"James Heinrich"
May 2004
exNorthern Ontario
D9E_{16} Posts 
Quote:
This is what I have in PHP: Code:
function ECMeffortExpected() { return array( 20 => array( 11000, 100), 25 => array( 50000, 280), 30 => array( 250000, 640), 35 => array( 1000000, 1580), 40 => array( 3000000, 4700), 45 => array( 11000000, 9700), 50 => array( 44000000, 17100), 55 => array(110000000, 46500), 60 => array(260000000, 112000), 65 => array(800000000, 360000), ); } function ECMeffort($curves) { // https://www.mersenneforum.org/showpost.php?p=573887&postcount=569 // https://github.com/sethtroisi/miscscripts/blob/main/mersenne/ecm_progress/test.py // Try to convert list of <B1, B2, number_of_curves> to t40, t45, t50 progress // original code by SethTro <https://www.mersenneforum.org/member.php?u=15870> 20210316 // translated to PHP by James Heinrich <james@mersenne.ca> 20210317 $completed = array(); foreach (ECMeffortExpected() as $digits => $B1curvesNeeded) { $complete = 0; $curvecount = 0; list($min_B1, $curves_needed) = $B1curvesNeeded; foreach ($curves as $curvedata) { list($B1, $B2, $count) = $curvedata; if (($B1 >= $min_B1) && ($B2 >= (20 * $min_B1))) { $complete += ($count / $curves_needed); $curvecount += $count; } } if ($complete > 0.001) { $completed[$digits] = array($complete, $curvecount); } } return $completed; } function ECMprobabilityMissedSmallerFactor($curves_complete, $needed_curves) { // https://www.mersenneforum.org/showpost.php?p=573952&postcount=574 // Probability of missed (smaller) factor is // (1  1/needed_curves) ^ curves_complete // which can be transformed to // math.exp(curve_count * math.log(1  1/curve_count)) // or approximated by // (1/e)^(complete) //return exp($curves_complete * log(1  (1 / $needed_curves))); return pow(1  (1 / $needed_curves), $curves_complete); } 

20210730, 01:39  #662  
"Seth"
Apr 2019
2·181 Posts 
Quote:
https://gist.github.com/sethtroisi/6...8d6372579f09bd The output is $complete, $FacMissed 

20210730, 15:21  #663 
"James Heinrich"
May 2004
exNorthern Ontario
2·3·7·83 Posts 
Thanks for the revised code. For now I have placed the old and new calculations sidebyside so you can compare and review them and see if the new numbers appear to make sense.
However with the new calculation method and lookup table I wasn't able to figure out how to replicate the B1 and Curves columns. I have left in the B1 calculated from the old lookup table but removed the Curves column. I don't understand this stuff quite well enough to be certain this makes sense, so please review and suggest changes. 
20210817, 06:14  #666 
"Seth"
Apr 2019
2·181 Posts 
The "Actual" numbers on https://www.mersenne.ca/manyfactors.php don't always match the numbers under the exponent, they seem to be updated several days after the exponent page updates but not always.
https://www.mersenne.ca/manyfactors...._max=202528981 shows actual 2^67 But https://www.mersenne.ca/exponent/202528981 shows actual 2^78 with no update in the last year https://www.mersenne.ca/manyfactors...._max=340011179 shows actual 2^72 but was updated 10 days ago.  Another confusion I have is on an exponent page (for example https://www.mersenne.ca/exponent/340011179) what is actual ghzday measuring? It doesn't seem to match the sum of TF results. 
20210817, 12:06  #667  
"James Heinrich"
May 2004
exNorthern Ontario
2×3×7×83 Posts 
Quote:
Quote:
Quote:
But your screenshot confuses me, I'm not sure how you saw (or still seeing?) nearzero effort, it should be about 90 GHd (attached). 

20210827, 02:46  #668 
"Seth"
Apr 2019
2·181 Posts 
I appreciate the last updated timestamp and the fact that the column updates as soon as I visit one of these numbers.
the actual TF value and GHzDays is fixed by loading the page. But there are still hundreds of exponents with the wrong "actual" TF level. Totally understand if you don't have time or energy to fix. The root cause is that I submitted a large number of NF TF results from bitranges I tested in 2019 that had a composite factor at that time manual results didn't mark these ranges at complete (not knowing if I had finished TF the entire range), now the server records these as NF. Often I had done additional TF past that bit level leaving a hole that I recently plugged. I ran this to find all the exponents I found a factor for and bulk query those 5K so the entries get updated. Code:
cat megaresults.txt  grep oP "(?<=M)[09]*(?= )"  sort n u  sed 's#^#https://www.mersenne.ca/exponent/#' > open.txt cat ../open.txt  while read f; do sleep 2; wget "${f}"; done; 
20210827, 10:01  #669 
"James Heinrich"
May 2004
exNorthern Ontario
2×3×7×83 Posts 

20210827, 10:18  #670 
"Seth"
Apr 2019
2×181 Posts 
I lovingly smoothed the server for a bit and I can't find any that are incorrect any more.

20210902, 05:15  #671 
"Seth"
Apr 2019
16A_{16} Posts 
I'm frequently looking at the TWOK (to below two thousand) report on https://www.gpu72.com/twok/
Each row links to James' wonderful TF status page e.g. https://www.mersenne.ca/status/tf/0/3/5/1460 Which in turn link to the list of unfactored exponents e.g. https://www.mersenne.org/report_fact...xp_hi=14619999 I was wondering if it would be possible when looking at a comparison with X days ago to link to the report_factors when I click on the factored column? e.g. https://www.mersenne.org/report_fact...xp_lo=14600000 
Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Small inconsistencies between mersenne.org and mersenne.ca factor databases  GP2  mersenne.ca  44  20160619 19:29 
mersenne.ca (ex mersennearies.sili.net)  LaurV  mersenne.ca  8  20131125 21:01 
GaussianMersenne & EisensteinMersenne primes  siegert81  Math  2  20110919 17:36 
Mersenne Wiki: Improving the mersenne primes web site by FOSS methods  optim  PrimeNet  13  20040709 13:51 