1

我正在尝试使用Exact_circular_kernel_2. 当我在 Arrangement 中插入一个 Circle 时,出现分段错误。
这是我的代码:

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;

代码在没有警告的情况下编译,无论我使用circle1还是circle2在第 3 行,结果都是一样的。

该程序在分段错误发生之前打印以下内容:

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似乎是有效的。

任何想法我做错了什么?

我知道我可以使用另一个内核,但我想比较不同内核的性能,所以不使用这个内核并不是一个真正的解决方案。

4

1 回答 1

3

我无法重现该问题。

下面是一个稍微干净一点的版本,但你的版本也适用于我的平台,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;
}
于 2012-12-19T09:58:36.323 回答