mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Aliquot Sequences

Reply
 
Thread Tools
Old 2016-10-13, 01:29   #353
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

2·32·197 Posts
Default

Let's see if I learned anything or if I'm way off track...

If i try to figure out how much ecm to do on the c144, (0.31*144=~45). This means I should go for a t45 threshold.

All the little B1 values equate to very tiny percentages, so if I start at:
2390@3e6
3860@11e6
2000@43e6

I get just over 220% of t45, or just about 1/3 of t50.

Have I gotten any of this correct?

Running ecm-toy.py with the 2000@43e6, I am told to run 300 curves at 110e6.

I currently am passing 830@110e6, which, if I have things correct, adds 23% to the t50 value, taking it over 50% of t50.

Does this signal that enough ECM has been done?

Am I correct at all in my calculations?

If the above is correct, should I move into the gnfs work, or would someone else like to run this one?
EdH is offline   Reply With Quote
Old 2016-10-13, 02:56   #354
VBCurtis
 
VBCurtis's Avatar
 
"Curtis"
Feb 2005
Riverside, CA

53×37 Posts
Default

I think your calcs are fine, and half a t50 is about 3*t45, so you're about triple the usual amount of ECM for a number this size. Head to GNFS!

Take it yourself, IMO; you've surely done enough work on the sequence recently to have the honor (?) of advancing it a couple more steps.
VBCurtis is offline   Reply With Quote
Old 2016-10-13, 03:29   #355
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

2·32·197 Posts
Default

Thanks!

I'll have one of my machines run a poly overnight while the rest keep running curves.

I'll look to see if anyone else chimes in by tomorrow morning and run with it, if not.
EdH is offline   Reply With Quote
Old 2016-10-13, 09:04   #356
henryzz
Just call me Henry
 
henryzz's Avatar
 
"David"
Sep 2007
Cambridge (GMT/BST)

5×19×61 Posts
Default

After something like a t35 is done I would then ask the script what it things. Maybe more of the curves should have been at 110e6 for this number.
henryzz is offline   Reply With Quote
Old 2016-10-13, 09:37   #357
pinhodecarlos
 
pinhodecarlos's Avatar
 
"Carlos Pinho"
Oct 2011
Milton Keynes, UK

10010111101012 Posts
Default

Ed, please let me know if you need support from myself for the next LA. Currently I am on a job which will finish in 43 hours and already queued up the next job but between these two I can squeeze a LA to help you out. PM me if you are interested. Take care.
pinhodecarlos is online now   Reply With Quote
Old 2016-10-13, 13:09   #358
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

2×32×197 Posts
Default

Quote:
Originally Posted by henryzz View Post
After something like a t35 is done I would then ask the script what it things. Maybe more of the curves should have been at 110e6 for this number.
My scripts from long ago have a fixed pattern, but I don't know where I got that pattern. I suppose to be most efficient I should make a specific pattern for each job. I have the easy capability to set the curves and values in the main script. I just don't know how to make up the table. Would entering a low value into the ecm-toy script tell me the best amounts for each in any rough manner? I don't have the true data for the ecm-toy gnfs table, but would that matter very much at the digit levels I can run?

Quote:
Originally Posted by pinhodecarlos View Post
Ed, please let me know if you need support from myself for the next LA. Currently I am on a job which will finish in 43 hours and already queued up the next job but between these two I can squeeze a LA to help you out. PM me if you are interested. Take care.
Thanks Carlos, I can do this one, if everyone is patient enough to wait for my ancient hardware.-Ed

Overnight my machines finished 2000@110e6. They were supposed to start 2000@430e6, but my script had a bug that prevented the 430e6 runs from being assigned. My calculations come up to 90% of a t50 now. I have a poly and will be swapping everything over to gnfs. How is sieving range normally determined? Is it by testing or is there a calculated start point? I was guessing for the others but my guess here would be, maybe 15e6.

Thanks for all the "learnin'" everyone.
EdH is offline   Reply With Quote
Old 2016-10-13, 15:01   #359
VBCurtis
 
VBCurtis's Avatar
 
"Curtis"
Feb 2005
Riverside, CA

53·37 Posts
Default

For GNFS jobs, I start sieving at 1/3rd of the factor base bound. The original factmsieve script started at 1/2 the bound, so somewhere in that range is reasonable.

If your tools don't supply a default factor base (factmsieve.py does) it doesn't matter very much (within a factor of two of "best", say) what factor-base bounds you pick; perhaps 30M for the bounds with sieving starting at 12M should work fine. 14e is probly best; if you want to play with a little test-sieving, try using 14e the regular way and also invoking with with "-J 12" in the flag list. The latter may be a bit faster, at the expense of less yield (and thus a larger range of Q to be sieved).
VBCurtis is offline   Reply With Quote
Old 2016-10-13, 15:30   #360
henryzz
Just call me Henry
 
