sampling from watson/bingham distribution
The current implementation of watson_convolution uses uniform sampling and weights samples based on the watson odf, it might be a good idea to have a function that directly samples from watson distribution. Its better to sample from Bingham and use the special case of that (k1=k2) to sample from watson.
Edited by Hossein Rafipoor