zero mean value L^2

Hi guys,

I was wondering how to implement L^2 polynomial spaces with zero mean value?(i.e. L^2_0(\Omega)) space

Thanks a lot!


Hello Rob,

one option would be to enforce the constraint via a Lagrange multiplier:

V =  L2(mesh, order = k)
N = NumberSpace(mesh)
X = FESpace([V,N])
(u,lam),(v,mu) = X.TnT()

and then add

u*mu + v*lam

to your bilinear form.

Besh wishes,

Hi Henry,

That works! Thanks a lot!

Best regards,