tsolver-lapacklu.hxx

Go to the documentation of this file.

#ifndef TSOLVER_LAPACK_LU_HXX
#define TSOLVER_LAPACK_LU_HXX

#include "tmatrix.hxx"

namespace numcxx
{
  template<typename T> 
  class TSolverLapackLU: public  TLinSolver<T>
  {
    const std::shared_ptr< TMatrix<T> > pMatrix;
    const std::shared_ptr< TMatrix<T> > pLU;
    const std::shared_ptr< TArray1<int>> pIPiv;
  public:
    TSolverLapackLU(const std::shared_ptr<TMatrix<T>> pMatrix);

    TSolverLapackLU(const TMatrix<T>& Matrix);

    static std::shared_ptr<TSolverLapackLU<T>> create(const std::shared_ptr<TMatrix<T>> pMatrix);

    void update();

    void update(const TMatrix<T>& Matrix);

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

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

    std::shared_ptr<TMatrix<T>> calculate_inverse();

    TSolverLapackLU() {};



    TMatrix<T> & LU(){ return *pLU;}
    TArray1<int> & IPiv(){ return *pIPiv;}
  };

  template<typename T> 
  inline std::ostream & operator << (std::ostream & s, TSolverLapackLU<T> &LU)
  {
    s << "LU:" << std::endl;
    s << LU.LU() << std::endl;
    s << "IPiv:" << std::endl;
    s << LU.IPiv() << std::endl;
  }

}
#include "tsolver-lapacklu.ixx"
#endif