fdScheme.hpp

#ifndef MARIAN_FDSCHEME_HPP
#define MARIAN_FDSCHEME_HPP

#include <vector>
#include <utils/smartPointer.hpp>
#include <FDM/tridiagonalSolver.hpp>
#include <FDM/boundaryConditions/boundaryCondition.hpp>

namespace marian {

  class FDScheme {
  public:
  
    virtual void setSolver(const SmartPointer<TridiagonalSolver>& solver) = 0; 
    virtual std::vector<double> solve(std::vector<double> f,
                      const std::vector<SmartPointer<BoundaryCondition> >& bcs,
                      const std::vector<double>& time_grid,
                      const TridiagonalOperator& L) = 0;
                      
    virtual 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) = 0;
                          
    virtual std::string info() const = 0;
    
    virtual FDScheme* clone() const = 0;
    
    virtual ~FDScheme(){};
  };

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


#endif /* MARIAN_FDSCHEME_HPP */