00001
00003 inline _zhematrix t(const _zhematrix& mat)
00004 {
00005 #ifdef CPPL_VERBOSE
00006 std::cerr << "# [MARK] t(const _zhematrix&)"
00007 << std::endl;
00008 #endif//CPPL_VERBOSE
00009
00010 #ifdef CPPL_DEBUG
00011 std::cerr << "[WARNING] t(const zhematrix&)" << 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 _zgematrix i(const _zhematrix& mat)
00022 {
00023 #ifdef CPPL_VERBOSE
00024 std::cerr << "# [MARK] i(const _zhematrix&)"
00025 << std::endl;
00026 #endif//CPPL_VERBOSE
00027
00028 zhematrix mat_cp;
00029 mat_cp.shallow_copy(mat);
00030
00031 zgematrix mat_inv(mat.N,mat.N);
00032 mat_inv.identity();
00033
00034 mat_cp.zhesv(mat_inv);
00035
00036 return _(mat_inv);
00037 }