AssemblerInterface.hpp

#pragma once

#include <type_traits>

#include <dune/common/dynmatrix.hh>
#include <dune/common/dynvector.hh>

#include <amdis/ContextGeometry.hpp>

namespace AMDiS
{
  template <class LC, class C>
  struct DefaultAssemblerTraits
  {
    using LocalContext = LC;
    using ElementContainer = C;
  };

  template <class Traits, class... Nodes>
  class AssemblerInterface
  {
    using LocalContext = typename Traits::LocalContext;
    using ContextType = Impl::ContextTypes<LocalContext>;

  public:
    using Element = typename ContextType::Entity;
    using Geometry = typename Element::Geometry;

    static constexpr int numNodes = sizeof...(Nodes);
    static_assert( numNodes == 1 || numNodes == 2,
      "VectorAssembler gets 1 Node, MatrixAssembler gets 2 Nodes!");

  public:
    virtual ~AssemblerInterface() = default;

    virtual void bind(Element const& element, Geometry const& geometry) = 0;

    virtual void unbind() = 0;

    virtual void assemble(typename Traits::LocalContext const& localContext,
                          Nodes const&... nodes,
                          typename Traits::ElementContainer& elementMatrixVector) = 0;
  };

} // end namespace AMDiS