![]() |
[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? |
[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. |
[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] |
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. |
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! |
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.
|
[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. |
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. |
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=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] |
[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.