This is now
#27374 <https://trac.sagemath.org/ticket/27374>
On Wednesday, February 27, 2019 at 9:53:35 AM UTC+1, Simon Brandhorst wrote:
>
> ---------------------------------------------------------------------
> 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.
