Changeset 4476 for branches/f4grobner
- Timestamp:
- 2016-06-15T08:42:32-07:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/f4grobner/5am-buchberger.lisp
r4470 r4476 111 111 ) 112 112 113 ;; Cyclic roots of degree 5 114 ;; poly_reduced_grobner([x+y+z+u+v,x*y+y*z+z*u+u*v+v*x,x*y*z+y*z*u+z*u*v+u*v*x+v*x*y,x*y*z*u+y*z*u*v+z*u*v*x+u*v*x*y+v*x*y*z,x*y*z*u*v-1],[u,v,x,y,z]); 115 (def-fixture buchberger-cyclic-5-context () 116 (let ((fl 117 (cdr (string->poly 118 "[x+y+z+u+v,x*y+y*z+z*u+u*v+v*x,x*y*z+y*z*u+z*u*v+u*v*x+v*x*y,x*y*z*u+y*z*u*v+z*u*v*x+u*v*x*y+v*x*y*z,x*y*z*u*v-1]" 119 '(u v x y z)))) 120 (gb (cdr (string->poly "[z+y+x+v+u, 121 8*z^12+2*y*z^11+979*z^7+231*y*z^6-55*y^2*z^5-987*z^2-233*y*z+55*y^2, 122 z^15+122*z^10-122*z^5-1, 123 1467*z^12+566*y*z^11+178981*z^7+69003*y*z^6+550*y^3*z^4-275*y^4*z^3 124 -1650*y^5*z^2-179073*z^2-550*y^6*z-69019*y*z-825*v*z-275*y^2 125 -275*v^2,(-z^11)-143*z^6-55*v*z^5+144*z+55*v, 126 (-346*z^12)-124*y*z^11-42218*z^7-15092*y*z^6-275*y^3*z^4+275*y^4*z^3 127 +440*y^5*z^2+42124*z^2+110*y^6*z+15106*y*z+275*v*z-275*v*y, 128 867*z^12+334*y*z^11+105776*z^7+40722*y*z^6+275*y^3*z^4-275*y^4*z^3 129 -1045*y^5*z^2-105873*z^2-330*y^6*z-40726*y*z+550*x*z-275*v*z 130 -550*y^2+275*x^2+275*v*x, 131 (-568*z^13)-232*y*z^12-69289*z^8-28336*y*z^7-550*y^4*z^4+550*y^5*z^3 132 +69307*z^3+275*y^6*z^2+28018*y*z^2-550*x*z^2+550*y^2*z 133 +550*x^2*z+275*x^3,(-z^11)-143*z^6-55*x*z^5+144*z+55*x, 134 1121*z^12+442*y*z^11+136763*z^7+53911*y*z^6+275*y^3*z^4-1210*y^5*z^2 135 -136674*z^2-440*y^6*z-53913*y*z+275*x*z-275*y^2-275*x*y, 136 1042*z^12+398*y*z^11+127116*z^7+48554*y*z^6-55*y^5*z^2-128103*z^2 137 -165*y^6*z-48787*y*z-55*y^7+55*y^2]" '(u v x y z))))) 138 (&body))) 113 139 140 141 (test buchberger-cyclic-5 142 "Buchberger algorithm - cyclic 5" 143 (with-fixture buchberger-cyclic-5-context () 144 (is-true (grobner-test gb fl)) 145 ;; NOTE: Cannot compare using SET-EXCLUSIVE-OR, as the Grobner basis 146 ;; FL was computed using "sugar" strategy and is different from 147 ;; the one obtained with BUCHBERGER, which uses the "normal" strategy. 148 (is (grobner-equal (buchberger fl) gb)) 149 ;;(is (every #'universal-equalp (parallel-buchberger fl) gb)) 150 ) 151 ) 114 152 115 153 (run! 'buchberger-suite)
Note:
See TracChangeset
for help on using the changeset viewer.