Function raw_to_bsp

tools/Boolean/BSP/BSPEngine.cpp:20–49

    typedef boost::shared_ptr<BSPlib::Bsp> BSPPtr;

    BSPPtr raw_to_bsp(const MatrixFr& vertices, const MatrixIr& faces) {
        const size_t num_vertices = vertices.rows();
        const size_t num_faces = faces.rows();
        assert(vertices.cols() == 3);
        assert(faces.cols() == 3);

        PLYmesh mesh;
        mesh.vertices.resize(num_vertices * 3);
        std::copy(vertices.data(), vertices.data() + num_vertices*3,
                mesh.vertices.begin());

        mesh.faces.resize(num_faces);
        for (size_t i=0; i<num_faces; i++) {
            mesh.faces[i] = {faces(i,0), faces(i,1), faces(i,2)};
        }
        BRep brep;
        Ply2Convex(mesh, brep);
        {
            auto r = BSPlib::Bsp::BuildTree(brep);
            auto convex_polys = r->GetBoundary();
            std::vector<BSPlib::MetricPolygon> metric_polys(convex_polys.size());
            std::transform(convex_polys.begin(), convex_polys.end(),
                    metric_polys.begin(), [](const BSPlib::ConvexPolygon& p) {
                    return BSPlib::Approximate(p); });
            PLYmesh tmp;
            Metric2Ply(metric_polys, tmp);
            Write("tmp.ply", tmp);
        }
        return BSPlib::Bsp::BuildTree(brep);
    }

    void bsp_to_raw(const BSPPtr bsp, MatrixFr& vertices, MatrixIr& faces) {
        const auto convex_polys = bsp->GetBoundary();