print ps = """ ( ( z1^4 z2^2 - w1 ), ( z1^5 z2^3 - w2 ), ( z1 - w3 ), ( z2 - w4 ) ) """ f = Ideal(r, ps) print "Ideal: " + str(f) print rg = f.GB() print "seq Output:", rg print pf = """ ( ( z1^37 z2^20 ) ) """ fp = Ideal(r, pf) print "Ideal: " + str(fp) print nf = fp.NF(rg) print "NFs: " + str(nf)
ps = """ ( ( y1 + y2 - 1 ), ( x1 - y1^2 - y1 - y2 ), ( x2 - y1 - y2^2 ), ( x3 - y1 y2 ) ) """ r = Ring(rs) print "Ring: " + str(r) i = Ideal(r, ps) print "Ideal: " + str(i) g = i.GB() print "seq GB:", g rsi = """ # polynomial ring: Rat(x1,x2,x3) G """ ri = Ring(rsi) print "Ring: " + str(ri) y = application.Ideal(g.pset).intersect(ri.ring) len = y.list.size() print "seq intersect y: ", y rs = """