mersenneforum.org

mersenneforum.org (https://www.mersenneforum.org/index.php)
-   Msieve (https://www.mersenneforum.org/forumdisplay.php?f=83)
-   -   Polynomial Request Thread (https://www.mersenneforum.org/showthread.php?t=18368)

henryzz 2017-12-12 11:24

[QUOTE=VBCurtis;473795]No, a 202-digit sequence is for shared-forum efforts, or RyanP.
I'm still working on tweaking CADO with GNFS inputs in 120-150 digit range, likely will be for the next month or so.[/QUOTE]

How close are you getting to ggnfs for 120 digits?

VBCurtis 2017-12-12 17:41

[QUOTE=henryzz;473806]How close are you getting to ggnfs for 120 digits?[/QUOTE]

What siever and LP would you use for 120 with GGNFS? I'll run GGNFS with your settings, and then CADO with the best I've found, and report timings from the same input number.

Dubslow 2017-12-12 17:49

[QUOTE=VBCurtis;473858]What siever and LP would you use for 120 with GGNFS? I'll run GGNFS with your settings, and then CADO with the best I've found, and report timings from the same input number.[/QUOTE]

I assume you mean msieve/ggnfs/msieve?

Edit: I'm doing a C119 right now for an aliquot. yafu uses a hardcoded cpu-time limit for polysearch, followed by [URL="https://github.com/MersenneForum/yafu/blob/master/factor/nfs/nfs.c#L955"]this table[/URL] to [URL="https://github.com/MersenneForum/yafu/blob/master/factor/nfs/nfs.c#L1076"]acquire[/URL] sieving settings.

Here's a paste of relevant output:

[code]elapsed time of 1563.2779 seconds exceeds 1462 second deadline; poly select done
searching for best poly...
new best score = 2.752000e-10, new best line = 1
new best score = 2.790000e-10, new best line = 2
new best score = 2.878000e-10, new best line = 5
new best score = 3.057000e-10, new best line = 10
new best score = 3.078000e-10, new best line = 23
new best score = 3.192000e-10, new best line = 43
new best score = 3.258000e-10, new best line = 91
new best score = 3.295000e-10, new best line = 359
new best score = 3.298000e-10, new best line = 452
nfs: creating ggnfs job parameters for input of size 119
best poly:
# norm 2.549784e-11 alpha -6.791000 e 3.298e-10 rroots 5
n: 14990310173172968149396076235298536707594665367836772041450877325008250952617611478911487775887067660890725454134245487
skew: 66311.79
c0: 20051580705519104784539047600
c1: 1430536546294158310640730
c2: -84327603391838242499
c3: -686423322375450
c4: 18811126798
c5: 61404
Y0: -47590727601634155322689
Y1: 798699592409
nfs: commencing algebraic side lattice sieving over range: 2685000 - 2690000
nfs: commencing algebraic side lattice sieving over range: 2660000 - 2665000
nfs: commencing algebraic side lattice sieving over range: 2670000 - 2675000
nfs: commencing algebraic side lattice sieving over range: 2655000 - 2660000
nfs: commencing algebraic side lattice sieving over range: 2665000 - 2670000
nfs: commencing algebraic side lattice sieving over range: 2650000 - 2655000
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2684999.
Warning: lowering FB_bound to 2659999.
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2669999.
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2649999.
nfs: commencing algebraic side lattice sieving over range: 2675000 - 2680000
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2664999.
nfs: commencing algebraic side lattice sieving over range: 2680000 - 2685000
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2674999.
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2679999.
gnfs-lasieve4I13e (with asm64): L1_BITS=15, SVN $Revision$
Warning: lowering FB_bound to 2654999.
FBsize 194818+0 (deg 5), 367899+0 (deg 1)
FBsize 195162+0 (deg 5), 367899+0 (deg 1)
FBsize 193702+0 (deg 5), 367899+0 (deg 1)
FBsize 196172+0 (deg 5), 367899+0 (deg 1)
FBsize 194448+0 (deg 5), 367899+0 (deg 1)
FBsize 195483+0 (deg 5), 367899+0 (deg 1)
FBsize 195837+0 (deg 5), 367899+0 (deg 1)
FBsize 194074+0 (deg 5), 367899+0 (deg 1)
total yield: 33366, q=2675003 (0.00730 sec/rel) ETA 0h00m)
321 Special q, 954 reduction iterations
reports: 53311356->11828017->10708235->2877270->2848476->2497604
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 33366
0/0 mpqs failures, 25319/14051 vain mpqs
milliseconds total: Sieve 114047 Sched 0 medsched 39082
TD 46807 (Init 2109, MPQS 17277) Sieve-Change 43505
TD side 0: init/small/medium/large/search: 1409 3404 2026 3271 4540
sieve: init/small/medium/large/search: 2832 31851 2588 23370 2634
TD side 1: init/small/medium/large/search: 998 3075 2021 3506 2916
sieve: init/small/medium/large/search: 1703 19122 2638 25369 1941
total yield: 34581, q=2690003 (0.00712 sec/rel)
324 Special q, 955 reduction iterations
reports: 54392316->12040493->10912857->2948978->2917727->2557623
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 34581
0/0 mpqs failures, 26155/14386 vain mpqs
milliseconds total: Sieve 115149 Sched 0 medsched 39470
TD 47758 (Init 2147, MPQS 17913) Sieve-Change 43931
TD side 0: init/small/medium/large/search: 1421 3444 2037 3281 4637
sieve: init/small/medium/large/search: 2830 32228 2604 23612 2668
TD side 1: init/small/medium/large/search: 1002 3120 2009 3517 2970
sieve: init/small/medium/large/search: 1722 19303 2657 25588 1936
total yield: 35479, q=2685017 (0.00721 sec/rel) ETA 0h00m)
338 Special q, 1011 reduction iterations
reports: 56590986->12508912->11329741->3064183->3032625->2659144
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 35479
0/0 mpqs failures, 27050/15080 vain mpqs
milliseconds total: Sieve 119501 Sched 0 medsched 40984
TD 49473 (Init 2222, MPQS 18485) Sieve-Change 45766
TD side 0: init/small/medium/large/search: 1474 3563 2127 3418 4798
sieve: init/small/medium/large/search: 2960 33468 2708 24456 2774
TD side 1: init/small/medium/large/search: 1040 3234 2097 3662 3081
sieve: init/small/medium/large/search: 1788 20011 2761 26551 2024
total yield: 36383, q=2670011 (0.00713 sec/rel) ETA 0h00m)
344 Special q, 1015 reduction iterations
reports: 57575451->12722329->11521886->3128889->3098230->2715891
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 36383
0/0 mpqs failures, 27691/15188 vain mpqs
milliseconds total: Sieve 121244 Sched 0 medsched 41642
TD 50250 (Init 2260, MPQS 18759) Sieve-Change 46325
TD side 0: init/small/medium/large/search: 1495 3636 2153 3462 4886
sieve: init/small/medium/large/search: 3005 33852 2738 24874 2819
TD side 1: init/small/medium/large/search: 1061 3291 2124 3711 3136
sieve: init/small/medium/large/search: 1813 20291 2797 27005 2048
total yield: 36864, q=2680003 (0.00719 sec/rel) ETA 0h00m)
354 Special q, 1031 reduction iterations
reports: 59003987->13071180->11841322->3190920->3158277->2769224
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 36864
0/0 mpqs failures, 28047/15714 vain mpqs
milliseconds total: Sieve 123704 Sched 0 medsched 42616
TD 51160 (Init 2311, MPQS 18999) Sieve-Change 47548
TD side 0: init/small/medium/large/search: 1538 3721 2198 3538 4965
sieve: init/small/medium/large/search: 3065 34562 2786 25392 2871
TD side 1: init/small/medium/large/search: 1095 3372 2175 3785 3183
sieve: init/small/medium/large/search: 1848 20681 2853 27540 2106
total yield: 38522, q=2665001 (0.00707 sec/rel) ETA 0h00m)
370 Special q, 1107 reduction iterations
reports: 61561333->13670471->12384952->3307745->3275916->2872052
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 38522
0/0 mpqs failures, 28793/16031 vain mpqs
milliseconds total: Sieve 126782 Sched 0 medsched 43904
TD 52669 (Init 2376, MPQS 19308) Sieve-Change 49003
TD side 0: init/small/medium/large/search: 1585 4068 2242 3600 5093
sieve: init/small/medium/large/search: 3153 35168 2839 26164 2945
TD side 1: init/small/medium/large/search: 1119 3664 2223 3854 3252
sieve: init/small/medium/large/search: 1904 21120 2902 28421 2166
total yield: 38594, q=2655017 (0.00709 sec/rel)
374 Special q, 1098 reduction iterations
reports: 61845210->13754046->12473419->3290171->3259284->2858343
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 38594
0/0 mpqs failures, 28803/16008 vain mpqs
milliseconds total: Sieve 127786 Sched 0 medsched 44292
TD 52219 (Init 2393, MPQS 19262) Sieve-Change 49313
TD side 0: init/small/medium/large/search: 1599 3782 2273 3635 5040
sieve: init/small/medium/large/search: 3161 35510 2863 26338 2959
TD side 1: init/small/medium/large/search: 1144 3431 2242 3893 3239
sieve: init/small/medium/large/search: 1894 21316 2916 28609 2222
total yield: 38992, q=2660003 (0.00701 sec/rel)
372 Special q, 1103 reduction iterations
reports: 61848418->13695251->12404411->3348310->3316274->2906721
Number of relations with k rational and l algebraic primes for (k,l)=:

