On Mon, Mar 16, 2009 at 3:42 PM, Roy Stogner <roystgnr@...> wrote:
> On Mon, 16 Mar 2009, Andrea Hawkins wrote:
>> In my code, I'm currently just generating my mesh using
>> MeshTools::Generation::build_square (mesh, x_elem, y_elem, min, max, min,
>> max, QUAD9);
>> In this routine, I see it setting boundary ids for the sides, but I don't
>> see anywhere it sets them (or would set them)
>> for the nodes. Where would be the appropriate place to set this? I'm just
>> wanting the entire boundary to be Dirichlet.
> I seem to recall that, when adding a boundary id for a side, the nodes
> for that side are given the same id, possibly overwriting any existing
> id of their own.
> Searching through the library, I can find absolutely no hint that any
> such behavior exists.
Agree, looks like we have definitely been lazy about adding nodes to
the boundary during build_cube. I think probably the min of all
possible boundary IDs of sides which contain the node might be a
logical choice. That shouldn't be too hard to add to build_cube, not
that I'm volunteering of course :-)
> I'd say that it might be a useful idea to at
> least add a method like "BoundaryInfo::sync_sides_to_nodes()"
Agree again, because I'm not sure the boundary info will be correctly
updated when MeshTools::renumber_nodes_and_elements() is called...
> at the
> moment I'm more concerned about how I hallucinated a memory of code
> that never existed...
Maybe you are thinking about the way that child elements' boundary IDs
are automatically gotten from their parents' boundary IDs w/o having
to add them to the BoundaryInfo object?