| 1 | (in-package :buchberger)
|
|---|
| 2 | (use-package :symbolic-polynomial)
|
|---|
| 3 |
|
|---|
| 4 | (proclaim '(special fl gb))
|
|---|
| 5 |
|
|---|
| 6 | (setf fl (cdr (string->poly "[x-3*u-3*u*v^2+u^3,y-3*v-3*u^2*v+v^3,z-3*u^2+3*v^2]" '(u v x y z)))
|
|---|
| 7 | gb (cdr (string->poly "[x-3*u*v^2+u^3-3*u,y+v^3-3*u^2*v-3*v,z+3*v^2-3*u^2,
|
|---|
| 8 | (-u*z)-3*x+6*u*v^2+9*u,(-v*z)+y-2*v^3-3*v,z^2+6*v^2*z-9*z-9*v*y+9*u*x,
|
|---|
| 9 | 4*u*v*z-3*u*y+3*v*x,2*u*z^2+6*x*z-18*u*z-9*u*v*y+9*v^2*x,
|
|---|
| 10 | (-8*u*z^3)-24*x*z^2+72*u*z^2-36*v^2*x*z+27*u*y^2-27*v*x*y,
|
|---|
| 11 | z^3+18*v^2*z^2-18*z^2-54*v*y*z+54*v^2*z+81*z+27*y^2-27*x^2,
|
|---|
| 12 | (-4*z^4)+48*z^3-108*v*y*z^2+108*z^2+135*y^2*z+324*v*y*z+108*x^2*z
|
|---|
| 13 | -1296*v^2*z-1944*z-243*v^2*y^2-648*y^2+243*v^2*x^2+648*x^2,
|
|---|
| 14 | 8*v*z^3-9*y*z^2+72*v*z^2+54*v^2*y*z-27*y*z-27*v*y^2+27*v*x^2,
|
|---|
| 15 | (-8*v*z^4)+12*y*z^3-96*v*z^3-216*v*z^2-135*v*y^2*z+324*y*z-27*v*x^2*z
|
|---|
| 16 | +81*y^3+81*v*y^2-81*x^2*y-81*v*x^2,
|
|---|
| 17 | (-64*v*z^6)+120*y*z^5-1152*v*z^5+288*y*z^4-5184*v*z^4-648*v*y^2*z^3
|
|---|
| 18 | -216*y*z^3+6912*v*z^3+81*y^3*z^2-9720*v*y^2*z^2
|
|---|
| 19 | -1539*x^2*y*z^2+31104*y*z^2+62208*v*z^2+8505*y^3*z
|
|---|
| 20 | +46656*v*y^2*z-8505*x^2*y*z-93312*y*z+729*v*y^4-23328*y^3
|
|---|
| 21 | -1458*v*x^2*y^2-23328*v*y^2+23328*x^2*y+729*v*x^4
|
|---|
| 22 | +23328*v*x^2,
|
|---|
| 23 | 8*z^6-72*z^5+648*v*y*z^4-648*z^4-945*y^2*z^3+5184*v*y*z^3-189*x^2*z^3
|
|---|
| 24 | +5832*z^3+972*y^2*z^2+17496*v*y*z^2-2430*x^2*z^2+8748*v*y^3*z
|
|---|
| 25 | -19683*y^2*z+2187*x^2*z-5103*y^4-4374*v*y^3+5832*x^2*y^2
|
|---|
| 26 | +4374*v*x^2*y-729*x^4,
|
|---|
| 27 | 8*z^7-48*z^6+648*v*y*z^5-864*z^5-945*y^2*z^4+5832*v*y*z^4-189*x^2*z^4
|
|---|
| 28 | +3888*z^4+81*y^2*z^3+17496*v*y*z^3-2997*x^2*z^3+17496*z^3
|
|---|
| 29 | +8748*v*y^3*z^2-16767*y^2*z^2+17496*v*y*z^2-5103*x^2*z^2
|
|---|
| 30 | -5103*y^4*z+5832*x^2*y^2*z-6561*y^2*z-729*x^4*z+6561*x^2*z
|
|---|
| 31 | -2187*y^4+4374*x^2*y^2-2187*x^4,
|
|---|
| 32 | 64*z^9-10368*z^7+1296*y^2*z^6-1296*x^2*z^6-34992*y^2*z^5-34992*x^2*z^5
|
|---|
| 33 | +419904*z^5+174960*y^2*z^4-174960*x^2*z^4-10935*y^4*z^3
|
|---|
| 34 | -56862*x^2*y^2*z^3+314928*y^2*z^3-10935*x^4*z^3+314928*x^2*z^3
|
|---|
| 35 | +118098*y^4*z^2-118098*x^4*z^2+59049*y^4*z-118098*x^2*y^2*z
|
|---|
| 36 | +59049*x^4*z+19683*y^6-59049*x^2*y^4+59049*x^4*y^2-19683*x^6]" '(u v x y z))))
|
|---|
| 37 |
|
|---|
| 38 |
|
|---|
| 39 |
|
|---|
| 40 | (grobner-test gb fl)
|
|---|
| 41 | (every #'universal-equalp (buchberger fl) gb)
|
|---|
| 42 |
|
|---|
