mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Operation Kibibit (https://www.mersenneforum.org/forumdisplay.php?f=97)
-   -   Call for volunteers: RSA896 (https://www.mersenneforum.org/showthread.php?t=17460)

WraithX 2013-02-09 07:38

[QUOTE=frmky;328653]Here's another 9.5 million. I've stopped it for now.
[url]https://www.dropbox.com/s/meef9g9dkq1j0pz/rsa896_4.dat.m.gz[/url][/QUOTE]

I have finished my 2 batches of 5M. I'll go ahead and take this batch. Jason, I'll upload all the results once I finish this batch.

poily 2013-02-09 21:09

Finished completely with batch 7. Best I got was 2.42e-19 of combined e.

firejuggler 2013-02-10 19:31

1 Attachment(s)
here are my 5000 best hit
only ran size optimisation...
will run root opti for the top 100
but from earlier run,
[code]
Sat Feb 09 12:13:30 2013 R0: -6203084730713792276611133070779722382751981
Sat Feb 09 12:13:30 2013 R1: 42260925725431380236248723
Sat Feb 09 12:13:30 2013 A0: 592294109350104075885455758432520983785765035911272504595840
Sat Feb 09 12:13:30 2013 A1: -44228198432626906079450591344684081042731028765034896
Sat Feb 09 12:13:30 2013 A2: -1266780839981054007483109448136912523638841412
Sat Feb 09 12:13:30 2013 A3: 25584881562241130646568639974012624680
Sat Feb 09 12:13:30 2013 A4: -32136804857417066806626443393
Sat Feb 09 12:13:30 2013 A5: -3165936771052409436964
Sat Feb 09 12:13:30 2013 A6: 7232286821760
Sat Feb 09 12:13:30 2013 skew 88161163.72, size 5.554e-020, alpha -12.726, combined = 2.853e-019 rroots = 4
[/code]
is my best

Dubslow 2013-02-10 21:13

Okay, the batch 3 top-5000 root opt is complete, and this was the not-very-good best:
[code]# norm 6.222097e-20 alpha -11.905801 e 2.809e-19 rroots 4
skew: 61769793.67
c0: 145232831029996093163552361863066803348899730697627094912576
c1: 1472498973799415233497204250419859073914693142439276
c2: -309381525988912811338453294676772024749774984
c3: -512903027606228113001250486017991403
c4: 32343323870153789851863490738
c5: 180328757325985585752
c6: 5966748627840
Y0: -6405170523828211176788244501877849165285655
Y1: 159602931489094517610579779[/code]

I'll start rootopt on the top-1000 of batch 6 momentarily, and will upload all results at once.

firejuggler 2013-02-10 22:23

best poly only slightly better
[code]
Sun Feb 10 23:04:13 2013 R0: -6121810884144972706146540517486139870590622
Sun Feb 10 23:04:13 2013 R1: 2435907600823901442520753
Sun Feb 10 23:04:13 2013 A0: -300111002775699634366748860938862248175411656644137966719225
Sun Feb 10 23:04:13 2013 A1: -17621784293264685805864325144713943857874136634393260
Sun Feb 10 23:04:13 2013 A2: 498095951660428193317679316044985385382946564
Sun Feb 10 23:04:13 2013 A3: 24445078002437030133657739422237219418
Sun Feb 10 23:04:13 2013 A4: -206436615132453702777897484267
Sun Feb 10 23:04:13 2013 A5: -2593745547217833521550
Sun Feb 10 23:04:13 2013 A6: 7827849469440
Sun Feb 10 23:04:13 2013 skew 67686502.19, size 5.627e-020, alpha -12.009, combined = 2.887e-019 rroots = 6
[/code]
2.853 to 2.887, 1% better

Dubslow 2013-02-11 14:45

Top 1000 of batch 6 revealed this better-but-not-great best:
[code]Sun Feb 10 21:03:13 2013 polynomial selection complete
Sun Feb 10 21:03:13 2013 R0: -6861489586646307495841370548381377833058191
Sun Feb 10 21:03:13 2013 R1: 3562982110672646277353362747
Sun Feb 10 21:03:13 2013 A0: 1397116357446289568248089240738277216190197418531172590400
Sun Feb 10 21:03:13 2013 A1: 2798522157134386626865301273289555923795160762285992
Sun Feb 10 21:03:13 2013 A2: -131505870094157041369095365530216420938114820
Sun Feb 10 21:03:13 2013 A3: -3167161386190221226627602128211193654
Sun Feb 10 21:03:13 2013 A4: 43173308215309686805724932635
Sun Feb 10 21:03:13 2013 A5: 618865219486539197262
Sun Feb 10 21:03:13 2013 A6: 3948311727360
Sun Feb 10 21:03:13 2013 skew 50359208.73, size 5.943e-20, alpha -11.761, combined = 3.105e-19 rroots = 4
Sun Feb 10 21:03:13 2013 elapsed time 05:18:06[/code]
Upload coming in some hours.

Dubslow 2013-02-11 23:56

[QUOTE=Dubslow;328966]
Upload coming in some hours.[/QUOTE]

I've finally gotten around to [URL="http://dubslow.tk/random/batches3and6.tar.gz"]it[/URL]. That's all results of batches 3 and 6 -- caution, it's a hefty 2.09 GiB.

jasonp 2013-02-14 21:52

[QUOTE=Dubslow;328882]Okay, the batch 3 top-5000 root opt is complete, and this was the not-very-good best:
[code]# norm 6.222097e-20 alpha -11.905801 e 2.809e-19 rroots 4
skew: 61769793.67
c0: 145232831029996093163552361863066803348899730697627094912576
c1: 1472498973799415233497204250419859073914693142439276
c2: -309381525988912811338453294676772024749774984
c3: -512903027606228113001250486017991403
c4: 32343323870153789851863490738
c5: 180328757325985585752
c6: 5966748627840
Y0: -6405170523828211176788244501877849165285655
Y1: 159602931489094517610579779[/code]
[/QUOTE]
Could you grep through your dataset for the size-optimized polynomial that generated this, then rerun stage 2 with

-ncr "stage2_norm=1e100"

to see if the best poly changes? If you take the stage 1 hit of this polynomial, then run the CADO size optimization and the CADO root sieve, you get a polynomial that is the 3rd-best across the entire batch, with an E value of 3.47e-19

I've also seen it go the other way, where Msieve produces a polynomial 20% better than the CADO tools, starting from the same stage 1 hit.

I wish this process was a more exact science :(

Dubslow 2013-02-15 16:53

Okay, I only found one size-opt hit where both leading coeffs were the same. The results aren't very encouraging, though per the last command it's entirely possible I screwed it up.

[code]bill@Gravemind ~/msieve ∰∂ grep "159602931489094517610579779" batch3.ms | grep "5966748627840" > thebigtest.ms
bill@Gravemind ~/msieve ∰∂ cat thebigtest.ms
5966748627840 292915719056610727512 36064024000366164152339517778 -84497976889235989825064880544623387 -4788238179361086188471687577637977225787162 4303872043414327223117808302910530617580272619100 81812578696380352446961680035616654504774205872670962790 159602931489094517610579779 -6405170523326284855312354542737323010776179 -3.25 7.963209e+32
bill@Gravemind ~/msieve ∰∂ ./msieve -s thebigtest -npr "stage2_norm=1e100" -v -v


Msieve v. 1.51 (SVN 839M)
Fri Feb 15 10:57:05 2013
random seeds: 3f19d3de de1d536a
factoring 412023436986659543855531365332575948179811699844327982845455626433876445565248426198098870423161841879261420247188869492560931776375033421130982397485150944909106910269861031862704114880866970564902903653658867433731720813104105190864254793282601391257624033946373269391 (270 digits)
no P-1/P+1/ECM available, skipping
commencing number field sieve (270-digit input)
R0: -6405170523828211176788244501877849165285655
R1: 159602931489094517610579779
A0: 145232831029996093163552361863066803348899730697627094912576
A1: 1472498973799415233497204250419859073914693142439276
A2: -309381525988912811338453294676772024749774984
A3: -512903027606228113001250486017991403
A4: 32343323870153789851863490738
A5: 180328757325985585752
A6: 5966748627840
skew 61769793.67, size 5.287e-20, alpha -11.906, combined = 2.809e-19 rroots = 4
elapsed time 00:00:02
bill@Gravemind ~/msieve ∰∂ ls thebigtest.p
ls: cannot access thebigtest.p: No such file or directory[/code]

jasonp 2013-02-15 17:40

make sure to delete the .fb file from previous runs; poly selection won't happen if a poly was previously selected :)

firejuggler 2013-02-15 18:06

Since I had a slightly better poly than Dubslow, could anyone run the root optimisation on CADO to see if it get any better?
ms line is
[code]
7827849469440 -1688698656210866333070 -412741298353303163926536130117 22088628597291494753673208243029819 2460026875023452236191886387968291153532877 -51597763639642042798524045753868032621520266693 -1903072627680759370872623350274028735027372531045504753 2435907600823901442520753 -6121810884098033236808830895701743366785951 -2.85 2.713411e+033
[/code]


All times are UTC. The time now is 18:06.

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