std::vector<double> vertices;
    std::vector<int> indices;
    int index_count = 0;
    int triangle_count = 0;
    for (TopExp_Explorer explorer(shape, TopAbs_FACE); explorer.More(); explorer.Next()) {
        TopoDS_Face face = TopoDS::Face(explorer.Current());
        TopAbs_Orientation face_orientation = face.Orientation();
        TopLoc_Location location;

        Handle(Poly_Triangulation) triangulation = BRep_Tool::Triangulation(face, location);

        if (!triangulation.IsNull()) {

            int nb_nodes = triangulation->NbNodes();
            for (int node_index = 1; node_index <= nb_nodes; node_index++) {
                gp_Pnt node = triangulation->Node(node_index);

                vertices.push_back(node.X()); 
                vertices.push_back(node.Y()); 
                vertices.push_back(node.Z());
            }

            int nb_triangles = triangulation->NbTriangles();
            for (int triangle_index = 1; triangle_index <= nb_triangles; triangle_index++) {
                int i0, i1, i2;
                Poly_Triangle triangle = triangulation->Triangle(triangle_index);
                triangle.Get(i0, i1, i2);
                indices.push_back(i0 - 1 + index_count);
                indices.push_back(i1 - 1 + index_count);
                indices.push_back(i2 - 1 + index_count);
            }

            index_count += nb_nodes;
            triangle_count += nb_triangles;
        }
    }