Changeset 1019

Timestamp:

2015-06-10T06:50:55-07:00 (10 years ago)

Author:

Marek Rychlik

Message:

* empty log message *

File:

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

branches/f4grobner/parse.lisp

@@ -255 +255 @@
                       (ring (ro-ring ring-and-order))
                       (n (length vars))
-                      (basis (monom-basis (length vars))))
+                      (basis (variable-basis ring (length vars))))
   "Evaluate an expression EXPR as polynomial by substituting operators
 + - * expt with corresponding polynomial operators and variables VARS
@@ -324 +324 @@
 
 ;; Return the standard basis of the monomials in n variables
-(defun monom-basis (n &aux
-                      (basis
-                       (copy-tree
-                        (make-list n :initial-element
-                                   (list (cons
-                                          (make-list
-                                           n
-                                           :initial-element 0)
-                                          1))))))
+(defun variable-basis (ring n &aux (basis (make-list n)))
   "Generate a list of monomials ((1 0 ... 0) (0 1 0 ... 0) ... (0 0 ... 1)
 which correspond to linear monomials X1, X2, ... XN."
   (dotimes (i n basis)
-    (setf (elt (caar (elt basis i)) i) 1)))
+    (setf (elt basis i) (make-variable ring n i))))