Rahul,
I was responsible for starting the discussion regarding Gmsh meshes and
subdomain_ids() earlier this week :) . I confirm that things work with
Gmsh meshes.
You could try the following code to the example ex3.C, right after you call
mesh.print_info(); It essentially does what Derek men
Rahul,
What type of boundary condition (Dirichlet or Neumann)? Are you
planning on enforcing Dirichlet conditions using a penalty (as all of
the libMesh examples do)?
In general in libMesh we don't operate on _nodes_ for boundary
conditions we operate on sides. The reason for this is
Hi:
I would like to know how to tag nodes with different boundary ids for
possibly different sets of boundary conditions and loop over the
subset of elements that contain nodes of a particular boundary id. I
would also be able to do this in parallel. Could someone kindly tell
me how to do this?
T