CPPLapack Code

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