mersenneforum.org  

Go Back   mersenneforum.org > Factoring Projects > Aliquot Sequences

Reply
 
Thread Tools
Old 2021-05-19, 16:05   #67
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

2A716 Posts
Default

Quote:
Originally Posted by EdH View Post
I've provided an interesting example, where in the case of the regina_file only containing up to 1e5, there is only one even sequence that ends with the prime 269. I wonder if there are any others within a larger set.
Thank you for this new program !
I did the same manipulation as you for the prime number 269, with regina_file until 5e6.
Result : 13 even sequences found !
Code:
Data available for sequences 2 through 5000000
Sequence endings - prime: 3876810, cycle: 78879, open: 1044310
Enter sequence (##/a/h/p/p##/q/u): 45987
45987 terminates with prime 269.
List all sequences that terminate with 269? (y/n/c/f): c
3684 sequences found.
Perform Advanced Filtering on these results? (y/n): y
Advanced Routine using previous results
Sequence filter a (mod b) (a b): 2
Sequence range (begin end): 
All, Prime, Open, Cylce (a/p/o/c): 
Range for number of iterations (min max): 
Range for number of relative minimums (min max): 
Range for number of relative maximums (min max): 
Range for number of parity changes (min max): 
Range for number of consecutive even abundant terms (min max): 
Range for number of consecutive even deficient terms (min max): 
Range for number of consecutive odd deficient terms (min max): 
Range for number of consecutive odd abundant terms (min max): 
Range for number of downdriver extracts (min max): 
Range for smallest quotient of two consecutive terms (min max): 
Range for largest quotient of two consecutive terms (min max): 
Range for arith mean of all quotients of 2 consequtive terms (min max): 
Range for geo mean of all quotients of 2 consequtive terms (min max): 
Range for arith mean of # of digits of all the mins (min max): 
Range for arith mean of # of digits of all the maxes (min max): 
List filtered sequences for 0 (mod 2) from 2 through 5000000? (y/n/c/f): y
2888      436140    571804    628892    745148    887364    1143602   1216800   1252450   2876416   3476924   3511876   3969836   
13 sequences found.
Enter sequence (##/a/h/p/p##/q/u):
garambois is offline   Reply With Quote
Old 2021-05-19, 17:33   #68
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

76458 Posts
Default

Thanks! I find it interesting that such a large gap exists between the first and second sequence in the final list.

I also see I've still got some spelling errors that I thought I had fixed. Let me know of any you may find.

More later. . .
EdH is offline   Reply With Quote
Old 2021-08-09, 08:58   #69
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

7×97 Posts
Default

After several weeks of uninterrupted calculation, I have the pleasure to put online the brand new "regina" file with all the sequences up to 18 millions !
This file has a size of 449.5 MB compressed and 2 GB uncompressed.
Edwin, your program works very well !
Here is exactly the same entry as in post #67, but with the new file :


Code:
Data available for sequences 2 through 18000000
Sequence endings - prime: 13623702, cycle: 260439, open: 4115858
Enter sequence (##/a/h/p/p##/q/u): 45987
45987 terminates with prime 269.
List all sequences that terminate with 269? (y/n/c/f): c
12693 sequences found.
Perform Advanced Filtering on these results? (y/n): y
Advanced Routine using previous results
Sequence filter a (mod b) (a b): 2
Sequence range (begin end): 
All, Prime, Open, Cylce (a/p/o/c): 
Range for number of iterations (min max): 
Range for number of relative minimums (min max): 
Range for number of relative maximums (min max): 
Range for number of parity changes (min max): 
Range for number of consecutive even abundant terms (min max): 
Range for number of consecutive even deficient terms (min max): 
Range for number of consecutive odd deficient terms (min max): 
Range for number of consecutive odd abundant terms (min max): 
Range for number of downdriver extracts (min max): 
Range for smallest quotient of two consecutive terms (min max): 
Range for largest quotient of two consecutive terms (min max): 
Range for arith mean of all quotients of 2 consequtive terms (min max): 
Range for geo mean of all quotients of 2 consequtive terms (min max): 
Range for arith mean of # of digits of all the mins (min max): 
Range for arith mean of # of digits of all the maxes (min max): 
List filtered sequences for 0 (mod 2) from 2 through 18000000? (y/n/c/f): y
2888      436140    571804    628892    745148    887364    1143602   1216800   1252450   2876416   3476924   3511876   3969836   7195662   7538382   7967862   8358570   8538174   8794806   8794818   8916102   8969544   9961242   9961254   10276902  11767566  11954418  12396420  12560778  13101980  13616100  14231820  14537380  14806494  15081678  15182268  15323166  15352182  15423458  15478542  17306076  
41 sequences found.
Enter sequence (##/a/h/p/p##/q/u):
41 sequences found VS 13 sequences found !

I let the program run.
It is likely that by the end of the year we will have 30 million sequences.
At the moment, the program is adding one million sequences every week, but gradually it will take 8, 9, then 10 days to add one million sequences.

I will do some data analysis with this big file later.

Please let me know if you notice any malfunction when you download this file and work with it !
garambois is offline   Reply With Quote
Old 2021-08-09, 12:49   #70
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

FA516 Posts
Default

This sounds great, Jean-Luc! However, the current program (seqinfo2) can only accept up to 20M. I will need to rewrite it to support 30M (maybe call it seqinfo3). I don't think simply increasing the elements will work, since I already ran into a barrier prior to 20M. I guess I'd better get working on it.
EdH is offline   Reply With Quote
Old 2021-08-09, 14:51   #71
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

7·97 Posts
Default

OK, thanks Ed !
But I don't think I'll stop at 30 million.
I'm going to let the program run for several months.
And within a year, we may even have 50 million aliquot sequences.
However, after a while, even with this great program, there will be the RAM barrier.
But with this perfectly successful 18 million test, I am still amazed at how little RAM is needed !

In the more distant future, I'm even thinking of modifying the "regina" file by scanning all its Open-End sequences on FactorDB and filling them with the more complete data.
Then we will have a great tool to do statistics on aliquot sequences !!!
garambois is offline   Reply With Quote
Old 2021-08-09, 15:05   #72
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

32·5·89 Posts
Default

Quote:
Originally Posted by garambois View Post
. . .
In the more distant future, I'm even thinking of modifying the "regina" file by scanning all its Open-End sequences on FactorDB and filling them with the more complete data.
Then we will have a great tool to do statistics on aliquot sequences !!!
This may match with my intent for the "U"pdate function of my program. I intended to have it gather updated info for any sequences that may have changed. However, I don't have all the details for calculating all the elements. I was expecting to simply have a separate file of all the changes that the program would read in at runtime, without altering the regina_file.
EdH is offline   Reply With Quote
Old 2021-08-09, 18:57   #73
garambois
 
garambois's Avatar
 
"Garambois Jean-Luc"
Oct 2011
France

2A716 Posts
Default

Quote:
Originally Posted by EdH View Post
This may match with my intent for the "U"pdate function of my program. I intended to have it gather updated info for any sequences that may have changed. However, I don't have all the details for calculating all the elements. I was expecting to simply have a separate file of all the changes that the program would read in at runtime, without altering the regina_file.
I seem to recall that your "u" function was not yet available, is that correct ?
Let me know if you manage to do any work in this direction, I'm very interested !
garambois is offline   Reply With Quote
Old 2021-08-09, 19:29   #74
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

400510 Posts
Default

Quote:
Originally Posted by garambois View Post
I seem to recall that your "u" function was not yet available, is that correct ?
Let me know if you manage to do any work in this direction, I'm very interested !
I haven't been working on that yet, because a true update file would need all of the calculations for the updated sequences. I was thinking of the update file at the time when I was only retrieving a few things like last prime, merges, etc. Now that we are looking at all the means, averages, etc., I will need all those formulas to correctly provide a useful update file and then, it will be a heavy task to implement it. Maybe later. . .
EdH is offline   Reply With Quote
Old 2021-09-01, 23:19   #75
Aillas
 
Aillas's Avatar
 
Oct 2002
France

22·3·13 Posts
Default

Hi Ed,

I've start looking at your seqInfo program. I've planned to release 2 updates (at least):

- First release: Remove the 20M limitation. This is almost done. The only limitation will be your computer memory.Plus some minor changes/optimizations.

- The second release could be a **huge** speed optimization. I plan to use multi threading and thread pool.
- Thread pool will be used to download sequences from FactorDB. By default 4 concurrent downloads could be done. This is the limitation I found on FactorDB with my account. If you have special account with less limitation, you could change this value.
- MultiThreading will be used when parsing the regina file. So if you have a Xeon Platinium 56 cores MT 112 HT, it could speed up a little the analysis :). Instead of reading the 20M array, line by line, we could split the job in a specified amount of threads.


Before releasing the second part, I will first experiment the thread pool / MT on AliMerge (when running alimerge 41 1 60 for ex, I will use the TP to download all the factordb sequences. No need to MT the computation as it's really fast)
If it works as expected, I will try it on seqinfo.


One more thing: Sharing source code as attachment in message is not really easy and difficult to found. So I've created a Github repository to manage both programs. It will be easier to track issues / features and maybe store binaries too. I will provide the links by the end of the week.
Aillas is offline   Reply With Quote
Old 2021-09-02, 02:28   #76
EdH
 
EdH's Avatar
 
"Ed Hall"
Dec 2009
Adirondack Mtns

32×5×89 Posts
Default

This all sounds great, Aillas!

I look forward to seeing what you come up with.I, too, am under the four simultaneous processes limitation with the db. I was looking at raising the limit from 20M as well, but haven't done anything with the program recently.

Thanks!
EdH is offline   Reply With Quote
Old 2021-09-03, 15:16   #77
Aillas
 
Aillas's Avatar
 
Oct 2002
France

22×3×13 Posts
Default

Hi,

I've now made the repositories public on my github

There are still plenty works to do, updating the README with more sections (Intro, how to build, how to use, ...)
I will upload later VC19 projects...

First try to share development...

PS: The AliSeq version on github has the fix for bad counter display when downloading the C80 file.
Aillas is offline   Reply With Quote
Reply

Thread Tools


Similar Threads
Thread Thread Starter Forum Replies Last Post
Aliquot sequence reservations schickel Aliquot Sequences 3608 2021-10-21 07:03
Extending an aliquot sequence backwards arbooker Aliquot Sequences 11 2021-09-10 05:28
A new tool to identify aliquot sequence margins and acquisitions garambois Aliquot Sequences 24 2021-02-25 23:31
Another Aliquot Sequence site schickel Aliquot Sequences 67 2012-01-20 17:53
Useful aliquot-sequence links 10metreh Aliquot Sequences 2 2009-07-31 17:43

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


Fri Oct 22 12:11:27 UTC 2021 up 91 days, 6:40, 1 user, load averages: 2.33, 1.65, 1.43

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.