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
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
|
set terminal cairolatex pdf color standalone
set output "harten-Pk-err-before-choc.tex"
a0 = 0.955006 # +/- 0.005645 (0.5911%)
b0 = -0.0492269 # +/- 0.01061 (21.56%)
log10_err_p0(log10_h) = a0*log10_h + b0
fit log10_err_p0(x) "harten-P0-err-before-choc.gdat" u (log10(2./$1)):(log10($3)) via a0,b0
err_p0(h) = 10.0**log10_err_p0(log10(h))
a1 = 1.73725 # +/- 0.06837 (3.935%)
b1 = -0.356635 # +/- 0.1398 (39.2%)
log10_err_p1(log10_h) = a1*log10_h + b1
fit log10_err_p1(x) "harten-P1d-err-before-choc.gdat" u (log10(2./$1)):(log10($3)) via a1,b1
err_p1(h) = 10.0**log10_err_p1(log10(h))
set size square 0.9
set colors classic
set log
set cbrange [1:2] # avoid error cb value is 0 in log scale when using fit
set cbrange [1:2]
set key bottom
set xrange [1e-3:1e-1]
set yrange [1e-5:1e-1]
graph_ratio = 2.0/4.0
set xtics (\
'[c]{$10^{-3}$}' 1e-3, \
'[c]{$10^{-2}$}' 1e-2, \
'[c]{$10^{-1}$}' 1e-1)
set ytics (\
'[r]{$10^{-5}$}' 1e-5, \
'[r]{$10^{-4}$}' 1e-4, \
'[r]{$10^{-3}$}' 1e-3, \
'[r]{$10^{-2}$}' 1e-2, \
'[r]{$10^{-1}$}' 1e-1)
set xlabel '[c]{$h$}'
set label '[l]{$\|u-u_h\|_{L^{\infty}(0,T;L^1)}$}' at graph 0.02, 0.95
set label '[l]{(a) $T=2/\pi$}' at graph 0.08, 0.85
# triangle a droite
slope_A = graph_ratio*a0
xA = 0.20
yA = 0.59
dxA = 0.10
dyA = dxA*slope_A
set label sprintf("[l]{\\scriptsize $%.2g$}",a0) at graph xA+dxA+0.02, yA+0.5*dyA right
set arrow from graph xA, yA to graph xA+dxA, yA nohead
set arrow from graph xA+dxA, yA to graph xA+dxA, yA+dyA nohead
set arrow from graph xA+dxA, yA+dyA to graph xA, yA nohead
# triangle a droite
slope_B = graph_ratio*a1
xB = 0.25
yB = 0.03
dxB = 0.10
dyB = dxB*slope_B
set label sprintf("[l]{\\scriptsize $%.2g$}",a1) at graph xB+dxB+0.02, yB+0.5*dyB right
set arrow from graph xB, yB to graph xB+dxB, yB nohead
set arrow from graph xB+dxB, yB to graph xB+dxB, yB+dyB nohead
set arrow from graph xB+dxB, yB+dyB to graph xB, yB nohead
plot \
"harten-P0-err-before-choc.gdat" \
using (2./$1):3 \
t '[r]{$k=0$}' \
w lp lt 1 lw 2 lc 0, \
"harten-P1d-err-before-choc.gdat" \
using (2./$1):3 \
t '[r]{$k=1$}' \
w lp lt 1 lw 2 lc 1 \
#pause -1
|
