VERB_code_2.3
_dsymatrix-calc.hpp
1 //=============================================================================
3 inline _dsymatrix t(const _dsymatrix& mat)
4 {
5 #ifdef CPPL_VERBOSE
6  std::cerr << "# [MARK] t(const _dsymatrix&)"
7  << std::endl;
8 #endif//CPPL_VERBOSE
9 
10 #ifdef CPPL_DEBUG
11  std::cerr << "[WARNING] t(const dsymatrix&)" << std::endl
12  << "This function call has no effect since the matrix is symmetric."
13  << std::endl;
14 #endif//CPPL_DEBUG
15 
16  return mat;
17 }
18 
19 //=============================================================================
21 inline _dgematrix i(const _dsymatrix& mat)
22 {
23 #ifdef CPPL_VERBOSE
24  std::cerr << "# [MARK] i(const _dsymatrix&)"
25  << std::endl;
26 #endif//CPPL_VERBOSE
27 
28  dsymatrix mat_cp;
29  mat_cp.shallow_copy(mat);
30 
31  dgematrix mat_inv(mat.N,mat.N);
32  mat_inv.identity();
33 
34  mat_cp.dsysv(mat_inv);
35 
36  return _(mat_inv);
37 }
long dsysv(dgematrix &)
Definition: dsymatrix-lapack.hpp:6
Real Double-precision General Dence Matrix Class.
Definition: dgematrix.hpp:3
Real Double-precision Symmetric Matrix Class [L-type (UPLO=L) Strage].
Definition: dsymatrix.hpp:3
(DO NOT USE) Smart-temporary Real Double-precision General Dence Matrix Class
Definition: _dgematrix.hpp:3
void identity()
Definition: dgematrix-misc.hpp:38
friend _dcovector t(const drovector &)
Definition: drovector-calc.hpp:3
void shallow_copy(const _dsymatrix &)
Definition: dsymatrix-misc.hpp:108
long N
matrix column or row size
Definition: _dsymatrix.hpp:7
(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class
Definition: _dsymatrix.hpp:3