Hi,
I'm working with XFEM and I would assume, that the boundary conditions should also hold for u_x (not only for u_std).
With
Code:
Vh = H1(mesh, order=1, dim=mesh.dim, dirichlet="left")
Vhx = XFESpace(Vh,ci)
VhG = FESpace([Vh,Vhx])
freedofs = VhG.FreeDofs()
print(freedofs)
I get
0: 0110111111111000111111111111111111111
but when clearing the dofs by hand
Code:
VhGx, VhGy = VhG.components
for w in mesh.vertices:
if abs(w.point[0])<0.00001:
for dof in VhG.GetDofNrs(w):
freedofs.Clear(dof)
freedofs.Clear(VhGx.ndof+dof)
print(freedofs)
also the dofs in the u_x part are also cleared:
0: 0110111111111000111111111110010011111
Does XFEM work like this or is this maybe a bug?
BR,
Nepomuk