Re: manipulating block diagonal matrices
- To: mathgroup at smc.vnet.net
- Subject: [mg37604] Re: manipulating block diagonal matrices
- From: Jens-Peer Kuska <kuska at informatik.uni-leipzig.de>
- Date: Wed, 6 Nov 2002 06:53:58 -0500 (EST)
- Organization: Universitaet Leipzig
- References: <aq876p$e6n$1@smc.vnet.net>
- Reply-to: kuska at informatik.uni-leipzig.de
- Sender: owner-wri-mathgroup at wolfram.com
Hi, no special blockdiagonal code but Developer`SparseLinearSolve[] may help you. Regards Jens "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 > +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++