From: David Knezevic <dave.knez@gm...>  20080814 10:46:35

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 difference since the Gauss rules for 2*order and 2*order+1 are the same, but it can make a difference on simplices, e.g. the order 4 and order 5 "Gauss" rules on triangles are different... 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? Cheers, Dave 