Scythe-1.0.3
template<matrix_order RO, matrix_style RS, typename T , matrix_order PO, matrix_style PS>
Matrix<T, RO, RS> scythe::cholesky ( const Matrix< T, PO, PS > &  A)

Cholesky decomposition of a symmetric positive-definite matrix.

This function performs Cholesky decomposition. That is, given a symmetric positive definite Matrix, $A$, cholesky() returns a lower triangular Matrix $L$ such that $A = LL^T$. This function is faster than lu_decomp() and, therefore, preferable in cases where one's Matrix is symmetric positive definite.

Parameters:
AThe 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:
scythe_alloc_error(Level 1)
scythe_dimension_error(Level 1)
scythe_null_error(Level 1)
scythe_type_error(Level 2)
scythe_alloc_error(Level 1)

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().