Frontal-Delaunay for quads with size callback
I have had success meshing in 2D with a size callback (set from the C++ API with gmsh::model::mesh::setSizeCallback
) + the Frontal-Delaunay algorithm, but when I try to use the Frontal-Delaunay for quads some of the surfaces in my model will not be meshed according to the size field. It appears to use the size field for most of the surfaces, and it's just a few, often exactly rectangular, which are "missed" and have no new nodes inserted. See below for an example, the curves are meshed according to the size field but the mesh size is not correct in the middle surface:
For comparison, this is what it looks like with the standard Frontal-Delaunay algorithm and same size field:
I'm wondering if this is a known restriction or bug, or if I need to specify something else before generating the mesh? Reading the paper it seems like this shouldn't be an issue, and many models are meshed just fine.
Thanks so much!