Scythe-1.0.3
|
template<matrix_order RO, matrix_style RS, typename T , matrix_order PO1, matrix_style PS1, matrix_order PO2, matrix_style PS2, matrix_order PO3, matrix_style PS3>
Solve for x via backward substitution, given a lower triangular matrix resulting from Cholesky decomposition. This function solves the system of equations via backward substitution. L is the lower triangular matrix generated by Cholesky decomposition such that . This function is intended for repeatedly solving systems of equations based on A. That is A stays constant while b varies.
References scythe::Matrix_base< ORDER, STYLE >::isColVector(), scythe::Matrix_base< ORDER, STYLE >::isNull(), scythe::Matrix_base< ORDER, STYLE >::isSquare(), scythe::Matrix_base< ORDER, STYLE >::rows(), and SCYTHE_CHECK_10. |