Issues with Defining Boundary Conditions for Stokes Problem Using GMSH Mesh Tags

[dokken]

@Mehran1

boundary_indices = fem.locate_dofs_topological(V, fdim, facet_markers.indices[facet_markers.values == boundary_numbers['inlet']])
bcs.append(fem.dirichletbc(velocity_bc_value, boundary_indices, V))
boundary_indices = fem.locate_dofs_topological(V, fdim, facet_markers.indices[facet_markers.values == boundary_numbers['walls']])
bcs.append(fem.dirichletbc(noslip_bc_value, boundary_indices, V))

boundary_indices = fem.locate_dofs_topological(Q, fdim, facet_markers.indices[facet_markers.values == boundary_numbers['outlet']])
bcs.append(fem.dirichletbc(pressure_bc_value, boundary_indices, Q))

These are wrong, see: Mixed formulation of stokes equations - #6 by dokken