zero mean value L^2

5 years 3 months ago #1774 by Presley

zero mean value L^2 was created by Presley

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!

Rob

5 years 3 months ago - 5 years 3 months ago #1775 by hvwahl

Replied by hvwahl on topic zero mean value L^2

Hello Rob,

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

Code:

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

and then add

Code:

u*mu + v*lam

to your bilinear form.

Besh wishes,
Henry

Last edit: 5 years 3 months ago by hvwahl.

5 years 3 months ago #1787 by Presley

Replied by Presley on topic zero mean value L^2

Hi Henry,

That works! Thanks a lot!

Best regards,
Rob

Time to create page: 0.108 seconds