fix: fix symbolic idxs interpolation for RODESolution

Author: AayushSabharwal

src/solutions/rode_solutions.jl

@@ -63,7 +63,7 @@ end
 
 function (sol::RODESolution)(t, ::Type{deriv} = Val{0}; idxs = nothing,
         continuity = :left) where {deriv}
-    sol.interp(t, idxs, deriv, sol.prob.p, continuity)
+    sol(t, deriv, idxs, continuity)
 end
 function (sol::RODESolution)(v, t, ::Type{deriv} = Val{0}; idxs = nothing,
         continuity = :left) where {deriv}

test/downstream/solution_interface.jl

@@ -62,6 +62,7 @@ sol = solve(sprob, ImplicitEM())
 @test_throws Exception sol[a]
 @test_throws Exception sol[noisy_population_model.a]
 @test_throws Exception sol[:a]
+@test_nowarn sol(0.5, idxs = noisy_population_model.s1)
 ### Tests on layered model (some things should not work). ###
 
 @parameters σ ρ β