From: Manav Bhatia <manav@u.washington.edu>  20051229 01:36:16

This is a critical point that you mentioned. I did go through the FE class in libMesh, and it appeared to me that it makes the assumtion of a 2D elem in a 2D space only, and like wise. For example, if I give FE a plane element with a nonzero zcoordinate, it seems to omit it for the calculation of shape functions, and derivatives. Hence, I always reduce my shell elems to a local 2D plane and then generate the matrices. However, the FE::inverse_map() method mentions something about NOT assuming that a lower dimensional elem lies in same dimensional elem itself. This makes it a little more confusing. Thanks Manav On Wed, 28 Dec 2005, Roy Stogner wrote: > On Wed, 28 Dec 2005, Manav Bhatia wrote: > >> I used the elem.contains_point() method to find if a point lies >> inside an element or not. For the following, I get true, which should >> otherwise be false. > > Does libMesh support lowerdimensional elements contained in > higherdimensional space? There's some support for "virtual" elements > generated as faces and edges of real elements, but in other cases I > think the code may assume that 2D elements are entirely in the XY > plane. > > I'll look through the code and see if there is such an assumption (and > see how hard it would be to remove); it would obviously be nice to > be able to use beam and shell elements in 3space. >  > Roy Stogner > 