#13458: Map to the Weierstrass form
--------------------------------------+-------------------------------------
Reporter: vbraun | Owner: AlexGhitza
Type: enhancement | Status: needs_review
Priority: major | Milestone: sage-5.11
Component: algebraic geometry | Resolution:
Keywords: | Work issues:
Report Upstream: N/A | Reviewers: Andrey Novoseltsev
Authors: Volker Braun | Merged in:
Dependencies: #13084 | Stopgaps:
--------------------------------------+-------------------------------------
Comment (by novoselt):
When I try
{{{
for P in ReflexivePolytopes(2):
E = CPRFanoToricVariety(P).anticanonical_hypersurface()
p = E.defining_polynomials()[0]
print WeierstrassForm(p, transformation=True)
}}}
it crashes on the 10th polytope with
{{{
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("Zm9yIFAgaW4gUmVmbGV4aXZlUG9seXRvcGVzKDIpOgogICAgRSA9IENQUkZhbm9Ub3JpY1ZhcmlldHkoUCkuYW50aWNhbm9uaWNhbF9oeXBlcnN1cmZhY2UoKQogICAgcCA9IEUuZGVmaW5pbmdfcG9seW5vbWlhbHMoKVswXQogICAgcHJpbnQgV2VpZXJzdHJhc3NGb3JtKHAsIHRyYW5zZm9ybWF0aW9uPVRydWUp"),globals())+"\\n");
execfile(os.path.abspath("___code___.py"))
File "", line 1, in <module>
File "/tmp/tmpZDkHaU/___code___.py", line 3, in <module>
exec compile(u'for P in ReflexivePolytopes(_sage_const_2 ):\n E =
CPRFanoToricVariety(P).anticanonical_hypersurface()\n p =
E.defining_polynomials()[_sage_const_0 ]\n print WeierstrassForm(p,
transformation=True)
File "", line 4, in <module>
File "lazy_import.pyx", line 313, in
sage.misc.lazy_import.LazyImport.__call__ (sage/misc/lazy_import.c:2475)
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/schemes/toric/weierstrass.py", line 479, in WeierstrassForm
return WeierstrassMap(polynomial, variables=variables)
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/schemes/toric/weierstrass_covering.py", line 255, in
WeierstrassMap
X,Y,Z = WeierstrassMap_P2(polynomial_aff, variables_aff)
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/schemes/toric/weierstrass_covering.py", line 339, in
WeierstrassMap_P2
H = cubic.Hessian()
File "cachefunc.pyx", line 1722, in
sage.misc.cachefunc.CachedMethodCallerNoArgs.__call__
(sage/misc/cachefunc.c:9112)
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/rings/invariant_theory.py", line 1615, in Hessian
return 1/F(216) * H.det()
File "matrix2.pyx", line 1257, in sage.matrix.matrix2.Matrix.det
(sage/matrix/matrix2.c:9875)
File "matrix_mpolynomial_dense.pyx", line 584, in
sage.matrix.matrix_mpolynomial_dense.Matrix_mpolynomial_dense.determinant
(sage/matrix/matrix_mpolynomial_dense.cpp:5877)
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/rings/polynomial/multi_polynomial_ring.py", line 456, in
__call__
return x.sage_poly(self)
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/interfaces/singular.py", line 1653, in sage_poly
self.name(),variable_str)).split(",")
File "/home/novoselt/sage-5.11.beta3/local/lib/python2.7/site-
packages/sage/interfaces/singular.py", line 590, in eval
raise SingularError('Singular error:\n%s'%s)
sage.interfaces.singular.SingularError: Singular error:
? `sage72` is not defined
? error occurred in or before STDIN line 159:
`string(coef(sage72,z0*z1*z2*z3));`
? wrong type declaration. type 'help string;'
}}}
Running it just for the 10th is OK, so looks more like a singular
interface issue, but may be worth investigation...
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/13458#comment:13>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" 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-trac.
For more options, visit https://groups.google.com/groups/opt_out.
