bump trajectories

Author: frankschae

test/weak_convergence/srk_weak_final.jl

@@ -187,7 +187,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(8e6)
+numtraj = Int(9e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -353,7 +353,7 @@ sim = test_convergence(dts,ensemble_prob,RI5(),
 println("RI5:", sim.𝒪est[:weak_final])
 
 
-numtraj = Int(5e6)
+numtraj = Int(6e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -379,7 +379,7 @@ sim = test_convergence(dts,ensemble_prob,RDI1WM(),
 println("RDI1WM:", sim.𝒪est[:weak_final])
 
 
-numtraj = Int(5e6)
+numtraj = Int(6e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 

test/weak_convergence/weak_strat.jl

@@ -127,6 +127,11 @@ sim = test_convergence(dts,ensemble_prob,RS1(),
 @test abs(sim.𝒪est[:weak_final]-2) < 0.4
 println("RS1:", sim.𝒪est[:weak_final])
 
+numtraj = Int(7e6)
+seed = 100
+Random.seed!(seed)
+seeds = rand(UInt, numtraj)
+
 sim = test_convergence(dts,ensemble_prob,RS2(),
     save_everystep=false,trajectories=numtraj,save_start=false,adaptive=false,
     weak_timeseries_errors=false,weak_dense_errors=false,
@@ -135,6 +140,11 @@ sim = test_convergence(dts,ensemble_prob,RS2(),
 @test abs(sim.𝒪est[:weak_final]-2) < 0.4 #order is 2.36
 println("RS2:", sim.𝒪est[:weak_final])
 
+numtraj = Int(5e6)
+seed = 100
+Random.seed!(seed)
+seeds = rand(UInt, numtraj)
+
 sim = test_convergence(dts,ensemble_prob,NON(),
     save_everystep=false,trajectories=numtraj,save_start=false,adaptive=false,
     weak_timeseries_errors=false,weak_dense_errors=false,