CGAL users discussion list

[nolvis <>]

    Mesh out_mesh; // surface mesh data structure
    std::vector<Alpha_shape_3::Facet > Facets;
    alpha_shape.get_alpha_shape_facets(std::back_inserter(Facets),
Alpha_shape_3::REGULAR);
    for (auto i = 0; i < Facets.size(); i++)
    {

      //checks for exterior cells
      if (alpha_shape.classify(Facets[i].first) != Alpha_shape_3::EXTERIOR)
{
        Facets[i] = alpha_shape.mirror_facet(Facets[i]);
      }

      CGAL_assertion(alpha_shape.classify(Facets[i].first) ==
Alpha_shape_3::EXTERIOR);

      // gets indices of alpha shape and gets consistent orientation
      int indices[3] = { (Facets[i].second + 1) % 4, (Facets[i].second + 2)
% 4, (Facets[i].second + 3) % 4 };
      if (Facets[i].second % 2 == 0) {
        std::swap(indices[0], indices[1]);
      }

      // adds data to cgal mesh
      for (auto j = 0; j < 3; ++j) {
        out_mesh.add_vertex(Facets[i].first->vertex(indices[j])->point());
      }
      auto v0 = static_cast<vertex_descriptor>(3 * i);
      auto v1 = static_cast<vertex_descriptor>(3 * i + 1);
      auto v2 = static_cast<vertex_descriptor>(3 * i + 2);
      out_mesh.add_face(v0, v1, v2);
    }



--
Sent from: http://cgal-discuss.949826.n4.nabble.com/