BoundaryCondition.hpp

#pragma once

#include <list>
#include <memory>

#include <amdis/typetree/TreeContainer.hpp>

namespace AMDiS
{

  template <class Mat, class Sol, class Rhs>
  class BoundaryCondition
  {
  public:
    BoundaryCondition() = default;


    virtual void init() = 0;


    virtual void apply(Mat& A, Sol& x, Rhs& b) = 0;
  };

  template <class Mat, class Sol, class Rhs>
  struct BCData
  {
    template <class RN, class CN>
    using type = std::list<std::shared_ptr<BoundaryCondition<Mat, Sol, Rhs>>>;
  };

  template <class Mat, class Sol, class Rhs, class RB, class CB>
  using BoundaryConditions
    = TreeMatrix<BCData<Mat, Sol, Rhs>::template type,
                 typename RB::LocalView::Tree, typename CB::LocalView::Tree>;

} // end namespace AMDiS