On Tue, 8 Nov 2011, Kirk, Benjamin (JSC-EG311) 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 self-consistent 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 FIN-S regression test results are still coming in under the 1e-15 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 ------------------------------------------------------------------------------ RSA(R) Conference 2012 Save $700 by Nov 18 Register now http://p.sf.net/sfu/rsa-sfdev2dev1 _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
