[MOAB-dev] Possible bug in 2-D quadratic elements

[Jason Kraftcheck]

Wollaber, Allan B. wrote:
> Jason (Tim),
> 
> The connectivity indices are now the correct length for all the elements
> that I've tried (HEX8/20/27, QUAD4/8/9, TRI3/6).  This is great!
> 
> I have a new question now, however: At one point in the mesh processing
> I must determine the local reference number of a facet (in a sideset
> with the NEUMANN_SET tag) with respect to the element it is bounding. 
> To do this, I use the MBCN_SideNumber routine.  This has worked well for
> every 3-D element that I've tried (hexes, tets) and for the linear 2-D
> elements (tri3, quad4), but not for the quadratic 2-D elements (tri6,
> quad8/9), even before the change that was just made to the repository. 
> For the quadratic 2-D elements, the local facet number is always
> returned as "1".  The following code snippet should help illuminate how
> I am doing this:
> 

The 'child_num_verts' argument to SideNumber should be the number of
corners, not the number of nodes.  The function works the same regardless of
the number of nodes in the element.

- jason