#17937: Cannot compute integral points of 0-dimensional Polyhedron
------------------------+----------------------------
   Reporter:  jdemeyer  |            Owner:
       Type:  defect    |           Status:  new
   Priority:  major     |        Milestone:  sage-6.6
  Component:  geometry  |         Keywords:
  Merged in:            |          Authors:
  Reviewers:            |  Report Upstream:  N/A
Work issues:            |           Branch:
     Commit:            |     Dependencies:
   Stopgaps:            |
------------------------+----------------------------
 A trivial case, but it should work anyway:
 {{{
 sage: P = Polyhedron([[]])
 sage: P
 A 0-dimensional polyhedron in ZZ^0 defined as the convex hull of 1 vertex
 sage: P.integral_points()
 ---------------------------------------------------------------------------
 IndexError                                Traceback (most recent call
 last)
 <ipython-input-13-2abbf2adff15> in <module>()
 ----> 1 P.integral_points()

 /usr/local/src/sage-git/local/lib/python2.7/site-
 packages/sage/geometry/polyhedron/base.pyc in integral_points(self,
 threshold)
    3826                 box_points<threshold:
    3827             from sage.geometry.integral_points import
 rectangular_box_points
 -> 3828             return rectangular_box_points(box_min, box_max, self)
    3829
    3830         # for more complicate polytopes, triangulate & use smith
 normal form

 /usr/local/src/sage-git/src/sage/geometry/integral_points.pyx in
 sage.geometry.integral_points.rectangular_box_points
 (build/cythonized/sage/geometry/integral_points.c:4761)()
     531     v = vector(ZZ, d)
     532     if not return_saturated:
 --> 533         for p in loop_over_rectangular_box_points(box_min,
 box_max, inequalities, d, count_only):
     534             #  v = vector(ZZ, orig_perm.action(p))   # too slow
     535             for i in range(0,d):

 /usr/local/src/sage-git/src/sage/geometry/integral_points.pyx in
 sage.geometry.integral_points.loop_over_rectangular_box_points
 (build/cythonized/sage/geometry/integral_points.c:5605)()
     580     while True:
     581         inequalities.prepare_inner_loop(p)
 --> 582         i_min = box_min[0]
     583         i_max = box_max[0]
     584         # Find the lower bound for the allowed region

 IndexError: list index out of range
 }}}

--
Ticket URL: <http://trac.sagemath.org/ticket/17937>
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/d/optout.

Reply via email to