Changeset 4458 for branches

Timestamp:

2016-06-13T20:41:57-07:00 (8 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

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

branches/f4grobner/test1.lisp

@@ -1 +1 @@
 (in-package :division)
 
+(proclaim '(special f fl))
+
+(setf f (string->poly "x^2" '(x y))
+      ;;y-sq (string->poly "y^2" '(x y))
+      fl (cdr (string->poly "[x+y,x-2*y]" '(x y)))
+      ;;quotients (cdr (string->poly "[x-y,0]" '(x y)))
+      )
 #+nil
-(let ((f (string->poly "x^2" '(x y)))
-      ;;(y-sq (string->poly "y^2" '(x y)))
-      (fl (cdr (string->poly "[x+y,x-2*y]" '(x y))))
-      ;;(quotients (cdr (string->poly "[x-y,0]" '(x y))))
-      )
-  (multiple-value-bind (f c division-count)
-      (step (normal-form f fl))
-    (format t "Quotient:        ~S~%" (->sexp f))
-    (format t "Coefficient:     ~S~%" c)
-    (format t "Division count:  ~S~%" division-count)))
+(multiple-value-bind (f c division-count)
+    (normal-form f fl)
+  (format t "Quotient:        ~S~%" (->sexp f))
+  (format t "Coefficient:     ~S~%" c)
+  (format t "Division count:  ~S~%" division-count))
 
-(let ((f (string->poly "x^2" '(x y)))
-      ;;(y-sq (string->poly "y^2" '(x y)))
-      (fl (cdr (string->poly "[x+y,x-2*y]" '(x y))))
-      ;;(quotients (cdr (string->poly "[x-y,0]" '(x y))))
-      )
-  (multiple-value-bind (a r c division-count)
-      (poly-pseudo-divide f fl)
-    (format t "Quotients:        ~S~%" (poly->string (cons :[ a) '(x y)))
-    (format t "Remainder:        ~S~%" (->sexp r '(x y)))
-    (format t "Coefficient:     ~S~%" c)
-    (format t "Division count:  ~S~%" division-count)))
+(multiple-value-bind (a r c division-count)
+    (poly-pseudo-divide f fl)
+  (format t "Quotients:        ~S~%" (poly->string (cons :[ a) '(x y)))
+  (format t "Remainder:        ~S~%" (->sexp r '(x y)))
+  (format t "Coefficient:     ~S~%" c)
+  (format t "Division count:  ~S~%" division-count))