Re: MATLAB COMPATIBILITY (sparse matrices)

[Ole Myren R|hne]

Sorry for the empty message:-|


>>>>> "Adler" == Andy Adler <address@hidden> writes:

    Adler> X=diag_inv(A,B) A is mxt where t+1 is number of
    Adler> non-zero superdiagonals B is mxn X is mxn

LAPACK handles band matrices:

> DPBTRF (l) - compute the Cholesky factorization of a real symmetric
> positive definite band matrix A
>
> DPBTRS (l) - solve a system of linear equations A*X = B with a
> symmetric positive definite band matrix A using the Cholesky
> factorization A = U**T*U or A = L*L**T computed by DPBTRF
>
> (and many more)

These routines should be well tested & optimized. As Octave already
uses LAPACK for linear algebra one might want to use LAPACK's band
matrix routines as well.

But what about matrices that are sparse but not narrow banded?

        Ole Myren Røhne