From: Roy Stogner <roystgnr@ic...>  20111109 02:16:46

On Tue, 8 Nov 2011, Kirk, Benjamin (JSCEG311) wrote: > None here. I was pretty much sequestered today  what's the gist? The short version is that it's a miracle I manage to dress myself and feed myself despite gross mental handicaps. The long version is that my "has_affine_map" optimization, which obtains the correct init_shape_functions result much faster on triangles and tets, obtains incorrect dphidxi etc. maps on rectangluar quads, rectangular prism hexes, and rectangular prisms. As far as I can tell, the only reason we managed to avoid noticing the problem for so long is that for most codes the errors can act as a selfconsistent weighting function in the weak formulation integrals, so e.g. manufactured solution tests still converge at the right order and exact solutions in the FEM space are still recovered. After the change all our FINS regression test results are still coming in under the 1e15 tolerance, even! I think Boyce only managed to catch things by using different order quadrature rules to integrate the matrix vs forcing function for an L2 projection.  Roy 