Changeset 1257 for branches

Timestamp:

2015-06-11T15:47:19-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

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

branches/f4grobner/division.lisp

@@ -74 +74 @@
 remainder, and P is the intermediate value which eventually becomes
 0."
+  (format t "Running loop-invariant check now.~%")
+  (format t "C=~A~%F=~A~%A=~A~%FL=~A~%R=~A~%P=~A~%" c f a fl r p)
   (flet ((p-add (x y) (poly-add ring-and-order x y))
          (p-sub (x y) (poly-sub ring-and-order x y))
          (p-mul (x y) (poly-mul ring-and-order x y)))
-    (let ((prod (inner-product a fl p-add p-mul p-zero)))
-      (poly-zerop
-       (p-sub
-        (scalar-times-poly ring c f)
-        (reduce #'p-add (list prod r p)))))))
+    (let* ((prod (inner-product a fl p-add p-mul p-zero))
+           (succeeded-p 
+            (poly-zerop
+             (p-sub
+              (scalar-times-poly ring c f)
+              (reduce #'p-add (list prod r p))))))
+      (unless succeeded-p (warn "Check failed"))
+      succeeded-p)))