How can I compute normally distributed 3D rotation matrices with Mathematica?
For 2D matrices I would sample a normal distributed angle and directly create a rotation matrix with:
normalRotation2[s_] := RotationMatrix[RandomVariate[NormalDistribution[0, s]]]
This has some kind of problems because the two tails of the distribution overlap at the far end due to the finite "size" of SO(2).
I am not quite sure what the correct definition of normal distribution would be here... If there isn't any, let's use some next best concept like in
Now the question is how to do this in 3D?
Edit: For my purposes it would be enough to have an isotropic normal distribution.