00001
00003 inline _dsymatrix t(const dsymatrix& mat)
00004 {
00005 #ifdef CPPL_VERBOSE
00006 std::cerr << "# [MARK] t(const dsymatrix&)"
00007 << std::endl;
00008 #endif//CPPL_VERBOSE
00009
00010 #ifdef CPPL_DEBUG
00011 std::cerr << "[WARNING] t(const dsymatrix&)" << std::endl
00012 << "This function call has no effect "
00013 << "since the matrix is symmetric." << std::endl;
00014 #endif//CPPL_DEBUG
00015
00016 dsymatrix newmat(mat);
00017 return _(newmat);
00018 }
00019
00020
00022 inline _dgematrix i(const dsymatrix& mat)
00023 {
00024 #ifdef CPPL_VERBOSE
00025 std::cerr << "# [MARK] i(const dsymatrix&)"
00026 << std::endl;
00027 #endif//CPPL_VERBOSE
00028
00029 dsymatrix mat_cp(mat);
00030 dgematrix mat_inv(mat.N,mat.N);
00031 mat_inv.identity();
00032 mat_cp.dsysv(mat_inv);
00033
00034 return _(mat_inv);
00035 }