tsolver-umfpack.hxx

Go to the documentation of this file.

#ifndef TSOLVER_UMFPACK_HXX
#define TSOLVER_UMFPACK_HXX

#include "tsparsematrix.hxx"

namespace numcxx
{

  template<typename T> 
  class TSolverUMFPACK: public  TLinSolver<T>
  {
    const std::shared_ptr< TSparseMatrix<T> > pMatrix;
        
    void * Symbolic;
        
    void * Numeric;

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

  public:
    TSolverUMFPACK(){};

    TSolverUMFPACK(const std::shared_ptr<TSparseMatrix<T>> pA);

    TSolverUMFPACK(const TSparseMatrix<T>&A);

    TSolverUMFPACK(const TSolverUMFPACK<T>& A)=delete;

    ~TSolverUMFPACK();

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

    void update();

    void update(const TSparseMatrix<T>&A);


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

    void solve( std::shared_ptr< TArray<T> > Sol,  const std::shared_ptr<TArray<T> > Rhs) {solve(*Sol,*Rhs);};

  };
}

#include "tsolver-umfpack.ixx"

#endif