CPPLapack Code

//=============================================================================
/*! double*dssmatrix operator */
inline _dssmatrix operator*(const double& d, const dssmatrix& mat)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator*(const double&, const dssmatrix&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  dssmatrix newmat(mat.N, mat.Data.size());
  for(std::vector<dcomponent>::const_iterator it=mat.Data.begin(); it!=mat.Data.end(); it++){
    newmat.fput(it->i, it->j, d*it->v);
  }
  return _(newmat);
}