Fwd: manipulate block matrices with maxima
- Subject: Fwd: manipulate block matrices with maxima
- From: Barton Willis
- Date: Sun, 2 May 2010 21:10:42 -0500
Maybe this will get you started:
(%i13) matrix_element_mult : "."$
(%i14) m : matrix([matrix([1,0],[0,-1]),matrix([0,1],[1,0])], [matrix([0,1],[1,0]),matrix([1,0],[0,-1])])$
(%i15) m-m;
(%o15) matrix([matrix([0,0],[0,0]),matrix([0,0],[0,0])],[matrix([0,0],[0,0]),matrix([0,0],[0,0])])
(%i16) 5 * m;
(%o16) matrix([matrix([5,0],[0,-5]),matrix([0,5],[5,0])],[matrix([0,5],[5,0]),matrix([5,0],[0,-5])])
(%i17) m.m;
(%o17) matrix([matrix([2,0],[0,2]),matrix([0,0],[0,0])],[matrix([0,0],[0,0]),matrix([2,0],[0,2])])
(%i18) get_lu_factors(lu_factor(m, 'noncommutingring));
(%o18) [matrix([matrix([1,0],[0,1]),matrix([0,0],[0,0])],[matrix([0,0],[0,0]),matrix([1,0],[0,1])]),matrix([matrix([1,0],[0,1]),matrix([0,0],[0,0])],[matrix([0,-1],[1,0]),matrix([1,0],[0,1])]),matrix([matrix([1,0],[0,-1]),matrix([0,1],[1,0])],[matrix([0,0],[0,0]),matrix([2,0],[0,-2])])]
(%i19) xreduce(".",%);
(%o19) matrix([matrix([1,0],[0,-1]),matrix([0,1],[1,0])],[matrix([0,1],[1,0]),matrix([1,0],[0,-1])])
--Barton
-----maxima-bounces at math.utexas.edu wrote: -----
>To:?maxima?<maxima at math.utexas.edu>
>From:?jmm?<mauricio.matera at gmail.com>
>Sent?by:?maxima-bounces at math.utexas.edu
>Date:?05/02/2010?08:33PM
>Subject:?[Maxima]?Fwd:?manipulate?block?matrices?with?maxima
>
>Hello,
>
>I?want?to?know?if?there?exist?a?package?to?manipulate?blocks?of?matrices,
>for?example,
>from?matrices?
>sx:matrix([0,1],[1,0]);
>sz:matrix([1,0],[0,-1]);
>
>create?a?matrix?with?shape
>
>
>[?sz??sx?]
>
>[?sx??sz?]
>
>matrix([1,0,0,1],[0,-1,1,0],[0,1,1,0],[1,0,0,-1])
>
>thanks
>
>
>Mauricio?Matera
>
>
>
>
>
>
>
>
>
>_______________________________________________
>Maxima?mailing?list
>Maxima at math.utexas.edu
>http://www.math.utexas.edu/mailman/listinfo/maxima