Default turn angle distribution about bactos.jl HOT 1 CLOSED

struct Arccos{T<:Real} <: ContinuousUnivariateDistribution
    a::T # lower limit
    b::T # upper limit
    Arccos{T}(a::T,b::T) where {T} = new{T}(a::T, b::T)
end

function Arccos(a::Real, b::Real; check_args::Bool=true)
    Distributions.@check_args Arccos a≤b -1≤a≤1 -1≤b≤1
    return Arccos{Float64}(Float64(a),Float64(b))
end # function

Arccos() = Arccos(-1,1)

function Base.rand(rng::AbstractRNG, d::Arccos)
    d₁ = Uniform(d.a, d.b)
    return acos(rand(rng, d₁))
end # function