JuliaMath / QuadGK.jl

adaptive 1d numerical Gauss–Kronrod integration in Julia
MIT License
268 stars 37 forks source link

unitful infinite bounds #64

Closed stevengj closed 1 year ago

stevengj commented 1 year ago

Reported by @sosock in https://github.com/PainterQubits/Unitful.jl/issues/597#issuecomment-1385075022:

julia> using Unitful, QuadGK

julia> quadgk(x->exp(-x/(1.0u"m")), 0.0u"m", Inf*u"m")
ERROR: DomainError with Inf m:
integrand produced NaN m in the interval (0.0 m, Inf m)
Stacktrace:
[...]