20120518, 06:43  #1 
Romulan Interpreter
Jun 2011
Thailand
2^{2}×7×317 Posts 
yafu trivia
Code:
==== sieving in progress ( 2 threads): 44444 relations needed ==== ==== Press ctrlc to abort and save state ==== 
20120519, 16:35  #2 
"Daniel Jackson"
May 2011
14285714285714285714
607 Posts 
What number did you factor?

20120520, 23:50  #3 
Basketry That Evening!
"Bunslow the Bold"
Jun 2011
40<A<43 89<O<88
3×2,399 Posts 
Hijacking this thread for some "trivia"
Is there a way to get the best poly found so far while poly select is running?

20120521, 02:24  #4 
Tribal Bullet
Oct 2004
DA1_{16} Posts 
Not without bruteforce parsing the results, unfortunately.

20120521, 03:01  #5 
Basketry That Evening!
"Bunslow the Bold"
Jun 2011
40<A<43 89<O<88
3·2,399 Posts 

20120521, 05:28  #6  
Basketry That Evening!
"Bunslow the Bold"
Jun 2011
40<A<43 89<O<88
1C1D_{16} Posts 
Quote:
Code:
bill@Gravemind:~/yafu∰∂ score.py nfs.dat.p  Best: 40024  score: 2.059e12 nfs.dat.1.p  Best: 129  score: 1.78e12 nfs.dat.2.p  Best: 562  score: 1.802e12 nfs.dat.3.p  Best: 383  score: 1.601e12 nfs.dat.4.p  Best: 707  score: 1.644e12 nfs.dat.5.p  Best: 4  score: 1.345e12 nfs.dat.6.p  Best: 121  score: 2.028e12 nfs.dat.7.p  Best: 828  score: 1.73e12 I think the best is the following: # norm 2.824831e15 alpha 8.666018 e 2.059e12 rroots 3 skew: 59896335.71 c0: 895897620029938237657596659356466915669760 c1: 9916866312411647273960145680089008 c2: 982422272849236864202218268 c3: 15047997790788080639 c4: 20444892482 c5: 3672 Y0: 6230395847942608324575436414951 Y1: 74116122267086123 This is the 40024th poly in nfs.dat.p; at ten lines per poly, this should be somewhere around line 400240. (Edit: It was actually at line 400253. Huh.) bill@Gravemind:~/yafu∰∂ Code:
#! /usr/bin/python3 dir = '/home/bill/yafu' # This is the absolute path to your poly select working directory files = ['nfs.dat.p', 'nfs.dat.1.p', 'nfs.dat.2.p', 'nfs.dat.3.p', 'nfs.dat.4.p', 'nfs.dat.5.p', 'nfs.dat.6.p', 'nfs.dat.7.p'] # This is the list of files with polys in the above directory # It's probably possible to automatically generate this, but I have no idea how standard Msieve is # about that sort of thing. # ################################################################################################### print() # Pet peeve import os, re os.chdir(dir) pattern = re.compile(r'(# norm .*? e ([09][.][09]{3}e1[123]).*?(?=#\Z))', flags=re.DOTALL) scores = [] meta = [] for file in files: with open(file, 'r') as data: mxsc = 0 # Max Murphye score found so far i = 0 # Rolling index best = 0 # Index corresponding to mxsc data = data.read() polys = pattern.findall(data) # This returns a list of matchtuples. match[0] gives the # first group, in this case the whole poly (see above regex) and match[1] gives the inner # group, which is the Murphy e score. while i < len(polys): score = float(polys[i][1]) #print(str(score)) # Debugging statement if score > mxsc: best = i mxsc = score i += 1 print(file+'  Best:', str(best), ' score:', str(mxsc)) scores.append(polys[best]) meta.append( (best, file) ) mxsc = 0 i = 0 best = 0 while i < len(scores): score = float(scores[i][1]) #print('Testing score:', scores[i][1]) if score > mxsc: best = i mxsc = score i += 1 print('\nI think the best is the following:') print(scores[best][0]) print('This is the '+str(meta[best][0])+'th poly in '+str(meta[best][1])+'; at ten lines per poly,\ this should be somewhere around line '+str(10*meta[best][0])+'.\n') As mentioned in the comments, how does Msieve organize the polyselect files? Are the polys always printed in the same format? On a mostly unrelated note, this is for a C158; the spreadsheet that fivemack maintains suggests that a score above 2e12 is above average, yet polyselect is only around 3540% through the amount of time it stated at the beginning. Is this a good poly? (I haven't trial sieved; I have no clue how to produce the other values in a .job file, like lp, q, etc. though I could make some decent guesses from that nice spreadsheet.) (In another segue, can someone explain what the various matrix and pmatrix columns of that spreadsheet mean?) Last fiddled with by Dubslow on 20120521 at 05:42 

20120521, 11:35  #7 
Tribal Bullet
Oct 2004
3·1,163 Posts 
Msieve saves every polynomial whose Evalue exceeds a bound, in the order they are found. It would be relatively straightforward to maintain a priority queue of the N best polynomials, and only print out a new one if it displaces another polynomial from the queue. But right now there are only two kinds of jobs: small jobs where the only thing that matters is the best score, and large jobs that require manual intervention no matter what.

