#11312: Speed up the computation of the Hilbert basis of a cone
-------------------------------------------+--------------------------------
Reporter: vbraun | Owner: mhampton
Type: enhancement | Status: needs_work
Priority: major | Milestone: sage-4.7.1
Component: geometry | Keywords:
Work_issues: non-full-dimensional errors | Upstream: N/A
Reviewer: Andrey Novoseltsev | Author: Volker Braun
Merged: | Dependencies:
-------------------------------------------+--------------------------------
Changes (by novoselt):
* status: needs_review => needs_work
* reviewer: => Andrey Novoseltsev
* work_issues: => non-full-dimensional errors
Comment:
The speed-up looks impressive! About 20x on trivial examples (no surprise
here since we are avoiding system calls to both cdd and PALP), on a
"complicated" 6-d example that I picked the new version worked for 15sec
and the old one didn't finish before my patience ran out ;-) Am I right
that this speed is partially due to taking into account the special
structure of the polytopes in which you are computing integral points?
There are, however, issues with non-full-dimensional cones::
{{{
sage: Cone([(1,1), (-1,1)], check=False).Hilbert_basis()
(N(1, 1), N(-1, 1), N(0, 1))
sage: Cone([(1,1,0), (-1,1,0)], check=False).Hilbert_basis()
(N(1, 1, 0), N(-1, 1, 0))
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/11312#comment:2>
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 post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
