CPPLapack Code

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

//=============================================================================
/*! -zhsmatrix operator */
inline _zhsmatrix operator-(const zhsmatrix& mat)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator-(const zhsmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  zhsmatrix newmat(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.N; i++){
    newmat.Line[i] =mat.Line[i];
  }
  
  return _(newmat);
}