Local mesh size of a facet element in 3D mesh

Hi all,

I am wondering if there is any simple way to get the local mesh size of a triangular facet element of a 3D mesh, or the simplest way is just to write some code to calculate it by the vertices?

Thanks in advance for any possible help!

the CoefficientFunction specialcf.mesh_size allows to evaluate the local mesh-size. You can interpolate it into a GridFunction having one value per facet:

gf_ms = GridFunction(FacetFESpace(mesh, order=0))
gf_ms.Set(specialcf.mesh_size, dual=True)
1 Like