CPPLapack Code

//=============================================================================
//! (DO NOT USE) Complex-double Class for Hermitian matrices
class zhecomplex : public comple
{
public:
  ///////////////////////////////////////////////
  /////////////////// objects ///////////////////
  ///////////////////////////////////////////////
  long i, j;
  comple& v;
  
  ///////////////////////////////////////////////
  ///////////////// constructors ////////////////
  ///////////////////////////////////////////////
  inline zhecomplex(const long&, const long&, comple&);
  
  ///////////////////////////////////////////////
  ////////////////// functions //////////////////
  ///////////////////////////////////////////////
  inline zhecomplex& operator=(const comple&);
};

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

//=============================================================================
/*! constructor */
inline zhecomplex::zhecomplex(const long& _i, const long& _j, comple& _v)
  : comple( _i < _j ? std::conj( _v ) : _v ), 
    v( _v )
{CPPL_VERBOSE_REPORT;
  i = _i;
  j = _j;
}

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

//=============================================================================
/*! operator= */
inline zhecomplex& zhecomplex::operator=(const comple& _v)
{CPPL_VERBOSE_REPORT;
#ifdef  CPPL_DEBUG
  if( i==j && std::fabs(_v.imag()) > DBL_MIN ){
    WARNING_REPORT;
    std::cerr << "Diagonal components of a hermitian matrix have to be real numbers." << std::endl
              << "Your input to the (" << i << "," << j << ") element was a complex number, " << _v << "." << std::endl;
  }
#endif//CPPL_DEBUG
  
  comple::operator=( _v );
  v = ( i < j ? std::conj( _v ) : _v );
  return *this;
}