f(t) = (1 + 2 rho cos(t - mu)) / (2 pi), 0 <= t < 2 pi, |rho| <= 1/2. Density, distribution function, quantile and draws by inversion of the morie Philox uniforms (identical in the Python arm).
f(t) = (1 + 2 rho cos(t - mu)) / (2 pi), 0 <= t < 2 pi, |rho| <= 1/2. Density, distribution function, quantile and draws by inversion of the morie Philox uniforms (identical in the Python arm).