I’m trying to explore some stabilized FE schemes involving the integrals of gradient jumps over the mesh facets. This works fine in 2D with intalledges. I know that the 3D analogue intallfaces exists but it does not seem to work. Here is the minimal example :
mesh3 Th = cube(n, n, n);
This gives the obviously wrong answer 0 (in version 4.2.1)
If I continue with
macro dN(u) (dx(u)*N.x+dy(u)*N.y+dz(u)*N.z) //
solve test(u,v)=intallfaces(Th)(jump(dN(u))*jump(dN(v)) )
then it gives
Assertion fail : (0)
line :202, in file problem.cpp
Could you please correct this issue?
Moreover, it would be useful to create a 3D analogue of nTonEdge and to check if jump and mean work properly in 3D.
Thanks in advance,