Changeset 4206

Timestamp:

2016-06-04T11:57:38-07:00 (8 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

• branches/f4grobner/division.lisp (modified) (2 diffs)

branches/f4grobner/division.lisp

@@ -180 +180 @@
 ;; method would be to test NORMAL-FORM using POLY-PSEUDO-DIVIDE.
 ;;
-(defun normal-form (f fl
-                    &optional
-                      (top-reduction-only $poly_top_reduction_only))
+(defun normal-form (f fl &optional (top-reduction-only $poly_top_reduction_only))
   #+grobner-check(when (null fl) (warn "normal-form: empty divisor list."))
   (do ((r (make-zero-for f))
-       (c 1)
+       (c (make-unit-for (leading-coefficient f)))
        (division-count 0))
       ((or (universal-zerop f)
@@ -194 +192 @@
          (when (universal-zerop r)
            (debug-cgb " ---> 0")))
-       (setf (poly-termlist f) (nreconc (poly-termlist r) (poly-termlist f)))
+       (setf f (add-to f r))
        (values f c division-count))
     (declare (fixnum division-count)