ALL the number fields created are being cached. Since you don't need them again, I would add "cache=False" to the NumberField() constructor and see how that goes.
It is also true that the finite fields created are cached, but taht does not look so serious since you are only constructin GF(p) for p up to 100, I think. John Cremona On Mar 4, 5:02 pm, Ben Linowitz <benjamin.linow...@gmail.com> wrote: > Marshall, > > Here is an indented version of the program. I did not think that I was > creating too many variables, as I intended the program to simply > overwrite their values during each loop, but would gladly delete them > during each iteration of the master loop if someone would let me know > what the proper command is. > > As for the total number of iterations, yes it is very big, but I never > intended the program to run until completion. Rather, I planned on > terminating the program after a few interesting examples were > outputted. The problem seems to be that the program expands its memory > footprint so rapidly that it cannot get through very many iterations. > > for c in range (1,5000): > for i in range (3,100): > if is_prime(i)==True: > c2=1 > c1=1 > while c1<30: > a4=c1*i > c1+=1 > while c2<30: > a7=c2*i > c2+=1 > a1=randint(3,100) > a2=randint(3,100) > a3=randint(3,100) > a5=randint(3,100) > a6=randint(3,100) > a8=randint(3,100) > a9=randint(3,100) > A= MatrixSpace(RationalField(),3)([[a1,a2,a3], > [a4,a5,a6],[a7,a8,a9]]) > f=A.charpoly() > if f.is_irreducible() == True: > K=NumberField(f,'t'); > if K.is_galois()==True: > if f.discriminant()%i != 0: > print > "======================================" > A > f > print "disc = > ",factor(f.discriminant()) > print "prime=",i > L = FiniteField(i,'q'); > P.<w>=L[] > f.factor_mod(i) > print "(2,1) is in prime:",0==a4%i > print "(3,1) is in prime:",0==a7%i > print "(3,2) is in prime:",0==a8%i > > Best, > Ben > > On Mar 4, 10:15 am, Marshall Hampton <hampto...@gmail.com> wrote: > > > If Ben or Justin could post a properly indented version I am willing > > to take a look. Maybe I'm not looking at it correctly, but I'm not > > sure what variables would be stored in large quantity. In the inner > > loop, the value of A is over-ridden, isn't it? What else would stick > > around? It seems possible that it might speed things up to declare > > > R = MatrixSpace(RationalField(),3) > > > outside the loop, and then do A= R([[23,4,17],[2,2,2],[3,4,5]]) inside > > the loop. But it looks like there are about 100*5000*30*30 = 450 > > million total iterations, so its not such a "little" program. > > > -Marshall > > > On Mar 3, 8:27 pm, Ben Linowitz <benjamin.linow...@gmail.com> wrote: > > > > Justin and Alasdair, > > > > Thanks for your help. > > > > I'm sorry that you had to reformat the code. I had it written up in a > > > text file and was just copying and pasting it into the terminal > > > running sage. I suppose that the formatting must have gotten screwed > > > up when I copied from the text file to post here. > > > > Anyway, I'm a complete novice when it comes to this. How do I delete > > > variables so that the program does not hog so much memory when set to > > > run for a long period of time? > > > > Thanks again, > > > Ben > > > > On Mar 3, 9:12 pm, "Justin C. Walker" <jus...@mac.com> wrote: > > > > > On Mar 3, 2010, at 14:09 , Ben Linowitz wrote: > > > > > > Thank you for your response Alex. Here is the code: > > > > > [snip] > > > > > The code isn't usable as-is (no indenting, wrapped comments), but if > > > > I've reformatted it correctly, using 'top' while this is running shows > > > > that the python process is continually expanding its memory footprint. > > > > > If this is the case, Alasdair may be correct, and deleting variables > > > > will help. > > > > > HTH > > > > > Justin > > > > > -- > > > > Justin C. Walker > > > > Curmudgeon at Large > > > > Director > > > > Institute for the Enhancement of the Director's Income > > > > -- > > > > Build a man a fire and he'll be warm > > > > for a night. > > > > Set a man on fire and he'll be warm > > > > for the rest of his life. > > > > > -- > > > > Justin C. Walker, Curmudgeon at Large > > > > Director > > > > Institute for the Enhancement of the Director's income > > > > ----------- > > > > -- > > > > They said it couldn't be done, but sometimes, > > > > it doesn't work out that way. > > > > - Casey Stengel > > > > -- -- To post to this group, send email to sage-support@googlegroups.com To unsubscribe from this group, send email to sage-support+unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-support URL: http://www.sagemath.org
