quad.h

 
#ifndef INCG9dc22ac9eb6645d3a27f60a7abcd52a4
#define INCG9dc22ac9eb6645d3a27f60a7abcd52a4
 
#include <lf/mesh/mesh.h>
 
namespace lf::mesh::hybrid2d {
 
class Point;
class Segment;
 
class Quadrilateral : public mesh::Entity {
  using size_type = mesh::Mesh::size_type;
 
 public:
  Quadrilateral() = default;
 
  Quadrilateral(const Quadrilateral&) = delete;
  Quadrilateral(Quadrilateral&&) noexcept = default;
  Quadrilateral& operator=(const Quadrilateral&) = delete;
  Quadrilateral& operator=(Quadrilateral&&) noexcept = default;
 
  explicit Quadrilateral(size_type index,
                         std::unique_ptr<geometry::Geometry>&& geometry,
                         const Point* corner0, const Point* corner1,
                         const Point* corner2, const Point* corner3,
                         const Segment* edge0, const Segment* edge1,
                         const Segment* edge2, const Segment* edge3);
 
  [[nodiscard]] unsigned Codim() const override { return 0; }
 
  [[nodiscard]] std::span<const Entity* const> SubEntities(
      unsigned rel_codim) const override;
 
  [[nodiscard]] std::span<const lf::mesh::Orientation> RelativeOrientations()
      const override {
    return edge_ori_;
  }
 
  [[nodiscard]] size_type index() const { return index_; }
 
  [[nodiscard]] const geometry::Geometry* Geometry() const override {
    return geometry_.get();
  }
  [[nodiscard]] base::RefEl RefEl() const override {
    return base::RefEl::kQuad();
  }
  [[nodiscard]] bool operator==(const mesh::Entity& rhs) const override {
    return this == &rhs;
  }
  ~Quadrilateral() override = default;
 
 private:
  size_type index_ = -1;  // zero-based index of this entity.
  std::unique_ptr<geometry::Geometry> geometry_;  // shape information
  std::array<const Point*, 4> nodes_{};           // nodes = corners of quad
  std::array<const Segment*, 4> edges_{};         // edges of quad
  std::array<lf::mesh::Orientation, 4>
      edge_ori_{};          // orientation of edges (set in constructor)
  Entity* this_ = nullptr;  // needed for SubEntity()
};
 
}  // namespace lf::mesh::hybrid2d
 
#endif  // INCG9dc22ac9eb6645d3a27f60a7abcd52a4