Hi,
for mesh refinement i need a good representation of curved boundaries. Refining the mesh, however does not give a better approximation of the geometry. In case that the curvature is known one can define a subdomain and use:
mesh.snap_boundary(SubDomain)
Now this fails in case that the boundary is not fully known before. A possible (very dirty) workaround would be to make a very dense BoundaryMesh, define a subdomain from this boundarymesh and use it for the vertex correction with snap_boundary. Is there also a more elegant solution for arbitrary geometries ? One way would be if it was possible to define a SubDomain from the geometry itself and not the mesh e.g. using any CAD format.
This is related to the Demo: demo_smoothing.py
Kind regards
Johann
