mersenneforum.org A Tool for Studying Aliquot Sequence Data
 Register FAQ Search Today's Posts Mark Forums Read

2021-05-13, 14:31   #56
EdH

"Ed Hall"
Dec 2009

3×5×251 Posts

Quote:
 Originally Posted by garambois Thanks Edwin for all your hard work ! I have some questions, remarks, answers, information... 1) Would it be possible to show us an example of an entry to achieve an "Advanced Search" ?
Code:
[n%2==0, a==0, b, c, d, e==0, f, g, h, i, j, 1.7<k<2.3, l, m]
and translating for the change of elements:
Code:
$./seqinfo2 regina_file_1e5 Data available for sequences 2 through 100000 Sequence endings - prime: 85544, cycle: 2225, open: 12230 Enter sequence (##/a/h/p/p##/q/u): a Advanced routine Sequence filter a (mod b) (a b): 2 Sequence range (begin end): All, Prime, Open, Cylce (a/p/o/c): o Range for number of relative minimums (min max): Range for number of relative maximums (min max): 0 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): 1.7 2.3 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 100000? (y/n/c/f): y . . . 17360 17544 17994 18006 18018 18138 18150 18490 18702 18786 18798 18888 19040 19410 19500 19560 19656 20204 20382 20640 21096 21198 21210 21234 21246 21258 21384 21582 21858 21870 22132 22392 . . . 842 sequences found. Enter sequence (##/a/h/p/p##/q/u): Of course, this was only with a regina_file that was limited to 1e5 and I showed only 32 of 842 sequences, to include 19560 which you referenced earlier. I entered "a" to call up the Advanced Routine. Then "2" for 0 (mod 2), "o" for open sequences, "0" for maximums, which became the range 0-0, and "1.7 2.3" for the arithmetic mean range. In this case I was inclusive of the values, whereas the original wasn't. I suppose "1.700001 2.29999" could have been used. Quote:  Originally Posted by garambois 2) After less than 3 weeks of computation, the new program came to more than 4,000,000. So I think it should reach 20 M sometime this summer. But then, the RAM needed to run the program should become too large for almost all computers, unless I am wrong ! For sequences up to 20 M, the file will be over 2 GB. And in python, in my experience, it is impossible to process such a file by storing it in RAM : it takes much more memory than the file size. But maybe it's different in C, I hope so ! 3) Just for information, I am thinking of letting the program run for several months and this new program could process all the sequences up to 50 M, or even more... if my old computer does not break down. I've placed a limit of 20M in the code at present, but I have no idea what the memory usage would actually be. Or, what happens when there's not enough RAM. Theoretically, I think it should move into swap space, but then what? Quote:  Originally Posted by garambois 4) There can be no other negative value than B in the table, that is for sure. Great! Thanks.  2021-05-14, 08:32 #57 garambois "Garambois Jean-Luc" Oct 2011 France 2×5×61 Posts Edwin, this new regina file analysis tool is very powerful ! Thank you very much for all this work ! I still have some suggestions to improve it a bit, if in the next times you work on this program again : 1) In the advanced search, would it be possible to add a query "Range for number of iterations (min max) :" ? 2) Unless I'm mistaken, there are only two types of searches I can't do with this wonderful program (The equivalent of ##p, but for open end sequences and for sequences that end in cycles) : I) Search for the number or list of sequences that merge with, for example, 276 or another open end sequence. II) Search for the number or list of sequences that arrive on a perfect number, and on another cycle. This is much more difficult, even impossible, because to enter for example a cycle of length 4, there are 4 possible entry doors and there are many different cycles 4 ! Maybe I will find other things to add to the program when I start to analyze the data very seriously in a few months, when the regina file will have 20 M sequences or more ! But for the moment, this is all I see... Congratulations, once again for the ease of use of this program and for the possibilities it offers ! 2021-05-14, 13:37 #58 EdH "Ed Hall" Dec 2009 Adirondack Mtns 376510 Posts Quote:  Originally Posted by garambois Edwin, this new regina file analysis tool is very powerful ! Thank you very much for all this work ! I still have some suggestions to improve it a bit, if in the next times you work on this program again : 1) In the advanced search, would it be possible to add a query "Range for number of iterations (min max) :" ? 2) Unless I'm mistaken, there are only two types of searches I can't do with this wonderful program (The equivalent of ##p, but for open end sequences and for sequences that end in cycles) : I) Search for the number or list of sequences that merge with, for example, 276 or another open end sequence. II) Search for the number or list of sequences that arrive on a perfect number, and on another cycle. This is much more difficult, even impossible, because to enter for example a cycle of length 4, there are 4 possible entry doors and there are many different cycles 4 ! Maybe I will find other things to add to the program when I start to analyze the data very seriously in a few months, when the regina file will have 20 M sequences or more ! But for the moment, this is all I see... Congratulations, once again for the ease of use of this program and for the possibilities it offers ! Those that merge with other sequences are shown in the basic search, but I'm guessing you would like to also filter those (and the others) through the Advanced Function?: Code: Enter sequence (##/a/h/p/p##/q/u): 70508 70508 is open ended. It merges with 276. Enter sequence (##/a/h/p/p##/q/u): 276 276 is open ended. List any sequences that merge with 276? (y/n/c/f): y . . . 77764 78718 79678 82604 82810 84094 85772 86158 87566 88210 90010 90772 91684 92158 93798 93932 94718 96098 96878 96912 97534 98008 99058 167 sequences found. Enter sequence (##/a/h/p/p##/q/u): 3142 3142 ends with the perfect number 496, Display all sequences that end with 496? (y/n/c/f): y 496 608 650 652 790 1294 1574 1778 2162 2582 3142 5158 12 sequences found. To perform cycle searches, you need to have wget and an Internet connection to reach factordb: Code: Enter sequence (##/a/h/p/p##/q/u): 83076 83076 ends in a cycle. Display cycle? (starts at entry point) (y/n/f): y 376736 381028 285778 152990 122410 97946 48976 45946 22976 22744 19916 17716 14316 19116 31704 47616 83328 177792 295488 629072 589786 294896 358336 418904 366556 274924 275444 243760 Display all sequences that end in this cycle? (y/n/c/f): y . . . 66312 67078 67128 67894 69542 70846 71612 75536 77506 77576 81898 82768 83076 83114 83328 84000 85438 88692 91886 94214 97946 85 sequences found. But, again, this does not include the Advance filtering, "yet." Last fiddled with by EdH on 2021-05-14 at 13:42 Reason: fixt spiiln errrer 2021-05-14, 16:07 #59 garambois "Garambois Jean-Luc" Oct 2011 France 2·5·61 Posts Quote:  Originally Posted by EdH but I'm guessing you would like to also filter those (and the others) through the Advanced Function ? YES ! Yes, but no problem if it's too complicated, especially when it comes to sequences that arrive into a cycle ! The possibilities of the program are really already fabulous now... And as said in my previous message, I think it would be very useful to add also the filter on the number of iterations of the sequence "Range for number of iterations (min max) :". 2021-05-14, 16:20 #60 EdH "Ed Hall" Dec 2009 Adirondack Mtns 3×5×251 Posts Quote:  Originally Posted by garambois . . . And as said in my previous message, I think it would be very useful to add also the filter on the number of iterations of the sequence "Range for number of iterations (min max) :". How strange. . . I added that and thought I posted it. I must have skipped hitting "Submit." Anyway, the basics were in the earlier source. A little copy/paste/editing and voilĂ . . . Attached Files  seqinfo2.cpp (36.8 KB, 4 views) 2021-05-14, 16:48 #61 garambois "Garambois Jean-Luc" Oct 2011 France 2·5·61 Posts Quote:  Originally Posted by EdH How strange. . . I added that and thought I posted it. I must have skipped hitting "Submit." Anyway, the basics were in the earlier source. A little copy/paste/editing and voilĂ . . . OK, many thanks. I will do some more testing in the next few days with a larger regina file. I'll let you know if I run into any problems, but everything should work fine. 2021-05-14, 17:56 #62 EdH "Ed Hall" Dec 2009 Adirondack Mtns 3×5×251 Posts Quote:  Originally Posted by garambois OK, many thanks. I will do some more testing in the next few days with a larger regina file. I'll let you know if I run into any problems, but everything should work fine. Sounds Good! Could you take a look at how much RAM is being used? If you run System Monitor, the Processes tab should show you the resting use of memory for seqinfo2. For the _1e5 file I have, my system is using 1.5 GB. 2021-05-14, 19:30 #63 garambois "Garambois Jean-Luc" Oct 2011 France 2·5·61 Posts Quote:  Originally Posted by EdH Sounds Good! Could you take a look at how much RAM is being used? If you run System Monitor, the Processes tab should show you the resting use of memory for seqinfo2. For the _1e5 file I have, my system is using 1.5 GB. This is truly amazing ! The program uses 1.5 GB if I use a 0.1M (1e5) regina file and it uses 1.7GB if I use a 3.5M (35e5) regina file ! I don't understand how this is possible ? 2021-05-14, 21:44 #64 EdH "Ed Hall" Dec 2009 Adirondack Mtns 3·5·251 Posts Quote:  Originally Posted by garambois This is truly amazing ! The program uses 1.5 GB if I use a 0.1M (1e5) regina file and it uses 1.7GB if I use a 3.5M (35e5) regina file ! I don't understand how this is possible ? I "think" it's because I'm initializing a huge structure (actually two), which allocates memory and then that memory is filled with the data. So, there is probably 1.5 GB set aside from the outset, and then it enlarges it as needed.  2021-05-15, 10:31 #65 garambois "Garambois Jean-Luc" Oct 2011 France 10011000102 Posts OK, thanks for the explanation. Let's wait and see in a few weeks when we're at 10M. 2021-05-18, 20:13 #66 EdH "Ed Hall" Dec 2009 Adirondack Mtns 72658 Posts I have added the capability to perform an Advanced Search on the results of a basic search. This doesn't cover everything that was mentioned, but it might present a way of doing most of the other things. Basically, the program now asks if you want to perform the Advanced Routine on the previous results. 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. Code: $ ./seqinfo2 regina_file_1e5
Data available for sequences 2 through 100000
Sequence endings - prime: 85544, cycle: 2225, open: 12230
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): y
269       595       1315      2093      2827      2888      4109      4267
6267      6357      6931      8451      9149      9307      9809      10435
12933     13133     13987     14333     15333     15931     16867     16913
17453     17469     17947     18789     19059     19733     20609     21305
22049     22891     25341     25709     28347     28673     29693     29813
32021     33627     36099     37133     38033     38165     39773     39955
40107     41451     42485     45675     45711     45987     47249     47587
47907     49015     49275     50573     51045     52109     52145     52347
54585     55289     57609     59013     59101     59309     59827     60653
61587     62609     63045     64469     65635     66555     66893     68695
69045     69547     70053     70253     71635     72833     72989     73853
76605     77393     78209     78453     79891     80813     81113     82459
83083     83309     85035     85073     86007     86609     87053     87809
87953     88013     89921     91533     92851     96031     97589     98635
112 sequences found.
Perform Advanced Filtering on these results? (y/n): y
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 100000? (y/n/c/f):
2888
1 sequence found.
Enter sequence (##/a/h/p/p##/q/u): q
Attached Files
 seqinfo2.cpp (38.6 KB, 2 views)

Last fiddled with by EdH on 2021-05-18 at 20:18 Reason: Helps to attach attachment. . .

 Similar Threads Thread Thread Starter Forum Replies Last Post schickel Aliquot Sequences 3487 2021-06-16 13:28 garambois Aliquot Sequences 24 2021-02-25 23:31 arbooker Aliquot Sequences 5 2020-11-07 15:58 schickel Aliquot Sequences 67 2012-01-20 17:53 10metreh Aliquot Sequences 2 2009-07-31 17:43

All times are UTC. The time now is 22:51.

Wed Jun 16 22:51:57 UTC 2021 up 19 days, 20:39, 0 users, load averages: 2.45, 2.36, 2.31