这个示例与《示例：基于广义地图的线性单元复合体》的相同

    1. #include <CGAL/Linear_cell_complex_for_combinatorial_map.h>
    2. #include <iostream>
    3. #include <algorithm>
    5. typedef CGAL::Linear_cell_complex_for_combinatorial_map<3> LCC_3;
    6. typedef LCC_3::Dart_handle                                 Dart_handle;
    7. typedef LCC_3::Point                                       Point;
    9. //一个仿函数，用于遍历一个体元的所有顶点
    10. // Functor used to display all the vertices of a given volume.
    11. template<class LCC>
    12. struct Display_vol_vertices : public CGAL::cpp98::unary_function<LCC, void>
    13. {
    14.     Display_vol_vertices(const LCC& alcc) :
    15.         lcc(alcc),
    16.         nb_volume(0)
    17.     {}
    19.     void operator() (typename LCC::Dart& d)
    20.     {
    21.         std::cout << "Volume " << ++nb_volume << " : ";
    22.         for (typename LCC::template One_dart_per_incident_cell_range<0, 3>::
    23.             const_iterator it = lcc.template one_dart_per_incident_cell<0, 3>
    24.             (lcc.dart_handle(d)).begin(),
    25.             itend = lcc.template one_dart_per_incident_cell<0, 3>
    26.             (lcc.dart_handle(d)).end();
    27.             it != itend; ++it)
    28.         {
    29.             std::cout << lcc.point(it) << "; ";
    30.         }
    31.         std::cout << std::endl;
    32.     }
    33. private:
    34.     const LCC& lcc;            //线性单元复合体（体模型）
    35.     unsigned int nb_volume;    //该体元对应的索引（第nb_volume+1个体元）
    36. };
    38. int main()
    39. {
    40.     LCC_3 lcc;
    42.     // Create two tetrahedra.
    43.     Dart_handle d1 = lcc.make_tetrahedron(Point(-1, 0, 0), Point(0, 2, 0),
    44.         Point(1, 0, 0), Point(1, 1, 2));
    45.     Dart_handle d2 = lcc.make_tetrahedron(Point(0, 2, -1),
    46.         Point(-1, 0, -1),
    47.         Point(1, 0, -1),
    48.         Point(1, 1, -3));
    50.     // Display all the vertices of the lcc by iterating on the
    51.     // Vertex_attribute container.
    52.     CGAL::IO::set_ascii_mode(std::cout);
    53.     std::cout << "Vertices: ";
    54.     for (LCC_3::Vertex_attribute_const_range::iterator
    55.         v = lcc.vertex_attributes().begin(),
    56.         vend = lcc.vertex_attributes().end();
    57.         v != vend; ++v)
    58.         std::cout << lcc.point_of_vertex_attribute(v) << "; ";
    59.     std::cout << std::endl;
    61.     // Display the vertices of each volume by iterating on darts.
    62.     std::for_each(lcc.one_dart_per_cell<3>().begin(),
    63.         lcc.one_dart_per_cell<3>().end(),
    64.         Display_vol_vertices<LCC_3>(lcc));
    66.     // 3-Sew the 2 tetrahedra along one facet
    67.     lcc.sew<3>(d1, d2);
    69.     // Display the vertices of each volume by iterating on darts.
    70.     std::for_each(lcc.one_dart_per_cell<3>().begin(),
    71.         lcc.one_dart_per_cell<3>().end(),
    72.         Display_vol_vertices<LCC_3>(lcc));
    74.     // Translate the second tetrahedra by a given vector
    75.     LCC_3::Vector v(3, 1, 1);
    76.     for (LCC_3::One_dart_per_incident_cell_range<0, 3>::iterator
    77.         it = lcc.one_dart_per_incident_cell<0, 3>(d2).begin(),
    78.         itend = lcc.one_dart_per_incident_cell<0, 3>(d2).end();
    79.         it != itend; ++it)
    80.     {
    81.         lcc.point(it) = LCC_3::Traits::Construct_translated_point_3()
    82.             (lcc.point(it), v);
    83.     }
    85.     // Display the vertices of each volume by iterating on darts.
    86.     std::for_each(lcc.one_dart_per_cell<3>().begin(),
    87.         lcc.one_dart_per_cell<3>().end(),
    88.         Display_vol_vertices<LCC_3>(lcc));
    90.     // We display the lcc characteristics.
    91.     std::cout << "LCC characteristics: ";
    92.     lcc.display_characteristics(std::cout) << ", valid=" << lcc.is_valid()
    93.         << std::endl;
    95.     return EXIT_SUCCESS;
    96. }