用我自己的Point类自定义CGAL内核

我想在CGAL约束delaunay三角测量中使用自定义Point类。但是，使用以下MyPoint类（它应该与CGAL :: Point_2的行为完全相同），我会遇到分段错误。如果我将MyKernel中的Point_2 typedef设置为CGAL :: Exact_predicates_inexact_constructions_kernel :: Point_2，它将完美工作。我究竟做错了什么？它出现segfaults在最后一行用我自己的Point类自定义CGAL内核

template<class P> 
struct MyPoint : public P { 
    MyPoint() : P() {} 
    MyPoint(const MyPoint& p) : P(p) {} 

    MyPoint(int x, int y) : P(x,y) {} 
    MyPoint(double x, double y) : P(x,y) {} 
}; 

struct MyKernel : CGAL::Exact_predicates_inexact_constructions_kernel { 
    typedef MyPoint<CGAL::Exact_predicates_inexact_constructions_kernel::Point_2> Point_2; 
}; 

typedef MyKernel K; 

typedef CGAL::Triangulation_vertex_base_2<K>      Vb; 
typedef CGAL::Constrained_triangulation_face_base_2<K>   Fb; 
typedef CGAL::Triangulation_data_structure_2<Vb,Fb>    TDS; 
typedef CGAL::Exact_predicates_tag        Itag; 
typedef CGAL::Constrained_Delaunay_triangulation_2<K, TDS, Itag> CDT; 
typedef CDT::Point   Point;

代码：改变内核

CDT cdt; 
Point cgal_p1; 
Point cgal_p2; 
cgal_p1 = Point(p1[1],p1[2]); 
cgal_p2 = Point(p2[1],p2[2]); 
cdt.insert_constraint(cgal_p1, 
         cgal_p2);

来源

2010-03-10 MHOOO

是棘手的。这里发生的事情是，从三角测量类看，与内核相比，唯一改变的是点类型，通过派生它。对于谓词函子，这可能是足够好的，但对于构造函子不适用，比如CDT所需的交集。

我认为你有两个选择：

写出ConstrainedTriangulationTraits_2您的类型充分性状接口。
使用Extensible Kernel机制应该做你想做的。

来源

2012-11-20 05:33:09

用我自己的Point类自定义CGAL内核

回答

相关问题