dolfin team mailing list archive

[Joachim Berdal Haga <jobh@xxxxxxxxx>]

In the following code snipped, the bolded line (bc_p_drained) has no effect.
If I change V_p to V.sub(1), it works. The line above it (bc_u_bedrock), on
the other hand, works with either V_u or V.sub(0).

Is this a bug, or is only the .sub() notation correct?

====

V_u = VectorFunctionSpace(mesh, "CG", 2)
V_p = FunctionSpace(mesh, "CG", 1)
V = V_u * V_p

...

bc_u_bedrock = DirichletBC(V_u, [0]*dim, bot_boundary)
*bc_p_drained = DirichletBC(V_p, 0, top_boundary)
*
bcs = [bc_p_drained, bc_u_bedrock]
problem = VariationalProblem(a, L, bcs)