henryzz's Avatar
 
"David"
Sep 2007
Cambridge (GMT/BST)

5·19·61 Posts
Default

Quote:
Originally Posted by EdH View Post
My scripts from long ago have a fixed pattern, but I don't know where I got that pattern. I suppose to be most efficient I should make a specific pattern for each job. I have the easy capability to set the curves and values in the main script. I just don't know how to make up the table. Would entering a low value into the ecm-toy script tell me the best amounts for each in any rough manner? I don't have the true data for the ecm-toy gnfs table, but would that matter very much at the digit levels I can run?

Thanks Carlos, I can do this one, if everyone is patient enough to wait for my ancient hardware.-Ed

Overnight my machines finished 2000@110e6. They were supposed to start 2000@430e6, but my script had a bug that prevented the 430e6 runs from being assigned. My calculations come up to 90% of a t50 now. I have a poly and will be swapping everything over to gnfs. How is sieving range normally determined? Is it by testing or is there a calculated start point? I was guessing for the others but my guess here would be, maybe 15e6.

Thanks for all the "learnin'" everyone.
It is difficult without the gnfs data matching the ecm data. It probably isn't too bad as long as you are within a factor of 2 of the best amount of ecm. These things usually have a fairly flat curve. There was the argument a while back whether it was better to tell aliqueit to do ecm to 0.25 digits or 0.33 digits. It didn't make that much difference.
You can run the script telling it you have run 0 curves (0@3e6 or something like that).
henryzz is offline   Reply With Quote
Old 2016-10-14, 03:08   #361
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

67328 Posts
Default

Quote:
Originally Posted by VBCurtis View Post
For GNFS jobs, I start sieving at 1/3rd of the factor base bound. The original factmsieve script started at 1/2 the bound, so somewhere in that range is reasonable.

If your tools don't supply a default factor base (factmsieve.py does) it doesn't matter very much (within a factor of two of "best", say) what factor-base bounds you pick; perhaps 30M for the bounds with sieving starting at 12M should work fine. 14e is probly best; if you want to play with a little test-sieving, try using 14e the regular way and also invoking with with "-J 12" in the flag list. The latter may be a bit faster, at the expense of less yield (and thus a larger range of Q to be sieved).
Unfortunately, I don't understand factor base. I am doing everything manually although I have some scripts I wrote to distribute sieving and ECMing among my machines. I have used factmsieve.py in the past and even modified a version to work with my cluster I had running then. I should probably resurrect that script...

I was out all day, so I didn't try the "-j 12" switch this time - maybe next.

Quote:
Originally Posted by henryzz View Post
...
You can run the script telling it you have run 0 curves (0@3e6 or something like that).
Excellent! This works great. I will do this and use the returned values to seed my ECM script from now on.

