Changeset 3064

Timestamp:

2015-06-21T18:08:06-07:00 (10 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

• branches/f4grobner/polynomial.lisp (modified) (1 diff)

branches/f4grobner/polynomial.lisp

@@ -301 +301 @@
         plist))
 
-#|
-
-(defun saturation-extension (ring f plist 
+(defun saturation-extension (plist 
                              &aux 
-                               (k (length plist))
-                               (d (monom-dimension (poly-lm (car plist))))
-                               f-x plist-x)
+                               (plist (poly-saturation-extension plist))
+                               (dim (poly-dimension (car plist))))
   "Calculate [F, U1*P1-1,U2*P2-1,...,UK*PK-1], where PLIST=[P1,P2,...,PK]."
-  (declare (type ring ring))
-  (setf f-x (poly-list-add-variables f k)
-        plist-x (mapcar #'(lambda (x)
-                            (setf (poly-termlist x)
-                                  (nconc (poly-termlist x)
-                                         (list (make-term :monom (make-monom :dimension d)
-                                                          :coeff (funcall (ring-uminus ring) 
-                                                                          (funcall (ring-unit ring)))))))
-                            x)
-                        (poly-standard-extension plist)))
-  (append f-x plist-x))
+  (flet ((subtract-1 (p) 
+           (insert-item (make-instance 'term :coeff -1 :dimension dim))))
+    (mapc #'subtract-1 plist))
+  (nconc F plist))