scalar_fe_space.h

#ifndef LF_USCALFE_SCALAR_FE_SPACE_H_
#define LF_USCALFE_SCALAR_FE_SPACE_H_
/***************************************************************************
 * LehrFEM++ - A simple C++ finite element libray for teaching
 * Developed from 2018 at the Seminar of Applied Mathematics of ETH Zurich,
 * lead developers Dr. R. Casagrande and Prof. R. Hiptmair
 ***************************************************************************/
 
#include <lf/assemble/assemble.h>
 
#include "scalar_reference_finite_element.h"
 
namespace lf::fe {
 
template <typename SCALAR>
class ScalarFESpace {
 public:
  using Scalar = SCALAR;
 
 protected:
  ScalarFESpace() = default;
  ScalarFESpace(const ScalarFESpace &) = default;
  ScalarFESpace(ScalarFESpace &&) noexcept = default;
  ScalarFESpace &operator=(const ScalarFESpace &) = default;
  ScalarFESpace &operator=(ScalarFESpace &&) noexcept = default;
 
 public:
  [[nodiscard]] virtual std::shared_ptr<const lf::mesh::Mesh> Mesh() const = 0;
  [[nodiscard]] virtual const lf::assemble::DofHandler &LocGlobMap() const = 0;
 
  [[nodiscard]] virtual ScalarReferenceFiniteElement<SCALAR> const *
  ShapeFunctionLayout(const lf::mesh::Entity &entity) const = 0;
 
  [[nodiscard]] virtual size_type NumRefShapeFunctions(
      const lf::mesh::Entity &entity) const = 0;
 
  virtual ~ScalarFESpace() = default;
};  // end class definition ScalarFESpace
 
template <typename SCALAR>
std::ostream &operator<<(std::ostream &o, const ScalarFESpace<SCALAR> &fes) {
  fes.PrintInfo(o, 0);
  return o;
}
 
}  // namespace lf::fe
 
#endif