sub_element_matrix_provider.h

#ifndef PROJECTS_DPG_SUB_ELEMENT_MATRIX_PROVIDER
#define PROJECTS_DPG_SUB_ELEMENT_MATRIX_PROVIDER
 
#include <lf/mesh/mesh.h>
#include <lf/uscalfe/uscalfe.h>
 
#include "dpg.h"
 
namespace projects::dpg {
 
template <typename SCALAR>
class SubElementMatrixProvider {
 public:
  using size_type = lf::uscalfe::size_type;
  using elem_mat_t = Eigen::Matrix<SCALAR, Eigen::Dynamic, Eigen::Dynamic>;
  using ElemMat = elem_mat_t;
 
  SubElementMatrixProvider() = default;
  virtual ~SubElementMatrixProvider() = default;
 
  SubElementMatrixProvider(const SubElementMatrixProvider&) = delete;
  SubElementMatrixProvider(SubElementMatrixProvider&&) noexcept = default;
  SubElementMatrixProvider& operator=(const SubElementMatrixProvider&) = delete;
  SubElementMatrixProvider& operator=(SubElementMatrixProvider&&) = delete;
 
  virtual bool isActive(const lf::mesh::Entity& /*cell*/) { return true; }
 
  virtual ElemMat Eval(const lf::mesh::Entity& cell) = 0;
 
  [[nodiscard]] virtual size_type TrialComponent() const = 0;
 
  [[nodiscard]] virtual size_type TestComponent() const = 0;
};
 
}  // namespace projects::dpg
 
#endif  // PROJECTS_DPG_SUB_ELEMENT_MATRIX_PROVIDER