-
Notifications
You must be signed in to change notification settings - Fork 11
/
burgers_test.jl
51 lines (46 loc) · 1.75 KB
/
burgers_test.jl
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
using Test, SummationByPartsOperators
for T in (Float32, Float64), split_form in (Val{true}(), Val{false}())
xmin = T(-1)
xmax = T(1)
N = 2^6
u0func = sinpi
tspan = (zero(T), one(T))
# Fourier
let D = fourier_derivative_operator(xmin, xmax, N)
Di = dissipation_operator(TadmorWaagan2012Standard(), D)
semidisc = BurgersPeriodicSemidiscretisation(D, Di, split_form)
println(devnull, semidisc)
ode = semidiscretise(u0func, semidisc, tspan)
du = similar(ode.u0)
semidisc(du, ode.u0, nothing, first(tspan))
end
# Periodic FD
for acc_order in 2:2:8
D = periodic_derivative_operator(1, acc_order, xmin, xmax, N)
Di = dissipation_operator(D)
semidisc = BurgersPeriodicSemidiscretisation(D, Di, split_form)
println(devnull, semidisc)
ode = semidiscretise(u0func, semidisc, tspan)
du = similar(ode.u0)
semidisc(du, ode.u0, nothing, first(tspan))
end
# Legendre
let D = legendre_derivative_operator(xmin, xmax, N)
Di = nothing
semidisc = BurgersNonperiodicSemidiscretisation(D, Di, split_form, zero, zero)
println(devnull, semidisc)
ode = semidiscretise(u0func, semidisc, tspan)
du = similar(ode.u0)
semidisc(du, ode.u0, nothing, first(tspan))
end
# SBP FD
for acc_order in 2:2:8
D = derivative_operator(MattssonSvärdNordström2004(), 1, acc_order, xmin, xmax, N)
Di = dissipation_operator(D)
semidisc = BurgersNonperiodicSemidiscretisation(D, Di, split_form, zero, zero)
println(devnull, semidisc)
ode = semidiscretise(u0func, semidisc, tspan)
du = similar(ode.u0)
semidisc(du, ode.u0, nothing, first(tspan))
end
end