Okay - this is basically the main time-consuming step.
At this point "rep" comes in as a list of 2 dimensional determinant 1
matrices whose entries are integers from the set {0,1,2,...,(q-1)}
And "Edge" is a list of integer tuples where each tuple is of the form
(a,b) with 0<= a,b <= (2n-1)
P = 0
from itertools import product
from itertools import izip
for X in product(rep,repeat = len (Edge)):
k = izip(Edge,X)
M = [[matrix(q^2, q^2, 0)]*(2*n) for _ in range(2*n)]
for ((a,b),Y) in k:
M [a][b] = Y
M [b][a] = Y.inverse()
Z = block_matrix(M)
P = P + Z.charpoly(x)
On Wednesday, May 13, 2015 at 7:29:47 PM UTC-5, William wrote:
>
> On Wed, May 13, 2015 at 5:20 PM, Phoenix <[email protected]
> <javascript:>> wrote:
> >
> >
> > So I have this SAGE code
>
> Post it. Please make every question you ask here as concrete and
> explicit as possible. At a bare minimum post enough so that the
> reader can reproduce exactly your problem with minimal effort.
> Thanks,
>
> William
>
> > which takes in two integers p and n and it
> > generates ~ (p^3)^(n^2) 0/1 matrices and does a sum of their
> characteristic
> > polynomials.
> >
> > I got answers for (p = 3, n = 2), (p =2 , n = 3) and (p=2, n=2)
> >
> > For any higher number the code runs for ~1hr and then it says,
> >
> > "
> >
> > Traceback (most recent call last):
> > File "<stdin>", line 1, in <module>
> > File "_sage_input_2.py", line 10, in <module>
> > exec compile(u'open("___code___.py","w").write("# -*- coding: utf-8
> > -*-\\n" +
> >
> _support_.preparse_worksheet_cell(base64.b64decode("UA=="),globals())+"\\n");
>
> > execfile(os.path.abspath("___code___.py"))
> > File "", line 1, in <module>
> >
> > File "/tmp/tmpPiPIoq/___code___.py", line 2, in <module>
> > exec compile(u'P
> > File "", line 1, in <module>
> >
> > "
> >
> >
> > Can someone help understand what is going on?
> >
> > Like any suggestion about how to go about it?
> >
> >
> > --
> > You received this message because you are subscribed to the Google
> Groups
> > "sage-support" group.
> > To unsubscribe from this group and stop receiving emails from it, send
> an
> > email to [email protected] <javascript:>.
> > To post to this group, send email to [email protected]
> <javascript:>.
> > Visit this group at http://groups.google.com/group/sage-support.
> > For more options, visit https://groups.google.com/d/optout.
>
>
>
> --
> William (http://wstein.org)
>
--
You received this message because you are subscribed to the Google Groups
"sage-support" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at http://groups.google.com/group/sage-support.
For more options, visit https://groups.google.com/d/optout.
