Hi all, In my continuing quadrature rule "audit", I came across another interesting one.
The current second-order (exact for quadratics) Gauss rule for triangles has points on the *boundaries* of the reference element. Roy pointed out to me a couple days ago that these types of rules are obviously not ideal for problems with discontinuous material coefficients along which elements have been aligned. Luckily, there is an alternate 3-point, second-order rule given by Flaherty with interior points. Unless there are objections I think I'm going to replace the current second-order rule with this one: points: (2/3, 1/6) (1/6, 2/3) (1/6, 1/6) weights: 1/6 for all -- John ------------------------------------------------------------------------- Check out the new SourceForge.net Marketplace. It's the best place to buy or sell services for just about anything Open Source. http://sourceforge.net/services/buy/index.php _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
