Update default_linsolve_structure.jl

Author: ChrisRackauckas

test/downstream/default_linsolve_structure.jl

@@ -6,23 +6,23 @@ jp_diag = Diagonal(zeros(2))
 fun = ODEFunction(f; jac = jac, jac_prototype = jp_diag)
 prob = ODEProblem(fun, ones(2), (1.0, 10.0))
 sol = solve(prob, Rosenbrock23())
-@test sol[end] ≈ [10.0, 10.0]
+@test sol.u[end] ≈ [10.0, 10.0]
 @test length(sol) < 60
 
 sol = solve(prob, Rosenbrock23(autodiff = false))
-@test sol[end] ≈ [10.0, 10.0]
+@test sol.u[end] ≈ [10.0, 10.0]
 @test length(sol) < 60
 
 jp = Tridiagonal(jp_diag)
 fun = ODEFunction(f; jac = jac, jac_prototype = jp)
 prob = ODEProblem(fun, ones(2), (1.0, 10.0))
 
 sol = solve(prob, Rosenbrock23())
-@test sol[end] ≈ [10.0, 10.0]
+@test sol.u[end] ≈ [10.0, 10.0]
 @test length(sol) < 60
 
 sol = solve(prob, Rosenbrock23(autodiff = false))
-@test sol[end] ≈ [10.0, 10.0]
+@test sol.u[end] ≈ [10.0, 10.0]
 @test length(sol) < 60
 
 #=
@@ -43,7 +43,7 @@ sol = solve(prob,Rosenbrock23())
     fun = ODEFunction(f; jac = jac, jac_prototype = jp)
     prob = ODEProblem(fun, ones(2), (1.0, 10.0))
     sol = solve(prob, Rosenbrock23(autodiff = false))
-    @test sol[end] ≈ [10.0, 10.0]
+    @test sol.u[end] ≈ [10.0, 10.0]
     @test length(sol) < 60
 end
 
@@ -52,6 +52,6 @@ end
     fun = ODEFunction(f; jac = jac, jac_prototype = jp)
     prob = ODEProblem(fun, ones(2), (1.0, 10.0))
     sol = solve(prob, Rosenbrock23(autodiff = false))
-    @test sol[end] ≈ [10.0, 10.0]
+    @test sol.u[end] ≈ [10.0, 10.0]
     @test length(sol) < 60
 end