//=============================================================================
//! Real Double-precision General Sparse Matrix Class
class dgsmatrix
{
private:
//// objects ////
long M; //!< matrix row size (NOT accessible)
long N; //!< matrix column size (NOT accessible)
std::vector<dcomponent> Data; //!< matrix data
std::vector< std::vector< std::pair<long,long> > > Rows; //!< array of vector to store the entry information of component for each row
std::vector< std::vector< std::pair<long,long> > > Cols; //!< array of vector to store the entry information of component for each column
public:
//// const references ////
long const& m; //!< matrix row size (readable)
long const& n; //!< matrix column size (readable)
std::vector<dcomponent> const& data; //!< reference of Data (readable)
std::vector< std::vector< std::pair<long,long> > > const& rows; //!< array of vector to store the entry information of component for each row (readable)
std::vector< std::vector< std::pair<long,long> > > const& cols; //!< array of vector to store the entry information of component for each column (readable)
//// constructor ////
inline dgsmatrix();
inline dgsmatrix(const dgsmatrix&);
inline dgsmatrix(const _dgsmatrix&);
inline dgsmatrix(const long&, const long&, const long=0);
inline dgsmatrix(const char*);
inline ~dgsmatrix(); //destructor
//// cast ////
inline _zgsmatrix to_zgsmatrix() const;
inline _dgematrix to_dgematrix() const;
//// io ////
inline double operator()(const long&, const long&) const;
inline double& operator()(const long&, const long&);
inline dcomponent component(const long&) const;
inline dcomponent& component(const long&);
inline dgsmatrix& put(const long&, const long&, const double&);
inline dgsmatrix& fput(const long&, const long&, const double&);
inline dgsmatrix& del(const long, const long); //<-- NOT (const long&)
inline dgsmatrix& fdel(const long); //<-- NOT (const long&)
inline friend std::ostream& operator<<(std::ostream&, const dgsmatrix&);
inline void write(const char*) const;
inline void read(const char*);
//// misc ////
inline void clear();
inline dgsmatrix& zero();
inline void chsign();
inline void copy(const dgsmatrix&);
inline void shallow_copy(const _dgsmatrix&);
inline dgsmatrix& resize(const long&, const long&, const long=0, const long=0);
inline void stretch(const long&, const long&);
inline bool isListed(const long&, const long&) const;
inline long number(const long&, const long&);
inline void diet(const double=DBL_MIN);
inline void checkup();
inline _drovector row(const long&) const;
inline _dcovector col(const long&) const;
inline friend void swap(dgsmatrix&, dgsmatrix&);
inline friend _dgsmatrix _(dgsmatrix&);
//// calc ////
inline friend _dgsmatrix t(const dgsmatrix&);
inline friend void idamax(long&, long&, const dgsmatrix&);
inline friend double damax(const dgsmatrix&);
//// unary ////
inline friend const dgsmatrix& operator+(const dgsmatrix&);
inline friend _dgsmatrix operator-(const dgsmatrix&);
//// operators with dgsmatrix ////
inline dgsmatrix& operator=(const dgsmatrix&);
inline dgsmatrix& operator+=(const dgsmatrix&);
inline dgsmatrix& operator-=(const dgsmatrix&);
inline dgsmatrix& operator*=(const dgsmatrix&);
//// operators with _dgsmatrix ////
inline dgsmatrix& operator=(const _dgsmatrix&);
inline dgsmatrix& operator+=(const _dgsmatrix&);
inline dgsmatrix& operator-=(const _dgsmatrix&);
inline dgsmatrix& operator*=(const _dgsmatrix&);
//// operators with double ////
inline dgsmatrix& operator*=(const double&);
inline dgsmatrix& operator/=(const double&);
//// + ////
inline friend _dgematrix operator+(const dgsmatrix&, const dgematrix&);
inline friend _dgematrix operator+(const dgsmatrix&, const _dgematrix&);
inline friend _dgematrix operator+(const dgematrix&, const dgsmatrix&);
inline friend _dgematrix operator+(const _dgematrix&, const dgsmatrix&);
inline friend _dgematrix operator+(const dgsmatrix&, const dgbmatrix&);
inline friend _dgematrix operator+(const dgsmatrix&, const _dgbmatrix&);
inline friend _dgematrix operator+(const dgbmatrix&, const dgsmatrix&);
inline friend _dgematrix operator+(const _dgbmatrix&, const dgsmatrix&);
inline friend _dgematrix operator+(const dgsmatrix&, const dsymatrix&);
inline friend _dgematrix operator+(const dgsmatrix&, const _dsymatrix&);
inline friend _dgematrix operator+(const dsymatrix&, const dgsmatrix&);
inline friend _dgematrix operator+(const _dsymatrix&, const dgsmatrix&);
inline friend _dgsmatrix operator+(const dgsmatrix&, const dgsmatrix&);
inline friend _dgsmatrix operator+(const dgsmatrix&, const _dgsmatrix&);
inline friend _dgsmatrix operator+(const _dgsmatrix&, const dgsmatrix&);
//// - ////
inline friend _dgematrix operator-(const dgsmatrix&, const dgematrix&);
inline friend _dgematrix operator-(const dgsmatrix&, const _dgematrix&);
inline friend _dgematrix operator-(const dgematrix&, const dgsmatrix&);
inline friend _dgematrix operator-(const _dgematrix&, const dgsmatrix&);
inline friend _dgematrix operator-(const dgsmatrix&, const dgbmatrix&);
inline friend _dgematrix operator-(const dgsmatrix&, const _dgbmatrix&);
inline friend _dgematrix operator-(const dgbmatrix&, const dgsmatrix&);
inline friend _dgematrix operator-(const _dgbmatrix&, const dgsmatrix&);
inline friend _dgematrix operator-(const dgsmatrix&, const dsymatrix&);
inline friend _dgematrix operator-(const dgsmatrix&, const _dsymatrix&);
inline friend _dgematrix operator-(const dsymatrix&, const dgsmatrix&);
inline friend _dgematrix operator-(const _dsymatrix&, const dgsmatrix&);
inline friend _dgsmatrix operator-(const dgsmatrix&, const dgsmatrix&);
inline friend _dgsmatrix operator-(const dgsmatrix&, const _dgsmatrix&);
inline friend _dgsmatrix operator-(const _dgsmatrix&, const dgsmatrix&);
//// * ////
inline friend _dgematrix operator*(const dgsmatrix&, const dgematrix&);
inline friend _dgematrix operator*(const dgsmatrix&, const _dgematrix&);
inline friend _dgematrix operator*(const dgematrix&, const dgsmatrix&);
inline friend _dgematrix operator*(const _dgematrix&, const dgsmatrix&);
inline friend _dgematrix operator*(const dgsmatrix&, const dgbmatrix&);
inline friend _dgematrix operator*(const dgsmatrix&, const _dgbmatrix&);
inline friend _dgematrix operator*(const dgbmatrix&, const dgsmatrix&);
inline friend _dgematrix operator*(const _dgbmatrix&, const dgsmatrix&);
inline friend _dgematrix operator*(const dsymatrix&, const dgsmatrix&);
inline friend _dgematrix operator*(const _dsymatrix&, const dgsmatrix&);
inline friend _dgematrix operator*(const dgsmatrix&, const dsymatrix&);
inline friend _dgematrix operator*(const dgsmatrix&, const _dsymatrix&);
inline friend _dgsmatrix operator*(const dgsmatrix&, const dgsmatrix&);
inline friend _dgsmatrix operator*(const dgsmatrix&, const _dgsmatrix&);
inline friend _dgsmatrix operator*(const _dgsmatrix&, const dgsmatrix&);
inline friend _dcovector operator*(const dgsmatrix&, const dcovector&);
inline friend _dcovector operator*(const dgsmatrix&, const _dcovector&);
inline friend _drovector operator*(const drovector&, const dgsmatrix&);
inline friend _drovector operator*(const _drovector&, const dgsmatrix&);
inline friend _dgsmatrix operator*(const dgsmatrix&, const double&);
inline friend _dgsmatrix operator*(const double&, const dgsmatrix&);
//// / ////
inline friend _dgsmatrix operator/(const dgsmatrix&, const double&);
//// friend classes ////
friend class dgematrix;
friend class _dgematrix;
friend class dgbmatrix;
friend class _dgbmatrix;
friend class dsymatrix;
friend class _dsymatrix;
friend class _dgsmatrix;
friend class dcovector;
friend class _dcovector;
friend class drovector;
friend class _drovector;
};