@@ -63,23 +63,23 @@ using OrdinaryDiffEqRosenbrock
63
63
using SciMLBase: NoInit
64
64
65
65
function main (;
66
- n = 10 , Plotter = nothing , verbose = false , tend = 1 ,
67
- unknown_storage = :sparse , assembly = :edgewise ,
68
- diffeq = false ,
69
- switchbc = false
66
+ n = 10 , Plotter = nothing , verbose = false , tend = 1 ,
67
+ unknown_storage = :sparse , assembly = :edgewise ,
68
+ diffeq = false ,
69
+ switchbc = false
70
70
)
71
71
h = 1.0 / convert (Float64, n)
72
72
X = collect (0.0 : h: 1.0 )
73
73
N = length (X)
74
74
75
- iCat= 1
76
- iBulk= 2
77
- inodeCat= 1
75
+ iCat = 1
76
+ iBulk = 2
77
+ inodeCat = 1
78
78
if switchbc
79
79
iCat, iBulk = iBulk, iCat
80
- inodeCat= N
80
+ inodeCat = N
81
81
end
82
-
82
+
83
83
grid = simplexgrid (X)
84
84
# # By default, \Gamma_1 at X[1] and \Gamma_2 is at X[end]
85
85
@@ -203,7 +203,7 @@ function main(;
203
203
)
204
204
end
205
205
206
- return tsol[iC, inodeCat, end ]
206
+ return tsol[iC, inodeCat, end ]
207
207
end
208
208
209
209
using Test
@@ -214,7 +214,7 @@ function runtests()
214
214
for assembly in (:edgewise , :cellwise )
215
215
for switchbc in (false , true )
216
216
@test isapprox (main (; unknown_storage, assembly, switchbc), testval; rtol = 1.0e-12 )
217
- @test isapprox (main (; diffeq= true , unknown_storage, assembly, switchbc), testvaldiffeq; rtol = 1.0e-12 )
217
+ @test isapprox (main (; diffeq = true , unknown_storage, assembly, switchbc), testvaldiffeq; rtol = 1.0e-12 )
218
218
end
219
219
end
220
220
end
0 commit comments