fn_chol.hpp
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021 template<typename T1>
00022 inline
00023 bool
00024 chol(Mat<typename T1::elem_type>& out, const Base<typename T1::elem_type,T1>& X)
00025 {
00026 arma_extra_debug_sigprint();
00027
00028 typedef typename T1::elem_type eT;
00029
00030 const unwrap_check<T1> tmp(X.get_ref(), out);
00031 const Mat<eT>& A = tmp.M;
00032
00033 arma_debug_check( (A.is_square() == false), "chol(): given matrix is not square");
00034
00035 return auxlib::chol(out, A);
00036 }
00037
00038
00039
00040 template<typename T1>
00041 inline
00042 Mat<typename T1::elem_type>
00043 chol(const Base<typename T1::elem_type,T1>& X)
00044 {
00045 arma_extra_debug_sigprint();
00046
00047 typedef typename T1::elem_type eT;
00048
00049 Mat<eT> out;
00050
00051 const bool ok = chol(out, X);
00052
00053 if(ok == false)
00054 {
00055 arma_print("chol(): matrix factorisation failed");
00056 }
00057
00058 return out;
00059 }
00060
00061
00062
00063