Method poke_holes

tools/Triangle/TriangleWrapper.cpp:424–456

}

void TriangleWrapper::poke_holes(Regions& regions) {
    std::list<size_t> interior_faces;
    for (auto& region : regions) {
        VectorF seed_p;
        bool seed_found = select_seed_point(region, seed_p);

        if (!seed_found) {
            // All faces are degenerated in this region.
            continue;
        }

        Float wind_num = compute_winding_number(seed_p, m_points, m_segments);
        if (fabs(wind_num) > 0.5) {
            interior_faces.splice(interior_faces.end(), region);
        }
    }

    MatrixIr faces(interior_faces.size(), 3);
    VectorI interior_regions(interior_faces.size());
    size_t count=0;
    for (const auto& f_index : interior_faces) {
        faces.row(count) = m_faces.row(f_index);
        interior_regions[count] = m_regions[f_index];
        count++;
    }

    IsolatedVertexRemoval remover(m_vertices, faces);
    remover.run();
    m_vertices = remover.get_vertices();
    m_faces = remover.get_faces();

    m_regions.swap(interior_regions);
}

bool TriangleWrapper::select_seed_point(const Region& region, VectorF& seed_p) {
    const Float degeneracy_tol = std::min(m_max_area * 0.01, 1e-6);