文件：linear_cell_complex_3_operations.cpp

    事先创建的三个六面体
    image.png image.png
    缝合结果
    image.pngimage.png
    插入的结果
    image.pngimage.png

    完整代码

    1. #include <CGAL/Linear_cell_complex_for_combinatorial_map.h>
    2. #include <CGAL/Linear_cell_complex_for_generalized_map.h>
    3. #include <vector>
    5. #ifdef _DEBUG
    6. #include <CGAL/draw_linear_cell_complex.h>
    7. #endif // _DEBUG
    10. typedef CGAL::Linear_cell_complex_for_combinatorial_map<3> LCC_3_cmap;
    11. typedef CGAL::Linear_cell_complex_for_generalized_map<3> LCC_3_gmap;
    13. template<typename LCC,
    14.     typename Map = typename LCC::Combinatorial_data_structure>
    15.     struct Alpha1
    16. {
    17.     static typename LCC::Dart_handle run(LCC&, typename LCC::Dart_handle dh)
    18.     {
    19.         return dh;
    20.     }
    21. };
    22. template<typename LCC>
    23. struct Alpha1<LCC, CGAL::Generalized_map_tag>
    24. {
    25.     static typename LCC::Dart_handle run(LCC& lcc, typename LCC::Dart_handle dh)
    26.     {
    27.         return lcc.template alpha<1>(dh);
    28.     }
    29. };
    31. template<typename LCC>
    32. void run_test()
    33. {
    34.     typedef typename LCC::Point Point;
    35.     typedef typename LCC::Dart_handle Dart_handle;
    37.     LCC lcc;
    39.     //创建三个六面体
    40.     Dart_handle dh1 = lcc.
    41.         make_hexahedron(Point(0, 0, 0), Point(1, 0, 0),
    42.             Point(1, 2, 0), Point(0, 2, 0),
    43.             Point(0, 3, 4), Point(0, 0, 4),
    44.             Point(6, 0, 4), Point(6, 3, 4));
    45.     Dart_handle dh2 = lcc.
    46.         make_hexahedron(Point(0, -5, 0), Point(2, -5, 0),
    47.             Point(2, -2, 0), Point(0, -2, 0),
    48.             Point(1, -1, 5), Point(1, -2, 5),
    49.             Point(5, -2, 5), Point(5, -2, 5));
    50.     Dart_handle dh3 = lcc.
    51.         make_hexahedron(Point(1, 0, 5), Point(0, 0, 6),
    52.             Point(0, 2, 5), Point(1, 2, 6),
    53.             Point(1, 3, 8), Point(0, 0, 8),
    54.             Point(5, 0, 9), Point(7, 3, 9));
    56. #ifdef _DEBUG
    57.     CGAL::draw(lcc);
    58. #endif // _DEBUG
    60.     //缝合
    61.     lcc.template sew<3>(dh1, lcc.other_orientation
    62.     (lcc.template opposite<2>
    63.         (lcc.next(lcc.next(lcc.template opposite<2>(dh2))))));
    64.     lcc.template sew<3>(lcc.template opposite<2>(lcc.next(dh1)),
    65.         lcc.other_orientation(lcc.template opposite<2>(lcc.previous(dh3))));
    67. #ifdef _DEBUG
    68.     CGAL::draw(lcc);
    69. #endif // _DEBUG
    71.     lcc.insert_cell_1_in_cell_2(lcc.next(dh1),
    72.         Alpha1<LCC>::run(lcc, lcc.previous(dh1)));
    73.     dh2 = lcc.template opposite<2>(lcc.next(lcc.next
    74.     (lcc.template opposite<2>(dh1))));
    75.     lcc.insert_cell_1_in_cell_2(dh2, Alpha1<LCC>::run
    76.     (lcc, lcc.next(lcc.next(dh2))));
    78.     std::vector<Dart_handle> path;
    79.     path.push_back(lcc.next(dh1));
    80.     path.push_back(lcc.next(lcc.template opposite<2>(lcc.previous(dh1))));
    81.     path.push_back(lcc.previous(dh2));
    82.     path.push_back(lcc.next(lcc.template opposite<2>(dh2)));
    83.     lcc.insert_cell_2_in_cell_3(path.begin(), path.end());
    85.     lcc.display_characteristics(std::cout) << ", valid="
    86.         << lcc.is_valid() << std::endl;
    88. #ifdef _DEBUG
    89.     CGAL::draw(lcc);
    90. #endif // _DEBUG
    91. }
    93. int main()
    94. {
    95.     //组合地图
    96.     run_test<LCC_3_cmap>();
    97.     //广义地图
    98.     //run_test<LCC_3_gmap>();
    100.     return EXIT_SUCCESS;
    101. }