//=============================================================================
/*! dssmatrix=_dssmatrix operator */
inline dssmatrix& dssmatrix::operator=(const _dssmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dssmatrix::operator=(const _dssmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
shallow_copy(mat);
return *this;
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//=============================================================================
/*! dssmatrix+=_dssmatrix operator */
inline dssmatrix& dssmatrix::operator+=(const _dssmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dssmatrix::operator+=(const _dssmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
#ifdef CPPL_DEBUG
if(N!=mat.N){
std::cerr << "[ERROR] dssmatrix::operator+=(const _dssmatrix&)"
<< std::endl
<< "These two matrises can not make a summation." << std::endl
<< "Your input was (" << N << "x" << N << ") += ("
<< mat.N << "x" << mat.N << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
for(long c=0; c<mat.VOL; c++){
(*this)(mat.Indx[c],mat.Jndx[c]) +=mat.Array[c];
}
mat.destroy();
return *this;
}
//=============================================================================
/*! dssmatrix-=_dssmatrix operator */
inline dssmatrix& dssmatrix::operator-=(const _dssmatrix& mat)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] dssmatrix::operator-=(const _dssmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
#ifdef CPPL_DEBUG
if(N!=mat.N){
std::cerr << "[ERROR] dssmatrix::operator-=(const _dssmatrix&)"
<< std::endl
<< "These two matrises can not make a sutraction." << std::endl
<< "Your input was (" << N << "x" << N << ") -= ("
<< mat.N << "x" << mat.N << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
for(long c=0; c<mat.VOL; c++){
(*this)(mat.Indx[c],mat.Jndx[c]) -=mat.Array[c];
}
mat.destroy();
return *this;
}
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
//=============================================================================
/*! dssmatrix+_dssmatrix operator */
inline _dssmatrix operator+(const dssmatrix& matA, const _dssmatrix& matB)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] operator+(const dssmatrix&, const _dssmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
#ifdef CPPL_DEBUG
if(matA.N!=matB.N){
std::cerr << "[ERROR] operator+(const dssmatrix&, const _dssmatrix&)"
<< std::endl
<< "These two matrises can not make a summation." << std::endl
<< "Your input was (" << matA.N << "x" << matA.N << ") + ("
<< matB.N << "x" << matB.N << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
dssmatrix newmat(matB);
for(long c=0; c<matA.VOL; c++){
newmat(matA.Indx[c], matA.Jndx[c]) +=matA.Array[c];
}
return _(newmat);
}
//=============================================================================
/*! dssmatrix-_dssmatrix operator */
inline _dssmatrix operator-(const dssmatrix& matA, const _dssmatrix& matB)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] operator-(const dssmatrix&, const _dssmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
#ifdef CPPL_DEBUG
if(matA.N!=matB.N){
std::cerr << "[ERROR] operator-(const dssmatrix&, const dssmatrix&)"
<< std::endl
<< "These two matrises can not make a subtraction." << std::endl
<< "Your input was (" << matA.N << "x" << matA.N << ") - ("
<< matB.N << "x" << matB.N << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
dssmatrix newmat(matB);
newmat.chsign();
//// add ////
for(long c=0; c<matA.VOL; c++){
newmat(matA.Indx[c], matA.Jndx[c]) +=matA.Array[c];
}
return _(newmat);
}
//=============================================================================
/*! dssmatrix*_dssmatrix operator */
/*
inline _dgsmatrix operator*(const dssmatrix& matA, const _dssmatrix& matB)
{
#ifdef CPPL_VERBOSE
std::cerr << "# [MARK] operator*(const dssmatrix&, const _dssmatrix&)"
<< std::endl;
#endif//CPPL_VERBOSE
#ifdef CPPL_DEBUG
if(matA.N!=matB.N){
std::cerr << "[ERROR] operator*(const dssmatrix&, const dssmatrix&)"
<< std::endl
<< "These two matrises can not make a product." << std::endl
<< "Your input was (" << matA.N << "x" << matA.N << ") * ("
<< matB.N << "x" << matB.N << ")." << std::endl;
exit(1);
}
#endif//CPPL_DEBUG
dssmatrix newmat( matA.N, 0 );
for(long c=0; c<matA.VOL; c++){
long k(matA.Jndx[c]);
std::vector< std::pair<long,long> >::iterator p;
for(p=matB.Line[k].begin(); p!=matB.Line[k].end(); p++){
newmat(matA.Indx[c],p->first) +=matA.Array[c]*matB.Array[p->second];
}
}
matB.destroy();
return _(newmat);
}
*/