Changeset 3694

Timestamp:

2015-09-06T06:37:59-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

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

branches/f4grobner/polynomial.lisp

@@ -467 +467 @@
   (scalar-divide-by object (poly-content object)))
 
+(defun saturation-extension (f plist 
+                             &aux 
+                               (k (length plist))
+                               (d (monom-dimension (poly-lm (car plist)))))
+  "Calculate [F', U1*P1-1,U2*P2-1,...,UK*PK-1], where F' has first K
+variables U1,U2,...,UK inserted as first K variables and
+PLIST=[P1,P2,...,PK]."
+  (setf f (poly-list-add-variables f k)
+        plist (standard-extension-1 plist))
+  (nconc f plist))
+
+
+(defun polysaturation-extension (f plist 
+                                 &aux (k (length plist))
+                                   (d (+ k (length (poly-lm (car plist))))))
+  "Calculate [F', U1*P1+U2*P2+...+UK*PK-1], where PLIST=[P1,P2,...,PK]."
+  (setf f (poly-list-add-variables f k)
+        plist (apply #'poly-append (poly-standard-extension plist))
+        (cdr (last (poly-termlist plist))) (list (make-term (make-monom d :initial-element 0)
+                                                            (funcall (ring-uminus ring) (funcall (ring-unit ring))))))
+  (append f (list plist)))
+
+
 (defun saturation-extension-1 (f p) 
   "Given family of polynomials F and a polynomial P, calculate [F', U*P-1], where