1
<?xml version="1.0" encoding="UTF-8" ?>
2
<!DOCTYPE testproblem SYSTEM "regressiontest.dtd">
5
<name>cv_mass_matrix_check_p1</name>
6
<owner userid="btollit"/>
8
<problem_definition length="short" nprocs="1">
10
../../bin/fluidity cv_mass_matrix_check_p1_1d.flml
11
../../bin/fluidity cv_mass_matrix_check_p1_2d.flml
12
../../bin/fluidity cv_mass_matrix_check_p1_3d.flml
14
<!-- One/two/three dimensional problem for using p1 elements to soley check that the CV mass matrix is correct for one element.-->
17
<variable name="element_volume_1d" language="python">
18
from fluidity_tools import stat_parser as stat
19
element_volume_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["ElementVolume"]["max"][-1]
21
<variable name="max_cv_mass_1d" language="python">
22
from fluidity_tools import stat_parser as stat
23
max_cv_mass_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["CVMassMatrix"]["max"][-1]
25
<variable name="min_cv_mass_1d" language="python">
26
from fluidity_tools import stat_parser as stat
27
min_cv_mass_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["CVMassMatrix"]["min"][-1]
29
<variable name="sum_cv_mass_1d" language="python">
30
from fluidity_tools import stat_parser as stat
31
sum_cv_mass_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["SumCVMassMatrix"]["max"][-1]
33
<variable name="max_fe_lumped_mass_1d" language="python">
34
from fluidity_tools import stat_parser as stat
35
max_fe_lumped_mass_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["FELumpedMassMatrix"]["max"][-1]
37
<variable name="min_fe_lumped_mass_1d" language="python">
38
from fluidity_tools import stat_parser as stat
39
min_fe_lumped_mass_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["FELumpedMassMatrix"]["min"][-1]
41
<variable name="sum_fe_lumped_mass_1d" language="python">
42
from fluidity_tools import stat_parser as stat
43
sum_fe_lumped_mass_1d = stat("cv_mass_matrix_check_p1_1d.stat")["fluid"]["SumFELumpedMassMatrix"]["max"][-1]
46
<variable name="element_volume_2d" language="python">
47
from fluidity_tools import stat_parser as stat
48
element_volume_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["ElementVolume"]["max"][-1]
50
<variable name="max_cv_mass_2d" language="python">
51
from fluidity_tools import stat_parser as stat
52
max_cv_mass_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["CVMassMatrix"]["max"][-1]
54
<variable name="min_cv_mass_2d" language="python">
55
from fluidity_tools import stat_parser as stat
56
min_cv_mass_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["CVMassMatrix"]["min"][-1]
58
<variable name="sum_cv_mass_2d" language="python">
59
from fluidity_tools import stat_parser as stat
60
sum_cv_mass_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["SumCVMassMatrix"]["max"][-1]
62
<variable name="max_fe_lumped_mass_2d" language="python">
63
from fluidity_tools import stat_parser as stat
64
max_fe_lumped_mass_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["FELumpedMassMatrix"]["max"][-1]
66
<variable name="min_fe_lumped_mass_2d" language="python">
67
from fluidity_tools import stat_parser as stat
68
min_fe_lumped_mass_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["FELumpedMassMatrix"]["min"][-1]
70
<variable name="sum_fe_lumped_mass_2d" language="python">
71
from fluidity_tools import stat_parser as stat
72
sum_fe_lumped_mass_2d = stat("cv_mass_matrix_check_p1_2d.stat")["fluid"]["SumFELumpedMassMatrix"]["max"][-1]
75
<variable name="element_volume_3d" language="python">
76
from fluidity_tools import stat_parser as stat
77
element_volume_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["ElementVolume"]["max"][-1]
79
<variable name="max_cv_mass_3d" language="python">
80
from fluidity_tools import stat_parser as stat
81
max_cv_mass_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["CVMassMatrix"]["max"][-1]
83
<variable name="min_cv_mass_3d" language="python">
84
from fluidity_tools import stat_parser as stat
85
min_cv_mass_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["CVMassMatrix"]["min"][-1]
87
<variable name="sum_cv_mass_3d" language="python">
88
from fluidity_tools import stat_parser as stat
89
sum_cv_mass_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["SumCVMassMatrix"]["max"][-1]
91
<variable name="max_fe_lumped_mass_3d" language="python">
92
from fluidity_tools import stat_parser as stat
93
max_fe_lumped_mass_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["FELumpedMassMatrix"]["max"][-1]
95
<variable name="min_fe_lumped_mass_3d" language="python">
96
from fluidity_tools import stat_parser as stat
97
min_fe_lumped_mass_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["FELumpedMassMatrix"]["min"][-1]
99
<variable name="sum_fe_lumped_mass_3d" language="python">
100
from fluidity_tools import stat_parser as stat
101
sum_fe_lumped_mass_3d = stat("cv_mass_matrix_check_p1_3d.stat")["fluid"]["SumFELumpedMassMatrix"]["max"][-1]
105
<test name="check that the sum of the CV mass matrix is the element volume for 1d to tolerance 1.0e-10 " language="python">
106
print "Difference of element volume to integral CV mass for 1d: ",abs(element_volume_1d - sum_cv_mass_1d)
107
assert abs(element_volume_1d - sum_cv_mass_1d) < 1e-10
109
<test name="check that the sum of the FE lumped mass matrix is the element volume for 1d to tolerance 1.0e-10 " language="python">
110
print "Difference of element volume to integral FE lumped mass for 1d: ",abs(element_volume_1d - sum_fe_lumped_mass_1d)
111
assert abs(element_volume_1d - sum_fe_lumped_mass_1d) < 1e-10
113
<test name="check that the max of the CV mass matrix is the element volume divided by 2 for 1d to tolerance 1.0e-10 " language="python">
114
print "Difference of element volume divided by 2 to max CV mass for 1d: ",abs((element_volume_1d/2.0) - max_cv_mass_1d)
115
assert abs((element_volume_1d/2.0) - max_cv_mass_1d) < 1e-10
117
<test name="check that the max of the FE lumped mass matrix is the element volume divided by 2 for 1d to tolerance 1.0e-10 " language="python">
118
print "Difference of element volume divided by 2 to max FE lumped mass for 1d: ",abs((element_volume_1d/2.0) - max_fe_lumped_mass_1d)
119
assert abs((element_volume_1d/2.0) - max_fe_lumped_mass_1d) < 1e-10
121
<test name="check that the max of the CV mass matrix is equal to the min of the CV mass matrix for 1d to tolerance 1.0e-10 " language="python">
122
print "Difference of max CV mass matrix to min CV mass matrix for 1d: ",abs(max_cv_mass_1d - min_cv_mass_1d)
123
assert abs(max_cv_mass_1d - min_cv_mass_1d) < 1e-10
125
<test name="check that the max of the FE lumped mass matrix is equal to the min of the FE lumped mass matrix for 1d to tolerance 1.0e-10 " language="python">
126
print "Difference of max FE lumped mass matrix to min FE lumped mass matrix for 1d: ",abs(max_fe_lumped_mass_1d - min_fe_lumped_mass_1d)
127
assert abs(max_fe_lumped_mass_1d - min_fe_lumped_mass_1d) < 1e-10
129
<test name="check that the max of the CV mass matrix is equal to the max of the FE lumped mass matrix for 1d to tolerance 1.0e-10 " language="python">
130
print "Difference of max CV mass matrix to max FE lumped mass matrix for 1d: ",abs(max_cv_mass_1d - max_fe_lumped_mass_1d)
131
assert abs(max_cv_mass_1d - max_fe_lumped_mass_1d) < 1e-10
134
<test name="check that the sum of the CV mass matrix is the element volume for 2d to tolerance 1.0e-10 " language="python">
135
print "Difference of element volume to integral CV mass for 2d: ",abs(element_volume_2d - sum_cv_mass_2d)
136
assert abs(element_volume_2d - sum_cv_mass_2d) < 1e-10
138
<test name="check that the sum of the FE lumped mass matrix is the element volume for 2d to tolerance 1.0e-10 " language="python">
139
print "Difference of element volume to integral FE lumped mass for 2d: ",abs(element_volume_2d - sum_fe_lumped_mass_2d)
140
assert abs(element_volume_2d - sum_fe_lumped_mass_2d) < 1e-10
142
<test name="check that the max of the CV mass matrix is the element volume divided by 3 for 2d to tolerance 1.0e-10 " language="python">
143
print "Difference of element volume divided by 3 to max CV mass for 2d: ",abs((element_volume_2d/3.0) - max_cv_mass_2d)
144
assert abs((element_volume_2d/3.0) - max_cv_mass_2d) < 1e-10
146
<test name="check that the max of the FE lumped mass matrix is the element volume divided by 3 for 2d to tolerance 1.0e-10 " language="python">
147
print "Difference of element volume divided by 3 to max FE lumped mass for 2d: ",abs((element_volume_2d/3.0) - max_fe_lumped_mass_2d)
148
assert abs((element_volume_2d/3.0) - max_fe_lumped_mass_2d) < 1e-10
150
<test name="check that the max of the CV mass matrix is equal to the min of the CV mass matrix for 2d to tolerance 1.0e-10 " language="python">
151
print "Difference of max CV mass matrix to min CV mass matrix for 2d: ",abs(max_cv_mass_2d - min_cv_mass_2d)
152
assert abs(max_cv_mass_2d - min_cv_mass_2d) < 1e-10
154
<test name="check that the max of the FE lumped mass matrix is equal to the min of the FE lumped mass matrix for 2d to tolerance 1.0e-10 " language="python">
155
print "Difference of max FE lumped mass matrix to min FE lumped mass matrix for 2d: ",abs(max_fe_lumped_mass_2d - min_fe_lumped_mass_2d)
156
assert abs(max_fe_lumped_mass_2d - min_fe_lumped_mass_2d) < 1e-10
158
<test name="check that the max of the CV mass matrix is equal to the max of the FE lumped mass matrix for 2d to tolerance 1.0e-10 " language="python">
159
print "Difference of max CV mass matrix to max FE lumped mass matrix for 2d: ",abs(max_cv_mass_2d - max_fe_lumped_mass_2d)
160
assert abs(max_cv_mass_2d - max_fe_lumped_mass_2d) < 1e-10
163
<test name="check that the sum of the CV mass matrix is the element volume for 3d to tolerance 1.0e-10 " language="python">
164
print "Difference of element volume to integral CV mass for 3d: ",abs(element_volume_3d - sum_cv_mass_3d)
165
assert abs(element_volume_3d - sum_cv_mass_3d) < 1e-10
167
<test name="check that the sum of the FE lumped mass matrix is the element volume for 3d to tolerance 1.0e-10 " language="python">
168
print "Difference of element volume to integral FE lumped mass for 3d: ",abs(element_volume_3d - sum_fe_lumped_mass_3d)
169
assert abs(element_volume_3d - sum_fe_lumped_mass_3d) < 1e-10
171
<test name="check that the max of the CV mass matrix is the element volume divided by 4 for 3d to tolerance 1.0e-10 " language="python">
172
print "Difference of element volume divided by 4 to max CV mass for 3d: ",abs((element_volume_3d/4.0) - max_cv_mass_3d)
173
assert abs((element_volume_3d/4.0) - max_cv_mass_3d) < 1e-10
175
<test name="check that the max of the FE lumped mass matrix is the element volume divided by 4 for 3d to tolerance 1.0e-10 " language="python">
176
print "Difference of element volume divided by 4 to max FE lumped mass for 3d: ",abs((element_volume_3d/4.0) - max_fe_lumped_mass_3d)
177
assert abs((element_volume_3d/4.0) - max_fe_lumped_mass_3d) < 1e-10
179
<test name="check that the max of the CV mass matrix is equal to the min of the CV mass matrix for 3d to tolerance 1.0e-10 " language="python">
180
print "Difference of max CV mass matrix to min CV mass matrix for 3d: ",abs(max_cv_mass_3d - min_cv_mass_3d)
181
assert abs(max_cv_mass_3d - min_cv_mass_3d) < 1e-10
183
<test name="check that the max of the FE lumped mass matrix is equal to the min of the FE lumped mass matrix for 3d to tolerance 1.0e-10 " language="python">
184
print "Difference of max FE lumped mass matrix to min FE lumped mass matrix for 3d: ",abs(max_fe_lumped_mass_3d - min_fe_lumped_mass_3d)
185
assert abs(max_fe_lumped_mass_3d - min_fe_lumped_mass_3d) < 1e-10
187
<test name="check that the max of the CV mass matrix is equal to the max of the FE lumped mass matrix for 3d to tolerance 1.0e-10 " language="python">
188
print "Difference of max CV mass matrix to max FE lumped mass matrix for 3d: ",abs(max_cv_mass_3d - max_fe_lumped_mass_3d)
189
assert abs(max_cv_mass_3d - max_fe_lumped_mass_3d) < 1e-10