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 4054 for branches


Ignore:
Timestamp:
2016-05-31T18:02:22-07:00 (8 years ago)
Author:
Marek Rychlik
Message:

* empty log message *

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/f4grobner/division.lisp

    r4052 r4054  
    6666(defun check-loop-invariant (c f a fl r p
    6767                             &aux
    68                                (p-zero (make-poly-zero))
     68                               (p-zero 0)
    6969                               (a (mapcar #'poly-reverse a))
    7070                               (r (poly-reverse r)))
     
    105105  (declare (type poly f) (list fl))
    106106  ;; Loop invariant: c*f=sum ai*fi+r+p, where p must eventually become 0
    107   (do ((r 0)
    108        (c 1)
    109        (a (make-list (length fl) :initial-element 0))
     107  (do ((r (make-zero-for f))
     108       (c (make-unit-for f))
     109       (a (make-list (length fl) :initial-element (make-zero-for f)))
    110110       (division-count 0)
    111111       (p f))
     
    203203                      (top-reduction-only $poly_top_reduction_only))
    204204  #+grobner-check(when (null fl) (warn "normal-form: empty divisor list."))
    205   (do ((r 0)
    206        (c 1)
     205  (do ((r (make-zero-for f))
     206       (c (make-zero-for f))
    207207       (division-count 0))
    208208      ((or (universal-zerop f)
     
    231231      (declare (ignore c))
    232232      (subtract
    233        (multiply cg (monom-times-poly mf f))
    234        (multiply cf (monom-times-poly mg g))))))
     233       (multiply cg (multiply mf f))
     234       (multiply cf (multiply mg g))))))
    235235
    236236(defun buchberger-criterion (g)
Note: See TracChangeset for help on using the changeset viewer.