sub_element_vector_provider.h

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