At the moment, you requiring the number of the 1D-quadrature points to be the same as the number of 1D-ansatz functions. This is not necessarily desired. Sometimes an under-integration has stabilizing effects or over-integration is used to decrease the integration error for non-polynomial right-hand sides.
In short: Make it an additional template argument.
At the moment, you requiring the number of the 1D-quadrature points to be the same as the number of 1D-ansatz functions. This is not necessarily desired. Sometimes an under-integration has stabilizing effects or over-integration is used to decrease the integration error for non-polynomial right-hand sides. In short: Make it an additional template argument.