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.

Ignore:
Timestamp:
2015-06-11T15:05:51-07:00 (9 years ago)
Author:
Marek Rychlik
Message:

* empty log message *

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/f4grobner/division.lisp

    r1241 r1242  
    6464            (scalar-times-poly-1 ring c1 (monom-times-poly m g))))
    6565
    66 (defun check-loop-invariant (ring-and-order c f0 a fl r p
     66(defun check-loop-invariant (ring-and-order c f a fl r p
    6767                             &aux
    6868                               (ring (ro-ring ring-and-order))
     
    7070  "Check loop invariant of division algorithms, when we divide a
    7171polynomial F by the list of polynomials FL. The invariant is the
    72 identity C*F0=SUM AI*FI+R+F, where F0 is the initial value of F, A is
     72identity C*F=SUM AI*FI+R+P, where F0 is the initial value of F, A is
    7373the list of partial quotients, R is the intermediate value of the
    74 remainder, and F is the intermediate value which eventually becomes
     74remainder, and P is the intermediate value which eventually becomes
    75750."
    76   (flet ((p-add (p q) (poly-add ring-and-order p q))
    77          (p-sub (p q) (poly-sub ring-and-order p q))
    78          (p-mul (p q) (poly-mul ring-and-order p q)))
     76  (flet ((p-add (x y) (poly-add ring-and-order x y))
     77         (p-sub (x y) (poly-sub ring-and-order x y))
     78         (p-mul (x y) (poly-mul ring-and-order x y)))
    7979    (poly-zerop
    8080     (p-sub
    81       (scalar-times-poly ring c f0)
     81      (scalar-times-poly ring c f)
    8282      (reduce #'p-add
    8383              (list (inner-product a fl p-add p-mul p-zero)
Note: See TracChangeset for help on using the changeset viewer.