//=============================================================================
//! (DO NOT USE) Smart-temporary Real Double-precision Symmetric Sparse Matrix Class
class _dssmatrix
{
public:
//// objects ////
mutable long N; //!< matrix row and column size
mutable std::vector<dcomponent> Data; //!< matrix data
mutable std::vector< std::vector< std::pair<long,long> > > Line; //!< vector of vector to store the entry information of component for each row and column
//// const references ////
long const& M; //!< matrix row and column size (readable)
//// constructor ////
inline _dssmatrix();
inline _dssmatrix(const _dssmatrix&);
inline ~_dssmatrix(); //destructor
//// cast ////
inline _zhsmatrix to_zhsmatrix() const;
inline _dgematrix to_dgematrix() const;
inline _dsymatrix to_dsymatrix() const;
inline _dgsmatrix to_dgsmatrix() const;
//// io ////
inline double operator()(const long&, const long&) const;//not return double&
inline friend std::ostream& operator<<(std::ostream&, const _dssmatrix&);
inline void write(const char *) const;
//// misc ////
inline void nullify() const;
inline void destroy() const;
//// calc ////
inline friend _dssmatrix t(const dssmatrix&);
inline friend void idamax(long&, long&, const dssmatrix&);
inline friend double damax(const dssmatrix&);
//// unary ////
inline friend const _dssmatrix& operator+(const _dssmatrix&);
inline friend _dssmatrix operator-(const _dssmatrix&);
//// + ////
inline friend _dssmatrix operator+(const dssmatrix&, const dssmatrix&);
inline friend _dssmatrix operator+(const dssmatrix&, const _dssmatrix&);
inline friend _dssmatrix operator+(const _dssmatrix&, const dssmatrix&);
inline friend _dssmatrix operator+(const _dssmatrix&, const _dssmatrix&);
inline friend _dgematrix operator+(const _dssmatrix&, const dgematrix&);
inline friend _dgematrix operator+(const _dssmatrix&, const _dgematrix&);
inline friend _dgematrix operator+(const _dgematrix&, const _dssmatrix&);
inline friend _dgematrix operator+(const dgematrix&, const _dssmatrix&);
inline friend _dgematrix operator+(const _dssmatrix&, const dgbmatrix&);
inline friend _dgematrix operator+(const _dssmatrix&, const _dgbmatrix&);
inline friend _dgematrix operator+(const _dgbmatrix&, const _dssmatrix&);
inline friend _dgematrix operator+(const dgbmatrix&, const _dssmatrix&);
inline friend _dgematrix operator+(const _dssmatrix&, const dsymatrix&);
inline friend _dgematrix operator+(const _dssmatrix&, const _dsymatrix&);
inline friend _dgematrix operator+(const _dsymatrix&, const _dssmatrix&);
inline friend _dgematrix operator+(const dsymatrix&, const _dssmatrix&);
//// - ////
inline friend _dssmatrix operator-(const dssmatrix&, const dssmatrix&);
inline friend _dssmatrix operator-(const dssmatrix&, const _dssmatrix&);
inline friend _dssmatrix operator-(const _dssmatrix&, const dssmatrix&);
inline friend _dssmatrix operator-(const _dssmatrix&, const _dssmatrix&);
inline friend _dgematrix operator-(const _dssmatrix&, const dgematrix&);
inline friend _dgematrix operator-(const _dssmatrix&, const _dgematrix&);
inline friend _dgematrix operator-(const _dgematrix&, const _dssmatrix&);
inline friend _dgematrix operator-(const dgematrix&, const _dssmatrix&);
inline friend _dgematrix operator-(const _dssmatrix&, const dgbmatrix&);
inline friend _dgematrix operator-(const _dssmatrix&, const _dgbmatrix&);
inline friend _dgematrix operator-(const _dgbmatrix&, const _dssmatrix&);
inline friend _dgematrix operator-(const dgbmatrix&, const _dssmatrix&);
inline friend _dgematrix operator-(const _dssmatrix&, const dsymatrix&);
inline friend _dgematrix operator-(const _dssmatrix&, const _dsymatrix&);
inline friend _dgematrix operator-(const _dsymatrix&, const _dssmatrix&);
inline friend _dgematrix operator-(const dsymatrix&, const _dssmatrix&);
//// * ////
inline friend _dgsmatrix operator*(const dssmatrix&, const dssmatrix&);
inline friend _dgsmatrix operator*(const dssmatrix&, const _dssmatrix&);
inline friend _dgsmatrix operator*(const _dssmatrix&, const dssmatrix&);
inline friend _dgsmatrix operator*(const _dssmatrix&, const _dssmatrix&);
inline friend _dgematrix operator*(const _dssmatrix&, const dgematrix&);
inline friend _dgematrix operator*(const _dssmatrix&, const _dgematrix&);
inline friend _dgematrix operator*(const _dgematrix&, const _dssmatrix&);
inline friend _dgematrix operator*(const dgematrix&, const _dssmatrix&);
inline friend _dgematrix operator*(const _dssmatrix&, const dgbmatrix&);
inline friend _dgematrix operator*(const _dssmatrix&, const _dgbmatrix&);
inline friend _dgematrix operator*(const _dgbmatrix&, const _dssmatrix&);
inline friend _dgematrix operator*(const dgbmatrix&, const _dssmatrix&);
inline friend _dgematrix operator*(const _dssmatrix&, const dsymatrix&);
inline friend _dgematrix operator*(const _dssmatrix&, const _dsymatrix&);
inline friend _dgematrix operator*(const _dsymatrix&, const _dssmatrix&);
inline friend _dgematrix operator*(const dsymatrix&, const _dssmatrix&);
inline friend _dcovector operator*(const dssmatrix&, const dcovector&);
inline friend _dcovector operator*(const dssmatrix&, const _dcovector&);
inline friend _drovector operator*(const drovector&, const _dssmatrix&);
inline friend _drovector operator*(const _drovector&, const _dssmatrix&);
inline friend _dssmatrix operator*(const dssmatrix&, const double&);
inline friend _dssmatrix operator*(const _dssmatrix&, const double&);
inline friend _dssmatrix operator*(const double&, const dssmatrix&);
inline friend _dssmatrix operator*(const double&, const _dssmatrix&);
//// / ////
inline friend _dssmatrix operator/(const dssmatrix&, const double&);
inline friend _dssmatrix operator/(const _dssmatrix&, const double&);
//// friend classes ////
friend class dgematrix;
friend class _dgematrix;
friend class dgbmatrix;
friend class _dgbmatrix;
friend class dsymatrix;
friend class _dsymatrix;
friend class dssmatrix;
friend class dcovector;
friend class _dcovector;
friend class drovector;
friend class _drovector;
};