Re: manipulating block diagonal matrices
- To: mathgroup at smc.vnet.net
- Subject: [mg37594] Re: manipulating block diagonal matrices
- From: "Lawrence A. Walker Jr." <lwalker701 at earthlink.net>
- Date: Wed, 6 Nov 2002 06:53:14 -0500 (EST)
- Organization: Morgan State University, COMSARE
- References: <aq876p$e6n$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
Hi David, I suggest the following. 1) The SparseLinearSolve package. 2) Finding C code on the web at http://cm.bell-labs.com/netlib/master/readme.html, LINPack or LaPack and then using Mathlink to link these to Mathematica. Hope this helps. Lawrence David E. Burmaster wrote: > Dear MathGroup, > > Can anyone suggest efficient algorithms for calculating either or both of > these? > > (i) the inverse and (ii) the determinant of a block diagonal matrix? > > Let blocki for i = 1, 2,..., n denote each of the square matrices along the > diagonal of a large block diagonal matrix. > > In my situation, the blocki matrices do not have a common size; a typical > blocki has dimensions in the range 6x6 to 20x20. > > In my situation, n is approx 1,000. > > In other words, the full block diagonal matrix has large dimensions -- but > it is sparse in a highly structured way. > > many thanks for your help > > dave > > > > > +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ > David E. Burmaster, Ph.D. > Alceon Corporation > POBox 382069 (new Box number effective 1 Sep 2001) > Harvard Square Station > Cambridge, MA 02238-2069 (new ZIP code effective 1 Sep 2001) > > Voice 617-864-4300 > > Web http://www.Alceon.com > Email deb at Alceon.com > +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ > > > -- Lawrence A. Walker Jr. http://www.kingshonor.com