Closed sathvikbhagavan closed 11 months ago
cc @lxvm
Seems to be a dispatch problem, since passing the argument as a function works:
julia> quadgk(A, 0.0, 10.0)
(51.76911234248546, 4.274582013885939)
julia> quadgk(t -> A(t), 0.0, 10.0)
(50.0, 0.0)
Okay, the problem is that LinearInterpolation
is a subtype of AbstractVector
, so it is calling the wrong method of QuadGK.evalrule
.
Right, I can rename evalrule(::AbstractVector,...)
to batchevalrule
to fix the issue
Should be fixed by #90. Will merge and tag when green.
Thanks!
An issue was caught in https://github.com/SciML/DataInterpolations.jl/pull/178 where tests started failing with QuadGK@2.9
I have made a MWE to illustrate the issue (using a simple LinearInterpolation from DataInterpolations):
It is a straight line with slope 1. So the integral should be 50.0.
with QuadGK@2.8
with QuadGK@2.9