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)
_______________________________________________ Mailing list: https://launchpad.net/~dolfin Post to : [email protected] Unsubscribe : https://launchpad.net/~dolfin More help : https://help.launchpad.net/ListHelp
