CPPLapack Code

//=============================================================================
/*! double*zgsmatrix operator */
inline _zgsmatrix operator*(const double& d, const zgsmatrix& mat)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator*(const double&, const zgsmatrix&)" << std::endl;
#endif//CPPL_VERBOSE
  
  zgsmatrix newmat(mat.M, mat.N, mat.CAP);
  for(long c=0; c<mat.VOL; c++){
    newmat.fput(mat.Indx[c], mat.Jndx[c], d*mat.Array[c]);
  }
  return _(newmat);
}