1 | initial version |
According to the documentation of
sage: I.variety?
we read:
Note that with "ring=RR" or "CC", computation is done numerically and potentially inaccurately; in particular, the number of points in the real variety may be miscomputed. With "ring=AA" or "QQbar", computation is done exactly (which may be much slower, of course).
Being miscomputed might also include "missing solutions"... This triangular decomposition also highly depends on the ideal and its groebner basis computation.
Currently, it seems indeed that getting a 0-dimensional variety out of an ideal is rather hazardous.
I am using Bertini to do such computations, eventually is could be interfaced in Sage.
Otherwise, maybe the interface to PHC could do the thing:
http://doc.sagemath.org/html/en/reference/interfaces/sage/interfaces/phc.html