Peter, The problem still occurs in version 6.2.beta1 (the version currently running in SMC). Unfortunately, I am not able to access a machine with/download to my current machine the most current development branch (I believe its at 6.2.beta5 right now). Sorry I can't be of more help at the moment.
On Monday, March 10, 2014 12:54:08 AM UTC-7, Peter Bruin wrote: > > Hello, > > This looks very much like a bug I fixed recently (Trac ticket #15446). > Could you check if the bug still exists in the latest development version? > > Peter > > > I'm trying to find asymptotic bounds for the number of polynomials of >> degree $n$ and height $\leq N$ with a specific Galois group (e.g. >> $Syl_2(S_n)$, or in the example below $S_n$). I build all monic polynomials >> with the given parameters, then ask pari for the order of their Galois >> groups using the command >> gp.polgalois(p)[1] >> The problem is, when both $n$ and $N$ are both still relatively small, >> some sort of data structure seems to be filling up, thereby causing a >> strange error. It appears that this structure is saving the commands I have >> passed to gp. The problem is cumulative, meaning this data structure is not >> cleared unless I quit my Sage session and restart. >> >> >> Perhaps the problem is best explained via an example. I can run the >> function I constructed for $n=4, 2 \leq N \leq 4$, but when I try to run it >> for $n=4, N=5$, the error appears: >> sage: E4(5) >> >> --------------------------------------------------------------------------- >> TypeError Traceback (most recent call >> last) >> <ipython-input-5-06eb9bccdf87> in <module>() >> ----> 1 E4(Integer(5)) >> >> <string> in E4(N) >> >> /opt/sage/local/lib/python2.7/site-packages/sage/interfaces/interface.py >> in __getitem__(self, n) >> 945 P = self._check_valid() >> 946 if not isinstance(n, tuple): >> --> 947 return P.new('%s[%s]'%(self._name, n)) >> 948 else: >> 949 return P.new('%s[%s]'%(self._name, str(n)[1:-1])) >> >> ............................................ >> >> >> /opt/sage/local/lib/python2.7/site-packages/sage/interfaces/gp.py in set( >> self, var, value) >> 509 out = self.eval(cmd) >> 510 if out.find('***') != -1: >> --> 511 raise TypeError, "Error executing code in >> GP:\nCODE:\n\t%s\nPARI/GP ERROR:\n%s"%(cmd, out) >> 512 >> 513 >> >> TypeError: Error executing code in GP: >> CODE: >> sage[65536]=sage[65535][1]; >> PARI/GP ERROR: >> *** at top-level: sage[65536]=sage[65535][1] >> *** ^--- >> *** _[_]: not a vector. >> The ellipses denote the error being traced through several different >> methods within the file >> >> /opt/sage/local/lib/python2.7/site-packages/sage/interfaces/gp.py >> If you would like to see the entire example session in which I produced >> the error message above, go to: >> http://sharetext.org/GlBk >> >> >> Noting that $65536=2^16$, it seems that this is quite a natural place for >> this data structure (denoted by 'sage[]' in the example) to be filling up. >> Unfortunately, this structure is essentially inaccessible from the Sage >> command line. >> I should also note that I attempted to fix the problem with the command >> gp.eval('allocatemem()') >> but this did not fix the problem. >> >> Any sort of explanation of what this data structure is, how I can clear >> it between calls to pari, or even some sort of black box workaround, would >> be greatly appreciated. >> > -- You received this message because you are subscribed to the Google Groups "sage-devel" group. To unsubscribe from this group and stop receiving emails from it, send an email to sage-devel+unsubscr...@googlegroups.com. To post to this group, send email to sage-devel@googlegroups.com. Visit this group at http://groups.google.com/group/sage-devel. For more options, visit https://groups.google.com/d/optout.