mersenneforum.org - Factoring for a publication

mersenneforum.org (https://www.mersenneforum.org/index.php)

- Factoring (https://www.mersenneforum.org/forumdisplay.php?f=19)

- - Factoring for a publication (https://www.mersenneforum.org/showthread.php?t=26875)

[QUOTE=Max0526;580802]This is what we call a nice split! c156=2*p78![/QUOTE]
Brilliant even!

[QUOTE=EdH;580810]I will have to play in that arena later today.

For now, all the smaller composites are factored and the c312 is nearing t50.[/QUOTE]My Bash script failed (miserably) to spin the tested SNFS polynomial on one of my machines. I'm trying the same machine on the latest GNFS requested poly (in the msieve thread) to see if that returns anything at least equal to the initial one.

As to the c312, it is at t50. Moving to the c326. . .

i'm jumping around but
(10,9)c204= 184374369527696042002774666493*P175

line 88 (5, -8) c209/snfs polys

[QUOTE=Max0526;580802]This is what we call a nice split! c156=2*p78! Done and done!
Thank you, fivemack!
Could you please send me that c156 GNFS poly? How did you select it? msieve GPU or CADO?
Could you also take a look at [url]https://mersenneforum.org/showpost.php?p=580745&postcount=228?[/url]

[B]EDIT:[/B]
fivemack, I hope you haven't started line 88 (5, -8) yet. If you haven't, could you please hold on for a bit, just to give me a chance to look at all available SNFS polys and their spins for that c209? It would be a perfect training ground for the definitely possible (just checked) polys 5-8 and the newly baked spinner.[/QUOTE]
Done! The spinner helped a lot.
The three best spun polys are below, all above 9e-13 (22% better). :stirpot:
All are in the format c4,c3,c2,c1,c0,Y1,Y0. All need test-sieving of course.
A full list is attached in a text file.
[code]
16, 24, 17, -9, 1, 9307684042037302713836537085227216320008466867809161, -3786905681649011959990757640861865401828382570683684
0.40317 9.62467489e-13

4, 12, 17, -18, 4, 9307684042037302713836537085227216320008466867809161, -7573811363298023919981515281723730803656765141367368
0.80637 9.55967455e-13

1, 6, 17, -36, 16, 9307684042037302713836537085227216320008466867809161, -15147622726596047839963030563447461607313530282734736
1.61274 9.09034083e-13
[/code]
The other curve is too weak (snfs214). Two best from the original 1-4 set:
[code]
(5, -8) c209 / snfs211 --> poly 1
n: 16565142879290504125073762823767688357919065144905054524206541804478927498295871968740614207813511789416877059764381461295498833321565662655267058786104869169924987607710560565878863482569594348028360005558029
# a = 41975122822309585932179048083335413789237187395320622/5839938684558745126126493478220245287305063414925575
Y0: -41975122822309585932179048083335413789237187395320622
Y1: 5839938684558745126126493478220245287305063414925575
# poly x^4 - 6*x^3 + 17*x^2 - 84*x + 196
c4: 1
c3: -6
c2: 17
c1: -84
c0: 196
skew: 4.63621
# E = 7.86184221e-13 <-- best in 1-4
---------------------------------------------
(5, -8) c209 / snfs211 --> poly 4
n: 16565142879290504125073762823767688357919065144905054524206541804478927498295871968740614207813511789416877059764381461295498833321565662655267058786104869169924987607710560565878863482569594348028360005558029
# a = 43070674852707955981472641819129110567338930886162219/9307684042037302713836537085227216320008466867809161
Y0: -43070674852707955981472641819129110567338930886162219
Y1: 9307684042037302713836537085227216320008466867809161
# poly x^4 - 6*x^3 + 17*x^2 - 84*x + 196
c4: 1
c3: -6
c2: 17
c1: -84
c0: 196
skew: 4.03775
# E = 7.79756343e-13 <-- second best in 1-4
[/code]

c146 from line 138 factored

[CODE]
Sat Jun 12 19:06:21 2021 p66 factor: 665143820458149897367011706930190969209165571709435291862102921201
Sat Jun 12 19:06:21 2021 p80 factor: 44893098614779250493405940343123948769333444856278197198665021119423641761101073
[/CODE]

Murphy e-score = 1.051e-11, 41.6M raw relations processed into a 3.3M matrix.

[QUOTE=Max0526;580762]I built the SNFS spinner and it works! :stirpot:
The two polys above easily spin up 9%, and above 1e-16.
[code]
1,0,2,-12,10,64721779846090870584687710177775610620232746158905507772251959674640,-635923536702326787817298118516458633781564513999305084724198320018881
2.80011 1.01256249e-16

2,2,1,4,8,126098023634219809899835084852688502789500914867461314423771553552960,-601881864613820064535786766952632109810808546510423702878545714894641
2.23217 1.01166588e-16

1,0,8,-96,160,32360889923045435292343855088887805310116373079452753886125979837320,-635923536702326787817298118516458633781564513999305084724198320018881
5.58749 1.00195264e-16
[/code]
If anybody wants to try spinning a poly, please let me know.

[B]@swishzzz[/B]
Is it possible to try with Python please?
1) feed in the text file with 7 comma separated values for c4,c3,c2,c1,c0,Y1,Y0 on each line, the test file is attached;
2) connect to [url]http://myfactorcollection.cownoise.com:8090/calculators.html[/url], the top form Optimal Skew;
3) send records through Optimal Skew form and receive the skew and E score for each;
4) pick the line with the best E score;
5) how long would it take to process 100 records?
[/QUOTE]

