Merge pull request #18 from SciML/convergence

diffusion delay free convergence test

Author: ChrisRackauckas

test/analyticless_convergence_tests.jl

@@ -9,7 +9,7 @@ begin
     end
     function hayes_modelg(du,u,h,p,t)
         τ,a,b,c,α,β,γ = p
-        du .= α.*u .+ β.*h(p,t-τ) .+ γ
+        du .= α.*u .+ γ
     end
     h(p,t) = (ones(1) .+ t);
     tspan = (0.,10.)
@@ -29,28 +29,35 @@ sim2 = analyticless_test_convergence(dts,prob,EM(),test_dt,trajectories=300, use
 sim2 = analyticless_test_convergence(dts,prob,LambaEM(),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-0.5) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,EulerHeun(),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,LambaEulerHeun(),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,RKMil(),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
-
-#=
 sim2 = analyticless_test_convergence(dts,prob,RKMil(interpretation=:Stratonovich),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
-=#
+sim2 = analyticless_test_convergence(dts,prob,WangLi3SMil_A(),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
+sim2 = analyticless_test_convergence(dts,prob,WangLi3SMil_B(),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
+sim2 = analyticless_test_convergence(dts,prob,WangLi3SMil_C(),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
+sim2 = analyticless_test_convergence(dts,prob,WangLi3SMil_D(),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
+sim2 = analyticless_test_convergence(dts,prob,WangLi3SMil_E(),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
+sim2 = analyticless_test_convergence(dts,prob,WangLi3SMil_F(),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
 
 # Test SROCK methods
 println("SROCK methods")
 prob.p .= pmul;
-#=
 sim2 = analyticless_test_convergence(dts,prob,SROCK1(),test_dt,trajectories=100, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,SROCK1(interpretation=:Stratonovich),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
-=#
 sim2 = analyticless_test_convergence(dts,prob,SROCKEM(),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,SROCKEM(strong_order_1=false),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-0.5) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,SKSROCK(),test_dt,trajectories=300, use_noise_grid=false)
@@ -65,23 +72,20 @@ sim2 = analyticless_test_convergence(dts,prob,ImplicitEM(),test_dt,trajectories=
 sim2 = analyticless_test_convergence(dts,prob,ImplicitEM(symplectic=true, theta = 1/2),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-0.5) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,ImplicitEulerHeun(),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,ImplicitEulerHeun(symplectic=true, theta = 1/2),test_dt,trajectories=300, use_noise_grid=false)
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
+sim2 = analyticless_test_convergence(dts,prob,ISSEM(),test_dt,trajectories=500, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-0.5) < 0.3
-sim2 = analyticless_test_convergence(dts,prob,ISSEM(),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
-sim2 = analyticless_test_convergence(dts,prob,ISSEM(symplectic=true, theta = 1/2),test_dt,trajectories=300, use_noise_grid=false)
+sim2 = analyticless_test_convergence(dts,prob,ISSEM(symplectic=true, theta = 1/2),test_dt,trajectories=500, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-0.5) < 0.3
-
-#=
 sim2 = analyticless_test_convergence(dts,prob,ImplicitRKMil(),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,ImplicitRKMil(symplectic=true, theta = 1/2),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,ImplicitRKMil(interpretation=:Stratonovich, symplectic = true, theta = 1/2),test_dt,trajectories=300, use_noise_grid=false)
 @test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,ISSEulerHeun(),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3
 sim2 = analyticless_test_convergence(dts,prob,ISSEulerHeun(symplectic=true, theta = 1/2),test_dt,trajectories=300, use_noise_grid=false)
-@test abs(sim2.𝒪est[:final]-0.5) < 0.3
-=#
+@test abs(sim2.𝒪est[:final]-1.0) < 0.3