@@ -34,8 +34,8 @@ for alg in algs
3434 adaptive = false , dt = 1.0f0 , callback = cb, merge_callbacks = true ,
3535 tstops = [2.40f0 ])
3636
37- @test norm (bench_sol (2.40f0 ) - sol[1 ](2.40f0 )) < 2e-3
38- @test norm (bench_sol. u - sol[1 ]. u) < 5e-3
37+ @test norm (bench_sol (2.40f0 ) - sol. u [1 ](2.40f0 )) < 2e-3
38+ @test norm (bench_sol. u - sol. u [1 ]. u) < 5e-3
3939
4040 # Test the truncation error due to floating point math, encountered when adjusting t for tstops
4141 local sol = solve (monteprob, alg, EnsembleGPUKernel (backend),
@@ -47,8 +47,8 @@ for alg in algs
4747 adaptive = false , dt = 0.01f0 , callback = cb, merge_callbacks = true ,
4848 tstops = [4.0f0 ])
4949
50- @test norm (bench_sol (4.0f0 ) - sol[1 ](4.0f0 )) < 2e-6
51- @test norm (bench_sol. u - sol[1 ]. u) < 3e-5
50+ @test norm (bench_sol (4.0f0 ) - sol. u [1 ](4.0f0 )) < 2e-6
51+ @test norm (bench_sol. u - sol. u [1 ]. u) < 3e-5
5252
5353 @info " Callback: CallbackSets"
5454
@@ -70,9 +70,9 @@ for alg in algs
7070 adaptive = false , dt = 1.0f0 , callback = cb, merge_callbacks = true ,
7171 tstops = [2.40f0 , 4.0f0 ])
7272
73- @test norm (bench_sol (2.40f0 ) - sol[1 ](2.40f0 )) < 2e-3
74- @test norm (bench_sol (4.0f0 ) - sol[1 ](4.0f0 )) < 3e-3
75- @test norm (bench_sol. u - sol[1 ]. u) < 7e-3
73+ @test norm (bench_sol (2.40f0 ) - sol. u [1 ](2.40f0 )) < 2e-3
74+ @test norm (bench_sol (4.0f0 ) - sol. u [1 ](4.0f0 )) < 3e-3
75+ @test norm (bench_sol. u - sol. u [1 ]. u) < 7e-3
7676
7777 @info " saveat and callbacks"
7878
@@ -85,9 +85,9 @@ for alg in algs
8585 adaptive = false , dt = 1.0f0 , callback = cb, merge_callbacks = true ,
8686 tstops = [2.40f0 , 4.0f0 ], saveat = [0.0f0 , 6.0f0 ])
8787
88- @test norm (bench_sol (2.40f0 ) - sol[1 ](2.40f0 )) < 1e-3
89- @test norm (bench_sol (6.0f0 ) - sol[1 ](6.0f0 )) < 3e-3
90- @test norm (bench_sol. u - sol[1 ]. u) < 3e-3
88+ @test norm (bench_sol (2.40f0 ) - sol. u [1 ](2.40f0 )) < 1e-3
89+ @test norm (bench_sol (6.0f0 ) - sol. u [1 ](6.0f0 )) < 3e-3
90+ @test norm (bench_sol. u - sol. u [1 ]. u) < 3e-3
9191
9292 @info " save_everystep and callbacks"
9393
@@ -100,9 +100,9 @@ for alg in algs
100100 adaptive = false , dt = 1.0f0 , callback = cb, merge_callbacks = true ,
101101 tstops = [2.40f0 , 4.0f0 ], save_everystep = false )
102102
103- @test norm (bench_sol (2.40f0 ) - sol[1 ](2.40f0 )) < 3e-5
104- @test norm (bench_sol (4.0f0 ) - sol[1 ](4.0f0 )) < 5e-5
105- @test norm (bench_sol. u - sol[1 ]. u) < 2e-4
103+ @test norm (bench_sol (2.40f0 ) - sol. u [1 ](2.40f0 )) < 3e-5
104+ @test norm (bench_sol (4.0f0 ) - sol. u [1 ](4.0f0 )) < 5e-5
105+ @test norm (bench_sol. u - sol. u [1 ]. u) < 2e-4
106106
107107 @info " Adaptive version"
108108
@@ -118,8 +118,8 @@ for alg in algs
118118 merge_callbacks = true ,
119119 tstops = [4.0f0 ])
120120
121- @test norm (bench_sol (4.0f0 ) - sol[1 ](4.0f0 )) < 5e-5
122- @test norm (bench_sol. u - sol[1 ]. u) < 2e-4
121+ @test norm (bench_sol (4.0f0 ) - sol. u [1 ](4.0f0 )) < 5e-5
122+ @test norm (bench_sol. u - sol. u [1 ]. u) < 2e-4
123123
124124 @info " Callback: CallbackSets"
125125
@@ -133,9 +133,9 @@ for alg in algs
133133 merge_callbacks = true ,
134134 tstops = [2.40f0 , 4.0f0 ])
135135
136- @test norm (bench_sol (2.40f0 ) - sol[1 ](2.40f0 )) < 6e-4
137- @test norm (bench_sol (4.0f0 ) - sol[1 ](4.0f0 )) < 1e-3
138- @test norm (bench_sol. u - sol[1 ]. u) < 3e-3
136+ @test norm (bench_sol (2.40f0 ) - sol. u [1 ](2.40f0 )) < 6e-4
137+ @test norm (bench_sol (4.0f0 ) - sol. u [1 ](4.0f0 )) < 1e-3
138+ @test norm (bench_sol. u - sol. u [1 ]. u) < 3e-3
139139
140140 @info " saveat and callbacks"
141141
@@ -151,9 +151,9 @@ for alg in algs
151151 tstops = [2.40f0 , 4.0f0 ], saveat = [0.0f0 , 6.0f0 ], reltol = 1.0f-7 ,
152152 abstol = 1.0f-7 )
153153
154- @test norm (bench_sol (2.40f0 ) - sol[1 ](2.40f0 )) < 7e-3
155- @test norm (bench_sol (6.0f0 ) - sol[1 ](6.0f0 )) < 2e-2
156- @test norm (bench_sol. u - sol[1 ]. u) < 2e-2
154+ @test norm (bench_sol (2.40f0 ) - sol. u [1 ](2.40f0 )) < 7e-3
155+ @test norm (bench_sol (6.0f0 ) - sol. u [1 ](6.0f0 )) < 2e-2
156+ @test norm (bench_sol. u - sol. u [1 ]. u) < 2e-2
157157
158158 @info " Unadaptive and Adaptive comparison"
159159
@@ -167,9 +167,9 @@ for alg in algs
167167 adaptive = true , dt = 1.0f0 , callback = cb, merge_callbacks = true ,
168168 tstops = [2.40f0 , 4.0f0 ], saveat = [0.0f0 , 4.0f0 ])
169169
170- @test norm (asol[1 ](2.40f0 ) - sol[1 ](2.40f0 )) < 3e-3
171- @test norm (asol[1 ](4.0f0 ) - sol[1 ](4.0f0 )) < 5e-3
172- @test norm (asol[1 ]. u - sol[1 ]. u) < 5e-3
170+ @test norm (asol. u [1 ](2.40f0 ) - sol. u [1 ](2.40f0 )) < 3e-3
171+ @test norm (asol. u [1 ](4.0f0 ) - sol. u [1 ](4.0f0 )) < 5e-3
172+ @test norm (asol. u [1 ]. u - sol. u [1 ]. u) < 5e-3
173173
174174 @info " Terminate callback"
175175
@@ -184,6 +184,6 @@ for alg in algs
184184 adaptive = false , dt = 1.0f0 , callback = cb, merge_callbacks = true ,
185185 tstops = [2.40f0 ])
186186
187- @test norm (bench_sol. t - sol[1 ]. t) < 2e-3
188- @test norm (bench_sol. u - sol[1 ]. u) < 5e-3
187+ @test norm (bench_sol. t - sol. u [1 ]. t) < 2e-3
188+ @test norm (bench_sol. u - sol. u [1 ]. u) < 5e-3
189189end
0 commit comments