Hello,

I have this question related to the two interesting specialcf for a very long time…

In particular, I want to check the specific normal direction of an edge in a triangular mesh.

For the mesh size function

h = specialcf.mesh_size,

I can plot it on a triangular mesh as a standard CoefficientFunction, and it looks like (should be) the area.

I guess, when restricting to the edges, it became the edge-length, right? But I can not verify this claim since I don’t know now to evaluate a Coefficient Function at an edge manually… Here is a related question concerning this issue: how to manually evaluate a function that is only defined on the edges? like the stabilization parameter for an HDG method, or its skeleton solution? The approach of MappedIntegrationPoints

mip = mesh(0.2,0.2)

f(mip)

seems not work for edge-wise evaluations…

Then, for the normal direction, the question is similar, I want to know n[0] and n[1], manually, at all the mesh skeleton…

Best regards,

Guosheng