dgematrix.hpp

Go to the documentation of this file.

//=============================================================================
//! double-precision general matrix class
class dgematrix
{
private:
  //// objects ////
  long M; //!< global matrix row size (NOT accessible)
  long N; //!< global matrix column size (NOT accessible)
  long Ml; //!< local matrix row size (NOT accessible)
  long Nl; //!< local matrix column size (NOT accessible)
  long Desc[9]; //!< matrix description (NOT accessible)
  double* Array; //!< 1D array to store matrix data (NOT accessible)
  
public:
  //// const reference ////
  long const& m; //!< global matrix row size (readable)
  long const& n; //!< global matrix column size (readable)
  long const& ml; //!< local matrix row size (readable)
  long const& nl; //!< local matrix column size (readable)
  long* const& desc; //!< matrix description (readable)
  double* const& array; //!< 1D array to store matrix data (readable)
  
  //// constructor ////
  inline dgematrix();
  inline dgematrix(const dgematrix&);
  inline dgematrix(const _dgematrix&);
  //inline dgematrix(const dgbmatrix&);
  //inline dgematrix(const _dgbmatrix&);
  inline dgematrix(const long&, const long&);
  inline ~dgematrix(); //destructor
  
  //// io ////
  inline CPPSL_double operator()(const long&, const long&);
  inline double operator()(const long&, const long&) const;
  inline void set(const long&, const long&, const double&) const;
  inline friend std::ostream& operator<<(std::ostream&, const dgematrix&);
  
  //// misc ////
  inline void clear();
  inline void zero();
  inline void identity();
  inline void chsign();
  inline void copy(const dgematrix&);
  inline void shallow_copy(const _dgematrix&);
  inline void resize(const long&, const long&);
  inline friend void swap(dgematrix&, dgematrix&);
  
  //// calc ////
  inline void transpose();
  inline friend _dgematrix t(const dgematrix&);
  //inline friend _dgematrix i(dgematrix);
  //inline friend void idamax(long&, long&, const dgematrix&);
  //inline friend double damax(const dgematrix&);

  //// lapack ////
  inline long dgesv(dgematrix&);
  inline long dgesv(dcovector&);
  //inline long dgels(dgematrix&);
  //inline long dgels(dcovector&);
  //inline long dgelss(dcovector&, dcovector&, long, double);
  //inline long dgeev(dcovector&, dcovector&, dgematrix&);
  //inline long dgeev(drovector&, drovector&, dgematrix&);
  //inline long dgesvd(dcovector&, dgematrix&, dgematrix&);
  
  //// unary ////
  inline friend const dgematrix& operator+(const dgematrix&);
  inline friend _dgematrix operator-(const dgematrix&);
  
  //// operators with dgematrix ////
  inline dgematrix& operator=(const dgematrix&);
  inline dgematrix& operator+=(const dgematrix&);
  inline dgematrix& operator-=(const dgematrix&);
  inline dgematrix& operator*=(const dgematrix&);
  //// operators with _dgematrix ////
  inline dgematrix& operator=(const _dgematrix&);
  inline dgematrix& operator+=(const _dgematrix&);
  inline dgematrix& operator-=(const _dgematrix&);
  inline dgematrix& operator*=(const _dgematrix&);
  //// operators with dgbmatrix ////
  //inline dgematrix& operator=(const dgbmatrix&);
  //inline dgematrix& operator+=(const dgbmatrix&);
  //inline dgematrix& operator-=(const dgbmatrix&);
  //inline dgematrix& operator*=(const dgbmatrix&);
  //// operators with _dgbmatrix ////
  //inline dgematrix& operator=(const _dgbmatrix&);
  //inline dgematrix& operator+=(const _dgbmatrix&);
  //inline dgematrix& operator-=(const _dgbmatrix&);
  //inline dgematrix& operator*=(const _dgbmatrix&);
  //// operators with double ////
  inline dgematrix& operator*=(const double&);
  inline dgematrix& operator/=(const double&);

  //// + ////
  inline friend _dgematrix operator+(const  dgematrix&, const  dgematrix&);
  inline friend _dgematrix operator+(const  dgematrix&, const _dgematrix&);
  inline friend _dgematrix operator+(const _dgematrix&, const  dgematrix&);
  
  //inline friend _dgematrix operator+(const  dgematrix&, const  dgbmatrix&);
  //inline friend _dgematrix operator+(const  dgematrix&, const _dgbmatrix&);
  //inline friend _dgematrix operator+(const  dgbmatrix&, const  dgematrix&);
  //inline friend _dgematrix operator+(const _dgbmatrix&, const  dgematrix&);
  
  //// - ////
  inline friend _dgematrix operator-(const  dgematrix&, const  dgematrix&);
  inline friend _dgematrix operator-(const  dgematrix&, const _dgematrix&);
  inline friend _dgematrix operator-(const _dgematrix&, const  dgematrix&);
  
  //inline friend _dgematrix operator-(const  dgematrix&, const  dgbmatrix&);
  //inline friend _dgematrix operator-(const  dgematrix&, const _dgbmatrix&);
  //inline friend _dgematrix operator-(const  dgbmatrix&, const  dgematrix&);
  //inline friend _dgematrix operator-(const _dgbmatrix&, const  dgematrix&);
  
  //// * ////
  inline friend _dgematrix operator*(const  dgematrix&, const  dgematrix&);
  inline friend _dgematrix operator*(const  dgematrix&, const _dgematrix&);
  inline friend _dgematrix operator*(const _dgematrix&, const  dgematrix&);
  
  //inline friend _dgematrix operator*(const  dgematrix&, const  dgbmatrix&);
  //inline friend _dgematrix operator*(const  dgematrix&, const _dgbmatrix&);  
  //inline friend _dgematrix operator*(const  dgbmatrix&, const  dgematrix&);
  //inline friend _dgematrix operator*(const _dgbmatrix&, const  dgematrix&);
  
  inline friend _dcovector operator*(const  dgematrix&, const  dcovector&);
  inline friend _dcovector operator*(const  dgematrix&, const _dcovector&);
  
  inline friend _drovector operator*(const  drovector&, const  dgematrix&);
  inline friend _drovector operator*(const _drovector&, const  dgematrix&);
  
  inline friend _dgematrix operator*(const  dgematrix&, const     double&);
  inline friend _dgematrix operator*(const     double&, const  dgematrix&);
  
  //// / ////
  inline friend _dgematrix operator/(const  dgematrix&, const     double&);

  //// friend classes ////
  friend class _dgematrix;
  //friend class  dgbmatrix;
  //friend class _dgbmatrix;
  friend class  dcovector;
  friend class _dcovector;
  friend class  drovector;
  friend class _drovector;
};

---