depwarn fixes

Author: ChrisRackauckas

test/adaptive/sde_weak_adaptive.jl

@@ -23,7 +23,7 @@ function reduction(u,batch,I)
 end
 
 function output_func(sol,i)
-  #h1(asinh(sol[end][1])),false
+  #h1(asinh(sol.u[end][1])),false
   h1.(asinh.(sol)),false
 end
 

test/complex_tests.jl

@@ -15,7 +15,7 @@ implicit_noautodiff = [SKenCarp(autodiff=false), ImplicitEulerHeun(autodiff=fals
 
   for alg in Iterators.flatten((explicit, implicit_noautodiff))
     sol = solve(prob, alg)
-    @test eltype(sol[end]) == ComplexF64
+    @test eltype(sol.u[end]) == ComplexF64
   end
 
   # currently broken
@@ -34,7 +34,7 @@ end
 
   for alg in Iterators.flatten((explicit, implicit_noautodiff))
     sol = solve(prob, alg)
-    @test eltype(sol[end]) == ComplexF64
+    @test eltype(sol.u[end]) == ComplexF64
   end
 
   # currently broken

test/gpu/sde_weak_adaptive_gpu.jl

@@ -25,7 +25,7 @@ function reduction(u,batch,I)
 end
 
 function output_func(sol,i)
-  #h1(asinh(sol[end][1])),false
+  #h1(asinh(sol.u[end][1])),false
   h1.(asinh.(sol)),false
 end
 

test/rode_linear_tests.jl

@@ -8,19 +8,19 @@ prob = RODEProblem(f,u0,tspan)
 sol = solve(prob,RandomEM(),dt=1/100, save_noise=true)
 prob = RODEProblem(f,u0,tspan, noise=NoiseWrapper(sol.W))
 sol2 = solve(prob,RandomHeun(),dt=1/100)
-@test abs(sol[end] - sol2[end]) < 0.1 * abs(sol[end])
+@test abs(sol.u[end] - sol2[end]) < 0.1 * abs(sol.u[end])
 sol3 = solve(prob,RandomTamedEM(),dt=1/100)
-@test abs(sol[end] - sol3[end]) < 0.1 * abs(sol[end])
+@test abs(sol.u[end] - sol3[end]) < 0.1 * abs(sol.u[end])
 
 f(du,u,p,t,W) = (du.=1.01u.+0.87u.*W)
 u0 = ones(4)
 prob = RODEProblem(f,u0,tspan)
 sol = solve(prob,RandomEM(),dt=1/100, save_noise=true)
 prob = RODEProblem(f,u0,tspan, noise=NoiseWrapper(sol.W))
 sol2 = solve(prob,RandomHeun(),dt=1/100)
-@test sum(abs,sol[end]-sol2[end]) < 0.1 * sum(abs, sol[end])
+@test sum(abs,sol.u[end]-sol2[end]) < 0.1 * sum(abs, sol.u[end])
 sol3 = solve(prob,RandomEM(),dt=1/100)
-@test sum(abs,sol[end]-sol3[end]) < 0.1 * sum(abs, sol[end])
+@test sum(abs,sol.u[end]-sol3[end]) < 0.1 * sum(abs, sol.u[end])
 
 f(u,p,t,W) = 2u*sin(W)
 u0 = 1.00
@@ -29,9 +29,9 @@ prob = RODEProblem{false}(f,u0,tspan)
 sol = solve(prob,RandomEM(),dt=1/100, save_noise=true)
 prob = RODEProblem{false}(f,u0,tspan, noise=NoiseWrapper(sol.W))
 sol2 = solve(prob,RandomHeun(),dt=1/100)
-@test abs(sol[end]-sol2[end]) < 0.1 * abs(sol[end])
+@test abs(sol.u[end]-sol2[end]) < 0.1 * abs(sol.u[end])
 sol3 = solve(prob,RandomTamedEM(),dt=1/100)
-@test abs(sol[end]-sol3[end]) < 0.1 * abs(sol[end])
+@test abs(sol.u[end]-sol3[end]) < 0.1 * abs(sol.u[end])
 
 function f(du,u,p,t,W)
   du[1] = 0.2u[1]*sin(W[1] - W[2])
@@ -43,9 +43,9 @@ prob = RODEProblem(f,u0,tspan)
 sol = solve(prob,RandomEM(),dt=1/100, save_noise=true)
 prob = RODEProblem(f,u0,tspan, noise=NoiseWrapper(sol.W))
 sol2 = solve(prob,RandomHeun(),dt=1/100)
-@test sum(abs,sol[end]-sol2[end]) < 0.1 * sum(abs, sol[end])
+@test sum(abs,sol.u[end]-sol2[end]) < 0.1 * sum(abs, sol.u[end])
 sol3 = solve(prob,RandomTamedEM(),dt=1/100)
-@test sum(abs,sol[end]-sol3[end]) < 0.1 * sum(abs, sol[end])
+@test sum(abs,sol.u[end]-sol3[end]) < 0.1 * sum(abs, sol.u[end])
 
 function f(du,u,p,t,W)
   du[1] = -0.2W[3]*u[1]*sin(W[1] - W[2])
@@ -57,6 +57,6 @@ prob = RODEProblem(f,u0,tspan,rand_prototype=zeros(3))
 sol = solve(prob,RandomEM(),dt=1/100, save_noise=true)
 prob = RODEProblem(f,u0,tspan, noise=NoiseWrapper(sol.W))
 sol2 = solve(prob,RandomHeun(),dt=1/100)
-@test sum(abs,sol[end]-sol2[end]) < 0.1 * sum(abs, sol[end])
+@test sum(abs,sol.u[end]-sol2[end]) < 0.1 * sum(abs, sol.u[end])
 sol3 = solve(prob,RandomTamedEM(),dt=1/100)
-@test sum(abs,sol[end]-sol3[end]) < 0.1 * sum(abs, sol3[end])
+@test sum(abs,sol.u[end]-sol3[end]) < 0.1 * sum(abs, sol3[end])

test/split_tests.jl

@@ -25,7 +25,7 @@ prob = SDEProblem{false}(f,σ,u0,(0.0,1.0),noise = NoiseWrapper(sol.W))
 
 sol2 = solve(prob,EM(),dt=1/10)
 
-@test sol[end][:] ≈ sol2[end][:]
+@test sol.u[end][:] ≈ sol2[end][:]
 
 ################################################################################
 

test/tdir_tests.jl

@@ -26,6 +26,6 @@ prob = SDEProblem(f,g,u₀,(0.0,1.0))
 sol =solve(prob,EulerHeun(),dt=0.01,save_noise=true)
 _sol = deepcopy(sol) # to make sure the plot is correct
 W3 = NoiseGrid(reverse!(_sol.t),reverse!(_sol.W))
-prob3 = SDEProblem(f,g,sol[end],(1.0,0.0),noise=W3)
+prob3 = SDEProblem(f,g,sol.u[end],(1.0,0.0),noise=W3)
 sol2 = solve(prob3,EulerHeun(),dt=0.01)
 @test sol.u ≈ reverse!(sol2.u) atol=1e-1

test/weak_convergence/PL1WM.jl

@@ -38,7 +38,7 @@ seeds = rand(UInt, numtraj)
 
 prob = SDEProblem(f,g,u₀,tspan,p)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h1(sol[end]),false),
+        output_func = (sol,i) -> (h1(sol.u[end]),false),
         prob_func = prob_func
         )
 
