> Actually I just had a thought that on quadratic isoparametric elements > the Jacobian will be a linear function, so is the idea of the 2*order+1 > to cover this case?
That was my original thought. Even on distorted bilinear quadrilaterals the Jacobian will be non-constant. This predates Roy's has_affine_map() stuff -- presumably you could use that information to make an intelligent decision about the order needed for simplices? The problem there is that the default quadrature order is set independent of an individual element, so you might only take advantage of this optimization when all of the elements are affine... Then again, the default is just that - a default. You can override it if you want, and in some cases, for some problems, there is certainly a performance benefit to under-integrating. This is possible in my compressible navier-stokes stuff so long as you lump the mass matrix. -Ben ------------------------------------------------------------------------- This SF.Net email is sponsored by the Moblin Your Move Developer's challenge Build the coolest Linux based applications with Moblin SDK & win great prizes Grand prize is a trip for two to an Open Source event anywhere in the world http://moblin-contest.org/redirect.php?banner_id=100&url=/ _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
