LatticePoints crash

[mcmillanbb]

The following example comes up in a calculation. However, it crashes consistently and quite severely.
In a jupyter notebook, the python kernel itself crashes even before printing its usual stack trace.

import PyNormaliz
lower = [[0,1,0,0,0],[1,0,0,0,0],[1,1,0,1,0],[1,1,1,0,0],[1,1,1,1,0]]                                 
upper = [[0,-1,0,0,1],[-1,0,0,0,0],[-1,-1,0,-1,2],[-1,-1,-1,0,2],[-1,-1,-1,-1,3]]
PyNormaliz.Cone(inhom_inequalities=lower+upper).LatticePoints()

Happy to provide more details if needed.

[mcmillanbb]

Ah, I just noticed the second normals of upper, lower are colinear, so the polyhedron is lower dimensional than expected. This seems to be important to reproducing.
The polyhedron is generated and consumed by other code, so a solution along lines of "manually restrict to that hyperplane" would not be ideal for me, unless there is some way to automatically do that.

[w-bruns]

There is definzely a bug. It makes Normaliz think that the polytope you define is a parallelotope what is not the case. At the time being, I can only suggest a workaround. It forbids the critical algorithm:

>>> import PyNormaliz
>>> lower = [[0,1,0,0,0],[1,0,0,0,0],[1,1,0,1,0],[1,1,1,0,0],[1,1,1,1,0]]                                 
>>> upper = [[0,-1,0,0,1],[-1,0,0,0,0],[-1,-1,0,-1,2],[-1,-1,-1,0,2],[-1,-1,-1,-1,3]]
>>> C = PyNormaliz.Cone(inhom_inequalities=lower+upper)
>>> C.Compute("LatticePoints", "NoProjection")
True
>>> C.LatticePoints();
[[0, 0, 0, 0, 1], [0, 0, 0, 1, 1], [0, 0, 0, 2, 1], [0, 0, 1, 0, 1], [0, 0, 1, 1, 1], [0, 0, 1, 2, 1], [0, 0, 2, 0, 1], [0, 0, 2, 1, 1], [0, 1, -1, 0, 1], [0, 1, -1, 1, 1], [0, 1, 0, -1, 1], [0, 1, 0, 0, 1], [0, 1, 0, 1, 1], [0, 1, 1, -1, 1], [0, 1, 1, 0, 1], [0, 1, 1, 1, 1]]
>>>