1. //自定义Point
    2. #include <CGAL/Linear_cell_complex_for_combinatorial_map.h>
    3. #include <iostream>
    4. #include <algorithm>
    6. template<typename K>
    7. struct mypoint : public K::Point_3
    8. {
    9.   typedef typename K::Point_3 Base;
    11.   mypoint() : mtype('a')
    12.   {}
    14.   mypoint(const Base& apoint) : Base(apoint),
    15.                                 mtype('b')
    16.   {}
    18.   mypoint(const mypoint& apoint) : Base(apoint),
    19.                                    mtype('c')
    20.   {}
    22.   mypoint(typename K::FT a, typename K::FT b, typename K::FT c) : Base(a,b,c),
    23.                                                                   mtype('d')
    24.   {}
    26.   char type() const
    27.   { return mtype; }
    29. private:
    30.   char mtype;
    31. };
    33. template<typename K>
    34. struct mytraits: public CGAL::Linear_cell_complex_traits<3, K>
    35. {
    36.   typedef mypoint<K> Point;
    37. };
    39. typedef mytraits<CGAL::Exact_predicates_inexact_constructions_kernel> Traits;
    41. typedef CGAL::Linear_cell_complex_for_combinatorial_map<3,3, Traits> LCC_3;
    42. typedef LCC_3::Dart_handle                     Dart_handle;
    43. typedef LCC_3::Point                           Point;
    44. typedef LCC_3::FT                              FT;
    46. //创建正方体
    47. Dart_handle make_iso_cuboid(LCC_3& lcc, const Point& basepoint, FT lg)
    48. {
    49.   return lcc.make_hexahedron(basepoint,
    50.                              Traits::Construct_translated_point()
    51.                              (basepoint,Traits::Vector(lg,0,0)),
    52.                              Traits::Construct_translated_point()
    53.                              (basepoint,Traits::Vector(lg,lg,0)),
    54.                              Traits::Construct_translated_point()
    55.                              (basepoint,Traits::Vector(0,lg,0)),
    56.                              Traits::Construct_translated_point()
    57.                              (basepoint,Traits::Vector(0,lg,lg)),
    58.                              Traits::Construct_translated_point()
    59.                              (basepoint,Traits::Vector(0,0,lg)),
    60.                              Traits::Construct_translated_point()
    61.                              (basepoint,Traits::Vector(lg,0,lg)),
    62.                              Traits::Construct_translated_point()
    63.                              (basepoint,Traits::Vector(lg,lg,lg)));
    64. }
    66. int main()
    67. {
    68.   LCC_3 lcc;
    70.   //创建两个正方体
    71.   // Create two iso_cuboids.
    72.   Dart_handle d1 = make_iso_cuboid(lcc, Point(-2, 0, 0), 1);
    73.   Dart_handle d2 = make_iso_cuboid(lcc, Point(0, 0, 0), 1);
    75.   //缝合两个正方体
    76.   // 3-Sew the two cubes along one facet.
    77.   lcc.sew<3>(lcc.beta(d1, 1, 1, 2), lcc.beta(d2, 2));
    79.   //将两个Cube相邻面做Barycentric三角剖分
    80.   // Barycentric triangulation of the facet between the two cubes.
    81.   lcc.insert_barycenter_in_cell<2>(lcc.beta(d2, 2));
    83.   // Display all the vertices of the map.
    84.   for (LCC_3::Vertex_attribute_range::iterator
    85.          it=lcc.vertex_attributes().begin(),
    86.          itend=lcc.vertex_attributes().end();
    87.        it!=itend; ++it)
    88.   {
    89.     std::cout<<"point: "<<lcc.point_of_vertex_attribute(it)
    90.              <<", "<<"type: "<<lcc.point_of_vertex_attribute(it).type()
    91.              <<std::endl;
    92.   }
    94.   return EXIT_SUCCESS;
    95. }