get the sizing correct

Author: ChrisRackauckas

test/odesystem.jl

@@ -466,8 +466,10 @@ sys = structural_simplify(sys)
 prob = ODEProblem(sys, [], (0, 1.0))
 sol = solve(prob, Tsit5())
 @test sol[2x[1] + 3x[3] + norm(x)] ≈
-      2sol[x[1]] + 3sol[x[3]] .+ vec(mapslices(norm, hcat(sol[x]...), dims = 2))
-@test sol[x + [y, 2y, 3y]] ≈ sol[x] .+ [sol[y], 2sol[y], 3sol[y]]
+      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]]
 
 # Mixed Difference Differential equations
 @parameters t a b c d