CPPLapack Code

//=============================================================================
/*! dsymatrix*=double operator */
inline dsymatrix& dsymatrix::operator*=(const double& d)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dsymatrix::operator*=(const double&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  for(long i=0; i<N; i++){
    for(long j=0; j<=i; j++){
      Darray[j][i] *=d;
    }
  }
  
  return *this;
}

//=============================================================================
/*! dsymatrix/=double operator */
inline dsymatrix& dsymatrix::operator/=(const double& d)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] dsymatrix::operator/=(const double&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  for(long i=0; i<N; i++){
    for(long j=0; j<=i; j++){
      Darray[j][i] /=d;
    }
  }
  return *this;
}

///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////
///////////////////////////////////////////////////////////////////////////////

//=============================================================================
/*! dsymatrix*double operator */
inline _dsymatrix operator*(const dsymatrix& mat, const double& d)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator*(const dsymatrix&, const double&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  dsymatrix newmat(mat.N);
  for(long i=0; i<mat.N; i++){
    for(long j=0; j<=i; j++){
      newmat.Darray[j][i] =mat.Darray[j][i]*d;
    }
  }
  
  return _(newmat);
}

//=============================================================================
/*! dsymatrix/double operator */
inline _dsymatrix operator/(const dsymatrix& mat, const double& d)
{
#ifdef  CPPL_VERBOSE
  std::cerr << "# [MARK] operator/(const dsymatrix&, const double&)"
            << std::endl;
#endif//CPPL_VERBOSE
  
  dsymatrix newmat(mat.N);
  for(long i=0; i<mat.N; i++){
    for(long j=0; j<=i; j++){
      newmat.Darray[j][i] =mat.Darray[j][i]/d;
    }
  }
  
  return _(newmat);
}