File: dirichlet_dg-cv-l2-2d-t-grid.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 (66 lines) | stat: -rw-r--r-- 1,824 bytes parent folder | download | duplicates (2)
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
set terminal cairolatex pdf color standalone
set output "dirichlet_dg-cv-l2-2d-t-grid.tex"

set size square 0.9
set colors classic
set log
set xrange [5e-3:0.5]
set yrange [1e-15:1]

set xtics (\
        "[r]{$10^{-3}$}" 1e-3, \
        "[r]{$10^{-2}$}" 1e-2, \
        "[r]{$10^{-1}$}" 1e-1, \
        "[r]{$1$}" 1 )
set ytics (\
        "[r]{$10^{-15}$}" 1e-15, \
        "[r]{$10^{-10}$}" 1e-10, \
        "[r]{$10^{-5}$}" 1e-5, \
        "[r]{$1$}" 1 )
set xlabel '[c]{\large $h$}'
set  label '[l]{$\left\|u - u_h\right\|_{0,\Omega}$}' at graph 0.05,0.92
set key bottom

graph_ratio = 2.0/15.0

# triangle a droite, pente +2
slope_A = graph_ratio*2.0
xA =  0.15
yA =  0.70
dxA = 0.10
dyA = dxA*slope_A
set label "[l]{\\scriptsize $2=k+1$}" 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

# placian-fixed-point-p=1,5-err-linf.plot:#triangle a droite, pente +5
slope_B = graph_ratio*5.0
xB =  0.44
yB =  0.34
dxB = 0.10
dyB = dxB*slope_B
set label "[l]{\\scriptsize $5=k+1$}" 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 \
 "dirichlet_dg-cv-2d-t-grid.gdat" \
    i 0 u (sqrt(2.)/$1):2 \
    t '[r]{$k=1$}' \
    w lp, \
 "dirichlet_dg-cv-2d-t-grid.gdat" \
    i 1 u (sqrt(2.)/$1):2 \
    t '[r]{$k=2$}' \
    w lp, \
 "dirichlet_dg-cv-2d-t-grid.gdat" \
    i 2 u (sqrt(2.)/$1):2\
    t '[r]{$k=3$}' \
    w lp, \
 "dirichlet_dg-cv-2d-t-grid.gdat" \
    i 3 u (sqrt(2.)/$1):2\
    t '[r]{$k=4$}' \
    w lp

#pause -1 "<retour>"