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

dmitrey openopt@ukr....
Tue Aug 28 02:25:12 CDT 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