Changeset 3682

Timestamp:

2015-09-06T05:42:59-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

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

branches/f4grobner/polynomial.lisp

@@ -444 +444 @@
             :initial-value (leading-coefficient self))))
 
+(defgeneric poly-primitive-part (self) 
+  (:documentation "Divide polynomial SELF by gcd of its
+coefficients. Return SELF.")
+  (:method ((self poly))
+    (if (universal-zerop object)
+        (values object 1)
+        (let ((c (poly-content object)))
+          (mapc #'(lambda (x)
+                    (setf (term-coeff x) (divide (term-coeff x) c)))
+                (poly-termlist p))))
+    self))
+
+
 #|
 
@@ -451 +464 @@
   (polysaturation-extension ring f (list p)))
 
-(defun poly-primitive-part (ring p)
-  "Divide polynomial P with integer coefficients by gcd of its
-coefficients and return the result."
-  (declare (type ring ring) (type poly p))
-  (if (poly-zerop p)
-      (values p 1)
-      (let ((c (poly-content ring p)))
-        (values (make-poly-from-termlist 
-                 (mapcar
-                  #'(lambda (x)
-                      (make-term :monom (term-monom x)
-                                 :coeff (funcall (ring-div ring) (term-coeff x) c)))
-                  (poly-termlist p))
-                 (poly-sugar p))
-                c))))