Benchmarks about fhist.jl HOT 5 CLOSED

function Hist1D(A::AbstractVector, r::AbstractRange{T}) where T <: AbstractFloat
    s = step(r)
    start = first(r)
    start2 = start + 0.5s
    stop = last(r)
    L = length(r) - 1
    counts = zeros(Int, L)
    @inbounds for idx in eachindex(A)
        # default
        i = A[idx]
        c = ifelse(i > stop, 0, 1)
        c = ifelse(i < start, 0, c)
        id = round(Int, (i - start2) / s) + 1
        counts[clamp(id, 1, L)] += c

        # push!() uses `searchsorted`
        # i = A[idx]
        # id = searchsortedlast(r, i)
        # c = ifelse(i > stop, 0, 1)
        # c = ifelse(i < start, 0, c)
        # counts[id] += c

    end
    return Hist1D(Histogram(r, counts))
end

@inline function calcbinindex(r::AbstractRange{<:Real}, x::Real)
    s = step(r)
    start2 = first(r) + 0.5s
    return round(Int, (x - start2) / s) + 1
end