bump trajectories

Author: frankschae

test/weak_convergence/srk_weak_final.jl

@@ -52,7 +52,7 @@ sim = test_convergence(dts,ensemble_prob,DRI1(),
 println("DRI1:", sim.𝒪est[:weak_final])
 
 
-numtraj = Int(4e6)
+numtraj = Int(1e7)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -88,7 +88,7 @@ sim = test_convergence(dts,ensemble_prob,RI5(),
 @test abs(sim.𝒪est[:weak_final]-2) < 0.3
 println("RI5:", sim.𝒪est[:weak_final])
 
-numtraj = Int(1e6)
+numtraj = Int(2e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -162,7 +162,7 @@ ensemble_prob = EnsembleProblem(prob;
         )
 
 
-numtraj = Int(3e6)
+numtraj = Int(5e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -184,7 +184,7 @@ sim = test_convergence(dts,ensemble_prob,DRI1NM(),
 @test abs(sim.𝒪est[:weak_final]-2) < 0.3
 println("DRI1NM:", sim.𝒪est[:weak_final])
 
-numtraj = Int(1e6)
+numtraj = Int(2e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -196,7 +196,7 @@ sim = test_convergence(dts,ensemble_prob,RI1(),
 @test abs(sim.𝒪est[:weak_final]-2) < 0.3
 println("RI1:", sim.𝒪est[:weak_final])
 
-numtraj = Int(5e5)
+numtraj = Int(1e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -209,7 +209,7 @@ sim = test_convergence(dts,ensemble_prob,RI3(),
 println("RI3:", sim.𝒪est[:weak_final])
 
 
-numtraj = Int(1e6)
+numtraj = Int(2e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -221,7 +221,7 @@ sim = test_convergence(dts,ensemble_prob,RI5(),
 @test abs(sim.𝒪est[:weak_final]-2) < 0.3
 println("RI5:", sim.𝒪est[:weak_final])
 
-numtraj = Int(1e6)
+numtraj = Int(2e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -301,7 +301,7 @@ ensemble_prob = EnsembleProblem(prob;
         prob_func = prob_func
         )
 
-numtraj = Int(4e6)
+numtraj = Int(5e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -360,7 +360,7 @@ sim = test_convergence(dts,ensemble_prob,RI6(),
 println("R6:", sim.𝒪est[:weak_final])
 
 
-numtraj = Int(1e6)
+numtraj = Int(2e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 

test/weak_convergence/weak_strat.jl

@@ -23,7 +23,7 @@ end
 
 @info "Scalar oop noise"
 
-numtraj = Int(5e6) # in the paper they use 1e9
+numtraj = Int(7e6) # in the paper they use 1e9
 u₀ = 0.1
 p = [1.5, 0.1]
 f(u,p,t) = p[1]*u

test/weak_convergence/weak_strat_non_diagonal.jl

@@ -47,7 +47,7 @@ ensemble_prob = EnsembleProblem(prob;
         prob_func = prob_func
         )
 
-numtraj = Int(2e7)
+numtraj = Int(3e7)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)