[FEM] Introduce higher accuracy quadrature rules

[th-skam]

For each Element, an option is added to use a higher-accuracy quadrature rule. The default behaviour of existing force fields is preserved. Lower-order rules are also included for completeness.

The purpose of this change is to enable more accurate error-norm computations in the V&V suite using the Method of Manufactured Solutions.

If the accuracy is insufficient, the integration introduces errors in the computed norms. This we observe either as the convergence rate approaching the expected one more slowly or as apparent superconvergence that counters the purpose of the test.

The degrees are:

Element	Up to degree	Default degree
Edge	3	1
Quad	3	2
Tri	2	1
Tetra	2	1
Hex	5	3

---

By submitting this pull request, I acknowledge that
I have read, understand, and agree SOFA Developer Certificate of Origin (DCO).

---

Reviewers will merge this pull-request only if

• it builds with SUCCESS for all platforms on the CI.
• it does not generate new warnings.
• it does not generate new unit test failures.
• it does not generate new scene test failures.
• it does not break API compatibility.
• it is more than 1 week old (or has fast-merge label).