main.cpp

#include <mesh.hpp>
#include <list>

// Forward declaration of some functions

void loop_through_half_edges_inside_a_face(HalfMesh::Mesh *mesh);

void loop_through_vertices_in_a_mesh(HalfMesh::Mesh *mesh);

void detect_boundary_edges(HalfMesh::Mesh *mesh);

void detect_boundary_half_edges(HalfMesh::Mesh *mesh);

int main() {
    HalfMesh::Mesh *new_mesh = new HalfMesh::Mesh();
    std::cout << "----> Reading Mesh (BEGIN)" << std::endl;
    new_mesh->read("bunny_merged.msh");
    std::cout << "----> Reading Mesh (END)" << std::endl;
    std::cout << "Multiply Connected : " << new_mesh->is_multiply_connected() << std::endl;
//    new_mesh->compute_number_of_holes();
//    std::cout << "Edges : " << new_mesh->get_edges().size() << std::endl;
//    std::cout << "Half Edges : " << new_mesh->get_half_edges().size() << std::endl;
//    new_mesh->save("test_1.msh");
//    new_mesh->save("test_1.bm");
//    new_mesh->save("test_1.vtk");
//    loop_through_vertices_in_a_mesh(new_mesh);
//    detect_boundary_edges(new_mesh);
    return 0;
}

// Create a simplical mesh
HalfMesh::Mesh *create_mesh() {
    HalfMesh::Mesh *new_mesh = new HalfMesh::Mesh;
    HalfMesh::Vertex *v1 = new HalfMesh::Vertex(0.0, 0.0, 0.0);
    HalfMesh::Vertex *v2 = new HalfMesh::Vertex(1.0, 0.0, 0.0);
    HalfMesh::Vertex *v3 = new HalfMesh::Vertex(0.0, 0.5, 0.0);
    HalfMesh::Vertex *v4 = new HalfMesh::Vertex(1.5, 0.5, 0.0);
    HalfMesh::Vertex *v5 = new HalfMesh::Vertex(2.5, 0.0, 0.0);
    new_mesh->add_vertex(v1);
    new_mesh->add_vertex(v2);
    new_mesh->add_vertex(v3);
    new_mesh->add_vertex(v4);
    new_mesh->add_vertex(v5);
    // add face
    new_mesh->add_face(v1, v2, v3);
    new_mesh->add_face(v2, v4, v3);
    new_mesh->add_face(v2, v5, v4);
    // Complete the mesh here
    new_mesh->complete_mesh();
    return new_mesh;
}

// Loop through half edges in a face
void loop_through_half_edges_inside_a_face(HalfMesh::Mesh *new_mesh) {
    for (unsigned int i = 0; i < new_mesh->get_faces().size(); ++i) {
        HalfMesh::Face *new_face = new_mesh->get_faces().at(i);
        std::cout << "One half edge for face F: " << new_face->handle() << " HE : "
                  << new_face->get_one_half_edge()->handle() << std::endl;
        unsigned int he_1, he_2, he_3, he_4;
        he_1 = new_face->get_one_half_edge()->handle();
        he_2 = new_mesh->get_next_half_edge(new_mesh->get_half_edge(he_1), new_face)->handle();
        he_3 = new_mesh->get_next_half_edge(new_mesh->get_half_edge(he_2), new_face)->handle();
        he_4 = new_mesh->get_next_half_edge(new_mesh->get_half_edge(he_3), new_face)->handle();
        std::cout << he_1 << "->" << he_2 << "->" << he_3 << "->" << he_4 << std::endl;
    }
}

// Loop through vertices
void loop_through_vertices_in_a_mesh(HalfMesh::Mesh *new_mesh) {
    for (unsigned int i = 0; i < new_mesh->get_vertices().size(); ++i) {
        HalfMesh::Vertex *new_vertex = new_mesh->get_vertices().at(i);
        std::cout << "Vertex " << i + 1 << " has " << new_vertex->get_incoming_half_edges().size()
                  << " incoming half edges and " << new_vertex->get_outgoing_half_edges().size()
                  << " outgoing half edges" << std::endl;
        HalfMesh::print_set(new_vertex->get_outgoing_half_edges());
        HalfMesh::print_set(new_vertex->get_incoming_half_edges());
    }
}

// Loop through edges and find boundaries
void detect_boundary_edges(HalfMesh::Mesh *new_mesh) {
    unsigned int count = 0;
    for (unsigned int i = 0; i < new_mesh->get_edges().size(); ++i) {
        HalfMesh::Edge *new_edge = new_mesh->get_edges().at(i);
        if (new_edge->is_boundary()) {
            count = count + 1;
        }
    }
    std::cout << "There is a total of " << count << " boundary edges in this mesh" << std::endl;
}

// Loop through half edges and find boundaries
void detect_boundary_half_edges(HalfMesh::Mesh *new_mesh) {
    unsigned int count = 0;
    for (unsigned int i = 0; i < new_mesh->get_half_edges().size(); ++i) {
        HalfMesh::HalfEdge *halfEdge = new_mesh->get_half_edges().at(i);
        if (halfEdge->is_boundary()) {
            count = count + 1;
        }
    }
    std::cout << "There is a total of " << count << " boundary half edges in this mesh" << std::endl;
}