Raise exception if quadrature rule does not fit the wavepackets
Unify quadrature initialization calls