[SciPy-user] solving Ax=b, known L: A = L L^T

dmitrey openopt at ukr.net
Tue Aug 28 03:25:12 EDT 2007


Hi all,
suppose I have a system Ax=b and I know (Cholesky) triangular L: A = L L^T
what's the simplest way to obtain x now?
(I meant python code)
is it solve(L^T, solve(L,b))?
or I should somehow inform solve() about L is triangular or, maybe, use 
other func?
Documentation of numpy.linalg.solve doesn't say anything.
Afaik MATLAB's '\' determ type of matrix automatically.
Regards, D.




More information about the SciPy-User mailing list