← Back to branch summary

~fluidity-core/fluidity/remove_tictoc

« back to all changes in this revision

Viewing changes to tests/mms_tracer_cv_3d_adv_excludemass/mms_tracer_cv_3d_adv_excludemass.xml

  • Committer: Jon Hill
  • Date: 2012-02-28 21:36:27 UTC
  • mfrom: (3880.1.59 fluidity)
  • Revision ID: jon.hill@imperial.ac.uk-20120228213627-75tqkd0za9nlcjrt
Merge from trunk

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/mms_tracer_cv_3d_adv_excludemass/mms_tracer_cv_3d_adv_excludemass.xml
 
1
<?xml version="1.0" encoding="UTF-8" ?>
 
2
<testproblem>
 
3
  <name>p1cv advection equation spatial convergence test</name>
 
4
  <owner userid="cwilson"/>
 
5
  <tags>flml</tags>
 
6
  <problem_definition length="medium" nprocs="1">
 
7
    <command_line>fluidity  MMS_A_cv.flml &gt; MMS_A.log; fluidity  MMS_B_cv.flml &gt; MMS_B.log; fluidity  MMS_C_cv.flml &gt; MMS_C.log; fluidity  MMS_D_cv.flml &gt; MMS_D.log</command_line>
 
8
  </problem_definition>
 
9
  <variables>
 
10
    <variable name="ab_convergence_gal_stat" language="python">
 
11
from fluidity_tools import stat_parser as stat
 
12
from math import log
 
13
 
 
14
a_error = stat("MMS_A_cv.stat")["Burgers"]["AbsoluteDifference"]["cv_l2norm"][-1]
 
15
b_error = stat("MMS_B_cv.stat")["Burgers"]["AbsoluteDifference"]["cv_l2norm"][-1]
 
16
 
 
17
a_error_inf = stat("MMS_A_cv.stat")["Burgers"]["AbsoluteDifference"]["max"][-1]
 
18
b_error_inf = stat("MMS_B_cv.stat")["Burgers"]["AbsoluteDifference"]["max"][-1]
 
19
 
 
20
print a_error
 
21
print b_error
 
22
 
 
23
print a_error_inf
 
24
print b_error_inf
 
25
 
 
26
ab_ratio = a_error / b_error
 
27
ab_ratio_inf = a_error_inf / b_error_inf
 
28
ab_convergence_gal_stat = [log(ab_ratio, 2), log(ab_ratio_inf, 2)]
 
29
    </variable>
 
30
    <variable name="bc_convergence_gal_stat" language="python">
 
31
from fluidity_tools import stat_parser as stat
 
32
from math import log
 
33
 
 
34
b_error = stat("MMS_B_cv.stat")["Burgers"]["AbsoluteDifference"]["cv_l2norm"][-1]
 
35
c_error = stat("MMS_C_cv.stat")["Burgers"]["AbsoluteDifference"]["cv_l2norm"][-1]
 
36
 
 
37
b_error_inf = stat("MMS_B_cv.stat")["Burgers"]["AbsoluteDifference"]["max"][-1]
 
38
c_error_inf = stat("MMS_C_cv.stat")["Burgers"]["AbsoluteDifference"]["max"][-1]
 
39
 
 
40
print b_error
 
41
print c_error
 
42
 
 
43
print b_error_inf
 
44
print c_error_inf
 
45
 
 
46
bc_ratio = b_error / c_error
 
47
bc_ratio_inf = b_error_inf / c_error_inf
 
48
bc_convergence_gal_stat = [log(bc_ratio, 2), log(bc_ratio_inf, 2)]
 
49
    </variable>
 
50
    <variable name="cd_convergence_gal_stat" language="python">
 
51
from fluidity_tools import stat_parser as stat
 
52
from math import log
 
53
 
 
54
c_error = stat("MMS_C_cv.stat")["Burgers"]["AbsoluteDifference"]["cv_l2norm"][-1]
 
55
d_error = stat("MMS_D_cv.stat")["Burgers"]["AbsoluteDifference"]["cv_l2norm"][-1]
 
56
 
 
57
c_error_inf = stat("MMS_C_cv.stat")["Burgers"]["AbsoluteDifference"]["max"][-1]
 
58
d_error_inf = stat("MMS_D_cv.stat")["Burgers"]["AbsoluteDifference"]["max"][-1]
 
59
 
 
60
print c_error
 
61
print d_error
 
62
 
 
63
print c_error_inf
 
64
print d_error_inf
 
65
 
 
66
cd_ratio = c_error / d_error
 
67
cd_ratio_inf = c_error_inf / d_error_inf
 
68
cd_convergence_gal_stat = [log(cd_ratio, 2), log(cd_ratio_inf, 2)]
 
69
    </variable>
 
70
    <variable name="a_finish_time" language="python">
 
71
from fluidity_tools import stat_parser as stat
 
