On Thu, 22 Aug 2013, Praveen C wrote: > No, not just on the shared face, but in the interior also.
In that case, it can't be done even in theory with most elements. Interiors of elements shouldn't overlap, and most element shape functions aren't always well-defined outside their own interior. In particular, it's easy to come up with valid meshes where the master->physical element transformation's Jacobian is invertible on the element interior but not outside the element; in that case just trying to map quadrature points from one master element to the other would necessarily fail. The XYZ bases theoretically don't need to use the master->physical element transformation, so you could still hypothetically make this work by avoiding that entirely. You'd first have the easy-but-tedious task of writing your own compatibly ordered FEXYZ reimplementation, since the FE::reinit API in libMesh can take a list of xi-eta-zeta points but not a list of x-y-z points. More challenging, you'd also have to invent your own custom quadrature rule, since ours are defined on the master elements and then transformed to physical space. --- Roy ------------------------------------------------------------------------------ Introducing Performance Central, a new site from SourceForge and AppDynamics. Performance Central is your source for news, insights, analysis and resources for efficient Application Performance Management. Visit us today! http://pubads.g.doubleclick.net/gampad/clk?id=48897511&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