Total yield: 38992
0/0 mpqs failures, 29233/16348 vain mpqs
milliseconds total: Sieve 127562 Sched 0 medsched 44100
TD 52642 (Init 2390, MPQS 19514) Sieve-Change 49022
TD side 0: init/small/medium/large/search: 1597 3856 2274 3629 5089
sieve: init/small/medium/large/search: 3162 35516 2865 26210 2962
TD side 1: init/small/medium/large/search: 1131 3470 2245 3901 3259
sieve: init/small/medium/large/search: 1909 21365 2924 28473 2177
nfs: found 292781 relations, need at least 10040346 (filtering ETA: 2h 59m), continuing with sieving ...[/code]

How useful these settings are for a C119 is of course debatable, and I don't know how it compares to e.g. fact_msieve.py (or whatever it's called), but it's a well-traveled benchmark at least.

Edit: Oops, sorry, didn't realize the previous paste doesn't include the actual ggnfs params, here they are:

[code]cat nfs.job
n: 14990310173172968149396076235298536707594665367836772041450877325008250952617611478911487775887067660890725454134245487
skew: 66311.79
c0: 20051580705519104784539047600
c1: 1430536546294158310640730
c2: -84327603391838242499
c3: -686423322375450
c4: 18811126798
c5: 61404
Y0: -47590727601634155322689
Y1: 798699592409
rlim: 5300000
alim: 5300000
lpbr: 27
lpba: 27
mfbr: 54
mfba: 54
rlambda: 2.5
alambda: 2.5[/code]

