A somewhat surprising request on X validated about the inverse cdf representation of a wrapped Cauchy distribution. I had not come across this distribution, but its density being
means that it is the superposition of shifted Cauchys on the unit circle (with nice complex representations). As such, it is easily simulated by re-shifting a Cauchy back to (-π,π), i.e. using the inverse transform
![\theta = [\gamma\tan(\pi U-\pi/2)+\pi]\ \text{mod}\,(2\pi) - \pi](https://s0.wp.com/latex.php?latex=%5Ctheta+%3D+%5B%5Cgamma%5Ctan%28%5Cpi+U-%5Cpi%2F2%29%2B%5Cpi%5D%5C+%5Ctext%7Bmod%7D%5C%2C%282%5Cpi%29+-+%5Cpi&bg=000000&fg=B0B0B0&s=0&c=20201002)
