[Maxima-commits] CVS: maxima/share/linearalgebra linalg-extra.lisp, 1.11, 1.12 test-linalg.mac, 1.15, 1.16

[Barton W. <wil...@us...>]

Update of /cvsroot/maxima/maxima/share/linearalgebra
In directory sc8-pr-cvs16.sourceforge.net:/tmp/cvs-serv21055/share/linearalgebra

Modified Files:
	linalg-extra.lisp test-linalg.mac 
Log Message:
o append tests for sylvester matrx 
o fix missing variialble in sylvester matrix

Index: linalg-extra.lisp
===================================================================
RCS file: /cvsroot/maxima/maxima/share/linearalgebra/linalg-extra.lisp,v
retrieving revision 1.11
retrieving revision 1.12
diff -u -d -r1.11 -r1.12
--- linalg-extra.lisp	29 Feb 2008 21:50:49 -0000	1.11
+++ linalg-extra.lisp	1 Mar 2008 04:07:00 -0000	1.12
@@ -142,7 +142,7 @@
 			(t '$pnz))))))))
 
 (defun $sylvester_matrix (p q z)
-  (let ((p-coeff nil) (q-coeff nil) (mat nil))
+  (let ((p-coeff nil) (q-coeff nil) (mat nil) (p-deg) (q-deg))
     (if (or (not ($polynomialp p `((mlist) ,z) `((lambda) ((mlist) c) (($freeof) ,z c))))
 	    (not ($polynomialp q `((mlist) ,z) `((lambda) ((mlist) c) (($freeof) ,z c)))))
 	(merror "The first two arguments to 'sylvester_matrix' must be polynomials"))

Index: test-linalg.mac
===================================================================
RCS file: /cvsroot/maxima/maxima/share/linearalgebra/test-linalg.mac,v
retrieving revision 1.15
retrieving revision 1.16
diff -u -d -r1.15 -r1.16
--- test-linalg.mac	19 Feb 2008 23:33:40 -0000	1.15
+++ test-linalg.mac	1 Mar 2008 04:07:00 -0000	1.16
@@ -634,5 +634,44 @@
 linalg_rank(matrix([0,1],[0,1]));
 1$
 
+sylvester_matrix(0, 0, x);
+matrix()$
+
+sylvester_matrix(1, 0, x);
+matrix()$
+
+sylvester_matrix(0, 1, x);
+matrix()$
+
+sylvester_matrix(1, 1, x);
+matrix()$
+
+sylvester_matrix(x, 1, x);
+matrix([1])$
+
+sylvester_matrix(1, x, x);
+matrix([1])$
+
+sylvester_matrix(x, x, x);
+matrix([1,0],[1,0])$
+
+sylvester_matrix(-x, x, x);
+matrix([-1,0],[1,0])$
+
+sylvester_matrix(a+b*x, x, x);
+matrix([b,a],[1,0])$
+
+sylvester_matrix(a+b*x, x^2, x);
+matrix([b,a,0],[0,b,a],[1,0,0])$
+
+
+
+
+
+
+
+
+
+