6 std::cerr <<
"# [MARK] dgematrix::operator=(const _dsymatrix&)"
32 std::cerr <<
"# [MARK] dgematrix::operator+=(const _dsymatrix&)"
37 if(N!=mat.
N || M!=mat.
N){
38 std::cerr <<
"[ERROR] dgematrix::operator+=(_dsymatrix&)" << std::endl
39 <<
"These two matrises can not make a summation." << std::endl
40 <<
"Your input was (" << M <<
"x" << N <<
") += ("
41 << mat.
N <<
"x" << mat.
N <<
")." << std::endl;
46 for(
long i=0;
i<M;
i++){
47 for(
long j=0; j<N; j++){
61 std::cerr <<
"# [MARK] dgematrix::operator-=(const _dsymatrix&)"
66 if(N!=mat.
N || M!=mat.
N){
67 std::cerr <<
"[ERROR] dgematrix::operator-=(_dsymatrix&)" << std::endl
68 <<
"These two matrises can not make a sutraction." << std::endl
69 <<
"Your input was (" << M <<
"x" << N <<
") -= ("
70 << mat.
N <<
"x" << mat.
N <<
")." << std::endl;
75 for(
long i=0;
i<M;
i++){
76 for(
long j=0; j<N; j++){
90 std::cerr <<
"# [MARK] dgematrix::operator*=(const _dsymatrix&)"
96 std::cerr <<
"[ERROR] dgematrix::operator*=(_dsymatrix&)" << std::endl
97 <<
"These two matrises can not make a product." << std::endl
98 <<
"Your input was (" << M <<
"x" << N <<
") *= ("
99 << mat.
N <<
"x" << mat.
N <<
")." << std::endl;
105 dsymm_(
'R',
'L', mat.
N, N, 1.0, mat.
Array, mat.
N,
106 Array, M, 0.0, newmat.Array, newmat.M );
122 std::cerr <<
"# [MARK] operator+(const dgematrix&, const _dsymatrix&)"
127 if(matA.N!=matB.
N || matA.M!=matB.
N){
128 std::cerr <<
"[ERROR] operator+(dgematrix&, _dgematrix&)" << std::endl
129 <<
"These two matrises can not make a summation." << std::endl
130 <<
"Your input was (" << matA.M <<
"x" << matA.N <<
") + ("
131 << matB.
N <<
"x" << matB.
N <<
")." << std::endl;
137 for(
long i=0; i<matA.M; i++){
138 for(
long j=0; j<matA.N; j++){
139 newmat(i,j) += matB(i,j);
152 std::cerr <<
"# [MARK] operator-(const dgematrix&, const _dsymatrix&)"
157 if(matA.N!=matB.
N || matA.M!=matB.
N){
158 std::cerr <<
"[ERROR] operator-(dgematrix&, _dgematrix&)" << std::endl
159 <<
"These two matrises can not make a subtraction." << std::endl
160 <<
"Your input was (" << matA.M <<
"x" << matA.N <<
") - ("
161 << matB.
N <<
"x" << matB.
N <<
")." << std::endl;
167 for(
long i=0; i<matA.M; i++){
168 for(
long j=0; j<matA.N; j++){
169 newmat(i,j) -= matB(i,j);
182 std::cerr <<
"# [MARK] operator*(const dgematrix&, const _dsymatrix&)"
188 std::cerr <<
"[ERROR] operator*(dgematrix&, _dgematrix&)" << std::endl
189 <<
"These two matrises can not make a product." << std::endl
190 <<
"Your input was (" << matA.M <<
"x" << matA.N <<
") * ("
191 << matB.
N <<
"x" << matB.
N <<
")." << std::endl;
197 dsymm_(
'R',
'L', matB.
N, matA.N, 1.0, matB.
Array, matB.
N,
198 matA.Array, matA.M, 0.0, newmat.Array, newmat.M );
void clear()
Definition: dgematrix-misc.hpp:3
double * Array
1D Array to store matrix data
Definition: _dsymatrix.hpp:8
void destroy() const
Definition: _dsymatrix-misc.hpp:3
friend void swap(dgematrix &, dgematrix &)
Definition: dgematrix-misc.hpp:154
Real Double-precision General Dence Matrix Class.
Definition: dgematrix.hpp:3
dgematrix & operator=(const dgematrix &)
Definition: dgematrix-dgematrix.hpp:3
dgematrix & operator*=(const dgematrix &)
Definition: dgematrix-dgematrix.hpp:68
friend _drovector operator-(const _drovector &)
Definition: _drovector-unary.hpp:15
(DO NOT USE) Smart-temporary Real Double-precision General Dence Matrix Class
Definition: _dgematrix.hpp:3
friend _drovector operator*(const drovector &, const dgematrix &)
Definition: drovector-dgematrix.hpp:3
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
void complete() const
Definition: _dsymatrix-misc.hpp:22
dgematrix & operator-=(const dgematrix &)
Definition: dgematrix-dgematrix.hpp:45
friend _dgematrix i(const dgematrix &)
Definition: dgematrix-calc.hpp:21
double ** Darray
array of pointers of column head addresses
Definition: _dsymatrix.hpp:9
double * Array
1D Array to store vector data
Definition: _drovector.hpp:8
dgematrix & operator+=(const dgematrix &)
Definition: dgematrix-dgematrix.hpp:22
friend const _drovector & operator+(const _drovector &)
Definition: _drovector-unary.hpp:3
double & operator()(const long &, const long &)
Definition: dgematrix-io.hpp:3