## [Maxima-commits] CVS: maxima/share/linearalgebra announcement.txt,1.1,1.2 linalg.mac,1.1,1.2 linalg.usage,1.1,1.2 matrixexp.usage,1.1,1.2

 [Maxima-commits] CVS: maxima/share/linearalgebra announcement.txt,1.1,1.2 linalg.mac,1.1,1.2 linalg.usage,1.1,1.2 matrixexp.usage,1.1,1.2 From: Robert Dodier - 2005-08-30 01:45:13 Update of /cvsroot/maxima/maxima/share/linearalgebra In directory sc8-pr-cvs1.sourceforge.net:/tmp/cvs-serv4600 Modified Files: announcement.txt linalg.mac linalg.usage matrixexp.usage Log Message: Change some names written in uppercase to lowercase. The only one which makes a functional difference is an instance of POLYNOMIALP -> polynomialp in linalg.mac. Index: announcement.txt =================================================================== RCS file: /cvsroot/maxima/maxima/share/linearalgebra/announcement.txt,v retrieving revision 1.1 retrieving revision 1.2 diff -u -d -r1.1 -r1.2 --- announcement.txt 29 Aug 2005 14:03:53 -0000 1.1 +++ announcement.txt 30 Aug 2005 01:45:03 -0000 1.2 @@ -11,7 +11,7 @@ (%o1) c:/maxima/linearalgebra/linalg.mac (%i2) declare(z,complex)\$ (%i3) conjugate([z,z^2,sqrt(z)]); -(%o3) [CONJUGATE(z),CONJUGATE(z)^2,CONJUGATE(sqrt(z))] +(%o3) [conjugate(z),conjugate(z)^2,conjugate(sqrt(z))] Off the branch cut of sqrt, conjugate and sqrt commute. Telling Maxima that z isn't on branch cut allows it to transform @@ -19,7 +19,7 @@ (%i4) assume(imagpart(z) > 0)\$ (%i5) conjugate(sqrt(z)); -(%o5) sqrt(CONJUGATE(z)) +(%o5) sqrt(conjugate(z)) (2) New functions for generating matrices: @@ -46,7 +46,7 @@ the matrix is numeric. (%i1) load("linalg"); -Warning - you are redefining the MACSYMA function RANK +Warning - you are redefining the Maxima function RANK (%o1) c:/maxima/linearalgebra/linalg.mac (%i2) lu_factor(matrix([a,b],[c,d])); (%o2) [matrix([a,b],[c/a,d-(b*c)/a]),[1,2]] @@ -84,4 +84,4 @@ As always, let me know how I can improve this code. -Barton Willis \ No newline at end of file +Barton Willis Index: linalg.mac =================================================================== RCS file: /cvsroot/maxima/maxima/share/linearalgebra/linalg.mac,v retrieving revision 1.1 retrieving revision 1.2 diff -u -d -r1.1 -r1.2 --- linalg.mac 29 Aug 2005 14:03:53 -0000 1.1 +++ linalg.mac 30 Aug 2005 01:45:03 -0000 1.2 @@ -287,7 +287,7 @@ funmake('matrix, makelist(makelist(vars[i]^j,j,0,n-1),i,1,n))); polytocompanion(p,x) := block([n], - if not POLYNOMIALP(p,[x], lambda([e], freeof(x,e))) then + if not polynomialp(p,[x], lambda([e], freeof(x,e))) then error("First argument to 'polytocompanion' must be a polynomial"), p : expand(p), n : hipow(p,x), Index: linalg.usage =================================================================== RCS file: /cvsroot/maxima/maxima/share/linearalgebra/linalg.usage,v retrieving revision 1.1 retrieving revision 1.2 diff -u -d -r1.1 -r1.2 --- linalg.usage 29 Aug 2005 14:03:53 -0000 1.1 +++ linalg.usage 30 Aug 2005 01:45:03 -0000 1.2 @@ -303,7 +303,7 @@ (%i1) batch("linalg.demo"); batching #pc:/maxima/linearalgebra/linalg.demo -(%i2) LOAD(linalg) +(%i2) load(linalg) (%o2) c:/maxima/linearalgebra/linalg.mac (%i3) m:matrix([1,2],[1,2]) (%o3) matrix([1,2],[1,2]) @@ -311,17 +311,17 @@ (%o4) span(matrix([1],[-1/2])) (%i5) columnspace(m) (%o5) span(matrix([1],[1])) -(%i6) ptriangularize(m-z*IDENT(2),z) +(%i6) ptriangularize(m-z*ident(2),z) (%o6) matrix([1,2-z],[0,3*z-z^2]) -(%i7) m:matrix([1,2,3],[4,5,6],[7,8,9])-z*IDENT(3) +(%i7) m:matrix([1,2,3],[4,5,6],[7,8,9])-z*ident(3) (%o7) matrix([1-z,2,3],[4,5-z,6],[7,8,9-z]) (%i8) mm:ptriangularize(m,z) (%o8) matrix([4,5-z,6],[0,66/49,-z^2/7+(102*z)/49+132/49],[0,0,(49*z^3)/264-(245*z^2)/88-(147*z)/44]) -(%i9) ALGEBRAIC:true +(%i9) algebraic:true (%o9) true -(%i10) TELLRAT(mm[3,3]) +(%i10) tellrat(mm[3,3]) (%o10) [z^3-15*z^2-18*z] -(%i11) mm:RATSIMP(mm) +(%i11) mm:ratsimp(mm) (%o11) matrix([4,5-z,6],[0,66/49,-(7*z^2-102*z-132)/49],[0,0,0]) (%i12) nullspace(mm) (%o12) span(matrix([1],[(z^2-14*z-16)/8],[-(z^2-18*z-12)/12])) @@ -329,7 +329,7 @@ (%o13) matrix([1,2,3,4],[5,6,7,8],[9,10,11,12],[13,14,15,16]) (%i14) column_space(m) (%o14) column_space(matrix([1,2,3,4],[5,6,7,8],[9,10,11,12],[13,14,15,16])) -(%i15) APPLY('orthogonal_complment,ARGS(nullspace(TRANSPOSE(m)))) +(%i15) apply('orthogonal_complment,args(nullspace(transpose(m)))) (%o15) orthogonal_complment(matrix([0],[1],[-2],[1]),matrix([1],[0],[-3],[2])) Index: matrixexp.usage =================================================================== RCS file: /cvsroot/maxima/maxima/share/linearalgebra/matrixexp.usage,v retrieving revision 1.1 retrieving revision 1.2 diff -u -d -r1.1 -r1.2 --- matrixexp.usage 29 Aug 2005 14:03:53 -0000 1.1 +++ matrixexp.usage 30 Aug 2005 01:45:03 -0000 1.2 @@ -65,7 +65,7 @@ (C3) matrixfun(lambda([x], f(x)),m); Attempt to differentiate with respect to a number: -1 -- an error. Quitting. To debug this try DEBUGMODE(TRUE);) +1 -- an error. Quitting. To debug this try debugmode(true); (C4) load("pdiff")\$ (C5) matrixfun(lambda([x], f(x)),m); @@ -77,4 +77,4 @@ Barton Willis wrote matexp. This code is covered by the GNU public -license. Please send bug reports to the Maxima mailing list. \ No newline at end of file +license. Please send bug reports to the Maxima mailing list.