From: Roy Stogner <roystgnr@ic...> - 2005-05-11 19:54:04
Currently, the projection-based constraints calculation makes the
following assumption: when the side of a child element and the face of
its parent overlap, the side numbers (i.e. the indices to pass to
build_side) of the overlapping sides are equal.
Is this true for every element? I've only tested a couple personally.
With every element except Tet10, a child whose sides overlap its
parent's sides will be an affine transformation of its parent, and so
giving it the same (local) node numbering seems natural. I'd rather
not pore over a dozen _embedding_matrix arrays to verify that, though,
if someone knows the answer offhand.
I think I've asked this before, but if so I can't find it in the
libmesh-devel archives. In any case I'm about to repeat this
assumption in some more code, so reassurance would be appreciated.