72
a_finish_time = stat("MMS_A_cv.stat")["ElapsedTime"]["value"][-1]
 
73
    </variable>
 
74
    <variable name="b_finish_time" language="python">
 
75
from fluidity_tools import stat_parser as stat
 
76
b_finish_time = stat("MMS_B_cv.stat")["ElapsedTime"]["value"][-1]
 
77
    </variable>
 
78
    <variable name="c_finish_time" language="python">
 
79
from fluidity_tools import stat_parser as stat
 
80
c_finish_time = stat("MMS_C_cv.stat")["ElapsedTime"]["value"][-1]
 
81
    </variable>
 
82
    <variable name="d_finish_time" language="python">
 
83
from fluidity_tools import stat_parser as stat
 
84
d_finish_time = stat("MMS_D_cv.stat")["ElapsedTime"]["value"][-1]
 
85
    </variable>
 
86
    <variable name="a_final_change" language="python">
 
87
import vtktools
 
88
from fluidity_tools import stat_parser as stat
 
89
vtu = vtktools.vtu("MMS_A_cv_1.vtu")
 
90
dt = stat("MMS_A_cv.stat")["dt"]["value"][-1]
 
91
a_final_change = max(abs(vtu.GetScalarField("NumericalSolution")-vtu.GetScalarField("OldNumericalSolution"))/dt)
 
92
    </variable>
 
93
    <variable name="b_final_change" language="python">
 
94
import vtktools
 
95
from fluidity_tools import stat_parser as stat
 
96
vtu = vtktools.vtu("MMS_B_cv_1.vtu")
 
97
dt = stat("MMS_B_cv.stat")["dt"]["value"][-1]
 
98
b_final_change = max(abs(vtu.GetScalarField("NumericalSolution")-vtu.GetScalarField("OldNumericalSolution"))/dt)
 
99
    </variable>
 
100
    <variable name="c_final_change" language="python">
 
101
import vtktools
 
102
from fluidity_tools import stat_parser as stat
 
103
vtu = vtktools.vtu("MMS_C_cv_1.vtu")
 
104
dt = stat("MMS_C_cv.stat")["dt"]["value"][-1]
 
105
c_final_change = max(abs(vtu.GetScalarField("NumericalSolution")-vtu.GetScalarField("OldNumericalSolution"))/dt)
 
106
    </variable>
 
107
    <variable name="d_final_change" language="python">
 
108
import vtktools
 
109
from fluidity_tools import stat_parser as stat
 
110
vtu = vtktools.vtu("MMS_D_cv_1.vtu")
 
111
dt = stat("MMS_D_cv.stat")["dt"]["value"][-1]
 
112
d_final_change = max(abs(vtu.GetScalarField("NumericalSolution")-vtu.GetScalarField("OldNumericalSolution"))/dt)
 
113
    </variable>
 
114
  </variables>
 
115
  <pass_tests>
 
116
    <test name="ab_convergence_stat: L2 order in [0.8,1.2]" language="python">
 
117
assert(abs(ab_convergence_gal_stat[0]-1.0) &lt; 0.2)
 
118
    </test>
 
119
    <test name="bc_convergence_stat: L2 order in [0.8,1.2]" language="python">
 
120
assert(abs(bc_convergence_gal_stat[0]-1.0) &lt; 0.2)
 
121
    </test>
 
122
    <test name="cd_convergence_stat: L2 order in [0.8,1.2]" language="python">
 
123
assert(abs(cd_convergence_gal_stat[0]-1.0) &lt; 0.2)
 
124
    </test>
 
125
    <test name="checking A finished in less than 3.0" language="python">
 
126
assert(a_finish_time-3.0 &lt; 1.E-10)
 
127
    </test>
 
128
    <test name="checking B finished in less than 2.0" language="python">
 
129
assert(b_finish_time-2.0 &lt; 1.E-10)
 
130
    </test>
 
131
    <test name="checking C finished in less than 1.5" language="python">
 
132
assert(c_finish_time-1.5 &lt; 1.E-10)
 
133
    </test>
 
134
    <test name="checking D finished in less than 1.5" language="python">
 
135
assert(d_finish_time-1.5 &lt; 1.E-10)
 
136
    </test>
 
137
    <test name="checking steady state was reached in A" language="python">
 
138
assert(a_final_change &lt; 1.E-10)
 
139
    </test>
 
140
    <test name="checking steady state was reached in B" language="python">
 
141
assert(b_final_change &lt; 1.E-10)
 
142
    </test>
 
143
    <test name="checking steady state was reached in C" language="python">
 
144
assert(c_final_change &lt; 1.E-10)
 
145
    </test>
 
146
    <test name="checking steady state was reached in D" language="python">
 
147
assert(d_final_change &lt; 1.E-10)
 
148
    </test>
 
149
  </pass_tests>
 
150
  <warn_tests>
 
151
  </warn_tests>
 
152
</testproblem>
 
153