@@ -67,7 +67,7 @@ p = [3//2,1//100]
 
 prob = SDEProblem(f1!,g1!,u₀,tspan,p)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h1(sol[end][1]),false),
+        output_func = (sol,i) -> (h1(sol.u[end][1]),false),
         prob_func = prob_func
         )
 
@@ -124,7 +124,7 @@ h2(z) = z
 
 prob = SDEProblem(f2!,g2!,u₀,tspan,p,noise_rate_prototype=zeros(4,4))
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h2(sol[end][1]),false),
+        output_func = (sol,i) -> (h2(sol.u[end][1]),false),
         prob_func = prob_func
         )
 
@@ -164,7 +164,7 @@ h3(z) = z^2 # == 1//10**exp(3//2*t) if h3(z) = z and  == 1//100**exp(301//100*t)
 
 prob = SDEProblem(f3!,g3!,u₀,tspan)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h3(sol[end][1]),false),
+        output_func = (sol,i) -> (h3(sol.u[end][1]),false),
         prob_func = prob_func
         )
 
@@ -200,7 +200,7 @@ seeds = rand(UInt, numtraj)
 
 prob = SDEProblem(fadd,gadd,u₀,tspan,p)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (sol[end],false),
+        output_func = (sol,i) -> (sol.u[end],false),
         prob_func = prob_func
         )
 
