6 std::cerr <<
"# [MARK] dgematrix::operator=(const dsymatrix&)"
11 for(
long i=0;
i<mat.N;
i++){
12 for(
long j=0; j<mat.N; j++){
29 std::cerr <<
"# [MARK] dgematrix::operator+=(const dsymatrix&)"
34 if(N!=mat.N || M!=mat.N){
35 std::cerr <<
"[ERROR] dgematrix::operator+=(dsymatrix&)" << std::endl
36 <<
"These two matrises can not make a summation." << std::endl
37 <<
"Your input was (" << M <<
"x" << N <<
") += ("
38 << mat.N <<
"x" << mat.N <<
")." << std::endl;
43 for(
long i=0;
i<M;
i++){
44 for(
long j=0; j<N; j++){
57 std::cerr <<
"# [MARK] dgematrix::operator-=(const dsymatrix&)"
62 if(N!=mat.N || M!=mat.N){
63 std::cerr <<
"[ERROR] dgematrix::operator-=(dsymatrix&)" << std::endl
64 <<
"These two matrises can not make a sutraction." << std::endl
65 <<
"Your input was (" << M <<
"x" << N <<
") -= ("
66 << mat.N <<
"x" << mat.N <<
")." << std::endl;
71 for(
long i=0;
i<M;
i++){
72 for(
long j=0; j<N; j++){
85 std::cerr <<
"# [MARK] dgematrix::operator*=(const dsymatrix&)"
91 std::cerr <<
"[ERROR] dgematrix::operator*=(dsymatrix&)" << std::endl
92 <<
"These two matrises can not make a product." << std::endl
93 <<
"Your input was (" << M <<
"x" << N <<
") *= ("
94 << mat.N <<
"x" << mat.N <<
")." << std::endl;
100 dsymm_(
'R',
'L', mat.N, N, 1.0, mat.Array, mat.N,
101 Array, M, 0.0, newmat.Array, newmat.M );
116 std::cerr <<
"# [MARK] operator+(const dgematrix&, const dsymatrix&)"
121 if(matA.N!=matB.N || matA.M!=matB.N){
122 std::cerr <<
"[ERROR] operator+(const dgematrix&, const dsymatrix&)"
124 <<
"These two matrises can not make a summation." << std::endl
125 <<
"Your input was (" << matA.M <<
"x" << matA.N <<
") + ("
126 << matB.N <<
"x" << matB.N <<
")." << std::endl;
132 for(
long i=0; i<matA.M; i++){
133 for(
long j=0; j<matA.N; j++){
134 newmat(i,j) += matB(i,j);
146 std::cerr <<
"# [MARK] operator-(const dgematrix&, const dsymatrix&)"
151 if(matA.N!=matB.N || matA.M!=matB.N){
152 std::cerr <<
"[ERROR] operator-(const dgematrix&, const dsymatrix&)"
154 <<
"These two matrises can not make a subtraction." << std::endl
155 <<
"Your input was (" << matA.M <<
"x" << matA.N <<
") - ("
156 << matB.N <<
"x" << matB.N <<
")." << std::endl;
162 for(
long i=0; i<matA.M; i++){
163 for(
long j=0; j<matA.N; j++){
164 newmat(i,j) -= matB(i,j);
176 std::cerr <<
"# [MARK] operator*(const dgematrix&, const dsymatrix&)"
182 std::cerr <<
"[ERROR] operator*(const dgematrix&, const dsymatrix&)"
184 <<
"These two matrises can not make a product." << std::endl
185 <<
"Your input was (" << matA.M <<
"x" << matA.N <<
") * ("
186 << matB.N <<
"x" << matB.N <<
")." << std::endl;
192 dsymm_(
'R',
'L', matB.N, matA.N, 1.0, matB.Array, matB.N,
193 matA.Array, matA.M, 0.0, newmat.Array, newmat.M );
void resize(const long &, const long &)
Definition: dgematrix-misc.hpp:126
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
Real Double-precision Symmetric Matrix Class [L-type (UPLO=L) Strage].
Definition: dsymatrix.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
dgematrix & operator-=(const dgematrix &)
Definition: dgematrix-dgematrix.hpp:45
friend _dgematrix i(const dgematrix &)
Definition: dgematrix-calc.hpp:21
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