File: mosolov-circle-Bi-0%2C2-n-0%2C5-error-u-linf.plot

package info (click to toggle)
rheolef 7.1-6
  • links: PTS, VCS
  • area: main
  • in suites: bullseye
  • size: 77,392 kB
  • sloc: cpp: 105,337; sh: 16,014; makefile: 5,293; python: 1,359; xml: 221; yacc: 218; javascript: 202; awk: 61; sed: 5
file content (64 lines) | stat: -rw-r--r-- 1,661 bytes parent folder | download | duplicates (5)
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
set terminal cairolatex pdf color standalone
set output "mosolov-circle-Bi-0,2-n-0,5-error-u-linf.tex"

set size square
set log xy
set colors classic
set xrange [1e-2:1e-1]
set yrange [1e-6:1]
graph_ratio_xy = 1./6.

set xlabel '[c]{\Large $h$}'
set  label '[l]{\Large $\|u-u_h\|_{0,\infty,\Omega}$}' at graph 0.03, 0.93
set xtics add ( \
	'[l]{$10^{-2}$}'	1e-2, \
	'[r]{$10^{-1}$}'	1e-1)
set ytics add ( \
	'[r]{$1$}'		1, \
	''			1e-1, \
	'[r]{$10^{-2}$}'	1e-2, \
	''			1e-3, \
	'[r]{$10^{-4}$}'	1e-4, \
	''			1e-5, \
	'[r]{$10^{-6}$}'	1e-6)

# triangle a droite
slope_A = graph_ratio_xy*2.0
xA =  0.35
yA =  0.37
dxA = 0.15
dyA = dxA*slope_A
set label "[l]{\\scriptsize $2=\min(k+1,1/n)$}" 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_xy*2.0
xB =  0.35
yB =  0.01
dxB = 0.15
dyB = dxB*slope_B
set label "[l]{\\scriptsize $2$}" 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 \
'mosolov-circle-Bi-0,2-n-0,5-error.gdat' \
  i 0 \
  u (1/$1):2 \
  t '[r]{$k=1$}' \
  w lp lw 4 lc 1, \
'mosolov-circle-Bi-0,2-n-0,5-error.gdat' \
  i 1 \
  u (1/$1):2 \
  t '[r]{$k=2$}' \
  w lp lw 4 lc rgb "#008800", \
'mosolov-circle-Bi-0,2-n-0,5-error.gdat' \
  i 2 \
  u (1/$1):2 \
  t '[r]{$k=3$}' \
  w lp lw 4 lc 3

#pause -1 "<return>"