On the good side, I was able to build a matrix with a little over 36M unique relations and I have LA running on my cluster. I should have the factors in the morning:
Code:
Thu Oct 13 22:32:44 2016  Msieve v. 1.53 (SVN 993)
Thu Oct 13 22:32:44 2016  random seeds: 78202805 bdd4a0ec
Thu Oct 13 22:32:44 2016  MPI process 0 of 3
Thu Oct 13 22:32:44 2016  factoring 131954905347588391934699304738827948133275634558635457624548533060467418907362122329509131668522067078677069743086471484822031174847391546400311 (144 digits)
Thu Oct 13 22:32:46 2016  searching for 15-digit factors
Thu Oct 13 22:32:47 2016  commencing number field sieve (144-digit input)
Thu Oct 13 22:32:47 2016  R0: -10853864674370357346890831166
Thu Oct 13 22:32:47 2016  R1: 2070175205206213
Thu Oct 13 22:32:47 2016  A0: 1567149971784611966940117631701000925
Thu Oct 13 22:32:47 2016  A1: 57830699450094690463301740533
Thu Oct 13 22:32:47 2016  A2: -399896828755837210648218
Thu Oct 13 22:32:47 2016  A3: -58194781578947499
Thu Oct 13 22:32:47 2016  A4: -6125301397
Thu Oct 13 22:32:47 2016  A5: 876
Thu Oct 13 22:32:47 2016  skew 5656895.57, size 6.878e-14, alpha -6.666, combined = 1.347e-11 rroots = 3
Thu Oct 13 22:32:47 2016  
Thu Oct 13 22:32:47 2016  commencing linear algebra
Thu Oct 13 22:32:47 2016  initialized process (0,0) of 3 x 1 grid
Thu Oct 13 22:32:48 2016  read 2789376 cycles
Thu Oct 13 22:32:57 2016  cycles contain 8822764 unique relations
Thu Oct 13 22:34:58 2016  read 8822764 relations
Thu Oct 13 22:35:10 2016  using 20 quadratic characters above 4294917295
Thu Oct 13 22:35:51 2016  building initial matrix
Thu Oct 13 22:37:43 2016  memory use: 1208.9 MB
Thu Oct 13 22:37:46 2016  read 2789376 cycles
Thu Oct 13 22:37:46 2016  matrix is 2789288 x 2789376 (849.1 MB) with weight 269812589 (96.73/col)
Thu Oct 13 22:37:46 2016  sparse part has weight 189118354 (67.80/col)
Thu Oct 13 22:38:18 2016  filtering completed in 2 passes
Thu Oct 13 22:38:19 2016  matrix is 2786778 x 2786925 (848.9 MB) with weight 269715636 (96.78/col)
Thu Oct 13 22:38:19 2016  sparse part has weight 189092868 (67.85/col)
Thu Oct 13 22:38:41 2016  matrix starts at (0, 0)
Thu Oct 13 22:38:41 2016  matrix is 929002 x 2786925 (372.8 MB) with weight 144909846 (52.00/col)
Thu Oct 13 22:38:41 2016  sparse part has weight 64287078 (23.07/col)
Thu Oct 13 22:38:41 2016  saving the first 48 matrix rows for later
Thu Oct 13 22:38:42 2016  matrix includes 64 packed rows
Thu Oct 13 22:38:42 2016  matrix is 928954 x 2786925 (346.7 MB) with weight 89762927 (32.21/col)
Thu Oct 13 22:38:42 2016  sparse part has weight 63022146 (22.61/col)
Thu Oct 13 22:38:42 2016  using block size 8192 and superblock size 294912 for processor cache size 3072 kB
Thu Oct 13 22:38:46 2016  commencing Lanczos iteration (2 threads)
Thu Oct 13 22:38:46 2016  memory use: 245.8 MB
Thu Oct 13 22:39:05 2016  linear algebra at 0.1%, ETA 9h 9m
Thu Oct 13 22:39:11 2016  checkpointing every 310000 dimensions
EdH is offline   Reply With Quote
Old 2016-10-14, 13:01   #362
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

67328 Posts
Default

Posted:
Code:
Fri Oct 14 08:43:00 2016  commencing square root phase
Fri Oct 14 08:43:00 2016  reading relations for dependency 1
Fri Oct 14 08:43:01 2016  read 1392884 cycles
Fri Oct 14 08:43:03 2016  cycles contain 4409660 unique relations
Fri Oct 14 08:43:50 2016  read 4409660 relations
Fri Oct 14 08:44:15 2016  multiplying 4409660 relations
Fri Oct 14 08:49:34 2016  multiply complete, coefficients have about 199.19 million bits
Fri Oct 14 08:49:35 2016  initial square root is modulo 14071643
Fri Oct 14 08:56:03 2016  sqrtTime: 783
Fri Oct 14 08:56:03 2016  p65 factor: 10034588955103289779660710855255598470263096964861235086632824637
Fri Oct 14 08:56:03 2016  p80 factor: 13150006037913501232049838541645294463485178998835505659359474149104392799427203
Fri Oct 14 08:56:03 2016  elapsed time 00:13:04
EdH is offline   Reply With Quote
Old 2016-10-14, 13:09   #363
unconnected
 
unconnected's Avatar
 
May 2009
Russia, Moscow

1001111000012 Posts
Default

200 digits!
unconnected is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Reserved for MF - Sequence 4788 schickel Aliquot Sequences 2934 2021-01-07 18:52
Reserved for MF - Sequence 3366 RichD Aliquot Sequences 463 2021-01-02 16:01
Reserved for MF - Sequence 276 kar_bon Aliquot Sequences 127 2020-12-17 10:05
Team Sieve #37: 3408:i1287 RichD Aliquot Sequences 14 2013-08-02 17:02
80M to 64 bits ... but not really reserved petrw1 Lone Mersenne Hunters 82 2010-01-11 01:57

All times are UTC. The time now is 13:39.

Mon Jan 25 13:39:34 UTC 2021 up 53 days, 9:50, 0 users, load averages: 2.42, 2.19, 2.09

Powered by vBulletin® Version 3.8.11
Copyright ©2000 - 2021, 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.