vector.hpp

Go to the documentation of this file.

#ifndef PARGEMSLR_VECTOR_H
#define PARGEMSLR_VECTOR_H
 
#include "../utils/utils.hpp"
#include "../utils/parallel.hpp"
#include "../utils/memory.hpp"
 
namespace pargemslr
{
   
   template <typename T>
   class VectorVirtualClass: public parallel_log
   {
   public:
      
      VectorVirtualClass();
      
      VectorVirtualClass(const VectorVirtualClass<T> &vec);
      
      VectorVirtualClass(VectorVirtualClass<T> &&vec);
      
      virtual int Clear();
      
      virtual ~VectorVirtualClass();
      
      virtual T&     operator[] (int i) = 0;
      
      virtual T* GetData() const = 0;
      
      virtual int GetDataLocation() const = 0;
      
      virtual int GetLengthLocal() const = 0;
      
      virtual int Fill(const T &v) = 0;
      
      virtual int MoveData( const int &location) = 0;
      
      virtual bool IsParallel() const;
      
      PrecisionEnum GetPrecision() const;
      
   };
   
   template<> struct PargemslrIsComplex<VectorVirtualClass<complexs> > : public std::true_type {};
   template<> struct PargemslrIsComplex<VectorVirtualClass<complexd> > : public std::true_type {};
   
   template <typename T>
   class VectorClass: public VectorVirtualClass<T>
   {
   public:
   
      VectorClass();
      
      VectorClass(const VectorClass<T> &vec);
      
      VectorClass(VectorClass<T> &&vec);
      
      virtual int Clear();
      
      virtual ~VectorClass();
      
      
      /* Operations - - - - - - - */
      
      virtual T&     operator[] (int i) = 0;
      
      virtual T* GetData() const = 0;
      
      virtual int GetDataLocation() const = 0;
      
      virtual int GetLengthLocal() const = 0;
      
      virtual long int GetLengthGlobal() const = 0;
      
      virtual long int GetStartGlobal() const = 0;
      
      virtual int Fill(const T &v) = 0;
      
      virtual int Rand() = 0;
      
      virtual int Scale( const T &alpha) = 0;
      
      virtual int Axpy( const T &alpha, const VectorClass<T> &x) = 0;
      
      virtual int Axpy( const T &alpha, const VectorClass<T> &x, VectorClass<T> &y) = 0;
      
      virtual int Axpy( const T &alpha, const VectorClass<T> &x, const T &beta, VectorClass<T> &y) = 0;
      
      virtual int Norm2( float &norm) const = 0;
      
      virtual int Norm2( double &norm) const = 0;
      
      virtual int NormInf( float &norm) = 0;
      
      virtual int NormInf( double &norm) = 0;
      
      virtual int Dot( const VectorClass<T> &y, T &t) const = 0;
 
#ifdef PARGEMSLR_CUDA 
#if (PARGEMSLR_CUDA_VERSION == 11)
 
      virtual cusparseDnVecDescr_t GetCusparseVec() const = 0;
#endif
#endif
 
      virtual int MoveData( const int &location) = 0;
      
   };
   
   typedef VectorClass<float>     vector_base_float;
   typedef VectorClass<double>    vector_base_double;
   typedef VectorClass<complexs>  vector_base_complexs;
   typedef VectorClass<complexd>  vector_base_complexd;
   template<> struct PargemslrIsComplex<VectorClass<complexs> > : public std::true_type {};
   template<> struct PargemslrIsComplex<VectorClass<complexd> > : public std::true_type {};
   
}
 
#endif