[SciPy-User] cholesky of sparse symmetric banded matrix
josef.pktd at gmail.com
josef.pktd at gmail.com
Tue Aug 21 19:16:29 EDT 2012
I would like to get a cholesky decomposition of a symmetric banded matrix
and multiply it with a dense array
I found
http://docs.scipy.org/doc/scipy/reference/generated/scipy.linalg.cholesky_banded.html
and this
http://docs.scipy.org/doc/scipy/reference/generated/scipy.sparse.dia_matrix.html
either
1) is there a way to do linear algebra (dot multiplication) directly in the
"upper diagonal ordered form"?
or
2) is there an efficient way to go from the "upper diagonal ordered form"
to a sparse diagonal matrix (or both ways)?
Is there code that uses this and that I can look at for the pattern?
-------
my problem is standard linear least squares, where I have an explicit
banded form for the (nobs, nobs) weighting matrix
X'WX and X'Wy
and I need a transformation X2 = W^(0.5) X and y2 = W^(0.5) y
so I get X2'X2 and X2'y2
(nobs: number of observations, prime is transpose)
My first example only has one upper and one lower off-diagonal, so working
with dense is wasteful.
Josef
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.scipy.org/pipermail/scipy-user/attachments/20120821/52e0272b/attachment.html>
More information about the SciPy-User
mailing list