我正在尝试使用Exact_circular_kernel_2
构建圆圈排列。当我在排列中插入圆时,出现分割错误。
这是我的代码:
CGAL::Circle_2< CGAL::Exact_circular_kernel_2 > circle1( CGAL::Point_2< CGAL::Exact_circular_kernel_2 >(1,1), 2 );
CGAL::Exact_circular_kernel_2::Circle_2 circle2 = circle1;
CGAL::Circular_arc_2< CGAL::Exact_circular_kernel_2 > arc( circle2 );
CGAL::Arr_circular_line_arc_traits_2< CGAL::Exact_circular_kernel_2 >::Curve_2 curve = arc;
CGAL::Arrangement_2< CGAL::Arr_circular_line_arc_traits_2< CGAL::Exact_circular_kernel_2 > > myArrangement;
cout<<curve<<endl;
cout<<myArrangement.is_valid()<<endl;
insert ( myArrangement, curve );
cout<<"done"<<endl;
代码编译时没有警告,无论我在第 3 行中使用 circle1
还是circle2
,结果都是一样的。
程序在分段错误发生之前打印以下内容:
1/1 1/1 2/1 1 0 1/1 -1/1 2/1 1 1/1 0 1/1 -1/1 2/1 1 1/1
1
因此,myArrangement
和curve
的值似乎是有效的。
知道我做错了什么吗?
我知道我可以使用另一个内核,但我想比较不同内核的性能,所以不使用此内核并不是真正的解决方案。
我无法重现该问题。
下面是一个稍微干净的版本,但你的版本也适用于我的平台,Ubuntu 12.04,g++ 4.6.3,CGAL 4.1(或至少接近4.1),你的是什么?
#include <iostream>
#include <CGAL/Exact_circular_kernel_2.h>
#include <CGAL/Arr_circular_line_arc_traits_2.h>
#include <CGAL/Arrangement_2.h>
typedef CGAL::Exact_circular_kernel_2 Kernel;
typedef Kernel::Point_2 Point_2;
typedef Kernel::Circle_2 Circle_2;
typedef CGAL::Arr_circular_line_arc_traits_2<Kernel> Traits;
typedef Traits::Curve_2 Curve_2;
typedef CGAL::Arrangement_2<Traits> Arrangement;
int main()
{
Circle_2 circle(Point_2(1,1), 2);
Curve_2 curve(circle);
Arrangement arr;
std::cout << curve << std::endl;
std::cout << arr.is_valid() << std::endl;
CGAL::insert(arr, curve);
std::cout << "done" << std::endl;
return 0;
}