VBCurtis 2017-12-12 19:10

13e/27LP, got it. I use the factmsieve python script myself, but I've hacked it so many times that my settings are likely nowhere near what normal people would use.

I'll pick a C120 off the aliquot "unreserved" list and compare timings as soon as the C146 from 4788 is complete.

richs 2018-01-28 20:13

Greetings,

Would anyone be willing to help me find a poly for the C158 of Aliquot sequence 829332 i3666?

[CODE]15990530642214633642900068032731441505142278352686412039745072782665004664553464803199860068369378806856559885009842001213044218648323149881127056386195916647[/CODE]

I am running a CPU search on this number using YAFU, but would appreciate any additional help since a C158 will take me three weeks or more of sieving on my laptop.

Thanks very much if you can help me!

VBCurtis 2018-01-29 05:16

Both my GPUs had unexpected CUDA updates destroy msieve stability, so I'll give this C158 a CADO run. The current version of CADO computes Murphy-E based on the alim/rlim/siever chosen, so I don't think they translate to msieve poly scores; you'll have to test-sieve mine against yours.

VBCurtis 2018-01-29 05:31

[QUOTE=VBCurtis;473870]13e/27LP, got it. I use the factmsieve python script myself, but I've hacked it so many times that my settings are likely nowhere near what normal people would use.

I'll pick a C120 off the aliquot "unreserved" list and compare timings as soon as the C146 from 4788 is complete.[/QUOTE]

I finally got around to this; turned out CADO has so many parameters to tweak that I did dozens of C100-C115 before applying learned knowledge to c120.
Stock CADO 2.3.0 on my system takes 93500 CPU-seconds to factor RSA-120. After I tweaked parameters, CADO spent 78500 CPU-seconds on the same factorization.
A factmsieve run with sieve threads set to 12 and LA threads set to 6 on my 6-core i7 spent 3500 thread-seconds on poly select, 80 wall-clock minutes sieving, 20 wall-clock minutes post-processing. Comparing wall-clock to CPU stats is iffy, but I timed CADO at a hair under 3 hrs start-to-finish.

CADO poly select is quite nice, and the sievers are fast. Filtering and LA seem to be the problem; CADO needs something like 40% more raw relations to build a matrix than msieve on similar size input/lim's/LP choice.

I believe that stitching together msieve post-processing with CADO poly and sieve programs would be faster than GGNFS. I have not yet explored how to export CADO relations in msieve-readable form, but if such a task were simple, I think I would choose CADO for any GNFS task.

VBCurtis 2018-01-29 14:39

