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
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
|
if test x$donotrun != xtrue; then
if gerris2D -DSCHEME=explicit bump.gfs; then :
else
exit 1
fi
if sed -e 's/SourceDiffusionExplicit Rhoe/#/' \
-e 's/# SourceDiffusion/ SourceDiffusion/' < bump.gfs | \
gerris2D -DSCHEME=implicit -; then :
else
exit 1
fi
fi
if gnuplot <<EOF ; then :
set term postscript eps color enhanced lw 2 20 solid
#
# time evolution of charge density at origin
#
set output 'profile.eps'
set xlabel 't'
set ylabel '{/Symbol r}_{e}'
a=0.05
perm=2.
K=1.
rhoinic=1./a/sqrt(2*pi)
rho(x)=(rhoinic*exp(-x*K/perm))
plot \
'timevol-explicit' u 1:3 t 'explicit', \
'timevol-implicit' u 1:3 t 'implicit', \
rho(x) t '{/Symbol r}_{e}(0,t)'
#
# time evolution of bump times 0, 2, 4 and 6
#
set output 'figure.eps'
set xlabel 'r'
set ylabel '{/Symbol r}_{e}'
time0=0
time2=2
time4=4
time6=6
rhoinic=1./a/sqrt(2*pi)
rho0(x)=(rhoinic*exp(-0.5*x*x/a/a)*exp(-time0*K/perm))
rho2(x)=(rhoinic*exp(-0.5*x*x/a/a)*exp(-time2*K/perm))
rho4(x)=(rhoinic*exp(-0.5*x*x/a/a)*exp(-time4*K/perm))
rho6(x)=(rhoinic*exp(-0.5*x*x/a/a)*exp(-time6*K/perm))
plot [0:0.2]\
't_0-explicit' u 1:2 t '' lt 1, \
't_0-implicit' u 1:2 t '' lt 1, \
't_2-explicit' u 1:2 t '' lt 2, \
't_2-implicit' u 1:2 t '' lt 2, \
't_4-explicit' u 1:2 t '' lt 3, \
't_4-implicit' u 1:2 t '' lt 3, \
't_6-explicit' u 1:2 t '' lt 4, \
't_6-implicit' u 1:2 t '' lt 4, \
rho0(x) t '{/Symbol r}_{e}(x,0)' lt 1,\
rho2(x) t '{/Symbol r}_{e}(x,2)' lt 2,\
rho4(x) t '{/Symbol r}_{e}(x,4)' lt 3,\
rho6(x) t '{/Symbol r}_{e}(x,6)' lt 4
#
# Evolution of error norms
#
set output 'error.eps'
set xlabel 't'
set ylabel 'Error'
plot 'norms-explicit.ref' u 3:9 w l t 'Max (ref Explicit)', \
'norms-implicit.ref' u 3:9 w l t 'Max (ref Implicit)', \
'norms-implicit' u 3:9 w p t '', 'norms-explicit' u 3:9 w p t '', \
'norms-explicit.ref' u 3:7 w l t 'L2 (ref Explicit)', \
'norms-implicit.ref' u 3:7 w l t 'L2 (ref Implicit)', \
'norms-implicit' u 3:7 w p t '', 'norms-explicit' u 3:7 w p t ''
EOF
else
exit 1
fi
if cat <<EOF | python ; then :
from check import *
if (Curve('norms-explicit',3,7) - Curve('norms-explicit.ref',3,7)).max() > 1e-5 or \
(Curve('norms-explicit',3,9) - Curve('norms-explicit.ref',3,9)).max() > 1e-5 or \
(Curve('norms-implicit',3,7) - Curve('norms-implicit.ref',3,7)).max() > 1e-5 or \
(Curve('norms-implicit',3,9) - Curve('norms-implicit.ref',3,9)).max() > 1e-5:
print (Curve('norms-explicit',3,7) - Curve('norms-explicit.ref',3,7)).max()
print (Curve('norms-explicit',3,9) - Curve('norms-explicit.ref',3,9)).max()
print (Curve('norms-explicit',3,7) - Curve('norms-explicit.ref',3,7)).max()
print (Curve('norms-explicit',3,9) - Curve('norms-explicit.ref',3,9)).max()
print (Curve('norms-implicit',3,7) - Curve('norms-implicit.ref',3,7)).max()
print (Curve('norms-implicit',3,9) - Curve('norms-implicit.ref',3,9)).max()
exit(1)
EOF
else
exit 1
fi
|
