Method subdivide

tools/MeshUtils/LoopSubdivision.cpp:8–21

using namespace PyMesh;

void LoopSubdivision::subdivide(
        MatrixFr vertices, MatrixIr faces, size_t num_iterations) {
    if (faces.cols() != 3) {
        throw NotImplementedError(
                "Only triangles mesh subdivision are supported!");
    }
    m_vertices = vertices;
    m_faces = faces;
    initialize_face_indices();

    for (size_t i=0; i<num_iterations; i++) {
        subdivide_once();
    }
}

void LoopSubdivision::subdivide_once() {
    extract_edges();