boundaryCondition.hpp

#ifndef MARIAN_BOUNDARYCONDITION_HPP
#define MARIAN_BOUNDARYCONDITION_HPP

#include <vector>
#include <FDM/tridiagonalOperator.hpp>

namespace marian {
  enum class BCSide {LOW, 
             UPP, 
             FREE 
  };

  class BoundaryCondition {    
  public:
    BoundaryCondition(){};

    virtual void beforeExplicitStep(TridiagonalOperator& L) = 0;

    virtual void afterExplicitStep(std::vector<double>& f, double t) = 0;

    virtual void beforeImplicitStep(TridiagonalOperator& L,
                    std::vector<double>& f,
                    double t) = 0;

    virtual void afterImplicitStep(std::vector<double>& f,double t) = 0;

    virtual std::string info() const = 0;
    
    virtual BoundaryCondition* clone() const = 0;

    virtual ~BoundaryCondition(){};

  private:

  };

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


#endif /* MARIAN_BOUNDARYCONDITION_HPP */