LUSolver.hpp

#ifndef MARIAN_LUSOLVER_HPP
#define MARIAN_LUSOLVER_HPP

#include <FDM/tridiagonalSolver.hpp>

namespace marian {
  class LUSolver : public DCTridiagonalSolver<LUSolver> {
  public:
    LUSolver(){};

    virtual std::vector<double> solve(const TridiagonalOperator& A,
                      const std::vector<double>& w) const override;

    ~LUSolver(){};
  };
  
}  // namespace marian



#endif /* MARIAN_LUSOLVER_HPP */