implicitScheme.hpp

#ifndef MARIAN_IMPLICITSCHEME_HPP
#define MARIAN_IMPLICITSCHEME_HPP

#include <FDM/schemes/fdScheme.hpp>
#include <FDM/tridiagonalSolver.hpp>

namespace marian {
  class ImplicitScheme : public DCFDScheme<ImplicitScheme> {
  public:
    ImplicitScheme(){};

    ImplicitScheme(SmartPointer<TridiagonalSolver> solver): solver_(solver) {};
    void setSolver(const SmartPointer<TridiagonalSolver>& solver) override {
      solver_ = solver;
    }
    std::vector<double> solve(std::vector<double> f,
                  const std::vector<SmartPointer<BoundaryCondition> >& bcs,
                  const std::vector<double>& time_grid,
                  const TridiagonalOperator& L)  override;
    std::vector<double> solveAndSave(std::vector<double> f,
                     const std::vector<SmartPointer<BoundaryCondition> >& bcs,
                     const std::vector<double>& spatial_grid,
                     const std::vector<double>& time_grid,
                     const TridiagonalOperator& L,
                     const std::string file_name) override;
    std::string info() const override {
      return "implicit";
    }
  private:
    SmartPointer<TridiagonalSolver> solver_; 
  };

} // namespace marian

#endif /* MARIAN_IMPLICITSCHEME_HPP */