bump trajectories

Author: frankschae

test/weak_convergence/srk_weak_final.jl

@@ -162,7 +162,7 @@ ensemble_prob = EnsembleProblem(prob;
         )
 
 
-numtraj = Int(5e6)
+numtraj = Int(6e6)
 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(2e6)
+numtraj = Int(3e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 
@@ -313,6 +313,9 @@ sim = test_convergence(dts,ensemble_prob,DRI1(),
 @test -(sim.𝒪est[:weak_final]-2) < 0.3 # order is 2.91
 println("DRI1:", sim.𝒪est[:weak_final])
 
+numtraj = Int(6e6)
+Random.seed!(seed)
+seeds = rand(UInt, numtraj)
 
 sim = test_convergence(dts,ensemble_prob,DRI1NM(),
     save_everystep=false,trajectories=numtraj,save_start=false,adaptive=false,
@@ -360,7 +363,7 @@ sim = test_convergence(dts,ensemble_prob,RI6(),
 println("R6:", sim.𝒪est[:weak_final])
 
 
-numtraj = Int(2e6)
+numtraj = Int(3e6)
 Random.seed!(seed)
 seeds = rand(UInt, numtraj)
 

test/weak_convergence/weak_strat.jl

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