tridiagonalSolver.hpp

#ifndef MARIAN_TRIDIAGONALSOLVER_HPP
#define MARIAN_TRIDIAGONALSOLVER_HPP

#include <FDM/tridiagonalOperator.hpp>

namespace marian {

  class TridiagonalSolver {
  public:
    TridiagonalSolver() {
    }

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

    virtual TridiagonalSolver* clone() const = 0;

    virtual ~TridiagonalSolver() {
    }
  };

  template<typename T>
  class DCTridiagonalSolver : public TridiagonalSolver {
  public:
    virtual TridiagonalSolver* clone() const {
      return new T(static_cast<const T&>(*this));
    }
  };
  
}  // namespace marian


#endif /* MARIAN_TRIDIAGONALSOLVER_HPP */