CGAL users discussion list

[魏华祎 <>]

Hi, everyone,

      I am doing some research about mean curvature flow problem. There is a square domain, and in the domain, there exists a poly-line interface, the vertices of which will move along the normal direction of the interface.  After the moving, I need to  adjust the mesh in the domain carefully to keep the mesh quality,  which is important for the Finite Element  computation.  So I need to have a better mesh sizing function to control the adjust process ( namely I can do some optimization, refine or coarsening) . But for the 2d mesh generation of CGAL,  I can't find any example about how to implement  a mesh sizing function. CGAL just has following criteria for the 2d mesh generation:

CGAL::Delaunay_mesh_criteria 2<CDT> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . page 2894

CGAL::Delaunay_mesh_size_criteria_2<CDT> . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . page 2896

which just give the lower bound for the minimal angle and the the upper bound of the longest edge. 

Does anybody can give an example or hints about how to implement mesh sizing function for the 2d mesh generation? For example, I give a size value for the domain vertices and the vertices of the poly-line interface, how to construct a mesh sizing function defined on the whole domain. By the way, all the edge of the domain and the poly-line interface are set as constrained edges.

Best 

Huayi