Cover indepsym indexing case

Author: sharanry

src/vector_of_array.jl

@@ -72,7 +72,11 @@ Base.@propagate_inbounds function Base.getindex(A::AbstractDiffEqArray{T, N},sym
   end
 
   if i === nothing
-    observed(A,sym,:)
+    if issymbollike(sym) && A.indepsym !== nothing && Symbol(sym) == A.indepsym
+      A.t
+    else
+      observed(A,sym,:)
+    end
   else
     Base.getindex.(A.u, i)
   end
@@ -85,7 +89,11 @@ Base.@propagate_inbounds function Base.getindex(A::AbstractDiffEqArray{T, N},sym
   end
 
   if i === nothing
-    observed(A,sym,args...)
+    if issymbollike(sym) && A.indepsym !== nothing && Symbol(sym) == A.indepsym
+      A.t[args...]
+    else
+      observed(A,sym,args...)
+    end
   else
     Base.getindex.(A.u, i, args...)
   end

test/downstream/symbol_indexing.jl

@@ -21,3 +21,5 @@ sol_new = DiffEqArray(
 )
 
 @test sol_new[RHS] ≈ (1 .- sol_new[x])./3.0
+@test sol_new[t] ≈ sol_new.t
+@test sol_new[t, 1:5] ≈ sol_new.t[1:5]