> 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.
Yeah, n
> 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
I was wondering why default_quadrature_order in fe_type returns
2*order+1? It says in the documentation that the idea is that the
default quadrature rule integrates the mass matrix exactly, so shouldn't
2*order be sufficient?
Using tensor product Gauss quadrature rules, the "+1" doesn't make a
