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 since the matrix is symmetric."
00013 << std::endl;
00014 #endif//CPPL_DEBUG
00015
00016 return mat;
00017 }
00018
00019
00021 inline _dgematrix i(const _dsymatrix& mat)
00022 {
00023 #ifdef CPPL_VERBOSE
00024 std::cerr << "# [MARK] i(const _dsymatrix&)"
00025 << std::endl;
00026 #endif//CPPL_VERBOSE
00027
00028 dsymatrix mat_cp;
00029 mat_cp.shallow_copy(mat);
00030
00031 dgematrix mat_inv(mat.N,mat.N);
00032 mat_inv.identity();
00033
00034 mat_cp.dsysv(mat_inv);
00035
00036 return _(mat_inv);
00037 }