CPPLapack
 All Classes Files Functions Variables Friends Pages
_zgematrix-zhematrix.hpp
Go to the documentation of this file.
1 //=============================================================================
2 /*! _zgematrix+zhematrix operator */
3 inline _zgematrix operator+(const _zgematrix& matA, const zhematrix& matB)
4 {CPPL_VERBOSE_REPORT;
5 #ifdef CPPL_DEBUG
6  if(matA.n!=matB.n || matA.m!=matB.n){
7  ERROR_REPORT;
8  std::cerr << "These two matrises can not make a summation." << std::endl
9  << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
10  exit(1);
11  }
12 #endif//CPPL_DEBUG
13 
14  for(CPPL_INT i=0; i<matB.n; i++){
15  for(CPPL_INT j=0; j<matB.n; j++){
16  matA(i,j) += matB(i,j);
17  }
18  }
19 
20  return matA;
21 }
22 
23 //=============================================================================
24 /*! _zgematrix-zhematrix operator */
25 inline _zgematrix operator-(const _zgematrix& matA, const zhematrix& matB)
26 {CPPL_VERBOSE_REPORT;
27 #ifdef CPPL_DEBUG
28  if(matA.n!=matB.n || matA.m!=matB.n){
29  ERROR_REPORT;
30  std::cerr << "These two matrises can not make a summation." << std::endl
31  << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.n << "x" << matB.n << ")." << std::endl;
32  exit(1);
33  }
34 #endif//CPPL_DEBUG
35 
36  for(CPPL_INT i=0; i<matB.n; i++){
37  for(CPPL_INT j=0; j<matB.n; j++){
38  matA(i,j) -= matB(i,j);
39  }
40  }
41 
42  return matA;
43 }
44 
45 //=============================================================================
46 /*! _zgematrix*zhematrix operator */
47 inline _zgematrix operator*(const _zgematrix& matA, const zhematrix& matB)
48 {CPPL_VERBOSE_REPORT;
49 #ifdef CPPL_DEBUG
50  if(matA.n!=matB.n){
51  ERROR_REPORT;
52  std::cerr << "These two matrises can not make a product." << std::endl
53  << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl;
54  exit(1);
55  }
56 #endif//CPPL_DEBUG
57 
58  zgematrix newmat( matA.m, matB.n );
59  char side ='R';
60  char uplo ='l';
61  comple alpha =comple(1.,0.);
62  comple beta =comple(0.,0.);
63 
64  zhemm_( &side, &uplo, &newmat.m, &newmat.n, &alpha, matB.array, &matB.n, matA.array, &matA.m, &beta, newmat.array, &newmat.m );
65 
66  matA.destroy();
67  return _(newmat);
68 }
CPPL_INT m
matrix row size
Definition: _zgematrix.hpp:9
_zgematrix operator+(const _zgematrix &matA, const zhematrix &matB)
_dgematrix i(const _dgbmatrix &mat)
CPPL_INT n
matrix column size
Definition: zgematrix.hpp:10
Complex Double-precision General Dence Matrix Class.
Definition: zgematrix.hpp:3
CPPL_INT n
matrix column size
Definition: zhematrix.hpp:11
comple * array
1D array to store matrix data
Definition: zgematrix.hpp:11
(DO NOT USE) Smart-temporary Complex Double-precision General Dence Matrix Class
Definition: _zgematrix.hpp:3
comple * array
1D array to store matrix data
Definition: zhematrix.hpp:12
void destroy() const
_zgematrix operator*(const _zgematrix &matA, const zhematrix &matB)
comple * array
1D array to store matrix data
Definition: _zgematrix.hpp:11
CPPL_INT m
matrix row size
Definition: zgematrix.hpp:9
Complex Double-precision Hermitian Matrix Class [l-type (UPLO=l) Strage].
Definition: zhematrix.hpp:4
_zgematrix operator-(const _zgematrix &matA, const zhematrix &matB)
_dcovector _(dcovector &vec)
CPPL_INT n
matrix column size
Definition: _zgematrix.hpp:10