Can try with this script (change .txt to .py and make sure you install the "requests" package). Less than 1 min for 100 records.

cownoise skew/E Python script

[QUOTE=swishzzz;580836]Can try with this script (change .txt to .py and make sure you install the "requests" package). Less than 1 min for 100 records.[/QUOTE]
Love it!
It's so easily scalable to c5 and c6!
[code]
# modify the following two values
fields = "c4,c3,c2,c1,c0,y1,y0"
entries = """
...
[/code]
It will be a great asset for everybody who does SNFS!
I hope you don't mind us advertising it a bit on your behalf.
Thank you so much, swishzzz!

The sieving goes on... 21,608,085/22,162,816

I was hoping the estimation was good, because I got a similar value of 17,000,000 with the function that yafu uses.

Well, as long as it won't take longer than a comparable GNFS run, I'm fine. :D The rels/q are already down to 1100 rels / 10000 q though.

[QUOTE=EdH;580824]My Bash script failed (miserably) to spin the tested SNFS polynomial on one of my machines. I'm trying the same machine on the latest GNFS requested poly (in the msieve thread) to see if that returns anything at least equal to the initial one.

As to the c312, it is at t50. Moving to the c326. . .[/QUOTE]My GNFS spin appears to be running properly on the number from another thread, so I guess I need to approach SNFS spinning in a different manner.

As for line 161, the c326 shed a small prime to become a c298 which has had t50 ECM. I'm curently ECMing the c328 a bit. . .

[QUOTE=EdH;580874]My GNFS spin appears to be running properly on the number from another thread, so I guess I need to approach SNFS spinning in a different manner.

As for line 161, the c326 shed a small prime to become a c298 which has had t50 ECM. I'm curently ECMing the c328 a bit. . .[/QUOTE]
EdH, please read through the explanation of what happens in CADO with SNFS polys:
[QUOTE=charybdis]
CADO's sizeopt and rootopt are designed for GNFS, so it makes sense that they would produce odd results when fed SNFS polynomials. SNFS polys have amazing size properties; this is why SNFS is so much faster than GNFS. The size properties of SNFS polynomials are so good that we don't usually care about the root properties at all (it's common to see positive alpha values for SNFS, indicating the algebraic norms are *less* likely to be smooth than random numbers of the same size!), and any attempt to run rootopt will likely ruin the size properties. CADO won't realise that the resulting polynomial is bad, because it doesn't know that it was given an SNFS polynomial!
[/QUOTE]