close Warning: Can't synchronize with repository "(default)" (The repository directory has changed, you should resynchronize the repository with: trac-admin $ENV repository resync '(default)'). Look in the Trac log for more information.

Changeset 4476 for branches


Ignore:
Timestamp:
2016-06-15T08:42:32-07:00 (8 years ago)
Author:
Marek Rychlik
Message:

Summary: Added Cyclic-5 example

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/f4grobner/5am-buchberger.lisp

    r4470 r4476  
    111111  )
    112112
     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)))
    113139
     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  )
    114152
    115153(run! 'buchberger-suite)
Note: See TracChangeset for help on using the changeset viewer.