sparce matrices

[Mike Miller]

It would be great if I could specify the non-zero elements of a big block-
diagonal matrix, call it X, say, then do things like chol(X)*Y without
having to allocate a lot of memory.  I'm dealing with situations where I
might have, say, 1,000 8x8 blocks.  That means there would be 64,000
non-zero elements in a matrix with 64,000,000 elements.

Can octave deal with this kind of thing efficiently?  I've been dealing
with it by making 1,000 distinct 8x8 matrices, doing 1,000 Cholesky
factorizations, multiplying the 1,000 resulting matrices by 1,000 pieces
of Y and then assembling the products into one big matrix.

Mike



-------------------------------------------------------------
Octave is freely available under the terms of the GNU GPL.

Octave's home on the web:  http://www.octave.org
How to fund new projects:  http://www.octave.org/funding.html
Subscription information:  http://www.octave.org/archive.html
-------------------------------------------------------------