pzshtmat.h

Go to the documentation of this file.

#ifndef SHTMATRIXHPP
#define SHTMATRIXHPP

#include <iostream>

template <class TObj>
class TPZGenMatrix {
                
protected:
                TObj *fMem;
                int64_t   fRows, fCols;
                
public:
                
                TPZGenMatrix ();
                TPZGenMatrix (const int64_t rows ,const int64_t columns);
                
                TPZGenMatrix (const TPZGenMatrix & A);
                
                ~TPZGenMatrix();
                
                void Print (const char *mess,std::ostream & out = std::cout) const;
                
                int64_t Rows() const {return fRows;}
                
                int64_t Cols() const {return fCols;}
                
                void Resize(const int64_t newrow,const int64_t newcol);
                
                //DEFINE OPERATORS
                
                TPZGenMatrix<TObj>& operator=(const TPZGenMatrix<TObj> & rval); //makes this.p = rval.p
                
                
                TObj & operator()(const int64_t row,const int64_t column = 0) const; //get value, set value
                
};

template <class TObj>
class TPZGenAMatrix : public TPZGenMatrix<TObj> {
                
                TPZGenAMatrix () : TPZGenMatrix<TObj>() {}                                                      //creates NULL matrix
                
                TPZGenAMatrix (const int64_t rows ,const int64_t columns) : TPZGenMatrix<TObj>(rows,columns) {} //creates a rowsXcolumns matrix
                
                TPZGenAMatrix operator+(const TPZGenAMatrix<TObj> & rval) const;
                
                TPZGenAMatrix operator+(const TObj x) const;
                
                //  friend TPZGenAMatrix operator+(const TObj x, const TPZGenAMatrix<TObj>& rval); // {return (rval+x);}
                
                TPZGenAMatrix operator-(const TPZGenAMatrix<TObj> & rval) const;
                
                TPZGenAMatrix operator-(const TObj rval) const;
                
                //  friend TPZGenAMatrix operator-(const TObj x, const TPZGenAMatrix<TObj>& rval); // {return ((-rval)+x);}
                
                TPZGenAMatrix operator-() const;
                
                TPZGenAMatrix operator*(const TPZGenAMatrix<TObj> & rval) const;
                
                TPZGenAMatrix operator*(const TObj rval) const;
                
                //  friend TPZGenAMatrix operator*(const TObj x, const TPZGenAMatrix<TObj>& rval); // {return (rval*x);}
                
                TPZGenAMatrix<TObj>& operator+=(const TPZGenAMatrix<TObj> & rval);
                
                TPZGenAMatrix<TObj>& operator+=(const TObj x);
                
                TPZGenAMatrix<TObj>& operator-=(const TPZGenAMatrix<TObj> & rval);
                
                TPZGenAMatrix<TObj>& operator-=(const TObj x);
                
                TPZGenAMatrix<TObj>& operator*=(const TObj x);
                
                
                TPZGenAMatrix<TObj> Transpose() const;
                
                
};

#endif