test fails on Julia 1.8 about integrals.jl HOT 8 CLOSED

julia > dlb1, dub1, dp1 = Zygote.gradient(testf, lb, ub, p)
(-0.9092974268256817, -0.27941549819892586, -0.4091435516769757)

using Integrals, Zygote
f(x, p) = sum(sin.(p * x))
lb = 1.0
ub = 3.0
p = 2.0
function testf(lb, ub, p)
    prob = IntegralProblem(f, lb, ub, p)
    sin(solve(prob, QuadGKJL(), reltol = 1e-3, abstol = 1e-3)[1])
end
dlb1, dub1, dp1 = Zygote.gradient(testf, lb, ub, p)

julia> p = [2.0]
1-element Vector{Float64}:
 2.0

julia> function testf(lb, ub, p)
           prob = IntegralProblem(f, lb, ub, p)
           sin(solve(prob, QuadGKJL(), reltol = 1e-3, abstol = 1e-3)[1])
       end
testf (generic function with 1 method)

julia> dlb1, dub1, dp1 = Zygote.gradient(testf, lb, ub, p)
ERROR: BoundsError: attempt to access 0-element Vector{Any} at index []
Stacktrace:
  [1] throw_boundserror(A::Vector{Any}, I::Tuple{})
    @ Base ./abstractarray.jl:744
  [2] checkbounds
    @ ./abstractarray.jl:709 [inlined]
  [3] _getindex
    @ ./abstractarray.jl:1326 [inlined]
  [4] getindex
    @ ./abstractarray.jl:1294 [inlined]
  [5] macro expansion
    @ ~/.julia/packages/Zygote/SuKWp/src/lib/lib.jl:314 [inlined]
  [6] (::Zygote.Jnew{Base.MethodList, Vector{Any}, false})(Δ::Nothing)
    @ Zygote ~/.julia/packages/Zygote/SuKWp/src/lib/lib.jl:308

julia> p = [2.0]
1-element Vector{Float64}:
 2.0

julia> Zygote.gradient(f, 30, p)
(-1.9048259608303126, [-28.57238941245469])

julia> p = 2.0
2.0

julia> Zygote.gradient(f, 30, p)
(-1.9048259608303126, -28.57238941245469)