tprecon-ilu.hxx

Go to the documentation of this file.

#ifndef TPRECON_ILU_HXX
#define TPRECON_ILU_HXX

#include "tsparsematrix.hxx"

namespace numcxx
{

  template<typename T> 
  class TPreconILU: public  TLinSolver<T>
  {

  public:
    std::shared_ptr< TArray1<T> > pInvDiag;
    std::shared_ptr< TArray1<int> > pDiagIdx;

    TPreconILU(const std::shared_ptr<TSparseMatrix<T>> pA);
    TPreconILU(TSparseMatrix<T> &A);

    ~TPreconILU(){};

    static std::shared_ptr<TPreconILU<T>> create(const std::shared_ptr<TSparseMatrix<T>> pA);

    void update(TSparseMatrix<T> &A);

    void solve(TArray<T> & Sol,  const TArray<T> & Rhs) const;

  private:
    std::shared_ptr<TArray1<T>> pA;
    std::shared_ptr<TArray1<int>> pIA;
    std::shared_ptr<TArray1<int>> pJA;
  };
}

#include "tprecon-ilu.ixx"

#endif