Method run_triangle

tools/Triangle/TriangleWrapper.cpp:263–422

}

void TriangleWrapper::run_triangle(const std::string& flags) {
    const size_t num_points = m_points.rows();
    const size_t dim = m_points.cols();
    const size_t num_segments = m_segments.rows();
    const size_t pt_per_segment = m_segments.cols();
    const size_t num_holes = m_holes.rows();
    const size_t num_triangles = m_triangles.rows();
    const size_t pt_per_triangle = m_triangles.cols();
    const size_t num_areas = m_max_areas.size();

    if (num_segments > 0 && pt_per_segment != 2) {
        throw RuntimeError("Segments must consist of 2 points!");
    }
    if (num_triangles > 0 && pt_per_triangle != 3) {
        throw RuntimeError("Triangles must consists of 3 points!");
    }
    if (dim != 2) {
        throw RuntimeError("Triangle only works in 2D domain.");
    }

    auto generate_empty_triangulateio = []() {
        triangulateio io;
        io.numberofpoints = 0;
        io.numberofpointattributes = 0;
        io.numberoftriangles = 0;
        io.numberoftriangleattributes = 0;
        io.numberofcorners = 0;
        io.numberofsegments = 0;
        io.numberofholes = 0;
        io.numberofregions = 0;
        io.numberofedges = 0;
        io.pointlist = nullptr;
        io.pointattributelist = nullptr;
        io.pointmarkerlist = nullptr;
        io.trianglelist = nullptr;
        io.triangleattributelist = nullptr;
        io.trianglearealist = nullptr;
        io.neighborlist = nullptr;
        io.segmentlist = nullptr;
        io.segmentmarkerlist = nullptr;
        io.edgelist = nullptr;
        io.edgemarkerlist = nullptr;
        io.holelist = nullptr;
        io.regionlist = nullptr;
        io.normlist = nullptr;
        return io;
    };

    triangulateio in       = generate_empty_triangulateio();
    triangulateio out      = generate_empty_triangulateio();
    triangulateio out_voro = generate_empty_triangulateio();

    in.numberofpoints          = num_points;
    in.numberofpointattributes = 0;
    in.pointlist               = new REAL[num_points * dim];
    std::copy(m_points.data(), m_points.data() + num_points * dim,
            in.pointlist);