Changeset 2785

Timestamp:

2015-06-20T20:43:13-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

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

branches/f4grobner/polynomial.lisp

@@ -192 +192 @@
 
 
+;; Multiplication of polynomials
+;; Non-destructive version
+(defun termlist-mul (ring-and-order p q
+                     &aux (ring (ro-ring ring-and-order)))
+  (declare (ring-and-order ring-and-order))
+  (cond ((or (endp p) (endp q)) nil)    ;p or q is 0 (represented by NIL)
+        ;; If p=p0+p1 and q=q0+q1 then pq=p0q0+p0q1+p1q
+        ((endp (cdr p))
+         (term-times-termlist ring (car p) q))
+        ((endp (cdr q))
+         (termlist-times-term ring p (car q)))
+        (t
+         (let ((head (term-mul-lst ring (termlist-lt p) (termlist-lt q)))
+               (tail (termlist-add ring-and-order
+                                   (term-times-termlist ring (car p) (cdr q))
+                                   (termlist-mul ring-and-order (cdr p) q))))
+           (cond ((null head) tail)
+                 ((null tail) head)
+                 (t (nconc head tail)))))))
+
+
+
+#|
+
 (defun poly-standard-extension (plist &aux (k (length plist)))
   "Calculate [U1*P1,U2*P2,...,UK*PK], where PLIST=[P1,P2,...,PK]
@@ -204 +228 @@
       (incf-power (nth i plist) i))))
 
-#|
+
 
 (defun saturation-extension (ring f plist