Changeset 4392 for branches/f4grobner

Timestamp:

2016-06-07T14:19:54-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

• branches/f4grobner/test5.lisp (modified) (1 diff)

branches/f4grobner/test5.lisp

@@ +1 @@
+
 ;;poly_grobner([x^3-x^2+2,y-2*x^2+1,z-3*x+5],[x,y,z]);
 ;;[x^3-x^2+2,y-2*x^2+1,z-3*x+5,2*z^2+20*z-9*y+41,(-z^3)-12*z^2-45*z-104]
 
 (in-package :5am-buchberger)
+(proclaim '(special fl ref-gb sp sp-rem))
 
+(setf grobner-debug::$poly_grobner_debug nil)
 
 (setf fl (cdr (string->poly "[x^3-x^2+2,y-2*x^2+1,z-3*x+5]" '(x y z)))
       ref-gb (cdr (string->poly "[x^3-x^2+2,y-2*x^2+1,z-3*x+5,2*z^2+20*z-9*y+41,(-z^3)-12*z^2-45*z-104]" '(x y z))))
 
-(format t "~S~%" (poly->string (cons :[ fl)))
-(format t "~S~%" (poly->string (cons :[ ref-gb)))
+(format t "~%")
+(format t "FL --> ~A~%" (poly->string (cons :[ fl)))
+(format t "REF-GB --> ~A~%" (poly->string (cons :[ ref-gb)))
 
-(setf sp-0-1 (s-polynomial (elt fl 0) (elt fl 1)))
-(format t "~S~%" (poly->string sp-0-1))
+(loop for i from 0 below (length fl)
+     do
+     (loop for j from (1+ i) below (length fl)
+          do
+          (setf sp (s-polynomial (elt fl i) (elt fl j)))
+          (format t "SPOLY(~A,~A) --> ~A~%" (poly->string (elt fl i)) (poly->string (elt fl j)) (poly->string sp))
+          (setf sp-rem (normal-form sp fl))
+          (format t "NORMAL-FORM(~A,~A) --> ~A~%" (poly->string sp-rem) (poly->string (cons :[ fl)) (poly->string sp-rem))
+          )
+)
 
-(setf sp-0-2 (s-polynomial (elt fl 0) (elt fl 2)))
-(format t "~S~%" (poly->string sp-0-2))
 
-(setf sp-1-2 (s-polynomial (elt fl 1) (elt fl 2)))
-(format t "~S~%" (poly->string sp-1-2))
-
-(setf grobner-debug::$poly_grobner_debug t)
 ;;(setf gb (buchberger fl))