20120522, 07:25  #8 
Basketry That Evening!
"Bunslow the Bold"
Jun 2011
40<A<43 89<O<88
3·2,399 Posts 
Hmm... it finished, and YAFU printed this:
Code:
best poly: # norm 1.404811e15 alpha 7.157730 e 1.2# norm 1.639528e15 alpha 7.621592 e 1.427e12 rroots 5 n: 34473140218538108783457767847160803005570403677190199323501471286696391707357183201882067234295472747361321855598229568051445043790754713420994112701235355501 skew: 36163026.47 c0: 742886448093238800788055271017068780151535 c1: 61941015689400432499049921540181313 c2: 2387468922892267418028395733 c3: 162188279674099406689 c4: 2825808162878 c5: 23256 Y0: 4307145920676934752087794328126 Y1: 171378602537761387 Code:
# norm 3.364711e15 alpha 6.617400 e 2.303e12 rroots 3 skew: 5288187.23 c0: 8352362303696618979154069233422527850 c1: 13760873120991388185420545758555 c2: 13016221317246226972433428 c3: 7337464421806796798 c4: 756156744888 c5: 20664 Y0: 4410149585920118231836181983101 Y1: 41380670996658991 Edit: There appear to be multiple places where this clobbering occurred. Code:
bill@Gravemind:~/yafu∰∂ grep "#.*#" nfs.dat.p n 1691148:# norm 1.576873e15 alpha 6.920# norm 1.500226e15 alpha 7.108843 e 1.360e12 rroots 3 1707128:# norm 1.395646e15 al# norm 1.524023e15 alpha 6.490944 e 1.370e12 rroots 5 1708992:# n# norm 1.427331e15 alpha 7.023741 e 1.338e12 rroots 1 1749323:# norm 1.614413e15 alpha 7.686080 e 1.417e1# norm 1.836832e15 alpha 7.081007 e 1.517e12 rroots 3 1749979:# norm 1.404811e15 alpha 7.157730 e 1.2# norm 1.639528e15 alpha 7.621592 e 1.427e12 rroots 5 1750109:# norm 1.605257e15 alpha 6.958131 e 1.386e12 rroots# norm 1.414861e15 alpha 7.297132 e 1.308e12 rroots 5 Edit2: I deleted the parts in red from .dat.p, deleted .fb and .job, and "resumed" the job; it's redoing around 20 minutes worth of PS, which isn't too bad. Hopefully it'll get the poly right this time. Last fiddled with by Dubslow on 20120522 at 07:44 
20120522, 09:49  #9 
Basketry That Evening!
"Bunslow the Bold"
Jun 2011
40<A<43 89<O<88
3·2,399 Posts 
Well the 20 minutes was a total lie; it just finished now. (2 of 8 threads took an hour, the rest have slowly trickled off until the last one died just now.)
Regardless, this time there were no errors in the poly file. Code:
best poly: # norm 3.364711e15 alpha 6.617400 e 2.303e12 rroots 3 n: 34473140218538108783457767847160803005570403677190199323501471286696391707357183201882067234295472747361321855598229568051445043790754713420994112701235355501 skew: 5288187.23 c0: 8352362303696618979154069233422527850 c1: 13760873120991388185420545758555 c2: 13016221317246226972433428 c3: 7337464421806796798 c4: 756156744888 c5: 20664 Y0: 4410149585920118231836181983101 Y1: 41380670996658991 Code:
bill@Gravemind:~/yafu∰∂ cat nfs.job n: 34473140218538108783457767847160803005570403677190199323501471286696391707357183201882067234295472747361321855598229568051445043790754713420994112701235355501 skew: 5288187.23 c0: 8352362303696618979154069233422527850 c1: 13760873120991388185420545758555 c2: 13016221317246226972433428 c3: 7337464421806796798 c4: 756156744888 c5: 20664 Y0: 4410149585920118231836181983101 Y1: 41380670996658991 rlim: 30000000 alim: 30000000 lpbr: 29 lpba: 29 mfbr: 58 mfba: 58 rlambda: 2.500000 alambda: 2.500000 
20120522, 12:50  #10 
"Ben"
Feb 2007
2·13·127 Posts 
I don't know  I've never seen yafu or msieve hose its own output before, but neither have I been looking. The poly file parser in yafu has been working just fine, as far as I know, so that's as far as I've investigated. I would think that yafu's parser would become confused by the "e 1.2#" you showed below, so I tend to think that doesn't happen normally...

20120522, 15:40  #11 
Sep 2008
Kansas
6106_{8} Posts 

Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Running YAFU via Aliqueit doesn't find yafu.ini  EdH  YAFU  8  20180314 17:22 
Random trivia questions  chappy  Puzzles  40  20150312 03:59 
yafu1.33  bsquared  YAFU  12  20121108 04:12 
yafu1.32.1  bsquared  YAFU  21  20120904 19:44 
Sports Trivia  cheesehead  Lounge  10  20101031 01:57 