6 if(matA.
n!=matB.
m || matA.
n!=matB.
n){
8 std::cerr <<
"These two matrises can not make a summation." << std::endl
9 <<
"Your input was (" << matA.
n <<
"x" << matA.
n <<
") + (" << matB.
m <<
"x" << matB.
n <<
")." << std::endl;
16 const std::vector<dcomponent>::const_iterator matB_data_end =matB.
data.end();
17 for(std::vector<dcomponent>::const_iterator it=matB.
data.begin(); it!=matB_data_end; it++){
18 newmat(it->i,it->j) += it->v;
31 if(matA.
n!=matB.
m || matA.
n!=matB.
n){
33 std::cerr <<
"These two matrises can not make a subtraction." << std::endl
34 <<
"Your input was (" << matA.
n <<
"x" << matA.
n <<
") - (" << matB.
m <<
"x" << matB.
n <<
")." << std::endl;
41 const std::vector<dcomponent>::const_iterator matB_data_end =matB.
data.end();
42 for(std::vector<dcomponent>::const_iterator it=matB.
data.begin(); it!=matB_data_end; it++){
43 newmat(it->i,it->j) -= it->v;
58 std::cerr <<
"These two matrises can not make a product." << std::endl
59 <<
"Your input was (" << matA.
n <<
"x" << matA.
n <<
") * (" << matB.
m <<
"x" << matB.
n <<
")." << std::endl;
67 const std::vector<dcomponent>::const_iterator matB_data_end =matB.
data.end();
68 for(std::vector<dcomponent>::const_iterator it=matB.
data.begin(); it!=matB_data_end; it++){
69 for(CPPL_INT
i=0;
i<matA.
n;
i++){
70 newmat(
i,it->j) += matB(
i,it->
i)*it->v;
CPPL_INT n
matrix column size
_dgematrix i(const _dgbmatrix &mat)
Real Double-precision General Dence Matrix Class.
(DO NOT USE) Smart-temporary Real Double-precision General Sparse Matrix Class
friend _dgematrix i(const _dgematrix &)
CPPL_INT m
matrix row size
(DO NOT USE) Smart-temporary Real Double-precision General Dence Matrix Class
_dgematrix to_dgematrix() const
std::vector< dcomponent > data
matrix data
(DO NOT USE) Smart-temporary Real Double-precision Symmetric Matrix Class
CPPL_INT n
matrix column size
_dgematrix operator*(const _dsymatrix &matA, const _dgsmatrix &matB)
_dcovector _(dcovector &vec)
_dgematrix operator-(const _dsymatrix &matA, const _dgsmatrix &matB)
_dgematrix operator+(const _dsymatrix &matA, const _dgsmatrix &matB)