Ask Your Question

Revision history [back]

click to hide/show revision 1
initial version

L is the inverse of A (using Cholesky Decomposition), not the actual Cholesky Decomposition. To confirm this do:

cout << LA << endl; // this is inv(A)A.

You'll get the identity matrix (or very close to it).

L is the inverse of A (using Cholesky Decomposition), not the actual Cholesky Decomposition. To confirm this do:

cout << LL * A << endl; // this is inv(A)inv(A) * A.

You'll get the identity matrix (or very close to it).