First poly from C158 of Aliquot sequence 829332 i3666:
[code]n: 159905306422146336429000680327314415051422783526864120397450727826650046645534648031998600683693788068565598850$
skew: 5541261.714
c0: 82473567788845624920896417672066836468
c1: 27423148032934194452519197807998
c2: -31312939954267820807544534
c3: -1488984073909057653
c4: 411147040016
c5: 20580
Y0: -3761868079831828394655862647321
Y1: 5935412090906132067467
# MurphyE (Bf=6.00e+07,Bg=3.50e+07,area=8.05e+15) = 2.05e-10
# f(x) = 20580*x^5+411147040016*x^4-1488984073909057653*x^3-31312939954267820807544534*x^2+27423148032934194452519$
# g(x) = 5935412090906132067467*x-3761868079831828394655862647321
[/code]
I split the poly-select run into two pieces, so that I could get two "best" polys; CADO by default only keeps one, and I haven't yet explored how to pick out the second-best. I'll have another one this afternoon for you.

VBCurtis 2018-01-29 19:40

Second poly from C158 of Aliquot sequence 829332 i3666:
[code]n: 159905306422146336429000680327314415051422783526864120397450727826650046645534648031998600683693788068565598850$
skew: 465492.372
c0: 87988987760957895905901763151718162
c1: -1040289470331592214995858814619
c2: -4539563586060685936221575
c3: 5583685685250082642
c4: 3041174965840
c5: 6360000
Y0: -2188067582950374776883171160937
Y1: 11826921105939311134567
# MurphyE (Bf=6.00e+07,Bg=3.50e+07,area=8.05e+15) = 2.13e-10
# f(x) = 6360000*x^5+3041174965840*x^4+5583685685250082642*x^3-4539563586060685936221575*x^2-104028947033159221499$
# g(x) = 11826921105939311134567*x-2188067582950374776883171160937
[/code]
This concludes my search; I spent around 625000 core-seconds on poly select.

Max0526 2018-01-29 23:24

[QUOTE=VBCurtis;478718]First poly from C158 of Aliquot sequence 829332 i3666:
[code]n: 159905306422146336429000680327314415051422783526864120397450727826650046645534648031998600683693788068565598850$
skew: 5541261.714
c0: 82473567788845624920896417672066836468
c1: 27423148032934194452519197807998
c2: -31312939954267820807544534
c3: -1488984073909057653
c4: 411147040016
c5: 20580
Y0: -3761868079831828394655862647321
Y1: 5935412090906132067467
# MurphyE (Bf=6.00e+07,Bg=3.50e+07,area=8.05e+15) = 2.05e-10
# f(x) = 20580*x^5+411147040016*x^4-1488984073909057653*x^3-31312939954267820807544534*x^2+27423148032934194452519$
# g(x) = 5935412090906132067467*x-3761868079831828394655862647321
[/code]
I split the poly-select run into two pieces, so that I could get two "best" polys; CADO by default only keeps one, and I haven't yet explored how to pick out the second-best. I'll have another one this afternoon for you.[/QUOTE]

Msieve E score:
[code]
c0: 82473567788845624920896417672066836468
c1: 27423148032934194452519197807998
c2: -31312939954267820807544534
c3: -1488984073909057653
c4: 411147040016
c5: 20580
Y0: -3761868079831828394655862647321
Y1: 5935412090906132067467
skew: 7872873.95142
E: 2.27160457e-12
[/code]

Max0526 2018-01-29 23:27

[QUOTE=VBCurtis;478736]Second poly from C158 of Aliquot sequence 829332 i3666:
[code]n: 159905306422146336429000680327314415051422783526864120397450727826650046645534648031998600683693788068565598850$
skew: 465492.372
c0: 87988987760957895905901763151718162
c1: -1040289470331592214995858814619
c2: -4539563586060685936221575
c3: 5583685685250082642
c4: 3041174965840
c5: 6360000
Y0: -2188067582950374776883171160937
Y1: 11826921105939311134567
# MurphyE (Bf=6.00e+07,Bg=3.50e+07,area=8.05e+15) = 2.13e-10
# f(x) = 6360000*x^5+3041174965840*x^4+5583685685250082642*x^3-4539563586060685936221575*x^2-104028947033159221499$
# g(x) = 11826921105939311134567*x-2188067582950374776883171160937
[/code]
This concludes my search; I spent around 625000 core-seconds on poly select.[/QUOTE]

Msieve E score:
[code]
c0: 87988987760957895905901763151718162
c1: -1040289470331592214995858814619
c2: -4539563586060685936221575
c3: 5583685685250082642
c4: 3041174965840
c5: 6360000
Y0: -2188067582950374776883171160937
Y1: 11826921105939311134567
skew: 563763.14396
E: 2.34112257e-12
[/code]


All times are UTC. The time now is 23:12.

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