I tried to set up u=g on the boundary for Stokes problem. I defined the DG FEspace as follows.
V = VectorL2(mesh,order=2, dirichlet ='left|right|top|bottom',dgjumps=True)
Q = L2(mesh,order=1)
X = FESpace([V,Q],dgjumps=True)
When I run the code, there is a error at the following line
[quote]gfu.components[0].Set(g, BND, definedon=mesh.Boundaries(“left|right|top|bottom”))
[/quote]
VectorL2FESpace does not have an evaluator for BND!
The code was attached.
Could you please show me how to set up the boundary condition for VectorL2FESpace or fix this error? Thank you so much.
Attachment: fem2.ipynb