@@ -233,7 +233,7 @@ gadd!(du,u,p,t) = @.(du = p[2])
 
 prob = SDEProblem(fadd!,gadd!,u₀,tspan,p)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (sol[end][1],false),
+        output_func = (sol,i) -> (sol.u[end][1],false),
         prob_func = prob_func
         )
 

test/weak_convergence/SIE_SME.jl

@@ -38,7 +38,7 @@ seeds = rand(UInt, numtraj)
 
 prob = SDEProblem(f,g,u₀,tspan,p)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h1(sol[end]),false),
+        output_func = (sol,i) -> (h1(sol.u[end]),false),
         prob_func = prob_func
         )
 
@@ -90,7 +90,7 @@ p = [3//2,1//100]
 
 prob = SDEProblem(f1!,g1!,u₀,tspan,p)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h1(sol[end][1]),false),
+        output_func = (sol,i) -> (h1(sol.u[end][1]),false),
         prob_func = prob_func
         )
 
@@ -155,7 +155,7 @@ h3(z) = z^2 # == 1//10**exp(3//2*t) if h3(z) = z and  == 1//100**exp(301//100*t)
 
 prob = SDEProblem(f3!,g3!,u₀,tspan)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h3(sol[end][1]),false),
+        output_func = (sol,i) -> (h3(sol.u[end][1]),false),
         prob_func = prob_func
         )
 

test/weak_convergence/W2Ito1.jl

@@ -38,7 +38,7 @@ seeds = rand(UInt, numtraj)
 
 prob = SDEProblem(f, g, u₀, tspan)
 ensemble_prob = EnsembleProblem(prob;
-    output_func=(sol, i) -> (h1(asinh(sol[end])), false),
+    output_func=(sol, i) -> (h1(asinh(sol.u[end])), false),
     prob_func=prob_func
 )
 
@@ -67,7 +67,7 @@ h2(z) = z^2 # == 1//10**exp(3//2*t) if h3(z) = z and  == 1//100**exp(301//100*t)
 
 prob = SDEProblem(f2, g2, u₀, tspan)
 ensemble_prob = EnsembleProblem(prob;
-    output_func=(sol, i) -> (h2(sol[end][1]), false),
+    output_func=(sol, i) -> (h2(sol.u[end][1]), false),
     prob_func=prob_func
 )
 
@@ -101,7 +101,7 @@ h3(z) = z^2 # but apply it only to u[1]
 
 prob = SDEProblem(f3, g3, u₀, tspan, noise_rate_prototype=zeros(2, 2))
 ensemble_prob = EnsembleProblem(prob;
-    output_func=(sol, i) -> (h3(sol[end][1]), false),
+    output_func=(sol, i) -> (h3(sol.u[end][1]), false),
     prob_func=prob_func
 )
 
@@ -136,7 +136,7 @@ h1(z) = z^3 - 6 * z^2 + 8 * z
 
 prob = SDEProblem(f1!, g1!, u₀, tspan)
 ensemble_prob = EnsembleProblem(prob;
-    output_func=(sol, i) -> (h1(asinh(sol[end][1])), false),
+    output_func=(sol, i) -> (h1(asinh(sol.u[end][1])), false),
     prob_func=prob_func
 )
 
@@ -171,7 +171,7 @@ h2(z) = z^2 # == 1//10**exp(3//2*t) if h3(z) = z and  == 1//100**exp(301//100*t)
 
 prob = SDEProblem(f2!, g2!, u₀, tspan)
 ensemble_prob = EnsembleProblem(prob;
-    output_func=(sol, i) -> (h2(sol[end][1]), false),
+    output_func=(sol, i) -> (h2(sol.u[end][1]), false),
     prob_func=prob_func
 )
 

test/weak_convergence/srk_weak_diagonal_final.jl

@@ -37,7 +37,7 @@ h3(z) = z^2 # == 1//10**exp(3//2*t) if h3(z) = z and  == 1//100**exp(301//100*t)
 
 prob = SDEProblem(f3!,g3!,u₀,tspan)
 ensemble_prob = EnsembleProblem(prob;
-        output_func = (sol,i) -> (h3(sol[end][1]),false),
+        output_func = (sol,i) -> (h3(sol.u[end][1]),false),
         prob_func = prob_func
         )