Go to the source code of this file.
Functions | |
_zgematrix | operator+ (const _zgematrix &matA, const zhematrix &matB) |
_zgematrix | operator- (const _zgematrix &matA, const zhematrix &matB) |
_zgematrix | operator* (const _zgematrix &matA, const zhematrix &matB) |
_zgematrix operator+ | ( | const _zgematrix & | matA, |
const zhematrix & | matB | ||
) | [inline] |
_zgematrix+zhematrix operator
Definition at line 3 of file _zgematrix-zssmatrix.hpp.
References zhematrix::array, _zgematrix::m, zhematrix::m, _zgematrix::n, and zhematrix::n.
{VERBOSE_REPORT; #ifdef CPPL_DEBUG if(matA.m!=matB.m || matA.n!=matB.n){ ERROR_REPORT; std::cerr << "These two matrises can not make a summation." << std::endl << "Your input was (" << matA.m << "x" << matA.n << ") + (" << matB.m << "x" << matB.n << ")." << std::endl; exit(1); } #endif//CPPL_DEBUG for(long c=0; c<matB.vol; c++){ matA(matB.indx[c],matB.jndx[c]) += matB.array[c]; } return matA; }
_zgematrix operator- | ( | const _zgematrix & | matA, |
const zhematrix & | matB | ||
) | [inline] |
_zgematrix-zhematrix operator
Definition at line 23 of file _zgematrix-zssmatrix.hpp.
References _zgematrix::array, zhematrix::array, i(), _zgematrix::m, zhematrix::m, _zgematrix::n, and zhematrix::n.
{VERBOSE_REPORT; #ifdef CPPL_DEBUG if(matA.m!=matB.m || matA.n!=matB.n){ ERROR_REPORT; std::cerr << "These two matrises can not make a subtraction." << std::endl << "Your input was (" << matA.m << "x" << matA.n << ") - (" << matB.n << "x" << matB.n << ")." << std::endl; exit(1); } #endif//CPPL_DEBUG //// change sign //// for(long i=0; i<matA.m*matA.n; i++){ matA.array[i]=-matA.array[i]; } //// add //// for(long c=0; c<matB.vol; c++){ matA(matB.indx[c],matB.jndx[c]) += matB.array[c]; } return matA; }
_zgematrix operator* | ( | const _zgematrix & | matA, |
const zhematrix & | matB | ||
) | [inline] |
_zgematrix*zhematrix operator
Definition at line 49 of file _zgematrix-zssmatrix.hpp.
References _(), zhematrix::array, _zgematrix::destroy(), i(), _zgematrix::m, zhematrix::m, _zgematrix::n, zhematrix::n, and zgematrix::zero().
{VERBOSE_REPORT; #ifdef CPPL_DEBUG if(matA.n!=matB.m){ ERROR_REPORT; std::cerr << "These two matrises can not make a product." << std::endl << "Your input was (" << matA.m << "x" << matA.n << ") * (" << matB.n << "x" << matB.n << ")." << std::endl; exit(1); } #endif//CPPL_DEBUG zgematrix newmat(matA.m, matB.n); newmat.zero(); for(long c=0; c<matB.vol; c++){ for(long i=0; i<matA.m; i++){ newmat(i,matB.jndx[c]) += matA(i,matB.indx[c])*matB.array[c]; } } matA.destroy(); return _(newmat); }