matrix.hpp

Go to the documentation of this file.

#ifndef PARGEMSLR_MATRIX_H
#define PARGEMSLR_MATRIX_H
 
#include "../utils/utils.hpp"
#include "../utils/parallel.hpp"
#include "../utils/memory.hpp"
#include "../vectors/vector.hpp"
 
namespace pargemslr
{
   
   template <typename T>
   class MatrixClass: public parallel_log
   {
   public:
      
      MatrixClass();
      
      MatrixClass(const MatrixClass<T> &mat);
      
      MatrixClass(MatrixClass<T> &&mat);
      
      virtual int Clear();
      
      virtual ~MatrixClass();
      
      virtual int GetDataLocation() const;
      
      virtual int GetNumRowsLocal() const = 0;
      
      virtual int GetNumColsLocal() const = 0;
      
      virtual long int  GetNumNonzeros() const;
      
      virtual int Eye();
      
      virtual int Fill(const T &v);
      
      virtual int Scale(const T &alpha);
      
      virtual int MatVec( char trans, const T &alpha, const VectorClass<T> &x, const T &beta, VectorClass<T> &y);
      
      virtual int MatVec( char trans, const T &alpha, const VectorClass<T> &x, const T &beta, const VectorClass<T> &y, VectorClass<T> &z);
      
      virtual int MoveData( const int &location);
      
      PrecisionEnum GetPrecision() const;
      
   };
   
   typedef MatrixClass<float>     matrix_base_float;
   typedef MatrixClass<double>    matrix_base_double;
   typedef MatrixClass<complexs>  matrix_base_complexs;
   typedef MatrixClass<complexd>  matrix_base_complexd;
   template<> struct PargemslrIsComplex<MatrixClass<complexs> > : public std::true_type {};
   template<> struct PargemslrIsComplex<MatrixClass<complexd> > : public std::true_type {};
   
}
 
#endif