Re: [Libmesh-devel] [Libmesh-users] avoid inverse maps in shape reinit on sides

2011-08-01 Thread John Peterson
On Mon, Aug 1, 2011 at 1:45 PM, Roy Stogner wrote: > > John, double-check me on this: > > I believe the PYRAMID5 in Lorenzo's patch is wrong - the Pyramid5 we > use, going by both my own hazy memory and by the ::shape() equations > in fe_lagrange_shape_3D.C, has points 0-3 in the zeta=0 plane, not

Re: [Libmesh-devel] [Libmesh-users] avoid inverse maps in shape reinit on sides

2011-08-01 Thread Roy Stogner
On Wed, 27 Jul 2011, Lorenzo Alessio Botti wrote: > here is a patch to avoid inverse maps on sides reinit. > I hope that the main developers will consider it and > have some time to give it a look since it might break a lot of code... Aha, so it might. You're adding a new method that can make D