On Wed, Apr 8, 2009 at 6:27 PM, David Knezevic <[email protected]> wrote: > You could do: > > AutoPtr<QBase> qrule = > fe_type.default_quadrature_rule(dim,extra_quad_order); > > where extra_quad_order is negative.
If you want to achieve "mass lumping" there are also Trapezoidal (QTrap) rules available. There's also Simpson's Rule (QSimpson) which may lump the quadratics...? I can't remember exactly how that works. -- John ------------------------------------------------------------------------------ This SF.net email is sponsored by: High Quality Requirements in a Collaborative Environment. Download a free trial of Rational Requirements Composer Now! http://p.sf.net/sfu/www-ibm-com _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
