CPPLapack Code

//=============================================================================
/*! +dgsmatrix operator */
inline const dgsmatrix& operator+(const dgsmatrix& mat)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator+(const dgsmatrix&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  return mat;
}

//=============================================================================
/*! -dgsmatrix operator */
inline _dgsmatrix operator-(const dgsmatrix& mat)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator-(const dgsmatrix&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  dgsmatrix newmat(mat.M,mat.N,mat.CAP);
  
  newmat.VOL=mat.VOL;
  for(long i=0; i<mat.VOL; i++){
    newmat.Array[i] =-mat.Array[i];
    newmat.Indx[i] =mat.Indx[i];
    newmat.Jndx[i] =mat.Jndx[i];
  }
  for(long i=0; i<mat.M; i++){ newmat.Row[i] =mat.Row[i]; }
  for(long j=0; j<mat.N; j++){ newmat.Col[j] =mat.Col[j]; }
  
  return _(newmat);
}