Function triangulate

tools/Wires/Inflator/PeriodicBoundaryRemesher.cpp:29–48

    const short max_axis_marker[3] = { 1,  2,  3};

    void triangulate(MatrixFr vertices, MatrixIr edges,
            MatrixFr& output_vertices, MatrixIr& output_faces, Float max_area) {
        assert(edges.rows() >= 3);
        MeshCleaner cleaner;
        cleaner.remove_isolated_vertices(vertices, edges);
        cleaner.remove_duplicated_vertices(vertices, edges, 1e-12);
        assert(vertices.rows() >= 3);

        TriangleWrapper triangle;
        triangle.set_points(vertices);
        triangle.set_segments(edges);
        triangle.set_max_area(max_area);
        triangle.set_split_boundary(false);
        triangle.set_auto_hole_detection(true);
        triangle.set_verbosity(0);
        triangle.run();

        output_vertices = triangle.get_vertices();
        output_faces = triangle.get_faces();
    }

    void reorientate_triangles(const MatrixFr& vertices, MatrixIr& faces,
            const VectorF& n) {