Distributions.jl
Distributions.jl copied to clipboard
Is there a reason `rand(TDist(nu))` does not obey the type of `nu`?
Hi, it seems that
nu = 3f0
rand(TDist(nu))
returns a Float64
instead of a Float32
. Is this intended behavior? I find it unusual because the rand
implementation tries to be type-stable half way there by using one(T)
but ends up calling randn(rng)
without specifying the type.