* contributions of the faces (which are the integrated square of the jumps
* times some factor) of each cell and take the square root.
*
- * The integration is done using a quadrature formula on the face. For linear
- * trial functions (FEQ1), QGauss (with two points) or even the QMidpoint rule
- * will suffice. For higher order elements, it is necessary to utilize higher
- * order quadrature formulae as well.
+ * The integration is done using a quadrature formula on the face
+ * provided by the caller of the estimate() functions declared by this
+ * class. For linear trial functions (FE_Q(1)), QGauss with two points
+ * or even the QMidpoint rule might actually suffice. For higher order
+ * elements, it is necessary to utilize higher order quadrature
+ * formulae with `fe.degree+1` Gauss points.
*
* We store the contribution of each face in a @p map, as provided by the C++
* standard library, with the iterator pointing to that face being the key