#include <vector>
#include "fmesh/fixed_size_face.hpp"
#include "fmesh/fracture_mesh.hpp"

using namespace fmesh;

struct point {
  double x;
  double y;
  double z;

  point() = default;
  point(double tx, double ty, double tz) : x{tx}, y{ty}, z{tz} {}
};

int main() {
  fracture_mesh<point, tri_face> mesh;

  // Add vertices
  std::vector<vertex_index> v_ids;
  v_ids.push_back(mesh.add_vertex(0.0, 0.0, 0.0));
  v_ids.push_back(mesh.add_vertex(1.0, 0.0, 0.0));
  v_ids.push_back(mesh.add_vertex(0.0, 0.0, 1.0));
  v_ids.push_back(mesh.add_vertex(1.0, 0.0, 1.0));
  v_ids.push_back(mesh.add_vertex(0.0, 0.0, 2.0));
  v_ids.push_back(mesh.add_vertex(2.0, 0.0, 0.5));

  // Add faces
  std::vector<face_index> f_ids;
  f_ids.push_back(mesh.add_face(v_ids[0], v_ids[1], v_ids[2]));
  f_ids.push_back(mesh.add_face(v_ids[1], v_ids[3], v_ids[2]));
  f_ids.push_back(mesh.add_face(v_ids[2], v_ids[3], v_ids[4]));
  f_ids.push_back(mesh.add_face(v_ids[1], v_ids[5], v_ids[3]));

  // Access to a point
  const auto& p1 = mesh.vertex(v_idx[0]);

  // Access to a face
  const auto& f1 = mesh.face(f_idx[0]);

  // Range-based for
  for (auto&& vi : mesh.vertices()) {
      const auto& p = mesh.vertex(vi);
      // ...
  }

  for (auto&& ei : mesh.edges()) {
      const auto& e = mesh.edge(ei);
      // ...
  }

  for (auto&& fi : mesh.faces()) {
      const auto& f = mesh.face(fi);
      // ...
  }

  // Invalidate vertices and faces
  mesh.invalidate(v_ids[0]);
  mesh.invalidate(f_ids[1]);

  // Remove invalid mesh entities
  mesh.remove_invalid_entities();
}