//=============================================================================
/*! _zgsmatrix*std::complex<double> operator */
inline _zgsmatrix operator*(const _zgsmatrix& mat, const std::complex<double>& d)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator*(const _zgsmatrix&, const std::complex<double>&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  zscal_(mat.VOL, d, mat.Array, 1);
  return mat;
}
//=============================================================================
/*! _zgsmatrix/std::complex<double> operator */
inline _zgsmatrix operator/(const _zgsmatrix& mat, const std::complex<double>& d)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator/(const _zgsmatrix&, const std::complex<double>&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  zscal_(mat.VOL, 1./d, mat.Array, 1);
  return mat;
}