tprecon-jacobi.hxx

Go to the documentation of this file.

#ifndef TPRECON_JACOBI_HXX
#define TPRECON_JACOBI_HXX

#include "tsparsematrix.hxx"

namespace numcxx
{

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

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

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

    ~TPreconJacobi(){};

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

    void update(TSparseMatrix<T> &A);

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

#include "tprecon-jacobi.ixx"

#endif