6 std::cerr <<
"# [MARK] operator+(const _zgematrix&, const zgematrix&)"
11 if(matA.
N!=matB.N || matA.
M!=matB.M){
12 std::cerr <<
"[ERROR] operator+(_zgematrix&, zgematrix&)" << std::endl
13 <<
"These two matrises can not make a summation." << std::endl
14 <<
"Your input was (" << matA.
M <<
"x" << matA.
N <<
") + ("
15 << matB.M <<
"x" << matB.N <<
")." << std::endl;
20 for(
long i=0; i<matA.
M*matA.
N; i++){ matA.
Array[i]+=matB.Array[i]; }
30 std::cerr <<
"# [MARK] operator-(const _zgematrix&, const zgematrix&)"
35 if(matA.
N!=matB.N || matA.
M!=matB.M){
36 std::cerr <<
"[ERROR] operator-(_zgematrix&, zgematrix&)" << std::endl
37 <<
"These two matrises can not make a subtraction." << std::endl
38 <<
"Your input was (" << matA.
M <<
"x" << matA.
N <<
") - ("
39 << matB.M <<
"x" << matB.N <<
")." << std::endl;
44 for(
long i=0; i<matA.
M*matA.
N; i++){ matA.
Array[i]-=matB.Array[i]; }
54 std::cerr <<
"# [MARK] operator*(const _zgematrix&, const zgematrix&)"
60 std::cerr <<
"[ERROR] operator*(_zgematrix&, zgematrix&)" << std::endl
61 <<
"These two matrises can not make a product." << std::endl
62 <<
"Your input was (" << matA.
M <<
"x" << matA.
N <<
") * ("
63 << matB.M <<
"x" << matB.N <<
")." << std::endl;
69 zgemm_(
'N',
'N', matA.
M, matB.N, matA.
N, std::complex<double>(1.0,0.0),
70 matA.
Array, matA.
M, matB.Array, matB.M,
71 std::complex<double>(0.0,0.0), newmat.array, matA.
M );
long N
matrix column size
Definition: _zgematrix.hpp:8
friend _zrovector operator-(const _zrovector &)
Definition: _zrovector-unary.hpp:15
Complex Double-precision General Dence Matrix Class.
Definition: zgematrix.hpp:3
(DO NOT USE) Smart-temporary Complex Double-precision General Dence Matrix Class
Definition: _zgematrix.hpp:3
void destroy() const
Definition: _zgematrix-misc.hpp:3
std::complex< double > * Array
1D Array to store matrix data
Definition: _zgematrix.hpp:9
friend _zrovector operator*(const zrovector &, const zgematrix &)
Definition: zrovector-zgematrix.hpp:3
long M
matrix row size
Definition: _zgematrix.hpp:7
friend const _zrovector & operator+(const _zrovector &)
Definition: _zrovector-unary.hpp:3