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
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
|
if test x$donotrun != xtrue; then
for solver in gerris hypre; do
if test x$solver = xhypre; then
gmodule=GModule
else
gmodule="# GModule"
fi
for level in 3 4 5 6 7 8 9; do
for cycle in 1 2 3 4 5 6 7 8 9 10; do
if ( sed "s/GModule/$gmodule/" < $1 | \
gerris2D -DLEVEL=$level -DCYCLE=$cycle -DSOLVER=$solver - ) ; then :
else
exit 1
fi
done
join time proj > res-$level-$solver
rm -f proj time
done
done
fi
rm -f error error-hypre error-gerris
for solver in gerris hypre; do
for level in 3 4 5 6 7 8; do
next=`expr $level + 1`
echo -n "$level " >> error-$solver
if gfscompare2D -C -c -v sim-$level-$solver sim-$next-$solver P 2>&1 | \
awk '{if ($1 == "total") print $4 " " $6 " " $8;}' >> error-$solver; then :
else
exit 1
fi
done
if echo "Save solution.eps { format = EPS line_width = 0.25}" | \
gfsview-batch2D sim-9-$solver solution.gfv; then :
else
exit 1
fi
if awk '
BEGIN { n = 0 }
{
l[n] = $1; n1[n] = $2; n2[n] = $3; ni[n++] = $4;
}
END {
for (i = 1; i < n; i++)
print l[i] " " log(n1[i-1]/n1[i])/log(2.) " " log(n2[i-1]/n2[i])/log(2.) " " log(ni[i-1]/ni[i])/log(2.);
}' < error-$solver > order-$solver; then :
else
exit 1
fi
done
if cat <<EOF | gnuplot ; then :
set term postscript eps color lw 3 solid 20
set output 'residual.eps'
set xlabel 'CPU time'
set ylabel 'Maximum residual'
set logscale y
plot 'res-7.ref' u 2:3 t 'ref' w lp, \
'res-7-gerris' u 2:3 t 'Gerris' w lp, \
'res-7-hypre' u 2:3 t 'Hypre' w lp
set output 'rate.eps'
set xlabel 'V-cycle'
set ylabel 'Cumulative residual reduction factor'
unset logscale
plot 'res-7.ref' u 1:4 t 'ref' w lp, \
'res-7-gerris' u 1:4 t 'Gerris' w lp, \
'res-7-hypre' u 1:4 t 'Hypre' w lp
set output 'error.eps'
set xlabel 'Level'
set ylabel 'Error norms'
set key
set logscale y
plot 'error.ref' u 1:2 t '1 (ref)' w lp, \
'error.ref' u 1:3 t '2 (ref)' w lp, \
'error.ref' u 1:4 t 'max (ref)' w lp, \
'error-gerris' u 1:2 t '1 Gerris' w lp, \
'error-gerris' u 1:3 t '2 Gerris' w lp, \
'error-gerris' u 1:4 t 'max Gerris' w lp, \
'error-hypre' u 1:2 t '1 Hypre' w lp, \
'error-hypre' u 1:3 t '2 Hypre' w lp, \
'error-hypre' u 1:4 t 'max Hypre' w lp
set output 'order.eps'
set xlabel 'Level'
set ylabel 'Order'
set key
unset logscale
set xtics 0,1
set ytics 0,1
set grid
plot [][0:3] 'order.ref' u 1:2 t '1 (ref)' w lp, \
'order.ref' u 1:3 t '2 (ref)' w lp, \
'order.ref' u 1:4 t 'max (ref)' w lp, \
'order-gerris' u 1:2 t '1 Gerris' w lp, \
'order-gerris' u 1:3 t '2 Gerris' w lp, \
'order-gerris' u 1:4 t 'max Gerris' w lp, \
'order-hypre' u 1:2 t '1 Hypre' w lp, \
'order-hypre' u 1:3 t '2 Hypre' w lp, \
'order-hypre' u 1:4 t 'max Hypre' w lp
EOF
else
exit 1
fi
if cat <<EOF | python ; then :
from check import *
from sys import *
if (Curve('res-7-gerris',1,3) - Curve('res-7.ref',1,3)).max() > 1e-8 or\
(Curve('res-7-hypre',1,3) - Curve('res-7.ref',1,7)).max() > 1e-8 or\
(Curve('error-gerris',1,4) - Curve('error.ref',1,4)).max() > 1e-6 or\
(Curve('error-hypre',1,4) - Curve('error.ref',1,8)).max() > 1e-6:
print (Curve('res-7-gerris',1,3) - Curve('res-7.ref',1,3)).max()
print (Curve('res-7-hypre',1,3) - Curve('res-7.ref',1,7)).max()
print (Curve('error-gerris',1,4) - Curve('error.ref',1,4)).max()
print (Curve('error-hypre',1,4) - Curve('error.ref',1,8)).max()
exit(1)
EOF
else
exit 1
fi
|
