py3: polyhedron pynormaliz backend does not handle large coefficients correctly

[videlec]

sage: x = polygen(ZZ)
sage: K.<sqrt2> = NumberField(x^2 - 2, embedding=AA(2).sqrt())
sage: P = Polyhedron(vertices=[(2**100, sqrt2), (sqrt2,1)], backend='normaliz')
Error in rational_handler: Python int too large to convert to C long
Error in nfelem_handler: unsupported operand type(s) for +=: 'NoneType' and 'sage.rings.number_field.number_field_element_quadratic.NumberFieldElement_quadratic'
Error in rational_handler: Python int too large to convert to C unsigned long
Error in rational_handler: Python int too large to convert to C unsigned long
Error in nfelem_handler: unsupported operand type(s) for *: 'NoneType' and 'sage.rings.number_field.number_field_element_quadratic.NumberFieldElement_quadratic'
Error in rational_handler: Python int too large to convert to C long
Error in rational_handler: Python int too large to convert to C long
Error in nfelem_handler: unsupported operand type(s) for *: 'NoneType' and 'sage.rings.number_field.number_field_element_quadratic.NumberFieldElement_quadratic'
Error in rational_handler: Python int too large to convert to C long
Error in rational_handler: Python int too large to convert to C long
Error in nfelem_handler: unsupported operand type(s) for *: 'NoneType' and 'sage.rings.number_field.number_field_element_quadratic.NumberFieldElement_quadratic'
sage: P.vertices()
(A vertex at (sqrt2, 1), A vertex at (0, sqrt2))

Even though PyNormaliz is perfectly fine working with large coefficients

sage: cone = PyNormaliz.NmzCone(cone=[],
....:          number_field=['a^2 - 2', 'a', '[1.414213562373095 +/- 2.99e-16]'],
....:          vertices=[[1267650600228229401496703205376r, [[0r, 1r], [1r, 1r]], 1r], [[[0r, 1r], [1r, 1r]], 1r, 1r]])
sage: PyNormaliz.NmzResult(cone, "VerticesOfPolyhedron")
[[[[0, 1], [1, 1]], [[1, 1], [0, 1]], [[1, 1], [0, 1]]],
 [[[1267650600228229401496703205376, 1], [0, 1]],
  [[0, 1], [1, 1]],
  [[1, 1], [0, 1]]]]
sage: PyNormaliz.NmzResult(cone, "SupportHyperplanes")
[[[[0, 1267650600228229401496703205376],
   [-1, 1267650600228229401496703205376]],
  [[1, 1], [0, 1]],
  [[0, 1], [0, 1]]],
 [[[0, 2], [1, 2]], [[-1, 1], [0, 1]], [[0, 1], [0, 1]]]]

Depends on #28321

CC: @mkoeppe @jplab

Component: geometry

Issue created by migration from https://trac.sagemath.org/ticket/28370

[videlec]

Dependencies: #28321

[videlec]

comment:2

I just think that #28321 is to blame (Python 3 issue)

[videlec]

Author: Vincent Delecroix

[jplab]

comment:3

You mean that it would be a duplicate? i.e. that the backend is fine?

[videlec]

comment:4

It is probably not a problem with the backend.

[videlec]

comment:5

I confirm that it is solved with #28321... however, the fact that errors are not raised in the constructor is a bug in itself. In the example I wrote in the ticket description, I obtain a polytope P which is simply wrong. It is to my mind much better to get an error rather than false results.

[videlec]

Changed author from Vincent Delecroix to none

[jplab]

comment:7

Replying to @videlec:

I confirm that it is solved with #28321... however, the fact that errors are not raised in the constructor is a bug in itself. In the example I wrote in the ticket description, I obtain a polytope P which is simply wrong. It is to my mind much better to get an error rather than false results.

I agree. Where to raise this error? This should become another ticket?

[videlec]

comment:8

Replying to @jplab:

Replying to @videlec:

I confirm that it is solved with #28321... however, the fact that errors are not raised in the constructor is a bug in itself. In the example I wrote in the ticket description, I obtain a polytope P which is simply wrong. It is to my mind much better to get an error rather than false results.

I agree. Where to raise this error? This should become another ticket?

• Another ticket: yes
• What to do: just remove the try/except