Method run

tools/TetGen/TetgenWrapper.cpp:15–111

using namespace PyMesh;

void TetgenWrapper::run() {
    tetgenio in, out, addin, bgmin;

    const size_t num_points = m_points.rows();
    const size_t num_triangles = m_triangles.rows();
    const size_t num_tets = m_tets.rows();
    const size_t num_point_markers = m_point_markers.size();
    const size_t num_point_weights = m_point_weights.size();
    const size_t num_triangle_markers = m_triangle_markers.size();

    in.firstnumber = 0;
    in.mesh_dim = 3;

    in.numberofpoints = m_points.rows();
    in.pointlist = new REAL[in.numberofpoints * 3];
    std::copy(m_points.data(), m_points.data() + m_points.size(),
            in.pointlist);

    if (num_point_markers > 0) {
        if (num_point_markers != num_points) {
            throw RuntimeError(
                    "Size of point marker must equal to the number of points");
        }
        in.pointmarkerlist = new int[num_points];
        std::copy(m_point_markers.data(), m_point_markers.data() + num_points,
                in.pointmarkerlist);
    }

    if (num_point_weights > 0) {
        if (num_point_weights != num_points) {
            throw RuntimeError(
                    "Size of point weights must equal to the nubmer of points");
        }
        in.numberofpointattributes = num_points;
        in.pointattributelist = new REAL[num_points];
        std::copy(m_point_weights.data(), m_point_weights.data() + num_points,
                in.pointattributelist);
    }

    if (num_triangles > 0) {
        if (m_triangles.cols() != 3) {
            throw NotImplementedError(
                    "Only triangles are supported as polygon input to tetgen");
        }
        in.numberoffacets = num_triangles;
        in.facetlist = new tetgenio::facet[num_triangles];
        for (size_t i=0; i<num_triangles; i++) {
            auto& f = in.facetlist[i];
            f.numberofpolygons = 1;
            f.polygonlist = new tetgenio::polygon[1];
            f.numberofholes = 0;
            f.holelist = NULL;
            auto& poly = f.polygonlist[0];
            poly.numberofvertices = 3;
            poly.vertexlist = new int[3];
            poly.vertexlist[0] = m_triangles(i, 0);
            poly.vertexlist[1] = m_triangles(i, 1);
            poly.vertexlist[2] = m_triangles(i, 2);