point.h

 
#ifndef INCG818709b0104548a7b5e6f47bdba89f69
#define INCG818709b0104548a7b5e6f47bdba89f69
 
#include <lf/mesh/mesh.h>
 
#include <iostream>
 
namespace lf::mesh::hybrid2d {
 
class Point : public mesh::Entity {
  using size_type = mesh::Mesh::size_type;
 
 public:
  Point() = default;
 
  Point(const Point&) = delete;
  Point(Point&&) noexcept = default;
  Point& operator=(const Point&) = delete;
  Point& operator=(Point&&) noexcept = default;
  explicit Point(size_type index,
                 std::unique_ptr<geometry::Geometry>&& geometry)
      : index_(index), geometry_(std::move(geometry)), this_(this) {
    // DIAGNOSTICS
    // std::cout << "hybrid2d::Point(" << index_ << ") " << std::endl;
    LF_VERIFY_MSG(geometry_, "Point must be supplied with a geometry");
    LF_VERIFY_MSG(geometry_->DimLocal() == 0,
                  "Geometry must be that of a point");
    LF_VERIFY_MSG(geometry_->RefEl() == base::RefEl::kPoint(),
                  "Geometry must fit point");
  }
 
  [[nodiscard]] unsigned Codim() const override { return 2; }
 
  [[nodiscard]] std::span<const Entity* const> SubEntities(
      unsigned rel_codim) const override {
    LF_ASSERT_MSG(rel_codim == 0, "A point has only codim = 0 sub-entities");
    return {&this_, 1};
  }
 
  [[nodiscard]] std::span<const lf::mesh::Orientation> RelativeOrientations()
      const override {
    LF_ASSERT_MSG(false, "A point has not sub-entities");
    return std::span<Orientation>();
  }
 
  [[nodiscard]] const geometry::Geometry* Geometry() const override {
    return geometry_.get();
  }
 
  [[nodiscard]] size_type index() const { return index_; }
 
  [[nodiscard]] base::RefEl RefEl() const override {
    return base::RefEl::kPoint();
  }
 
  [[nodiscard]] bool operator==(const mesh::Entity& rhs) const override {
    return this == &rhs;
  }
 
  ~Point() override = default;
 
 private:
  size_type index_ = -1;  // zero-based index of this entity.
  std::unique_ptr<geometry::Geometry> geometry_ = nullptr;  // shape information
  static constexpr std::array<lf::mesh::Orientation, 1> dummy_or_{
      lf::mesh::Orientation::positive};
  Entity* this_ = nullptr;  // needed for SubEntity()
};
 
}  // namespace lf::mesh::hybrid2d
 
#endif  // INCG818709b0104548a7b5e6f47bdba89f69