I would like to try and write some code for approximating the curvature of the boundary of the mesh. The current solution proposed in the function “curvature” does not work too well and I would like to try other variants.

The normals computed in FreeFEM are normals to the edges of the boundary or normals approximated at the vertices on the boundary?