hcubature not type stable

[sharanry]

Related: SciML/Integrals.jl#65

using HCubature, Distributions, Test

μ = [0.00 , 0.00]
Σ = [0.4 0.0 ; 0.00 0.4]
d = MvNormal(μ , Σ)
f = (x) -> pdf(d, x)
lb, ub = [-1 , -1] , [1 , 1]
@inferred hcubature(f, lb, ub;
            maxevals=typemax(Int))

ERROR: return type Tuple{Float64, Float64} does not match inferred return type Tuple{Any, Any}
Stacktrace:
 [1] error(s::String)
   @ Base ./error.jl:33
 [2] top-level scope
   @ REPL[7]:1

[sharanry]

Calling with StaticArrays is sufficient to make the function call type stable.

[stevengj]

Yes, or a tuple should work too.

Converting from a generic Array to an SVector cannot be type-stable because the length is not part of the type.