Sage Folks,
I am running into what I think is a bug in a Sage calculation of a
zero-dimensional variety. Attached is a test case that illustrates
the problem. The test case shows that the variety contains a couple
extra points that are not really killed by the corresponding ideal.
Thanks in advance for helping me with this - and if I'm doing
something wrong, please let me know.
Jeff Stroomer
_______________________________________________________________________
#! /usr/bin/env sage
import sys, os
from sage.all import *
R = PolynomialRing(QQ, 4, names = 'u,z,y,x', order = 'lex')
u,z,y,x = R.gens()
G = [
x^2 + 6/91*z + 6/13*y - 409/91*x + 66/91,
y*x - 6/91*z - 45/13*y - 410/91*x + 1572/91,
y^2 + 66/91*z - 129/13*y - 222/91*x + 2000/91,
z*x - 277/91*z - 4/13*y - 607/91*x + 2049/91,
z*y - 53/7*z - 5*y - 12/7*x + 313/7,
z^2 - 849/91*z - 30/13*y - 48/91*x + 2582/91,
u - 17/91*z - 4/13*y - 9/91*x - 187/91,
]
I = ideal(G)
V = I.variety()
for v in V:
print v
for g in G:
print ' %s(%s): %s' % (g, v, g.subs(v))
--
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-support
URL: http://www.sagemath.org
