CGAL users discussion list

["gre6" <>]

hi, all.
I want to get the convex hull of a set of points. The number of the Points is about 10K. My code is like below:
--------------------------------------
...
typedef CGAL::Cartesian<double>                    K;
typedef K::Point_3                             Point_3;
typedef CGAL::Polyhedron_3< K>                 Polyhedron;
...
Polyhedron P; // define polyhedron to hold convex hull
convex_hull_incremental_3( pset.begin(), pset.end(), P); 
...
--------------------------------------------
It is compiled ok, but when run the programe it abort and output message like:

CGAL::Polyhedron_incremental_builder_3<HDS>::
lookup_halfedge(): input error: facet 1857 shares a halfedge from vertex 942 to vertex 881 with facet 1738.
CGAL error: assertion violation!
Expr: check_protocoll == 0
File: /usr/include/CGAL/Polyhedron_incremental_builder_3.h
Line: 198
Explanation:
terminate called after throwing an instance of 'CGAL::Assertion_exception'
  what():  CGAL ERROR: assertion violation!
Expr: check_protocoll == 0
File: /usr/include/CGAL/Polyhedron_incremental_builder_3.h
Line: 198
Aborted


When I reduced the number of points to 4k, I get the right result. I can't understand why it happen. Is there something wrong with the Kernel? 
Thank you for your attention!

===============================================
快来和我一起享受TOM免费邮箱吧！ 看看除了1.5G，还有什么？

看看女人都在聊什么? 美图尽赏 聊天世界

明星金曲免费送(http://mm.tom.com/ivr/)：周杰伦 林俊杰 庞龙 张惠妹

劲爆歌曲尽情点(http://mm.tom.com/ivr/)：霍元甲 吉祥三宝 人质 曹操

炫酷彩铃免费送(http://mm.tom.com/cailing/)：周杰伦帮你接电话 麻烦女朋友 七里香 小城故事
===============================================