mersenneforum.org Aliquot genealogy project
 User Name Remember Me? Password
 Register FAQ Search Today's Posts Mark Forums Read

2011-12-31, 22:28   #23
science_man_88

"Forget I exist"
Jul 2009
Dumbassville

836910 Posts

Quote:
 Originally Posted by EdH Now, that, I like! You guys are coming up with some good stuff. The clickable nodes looks good as does the familyecho.com, but, yes, it would entail a whole lot of work. Of course, an automated version is what is needed. Here's what I envision - let me know what you guys think: 1. create a program to query the db for all merges - this would be a long term project 2. have the program write the graphviz code - let's try to be less manually intensive** 3. periodically run the program against all OE sequences and update the graphviz file as needed 4. or, have the program provide whatever data is needed to update familyecho or some other choice of tree type **@firejuggler: did you automate the file or manually enter the data?
Quote:
 digraph G { graph [URL="default.html"] a [URL="a.html"] b [URL="b.html"] c [URL="c.html"] a -> b -> c a -> c }
the bold looks easiest to write but Pari's been my only thing for so long, (so basically no useful language) the rest depend on searching and recording it. I think the way we'd detect it by site normally is:

1) check the sequence lines above i0 for sequences to check
2) check those sequences with a last line search
3)check if it's prime or a cycle if so the value is a termination not a merge ?

the problem is I still don't understand what to do from there to write the -> part because the example they give has multiple pathways.

2011-12-31, 23:13   #24
firejuggler

Apr 2010
Over the rainbow

1001101001002 Posts

updated my file, there should be no mistake anymore. (added cycle, A->A, added amicable , A->B->A)
Attached Files
 graph1.txt (5.7 KB, 104 views)

 2012-01-01, 01:11 #25 science_man_88     "Forget I exist" Jul 2009 Dumbassville 8,369 Posts Code: aligen7(w,s)->print("digraph G {");forprime(z=w,s,print(eval(z)""eval(Str(p))"[shape=plaintext];");for(x=1,#ali(z),print(ali(z)[x]"[shape=plaintext];");print(ali(z)[x]"->"Str(z"p"));for(y=1,#ali(ali(z)[x]),print(ali(ali(z)[x])[y]"[shape=plaintext];");print(ali(ali(z)[x])[y]"->"ali(z)[x]))));print("}") is the best I can do without total modification of the aligen script ( and yes I know it's not in order, I guess I could break it apart a bit and do a forprime loop before to get the first few lines in order.) Last fiddled with by science_man_88 on 2012-01-01 at 01:12
2012-01-01, 05:47   #26
firejuggler

Apr 2010
Over the rainbow

9A416 Posts

another update, done manually, again...
added a few more, the WIP are not done yet (still have to look for the preceding term)
It seems that there is a lot (as observed for now) of sequence that end in 43
naming convention :
prime are named "Xp" in plain text
odd are in plain text
even are default
Attached Files
 graph1.txt (14.0 KB, 117 views)

2012-01-01, 08:28   #27
schickel

"Frank <^>"
Dec 2004
CDP Janesville

2×1,049 Posts

Quote:
 Originally Posted by firejuggler another update, done manually, again... added a few more, the WIP are not done yet (still have to look for the preceding term) It seems that there is a lot (as observed for now) of sequence that end in 43 naming convention : prime are named "Xp" in plain text odd are in plain text even are default
Wow! that's starting to get intricate. It's nice how it does draw the cycles. Reading the documentation, it said you had to change some of the internal edges on the cycles, since DOT won't render them, but I guess it does some magic to make it work.

I found a couple of duplicated links in your structure: 25->45, 455->217, 329->55, 185->82, and 377->82. Nothing major, it just results in a double-arrow.

It's taking a noticable time to reder now: 30 secs or so. Also, the resulting picutre as a PNG is 16560x3227 pixels and weighs in at 1.5MB....geez, and the jpg is a whopping 2.1MB.

Guess we'll have to use the same people that print the Mersenne Prime posters to print it out.....

2012-01-01, 14:42   #28
Mini-Geek
Account Deleted

"Tim Sorbera"
Aug 2006
San Antonio, TX USA

10000101010112 Posts

Quote:
 Originally Posted by schickel It's taking a noticable time to reder now: 30 secs or so. Also, the resulting picutre as a PNG is 16560x3227 pixels and weighs in at 1.5MB....geez, and the jpg is a whopping 2.1MB.
I was able to convert it to PNG in about 5 seconds, and to SVG in under a second. The SVG is 481 KB (attached). I downloaded graphviz and used the command "dot -ograph1.png -Tpng graph1.txt" (and with "png" replaced by "svg"). How did you convert it that it took 30 seconds? (I suspect it was different software, because my PNG was 16455x3227 pixels)
Quote:
 Originally Posted by schickel Guess we'll have to use the same people that print the Mersenne Prime posters to print it out.....
