CPPLapack
 All Classes Files Functions Variables Friends Pages
_zgsmatrix-_zgematrix.hpp
Go to the documentation of this file.
1 //=============================================================================
2 /*! _zgsmatrix+_zgematrix operator */
3 inline _zgematrix operator+(const _zgsmatrix& matA, const _zgematrix& matB)
4 {CPPL_VERBOSE_REPORT;
5 #ifdef CPPL_DEBUG
6  if(matA.m!=matB.m || matA.n!=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.m << "x" << matB.n << ")." << std::endl;
10  exit(1);
11  }
12 #endif//CPPL_DEBUG
13 
14  const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
15  for(std::vector<zcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
16  matB(it->i,it->j) += it->v;
17  }
18 
19  matA.destroy();
20  return matB;
21 }
22 
23 //=============================================================================
24 /*! _zgsmatrix-_zgematrix operator */
25 inline _zgematrix operator-(const _zgsmatrix& matA, const _zgematrix& matB)
26 {CPPL_VERBOSE_REPORT;
27 #ifdef CPPL_DEBUG
28  if(matA.m!=matB.m || matA.n!=matB.n){
29  ERROR_REPORT;
30  std::cerr << "These two matrises can not make a subtraction." << 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  //// change sign ////
37  for(int i=0; i<matB.m*matB.n; i++){
38  matB.array[i] = -matB.array[i];
39  }
40 
41  //// add ////
42  const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
43  for(std::vector<zcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
44  matB(it->i,it->j) += it->v;
45  }
46 
47  matA.destroy();
48  return matB;
49 }
50 
51 //=============================================================================
52 /*! _zgsmatrix*_zgematrix operator */
53 inline _zgematrix operator*(const _zgsmatrix& matA, const _zgematrix& matB)
54 {CPPL_VERBOSE_REPORT;
55 #ifdef CPPL_DEBUG
56  if(matA.n!=matB.m){
57  ERROR_REPORT;
58  std::cerr << "These two matrises can not make a product." << std::endl
59  << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl;
60  exit(1);
61  }
62 #endif//CPPL_DEBUG
63 
64  zgematrix newmat(matA.m, matB.n);
65  newmat.zero();
66 
67  const std::vector<zcomponent>::const_iterator matA_data_end =matA.data.end();
68  for(std::vector<zcomponent>::const_iterator it=matA.data.begin(); it!=matA_data_end; it++){
69  for(CPPL_INT j=0; j<matB.n; j++){
70  newmat(it->i,j) += it->v *matB(it->j,j);
71  }
72  }
73 
74  matA.destroy();
75  matB.destroy();
76  return _(newmat);
77 }
CPPL_INT m
matrix row size
Definition: _zgematrix.hpp:9
_dgematrix i(const _dgbmatrix &mat)
CPPL_INT n
matrix column size
Definition: _zgsmatrix.hpp:10
(DO NOT USE) Smart-temporary Real Double-precision General Sparse Matrix Class
Definition: _zgsmatrix.hpp:3
Complex Double-precision General Dence Matrix Class.
Definition: zgematrix.hpp:3
zgematrix & zero()
(DO NOT USE) Smart-temporary Complex Double-precision General Dence Matrix Class
Definition: _zgematrix.hpp:3
void destroy() const
CPPL_INT m
matrix row size
Definition: _zgsmatrix.hpp:9
void destroy() const
comple * array
1D array to store matrix data
Definition: _zgematrix.hpp:11
_zgematrix operator+(const _zgsmatrix &matA, const _zgematrix &matB)
_zgematrix operator*(const _zgsmatrix &matA, const _zgematrix &matB)
_dcovector _(dcovector &vec)
CPPL_INT n
matrix column size
Definition: _zgematrix.hpp:10
std::vector< zcomponent > data
matrix data
Definition: _zgsmatrix.hpp:11
_zgematrix operator-(const _zgsmatrix &matA, const _zgematrix &matB)