Default quadrature rule

General Discussion
1

Hello,

I have a question about the default choice of quadrature in FreeFEM. Does FreeFEM choose the default quadrature rule to be of highest degree that is exact for the corresponding finite element function? Or do we need to change it by hand looking at the documentation for P2, P3 elements?

Thank you, Beniamin

2

You need to set it appropriately. E.g., if you compute int2d(Th)(u*v) with u and v of type P4, you’ll get incorrect results, unless you use the appropriate qforder.

3

Thank you for your answer! I think I got it now.