mesh.h

 
#ifndef __62731052ee4a4a2d9f256c2caac43835
#define __62731052ee4a4a2d9f256c2caac43835
 
#include <lf/mesh/mesh.h>
#include <spdlog/logger.h>
 
#include "lf/base/base.h"
#include "lf/mesh/utils/print_info.h"
#include "point.h"
#include "quad.h"
#include "segment.h"
#include "triangle.h"
 
namespace lf::mesh::hybrid2d {
 
using size_type = lf::base::size_type;
using dim_t = lf::base::dim_t;
using sub_idx_t = lf::base::sub_idx_t;
using glb_idx_t = lf::base::glb_idx_t;
const unsigned int idx_nil = lf::base::kIdxNil;
 
class MeshFactory;
 
class Mesh : public mesh::Mesh {
 public:
  [[nodiscard]] unsigned DimMesh() const override { return 2; }
  [[nodiscard]] unsigned DimWorld() const override { return dim_world_; }
 
  [[nodiscard]] nonstd::span<const Entity* const> Entities(
      unsigned codim) const override;
  [[nodiscard]] size_type NumEntities(unsigned codim) const override;
  [[nodiscard]] size_type NumEntities(
      lf::base::RefEl ref_el_type) const override;
  [[nodiscard]] size_type Index(const Entity& e) const override;
  [[nodiscard]] const mesh::Entity* EntityByIndex(
      dim_t codim, glb_idx_t index) const override;
  [[nodiscard]] bool Contains(const mesh::Entity& e) const override;
 
 private:
  dim_t dim_world_{};
  std::vector<hybrid2d::Point> points_;
  std::vector<hybrid2d::Segment> segments_;
  std::vector<hybrid2d::Triangle> trias_;
  std::vector<hybrid2d::Quadrilateral> quads_;
 
  std::array<std::vector<const mesh::Entity*>, 3> entity_pointers_;
 
  using GeometryPtr = std::unique_ptr<geometry::Geometry>;
  using NodeCoordList = std::vector<GeometryPtr>;
  using EdgeList =
      std::vector<std::pair<std::array<size_type, 2>, GeometryPtr>>;
  using CellList =
      std::vector<std::pair<std::array<size_type, 4>, GeometryPtr>>;
 
  Mesh(dim_t dim_world, NodeCoordList nodes, EdgeList edges, CellList cells,
       bool check_completeness);
 
  friend class MeshFactory;
 
 public:
  static std::shared_ptr<spdlog::logger>& Logger() {
    static auto logger = base::InitLogger("lf::mesh::hybrid2d::Mesh::Logger");
    return logger;
  }
};
 
}  // namespace lf::mesh::hybrid2d
 
#endif  // __62731052ee4a4a2d9f256c2caac43835