image.png

    1. #include <CGAL/Linear_cell_complex_for_combinatorial_map.h>
    2. #include <CGAL/draw_linear_cell_complex.h>
    4. typedef CGAL::Linear_cell_complex_for_combinatorial_map<3> LCC;
    5. typedef LCC::Dart_handle Dart_handle;
    6. typedef LCC::Point Point;
    8. int main()
    9. {
    10.     //创建一个线性单元复合体（即CGAL中的体模型）
    11.     LCC lcc;
    13.     //创建两个六面体，返回Dart句柄（即仿指针）
    14.     Dart_handle dh1 =
    15.         lcc.make_hexahedron(Point(0, 0, 0), Point(5, 0, 0),
    16.             Point(5, 5, 0), Point(0, 5, 0),
    17.             Point(0, 5, 4), Point(0, 0, 4),
    18.             Point(5, 0, 4), Point(5, 5, 4)); //8个点
    19.     Dart_handle dh2 =
    20.         lcc.make_hexahedron(Point(5, 0, 0), Point(10, 0, 0),
    21.             Point(10, 5, 0), Point(5, 5, 0),
    22.             Point(5, 5, 4), Point(5, 0, 4),
    23.             Point(10, 0, 4), Point(10, 5, 4));
    25.     //缝合
    26.     lcc.sew<3>(lcc.beta(dh1, 1, 1, 2), lcc.beta(dh2, 2));
    27.     //问题：beta是做什么的？
    29.     //体模型是否有效
    30.     lcc.display_characteristics(std::cout) << ", valid="
    31.         << lcc.is_valid() << std::endl;
    32.     //绘制
    33.     CGAL::draw(lcc);
    35.     return EXIT_SUCCESS;
    36. }