At 300 DPI, the image would be about 54.85x10.76 inches (139x27 cm). If printed on normal 8.5x11 inch or A4 paper, it would take about 7 sheets.
Attached Files
 graph1.zip (70.8 KB, 87 views)

Last fiddled with by Mini-Geek on 2012-01-01 at 14:52

2012-01-01, 14:57   #29
science_man_88

"Forget I exist"
Jul 2009
Dumbassville

100000101100012 Posts

Quote:
 Originally Posted by Mini-Geek I was able to convert it to PNG in about 5 seconds, and to SVG in under a second. The SVG is 481 KB (attached). I downloaded graphviz and used the command "dot -ograph1.png -Tpng graph1.txt" (and with "png" replaced by "svg"). How did you convert it that it took 30 seconds? (I suspect it was different software, because my PNG was 16455x3227 pixels) At 300 DPI, the image would be about 54.85x10.76 inches (139x27 cm). If printed on normal 8.5x11 inch or A4 paper, it would take about 7 sheets.
doesn't look very clickable but still a lot of relations. does this count :

Code:
(23:55)>ali(923774)
%2 = [853356555529]

Last fiddled with by science_man_88 on 2012-01-01 at 14:59

2012-01-01, 15:35   #30
science_man_88

"Forget I exist"
Jul 2009
Dumbassville

100000101100012 Posts

Quote:
 Originally Posted by science_man_88 doesn't look very clickable but still a lot of relations. does this count : Code: (23:55)>ali(923774) %2 = [853356555529]
made a new script to check for and only report under 1000000.

Code:
 (n)->a=n-1;b=[];for(x=0,floor(.5*a),if(x==0,if(sigma(a^2)-a^2==n,b=concat(b,a^2)),if(sigma((x)*(a-x))-((x)*(a-x))==n,b=concat(b,(x)*(a-x)))));b=vecsort(b,,8);for(g=1,#b,if(b[g]<1000000,b[g]=b[g],b[g]=0));
b=vecsort(b,,8)

2012-01-01, 15:43   #31
schickel

"Frank <^>"
Dec 2004
CDP Janesville

2·1,049 Posts

Quote:
 Originally Posted by Mini-Geek I was able to convert it to PNG in about 5 seconds, and to SVG in under a second. The SVG is 481 KB (attached). I downloaded graphviz and used the command "dot -ograph1.png -Tpng graph1.txt" (and with "png" replaced by "svg"). How did you convert it that it took 30 seconds? (I suspect it was different software, because my PNG was 16455x3227 pixels) At 300 DPI, the image would be about 54.85x10.76 inches (139x27 cm). If printed on normal 8.5x11 inch or A4 paper, it would take about 7 sheets.
I was using gvedit which gets installed with the Windows package. It was helpful getting started having a GUI frontend. I assume that it's doing some post processing internally.

I hadn't considered SVG, since I'm still on IE8 and it needs a plugin to render SVG files....good to know that a portable format keeps the size down.

As far as printing, I was anticipating that as more gets added, the size is going to start to scale up, but the size info on the current set is encouraging....

 2012-01-01, 18:44 #32 henryzz Just call me Henry     "David" Sep 2007 Cambridge (GMT/BST) 5,743 Posts Nice graph firejuggler. If we can get this automated this would be brilliant. Currently there is a mistake around 341/343.
2012-01-01, 19:05   #33
science_man_88

"Forget I exist"
Jul 2009
Dumbassville

202618 Posts

Quote:
 Originally Posted by henryzz Nice graph firejuggler. If we can get this automated this would be brilliant. Currently there is a mistake around 341/343.
if I read it correct I see multiple errors concerning 341:

Quote:
 ... 1001->341 ...
Code:
(15:02)>sigma(1001)-1001
%22 = 343

Quote:
 ... 341->57 ...
Code:
(15:02)>sigma(341)-341
%23 = 43
Quote:
 341->82
Code:
(15:02)>sigma(341)-341
%23 = 43

 Similar Threads Thread Thread Starter Forum Replies Last Post firejuggler Aliquot Sequences 26 2012-01-19 08:15 schickel Aliquot Sequences 307 2011-10-28 01:29 JohnFullspeed Aliquot Sequences 18 2011-08-20 21:11 schickel Aliquot Sequences 29 2011-08-12 17:45 Andi47 Aliquot Sequences 3 2009-03-08 10:18

All times are UTC. The time now is 15:36.

Tue Nov 24 15:36:59 UTC 2020 up 75 days, 12:47, 4 users, load averages: 2.01, 1.60, 1.66

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.