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.
