20210519, 16:05  #67  
"Garambois JeanLuc"
Oct 2011
France
3^{2}×79 Posts 
Quote:
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): 

20210519, 17:33  #68 
"Ed Hall"
Dec 2009
Adirondack Mtns
2^{5}×131 Posts 
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. . . 
20210809, 08:58  #69 
"Garambois JeanLuc"
Oct 2011
France
3^{2}×79 Posts 
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): 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 ! 
20210809, 12:49  #70 
"Ed Hall"
Dec 2009
Adirondack Mtns
1060_{16} Posts 
This sounds great, JeanLuc! 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.

20210809, 14:51  #71 
"Garambois JeanLuc"
Oct 2011
France
3^{2}·79 Posts 
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 OpenEnd sequences on FactorDB and filling them with the more complete data. Then we will have a great tool to do statistics on aliquot sequences !!! 
20210809, 15:05  #72 
"Ed Hall"
Dec 2009
Adirondack Mtns
2^{5}·131 Posts 
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.

20210809, 18:57  #73  
"Garambois JeanLuc"
Oct 2011
France
3^{2}×79 Posts 
Quote:
Let me know if you manage to do any work in this direction, I'm very interested ! 

20210809, 19:29  #74 
"Ed Hall"
Dec 2009
Adirondack Mtns
10140_{8} Posts 
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. . .

20210901, 23:19  #75 
Oct 2002
France
3·53 Posts 
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. 
20210902, 02:28  #76 
"Ed Hall"
Dec 2009
Adirondack Mtns
1000001100000_{2} Posts 
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! 
20210903, 15:16  #77 
Oct 2002
France
9F_{16} Posts 
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. 
Thread Tools  
Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Aliquot sequence reservations  schickel  Aliquot Sequences  3641  20220121 13:47 
Extending an aliquot sequence backwards  arbooker  Aliquot Sequences  13  20220113 08:04 
Useful aliquotsequence links  10metreh  Aliquot Sequences  4  20211028 22:17 
A new tool to identify aliquot sequence margins and acquisitions  garambois  Aliquot Sequences  24  20210225 23:31 
Another Aliquot Sequence site  schickel  Aliquot Sequences  67  20120120 17:53 