|
template<matrix_order RO, matrix_style RS, typename T , matrix_order PO, matrix_style PS>
Cholesky decomposition of a symmetric positive-definite matrix.
This function performs Cholesky decomposition. That is, given a symmetric positive definite Matrix, , cholesky() returns a lower triangular Matrix such that . This function is faster than lu_decomp() and, therefore, preferable in cases where one's Matrix is symmetric positive definite.
- Parameters:
-
A | The symmetric positive definite Matrix to decompose. |
- See also:
- chol_solve(const Matrix<T,PO1,PS1> &, const Matrix<T,PO2,PS2> &)
-
chol_solve(const Matrix<T,PO1,PS1> &, const Matrix<T,PO2,PS2> &, const Matrix<T,PO3,PS3> &)
-
lu_decomp(Matrix<T,PO1,PS1>, Matrix<T,PO2,Concrete>&, Matrix<T,PO3,Concrete>&, Matrix<unsigned int, PO4, Concrete>&)
- Exceptions:
-
References scythe::Matrix_base< ORDER, STYLE >::cols(), scythe::Matrix_base< ORDER, STYLE >::isNull(), scythe::Matrix_base< ORDER, STYLE >::isSquare(), scythe::Matrix_base< ORDER, STYLE >::rows(), SCYTHE_CHECK_10, and sqrt().
Referenced by scythe::rng< mersenne >::rmvnorm().
|