JeschkeLab/DeerLab

Incorrect normalization factor for 3D Rice distribution

stestoll opened this issue · 0 comments

The normalization of the Rice function in

def _multirice3dfun(r,nu,sig):

is incorrect.

The first factor in

 P = nu**(n/2-1)/s2*r**(n/2)*np.exp(-(r**2+nu**2)/(2*s2)+nu*r/s2)*I_scaled

should be

 P = nu**(1-n/2)/s2*r**(n/2)*np.exp(-(r**2+nu**2)/(2*s2)+nu*r/s2)*I_scaled

Overall, this error has no effect, since P is numerically normalized afterwards.