我正在使用CGAL在三角形中生成2D非结构化网格。代码如下所示。
typedef CGAL::Constrained_Delaunay_triangulation_2<K, Tds> CDT;
typedef CDT::Vertex_handle Vertex_handle;
CDT cdt;
Vertex_handle va = cdt.insert(Point(-4,0));
Vertex_handle vb = cdt.insert(Point(0,-1));
Vertex_handle vc = cdt.insert(Point(4,0));
cdt.insert_constraint(va, vb);
cdt.insert_constraint(vb, vc);
cdt.insert_constraint(vc, va);
CGAL::refine_Delaunay_mesh_2(cdt, Criteria(0.125, 0.5));
乍一看,网状物很好。但是,当我放大图形时,我发现还有很多额外的边缘:
我的代码出了什么问题?非常感谢。