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
|
<?xml version='1.0' encoding='utf-8'?>
<testproblem>
<name>Tracer conservation when moving the mesh</name>
<owner userid="piggott"/>
<tags>flml</tags>
<problem_definition length="medium" nprocs="1">
<command_line>fluidity -v2 -l sloshing_tank_tracer.flml</command_line>
</problem_definition>
<variables>
<variable name="area_cons_error" language="python">from fluidity_tools import stat_parser as stat
area_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Density"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Density"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Density"]["integral"][0]</variable>
<variable name="tracercg_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracercg_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerCG"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerCG"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerCG"]["integral"][0]</variable>
<variable name="tracercgl_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracercgl_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerCGlumped"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerCGlumped"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerCGlumped"]["integral"][0]</variable>
<variable name="tracercgibp_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracercgibp_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerCGibp"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerCGibp"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerCGibp"]["integral"][0]</variable>
<variable name="tracerdgibp1_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracerdgibp1_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp1"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp1"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp1"]["integral"][0]</variable>
<variable name="tracerdgibp2_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracerdgibp2_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp2"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp2"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp2"]["integral"][0]</variable>
<variable name="tracerdgibp1icbp_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracerdgibp1icbp_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp1icbp"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp1icbp"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp1icbp"]["integral"][0]</variable>
<variable name="tracerdgibp2icbp_cons_error" language="python">from fluidity_tools import stat_parser as stat
tracerdgibp2icbp_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp2icbp"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp2icbp"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerDGibp2icbp"]["integral"][0]</variable>
<variable name="tracercv_cons_error" language="python">
from fluidity_tools import stat_parser as stat
tracercv_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["TracerCV"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["TracerCV"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["TracerCV"]["integral"][0]</variable>
<variable name="beta0tracercg_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracercg_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCG"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCG"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCG"]["integral"][0]</variable>
<variable name="beta0tracercgl_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracercgl_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCGlumped"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCGlumped"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCGlumped"]["integral"][0]</variable>
<variable name="beta0tracercgibp_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracercgibp_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCGibp"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCGibp"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCGibp"]["integral"][0]</variable>
<variable name="beta0tracerdgibp1_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracerdgibp1_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp1"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp1"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp1"]["integral"][0]</variable>
<variable name="beta0tracerdgibp2_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracerdgibp2_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp2"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp2"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp2"]["integral"][0]</variable>
<variable name="beta0tracerdgibp1icbp_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracerdgibp1icbp_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp1icbp"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp1icbp"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp1icbp"]["integral"][0]</variable>
<variable name="beta0tracerdgibp2icbp_cons_error" language="python">from fluidity_tools import stat_parser as stat
beta0tracerdgibp2icbp_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp2icbp"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp2icbp"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerDGibp2icbp"]["integral"][0]</variable>
<variable name="beta0tracercv_cons_error" language="python">
from fluidity_tools import stat_parser as stat
beta0tracercv_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["integral"][0]</variable>
<variable name="beta0tracerhyperc_cons_error" language="python">
from fluidity_tools import stat_parser as stat
beta0tracerhyperc_cons_error = max(abs(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["integral"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["integral"][0]))*100/stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["integral"][0]</variable>
<variable name="beta0tracercv_max_diff_max" language="python">
from fluidity_tools import stat_parser as stat
beta0tracercv_max_diff_max = max(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["max"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["max"][0])</variable>
<variable name="beta0tracercv_max_diff_min" language="python">
from fluidity_tools import stat_parser as stat
beta0tracercv_max_diff_min = min(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["max"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerCV"]["max"][0])</variable>
<variable name="beta0tracerhyperc_max_diff_max" language="python">
from fluidity_tools import stat_parser as stat
beta0tracerhyperc_max_diff_max = max(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["max"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["max"][0])</variable>
<variable name="beta0tracerhyperc_max_diff_min" language="python">
from fluidity_tools import stat_parser as stat
beta0tracerhyperc_max_diff_min = min(stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["max"]-stat("sloshing_tank_tracer.stat")["water"]["Beta0TracerHyperC"]["max"][0])</variable>
<variable name="momucg_cons_error" language="python">
from fluidity_tools import stat_parser as stat
statf = stat("sloshing_tank_tracer.stat")["water"]
initial = statf["Velocity%1"]["integral"][0] - statf["Pressure"]["surface_integral%LeftSide"][0] + statf["Pressure"]["surface_integral%RightSide"][0]
all = statf["Velocity%1"]["integral"] - statf["Pressure"]["surface_integral%LeftSide"] + statf["Pressure"]["surface_integral%RightSide"]
momucg_cons_error = max(abs(all-initial))*100/initial</variable>
<variable name="momvcg_cons_error" language="python">
from fluidity_tools import stat_parser as stat
statf = stat("sloshing_tank_tracer.stat")["water"]
initial = statf["Velocity%2"]["integral"][0] - statf["Pressure"]["surface_integral%Bottom"][0]
all = statf["Velocity%2"]["integral"] - statf["Pressure"]["surface_integral%Bottom"]
momvcg_cons_error = max(abs(all-initial))*100/initial</variable>
</variables>
<pass_tests>
<test name="area conservation error" language="python">
assert(area_cons_error < 1.E-7)</test>
<test name="cg conservation error" language="python">
assert(tracercg_cons_error < 1.E-5)</test>
<test name="cgl conservation error" language="python">
assert(tracercgl_cons_error < 1.E-7)</test>
<test name="cg ibp conservation error" language="python">
assert(tracercgibp_cons_error < 1.E-5)</test>
<test name="dg ibp1 conservation error" language="python">
assert(tracerdgibp1_cons_error < 1.E-8)</test>
<test name="dg ibp2 conservation error" language="python">
assert(tracerdgibp2_cons_error < 1.E-8)</test>
<test name="dg ibp1icbp conservation error" language="python">
assert(tracerdgibp1icbp_cons_error < 1.E-8)</test>
<test name="dg ibp2icbp conservation error" language="python">
assert(tracerdgibp2icbp_cons_error < 1.E-8)</test>
<test name="cv conservation error" language="python">
assert(tracercv_cons_error < 1.E-6)</test>
<test name="beta0 cg conservation error" language="python">
assert(beta0tracercg_cons_error > 1.E-5)</test>
<test name="beta0 cgl conservation error" language="python">
assert(beta0tracercgl_cons_error > 1.E-7)</test>
<test name="beta0 cg ibp conservation error" language="python">
assert(beta0tracercgibp_cons_error > 1.E-5)</test>
<test name="beta0 dg ibp1 conservation error" language="python">
assert(beta0tracerdgibp1_cons_error > 1.E-8)</test>
<test name="beta0 dg ibp2 conservation error" language="python">
assert(beta0tracerdgibp2_cons_error > 1.E-8)</test>
<test name="beta0 dg ibp1icbp conservation error" language="python">
assert(beta0tracerdgibp1icbp_cons_error > 1.E-8)</test>
<test name="beta0 dg ibp2icbp conservation error" language="python">
assert(beta0tracerdgibp2icbp_cons_error > 1.E-8)</test>
<test name="beta0 cv conservation error" language="python">
assert(beta0tracercv_cons_error > 1.E-6)</test>
<!-- <test name="momu cg conservation error" language="python">
assert(abs(momucg_cons_error) < 1.E-6)</test>
<test name="momv cg conservation error" language="python">
assert(abs(momvcg_cons_error) < 1.E-6)</test>-->
</pass_tests>
<warn_tests/>
</testproblem>
|