From: KIRK, BENJAMIN (JSCEG) (NASA) <benjamin.kirk1@na...>  20050511 21:23:11

That is true. When the intersection of a parent side and child side is not empty then the two sides are the same number. This assumption is used in the BoundaryInfo class. If you have a child and you want to know if side s has a BC associated with it then you ask the parent if they have a BC for side s, and so on recursively. That way the BoundaryInfo class is not affected by refinement. The only element where this was awkward is the Pyramid. Naturally, you might want the square base to be side 0, but that won't work because the 4 triangular faces need to be sides 0,1,2,3 so this assumption will work for the tetrahedral children. So, the square base is side 4. Ben Original Message From: libmeshdeveladmin@... [mailto:libmeshdeveladmin@...] On Behalf Of Roy Stogner Sent: Wednesday, May 11, 2005 2:54 PM To: libmeshdevel@... Subject: [Libmeshdevel] Refinement assumptions Currently, the projectionbased 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 libmeshdevel archives. In any case I'm about to repeat this assumption in some more code, so reassurance would be appreciated.  Roy  This SF.Net email is sponsored by Oracle Space Sweepstakes Want to be the first software developer in space? Enter now for the Oracle Space Sweepstakes! http://ads.osdn.com/?ad_id=7393&alloc_id=16281&op=click _______________________________________________ Libmeshdevel mailing list Libmeshdevel@... https://lists.sourceforge.net/lists/listinfo/libmeshdevel 