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 4206


Ignore:
Timestamp:
2016-06-04T11:57:38-07:00 (9 years ago)
Author:
Marek Rychlik
Message:

* empty log message *

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/f4grobner/division.lisp

    r4200 r4206  
    180180;; method would be to test NORMAL-FORM using POLY-PSEUDO-DIVIDE.
    181181;;
    182 (defun normal-form (f fl
    183                     &optional
    184                       (top-reduction-only $poly_top_reduction_only))
     182(defun normal-form (f fl &optional (top-reduction-only $poly_top_reduction_only))
    185183  #+grobner-check(when (null fl) (warn "normal-form: empty divisor list."))
    186184  (do ((r (make-zero-for f))
    187        (c 1)
     185       (c (make-unit-for (leading-coefficient f)))
    188186       (division-count 0))
    189187      ((or (universal-zerop f)
     
    194192         (when (universal-zerop r)
    195193           (debug-cgb " ---> 0")))
    196        (setf (poly-termlist f) (nreconc (poly-termlist r) (poly-termlist f)))
     194       (setf f (add-to f r))
    197195       (values f c division-count))
    198196    (declare (fixnum division-count)
Note: See TracChangeset for help on using the changeset viewer.