Implement option to rotate quadrature points

[ablech]

Applying random rotations to the quadrature points can be used to estimate the quadrature error and assess the robustness of the integration. This is particularly relevant for integrands with strong anisotropy.

For convenience, a corresponding parameter could be added to the S2 class.
Extension to the S3 and SO3 classes as well as other geometries should also be possible.