The surface mixed FEM now works as expected. I want to hybridize the surface mixed FEM, but the discontinuous flag seems not working properly for surface HDiv: I am expecting with the discontinuous flag, the space will be co-normal discontinuous across surface edges, but this is not the case in the attached code
Code:
from ngsolve.meshes import MakeStructured3DMesh
from ngsolve import *
mesh = MakeStructured3DMesh(nx=1, ny=1, nz=1)
VC = Compress(HDivSurface(mesh, order=0, RT=True, discontinuous=False, definedon=mesh.Boundaries("top")))
print(VC.ndof, " free: ", VC.FreeDofs())
VD = Compress(HDivSurface(mesh, order=0, RT=True, discontinuous=True, definedon=mesh.Boundaries("top")))
print(VD.ndof, " free: ", VD.FreeDofs())
