---------------------------------------------------------------------
SystemError Traceback (most recent call last)
<ipython-input-6-7b242311279b> in <module>()
----> 1 classify_ord_pe(L,Integer(2),Integer(2),"results/order4.txt","w")
/home/simon/.sage/temp/k3/8598/K3_aut_classification.sageqIg9CK.py in
classify_ord_pe(L, p, e, file_name, rw)
252 print(" ")
253 cofix = cofixed[k].twist(-_sage_const_1 )
--> 254 for Aa in k3_prime_power(fix.genus(), p, e):
255 A, a, Oa = Aa
256 actsg = MaximalK3surfaceAut(A, cofix, a, Oa)
/home/simon/.sage/temp/k3/8598/prime_power.sage5OzLUx.py in
k3_prime_power(genus, prime, e)
138 signatures += [[ranks_E[k]]*(weights[k]//_sage_const_2
) for k in range(_sage_const_1 ,n)]
139 signatures[-_sage_const_1 ][_sage_const_0 ] -= _sage_const_1
--> 140 for act in
prime_power_actions(genus,prime,ranks,signatures):
141 yield act
142
/home/simon/.sage/temp/k3/8598/prime_power.sage5OzLUx.py in
prime_power_actions(genus, p, ranks, signatures)
276 # recurse
277 # print(R)
--> 278 for N in prime_power_actions(R, p,
ranks[:-_sage_const_1 ], R_signatures):
279 N, fN, GN = N
280 ext = extensions(M, fM, N, fN, GM, GN,
glue_order, p)
/home/simon/.sage/temp/k3/8598/prime_power.sage5OzLUx.py in
prime_power_actions(genus, p, ranks, signatures)
240 if G.is_even() and not M.is_even():
241 continue
--> 242 GM = Oq_equiv(M, fM, p**e)
243 DM = M.discriminant_group()
244 M_max_glue_group = (M.span((fM**(p**(e-_sage_const_1 ))
- fM**_sage_const_0 ).inverse())
/home/simon/.sage/temp/k3/8598/prime_power.sage5OzLUx.py in Oq_equiv(M, fM,
order)
514 """
515 if order==_sage_const_1 or order==_sage_const_2 :
--> 516 return M.image_in_Oq()
517 elif prod(M.signature_pair()) == _sage_const_0 :
518 #
This looks like a bug to me but I am not enough of an expert to decide.
/home/simon/sage/local/lib/python2.7/site-packages/sage/modules/free_quadratic_module_integer_symmetric.pyc
in image_in_Oq(self)
1443 generated by 5 elements
1444 """
-> 1445 Oq = self.discriminant_group().orthogonal_group()
1446 sig = self.signature_pair()
1447 if sig[0]*sig[1]==0 or self.rank()==2:
/home/simon/sage/local/lib/python2.7/site-packages/sage/misc/cachefunc.pyx
in sage.misc.cachefunc.CachedMethodCaller.__call__
(build/cythonized/sage/misc/cachefunc.c:10232)()
1948 return cache[k]
1949 except KeyError:
-> 1950 w = self._instance_call(*args, **kwds)
1951 cache[k] = w
1952 return w
/home/simon/sage/local/lib/python2.7/site-packages/sage/misc/cachefunc.pyx
in sage.misc.cachefunc.CachedMethodCaller._instance_call
(build/cythonized/sage/misc/cachefunc.c:9717)()
1824 True
1825 """
-> 1826 return self.f(self._instance, *args, **kwds)
1827
1828 cdef fix_args_kwds(self, tuple args, dict kwds):
/home/simon/sage/local/lib/python2.7/site-packages/sage/modules/torsion_quadratic_module.pyc
in orthogonal_group(self, gens, check)
1196 gens = [matrix(g) for g in _isom_fqf(self)]
1197 ambient = AbelianGroupGap(self.invariants()).aut()
-> 1198 gens = [ambient(g) for g in gens]
1199 gens = tuple(g for g in gens if g != ambient.one())
1200 return FqfOrthogonalGroup(ambient, gens, self, check=check)
/home/simon/sage/local/lib/python2.7/site-packages/sage/structure/parent.pyx
in sage.structure.parent.Parent.__call__
(build/cythonized/sage/structure/parent.c:9171)()
897 if mor is not None:
898 if no_extra_args:
--> 899 return mor._call_(x)
900 else:
901 return mor._call_with_args(x, args, kwds)
/home/simon/sage/local/lib/python2.7/site-packages/sage/structure/coerce_maps.pyx
in sage.structure.coerce_maps.DefaultConvertMap_unique._call_
(build/cythonized/sage/structure/coerce_maps.c:4547)()
160 print(type(C), C)
161 print(type(C._element_constructor),
C._element_constructor)
--> 162 raise
163
164 cpdef Element _call_with_args(self, x, args=(), kwds={}):
/home/simon/sage/local/lib/python2.7/site-packages/sage/structure/coerce_maps.pyx
in sage.structure.coerce_maps.DefaultConvertMap_unique._call_
(build/cythonized/sage/structure/coerce_maps.c:4439)()
155 cdef Parent C = self._codomain
156 try:
--> 157 return C._element_constructor(x)
158 except Exception:
159 if print_warnings:
/home/simon/sage/local/lib/python2.7/site-packages/sage/groups/abelian_gps/abelian_aut.pyc
in _element_constructor_(self, x, check)
289 if x in self._covering_matrix_ring:
290 dom = self._domain
--> 291 images = [dom(row).gap() for row in x.rows()]
292 x = dom.gap().GroupHomomorphismByImages(dom.gap(),
images)
293 from sage.modules.fg_pid.fgp_morphism import FGP_Morphism
/home/simon/sage/local/lib/python2.7/site-packages/sage/structure/parent.pyx
in sage.structure.parent.Parent.__call__
(build/cythonized/sage/structure/parent.c:9171)()
897 if mor is not None:
898 if no_extra_args:
--> 899 return mor._call_(x)
900 else:
901 return mor._call_with_args(x, args, kwds)
/home/simon/sage/local/lib/python2.7/site-packages/sage/structure/coerce_maps.pyx
in sage.structure.coerce_maps.DefaultConvertMap_unique._call_
(build/cythonized/sage/structure/coerce_maps.c:4547)()
160 print(type(C), C)
161 print(type(C._element_constructor),
C._element_constructor)
--> 162 raise
163
164 cpdef Element _call_with_args(self, x, args=(), kwds={}):
/home/simon/sage/local/lib/python2.7/site-packages/sage/structure/coerce_maps.pyx
in sage.structure.coerce_maps.DefaultConvertMap_unique._call_
(build/cythonized/sage/structure/coerce_maps.c:4439)()
155 cdef Parent C = self._codomain
156 try:
--> 157 return C._element_constructor(x)
158 except Exception:
159 if print_warnings:
/home/simon/sage/local/lib/python2.7/site-packages/sage/groups/abelian_gps/abelian_group_gap.pyc
in _element_constructor_(self, x, check)
357 x = gens_gap[0]**0
358 for i in range(len(exp)):
--> 359 x *= gens_gap[i]**exp[i]
360 x = x.gap()
361 return self.element_class(self, x, check=check)
/home/simon/sage/local/lib/python2.7/site-packages/sage/groups/libgap_wrapper.pyx
in sage.groups.libgap_wrapper.ElementLibGAP.__pow__
(build/cythonized/sage/groups/libgap_wrapper.c:7268)()
659 raise TypeError("exponent must be an integer")
660 P = self.parent()
--> 661 return P.element_class(P, self.gap() ** n)
662
663 def __invert__(self):
/home/simon/sage/local/lib/python2.7/site-packages/sage/groups/abelian_gps/abelian_group_gap.pyc
in __init__(self, parent, x, check)
71 sage: TestSuite(g).run()
72 """
---> 73 if check and x not in parent.gap():
74 raise ValueError("%s is not in the group %s" % (x,
parent))
75 ElementLibGAP.__init__(self, parent, x)
/home/simon/sage/local/lib/python2.7/site-packages/sage/libs/gap/element.pyx
in sage.libs.gap.element.GapElement.__contains__
(build/cythonized/sage/libs/gap/element.c:7087)()
581 from sage.libs.gap.libgap import libgap
582 GAP_IN = libgap.eval(r'\in')
--> 583 return GAP_IN(other, self).sage()
584
585 cpdef _type_number(self):
/home/simon/sage/local/lib/python2.7/site-packages/sage/libs/gap/element.pyx
in sage.libs.gap.element.GapElement_Function.__call__
(build/cythonized/sage/libs/gap/element.c:18667)()
2408 a = [x if isinstance(x, GapElement) else libgap(x) for
x in args]
2409
-> 2410 sig_on()
2411 try:
2412 GAP_Enter()
SystemError: calling remove_from_pari_stack() inside sig_on()
--
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 https://groups.google.com/group/sage-devel.
For more options, visit https://groups.google.com/d/optout.
