Changeset 160 for branches/f4grobner

Timestamp:

2015-06-05T13:36:44-07:00 (9 years ago)

Author:

Marek Rychlik

Message:

 

Location:

branches/f4grobner

Files:

• f4-matrix.lisp (modified) (1 diff)
• mx-parsing.lisp (deleted)

branches/f4grobner/f4-matrix.lisp

@@ -23 +23 @@
 
 (macsyma-module f4-maxima)
+
+(defun f4-grobner-op (ring c1 c2 m f g)
+  "Returns C2*F-C1*M*G, where F and G are polynomials M is a monomial.
+Assume that the leading terms will cancel."
+  #+grobner-check(funcall (ring-zerop ring)
+                          (funcall (ring-sub ring)
+                                   (funcall (ring-mul ring) c2 (poly-lc f))
+                                   (funcall (ring-mul ring) c1 (poly-lc g))))
+  #+grobner-check(monom-equal-p (poly-lm f) (monom-mul m (poly-lm g)))
+  (poly-sub ring
+            (scalar-times-poly ring c2 f)
+            (scalar-times-poly ring c1 (monom-times-poly m g))))
 
 ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;