Cauchy random number with specified location and scale parameters. The Cauchy distribution has heavier tails than the normal distribution, which can be useful for occasional large jumps in simulated annealing.
| Type | Intent | Optional | Attributes | Name | ||
|---|---|---|---|---|---|---|
| real(kind=wp), | intent(in) | :: | location |
location parameter (median) |
||
| real(kind=wp), | intent(in) | :: | scale |
scale parameter |
function cauchy(location, scale) real(wp),intent(in) :: location !! location parameter (median) real(wp),intent(in) :: scale !! scale parameter real(wp) :: cauchy real(wp) :: u u = uniform_random_number() cauchy = location + scale * tan(pi * (u - 0.5_wp)) end function cauchy