Random sample from LogNormal does not return values of the same type as its parameters

[itsdfish]

Hello,

I noticed that rand does not return the same type as the parameters in LogNormal. For example,

julia> rand(Normal(0f0, 1f0)) |> typeof
Float32

julia> rand(Gamma(1f0, 1f0)) |> typeof
Float32

julia> rand(LogNormal(0f0, 1f0)) |> typeof
Float64

The types for Gamma and Normal agree, but rand returns Float64 for LogNormal.

Version Info

[31c24e10] Distributions v0.25.118
Julia 1.11.5

[devmotion]

Closing as a duplicate as this issue has been discussed extensively in other issues (see a list of references in #1905). Would be fixed by #1905.