Changeset 3087

Timestamp:

2015-06-21T20:50:08-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

• branches/f4grobner/polynomial.lisp (modified) (3 diffs)

branches/f4grobner/polynomial.lisp

@@ -309 +309 @@
         (t (monom-dimension (leading-term poly)))))
 
-(defun mod-standard-extension (plist 
-                               &aux 
-                                 (k (length plist))
-                                 (plist (poly-standard-extension plist))
-                                 (nvars (poly-dimension (car plist))))
+(defun standard-extension-1 (plist 
+                             &aux 
+                               (k (length plist))
+                               (plist (poly-standard-extension plist))
+                               (nvars (poly-dimension (car plist))))
   "Calculate [U1*P1-1,U2*P2-1,...,UK*PK-1], where PLIST=[P1,P2,...,PK].
-Firstly, new K variables U1, U2, ..., UK, are inserted at the
-beginning of each polynomial.  Subsequently, P1, P2, ..., PK are
-destructively modified tantamount to replacing PI with UI*PI-1."
+Firstly, new K variables U1, U2, ..., UK, are inserted into each
+polynomial.  Subsequently, P1, P2, ..., PK are destructively modified
+tantamount to replacing PI with UI*PI-1."
   (flet ((subtract-1 (p) 
            (append-item p (make-instance 'term :coeff -1 :dimension (+ k nvars)))))
@@ -324 +324 @@
 
 
-(defun saturation-polynomial (F plist 
+(defun standard-sum (F plist 
                                  &aux 
                                    (k (length plist))
@@ -332 +332 @@
                                    ;; Set PLIST to [U1*P1,U2*P2,...,UK*PK]
                                    (plist (apply #'nconc (poly-standard-extension plist))))
-  "Calculate [F,U1*P1+U2*P2+...+UK*PK-1], where PLIST=[P1,P2,...,PK]."
-  ;; Add -1 as the last term
+  "Calculate the polynomial U1*P1+U2*P2+...+UK*PK-1, where PLIST=[P1,P2,...,PK].
+Firstly, new K variables, U1, U2, ..., UK, are inserted into each
+polynomial.  Subsequently, P1, P2, ..., PK are destructively modified
+tantamount to replacing PI with UI*PI, and the resulting polynomials
+are added.  It should be noted that the term order is not affected,
+which is equivalent to using a lexicographic order."
   (declare (type ring ring))
   (setf (cdr (last (poly-termlist plist)))
+        ;; Add -1 as the last term
         (list (make-term :monom (make-monom :dimension d)
                          :coeff (funcall (ring-uminus ring) (funcall (ring-unit ring))))))