dig in an actually do it

Author: ChrisRackauckas

test/odesystem.jl

@@ -467,9 +467,10 @@ prob = ODEProblem(sys, [], (0, 1.0))
 sol = solve(prob, Tsit5())
 @test sol[2x[1] + 3x[3] + norm(x)] ≈
       2sol[x[1]] + 3sol[x[3]] + sol[norm(x)]
-@test sol[x .+ [y, 2y, 3y]] ≈ [sol[x] .+ sol[y]
-                              sol[x] .+ 2sol[y]
-                              sol[x] .+ 3sol[y]]
+@test sol[x .+ [y, 2y, 3y]] ≈ map((x...)->[x...],
+                               map((x,y)->x[1].+y,sol[x],sol[y]),
+                               map((x,y)->x[2].+2y,sol[x],sol[y]),
+                               map((x,y)->x[3].+3y,sol[x],sol[y]))
 
 # Mixed Difference Differential equations
 @parameters t a b c d