mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Factoring

Reply
 
Thread Tools
Old 2020-04-15, 10:48   #1
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

2·52·127 Posts
Default Yields with 16e

This is a 16e factorisation of a GNFS-200.0 number (and about a third complete)

I am a little surprised that, with alim=rlim=400M, I'm seeing the yield increasing monotonically as Q goes down all the way from 300M to 200M ... I've generally expected the peak yield to be roughly at Q=lim

Maybe I should start this kind of large job by sieving 1M slots at 20M intervals then filling in the gaps. Much narrower sieve regions will just pick up Poisson noise rather than a 10% yield difference.

Actually the noise is much worse than Poisson, because there's a roughly Poisson process at each prime ideal plus non-uniform distributions of the prime ideals; if I fit a linear trend to the yield and look at \left| \frac{Y-f(X)}{\sqrt{f(X)}} \right| then the median is about seven.
Attached Thumbnails
Click image for larger version

Name:	Screen Shot 2020-04-15 at 11.30.53.png
Views:	44
Size:	28.6 KB
ID:	22032  

Last fiddled with by fivemack on 2020-04-15 at 10:49
fivemack is offline   Reply With Quote
Old 2020-04-19, 13:45   #2
fivemack
(loop (#_fork))
 
fivemack's Avatar
 
Feb 2006
Cambridge, England

2×52×127 Posts
Default

Quote:
Originally Posted by fivemack View Post
This is a 16e factorisation of a GNFS-200.0 number (and about a third complete)

I am a little surprised that, with alim=rlim=400M, I'm seeing the yield increasing monotonically as Q goes down all the way from 300M to 200M ... I've generally expected the peak yield to be roughly at Q=lim

Maybe I should start this kind of large job by sieving 1M slots at 20M intervals then filling in the gaps. Much narrower sieve regions will just pick up Poisson noise rather than a 10% yield difference.

Actually the noise is much worse than Poisson, because there's a roughly Poisson process at each prime ideal plus non-uniform distributions of the prime ideals; if I fit a linear trend to the yield and look at \left| \frac{Y-f(X)}{\sqrt{f(X)}} \right| then the median is about seven.
Dividing out by the counts of prime ideals (which are just length(polrootsmod(f,p)) ) gives something with much nicer statistical properties - the linear fit has an r^2 of 0.998 rather than 0.957, and you immediately see the <1% outliers caused by a segfault in the middle of some of the subintervals. This suggests that 100kQ slots with prime-ideal correction would be fine for yield estimation; in fact, looking at the first subinterval of each block, 10kQ slots at 10M intervals look to produce a pretty close fit.

Last fiddled with by fivemack on 2020-04-19 at 14:08
fivemack is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Zero/low yields in KF siever @~67M (badsched bug) Batalov Factoring 38 2015-02-09 07:57
M(57885161) yields no primitive trinomial over GF(2) ewmayer Computer Science & Computational Number Theory 5 2013-06-06 02:50

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

Mon Jul 13 18:18:08 UTC 2020 up 110 days, 15:51, 1 user, load averages: 2.26, 1.82, 1.70

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

This forum has received and complied with 0 (zero) government requests for information.

Permission is granted to copy, distribute and/or modify this document under the terms of the GNU Free Documentation License, Version 1.2 or any later version published by the Free Software Foundation.
A copy of the license is included in the FAQ.