fix bug in variety() for MPolynomialIdeals

malb commented 15 years ago

sage: B.<x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30>=BooleanPolynomialRing(30,order='lex')
sage: I1=ideal([x13*x15 + x1*x7*x13*x25 + x9*x22 + x17 + x18 + x3 + x13,1 + x6*x21 + x3,x2*x18 + 1 + x16*x25,x15 + x9 + x1*x10*x20,x23 + x9*x21*x23*x27 + x25 + x7,x13 + x1*x4*x15 + x6*x24 + 1 + x26*x27 + x3,1 + 1,1 + x13 + x22*x28 + x1*x29,x8 + x11*x29 + x15*x18 + x5*x11,x3*x9 + 1 + 1,x8,x14 + x14 + x6*x16*x28 + x6 + x13 + x23,x1*x15*x17 + 1 + x7 + x8 + x3 + x3*x15 + x12*x13*x25 + x3,x26,x4*x7*x8 + x3*x7 + x15*x17 + x15 + 1,x3*x17*x26*x28 + x2*x14*x20*x21 + x5*x18*x28 + x30,x5*x6 + x12 + x28 + x3*x25 + x12,1 + x7 + x16*x23 + 1 + 1,x1*x12*x24 + x4*x22*x26 + x6,x19*x21 + x2*x5*x9 + 1 + x21 + x6,x30 + x13,x17 + x10*x11 + x3 + x10 + x7*x8*x18,x14 + x7*x8*x11 + x5*x21 + x2*x15*x19 + x5*x7,x20*x21 + x15*x18*x19 + x9*x10,x2*x8*x10*x19*x22 + x8 + x14 + x16*x24 + x16*x21 + x29,x9*x25 + x10*x17*x18 + x2*x19 + x22*x25 + x6*x22*x25,x25 + x6 + x8*x26 + x10*x23 + x29 + x12 + x30 + x4*x13,x22 + x24 + x26 + x2 + x30 + x6*x10*x16 + x20,0,x28 + x8 + x9 + x2 + x15*x18*x29 + x16])
sage: G=I1.groebner_basis(faugere=False,linear_algebra_in_last_block=False)
sage: GG=Ideal(G)
sage: R.<x1,x2,x3,x4,x5,x6,x7,x8,x9,x10,x11,x12,x13,x14,x15,x16,x17,x18,x19,x20,x21,x22,x23,x24,x25,x26,x27,x28,x29,x30>=PolynomialRing(GF(2),30,order='lex')
sage: GB = ideal([R(f) for f in G]) + ideal(sage.rings.ideal.FieldIdeal(R))
sage: GB.variety()

this used to raise an exception, now it does what it is supposed to do.

Component: commutative algebra

Issue created by migration from https://trac.sagemath.org/ticket/5486

johnperry-math commented 15 years ago

comment:1

Attachment: variety_bugfix.patch.gz

doctests pass; changes make sense

85eec1a4-3d04-4b4d-b711-d4db03337c41 commented 15 years ago

comment:2

Merged in Sage 3.4.1.alpha0.

Cheers,

Michael

sagemath / sage

fix bug in variety() for MPolynomialIdeals #5486