idr.hpp

Go to the documentation of this file.

#ifndef PARALUTION_KRYLOV_IDR_HPP_
#define PARALUTION_KRYLOV_IDR_HPP_

#include "../solver.hpp"

#include <vector>

namespace paralution {

template <class OperatorType, class VectorType, typename ValueType>
class IDR : public IterativeLinearSolver<OperatorType, VectorType, ValueType> {

public:

  IDR();
  virtual ~IDR();

  virtual void Print(void) const;

  virtual void Build(void);
  virtual void Clear(void);

  virtual void SetShadowSpace(const int s);

protected:

  virtual void SolveNonPrecond_(const VectorType &rhs,
                                VectorType *x);
  virtual void SolvePrecond_(const VectorType &rhs,
                             VectorType *x);

  virtual void PrintStart_(void) const;
  virtual void PrintEnd_(void) const;

  virtual void MoveToHostLocalData_(void);
  virtual void MoveToAcceleratorLocalData_(void);

private:

  int s_;

  ValueType kappa_;

  ValueType *fhost_, *Mhost_;

  VectorType r_, v_, z_, t_;
  VectorType **g_, **u_, **P_;

};


}

#endif // PARALUTION_KRYLOV_IDR_HPP_