-
Notifications
You must be signed in to change notification settings - Fork 118
Expand file tree
/
Copy pathgauss_point_buoyancy.xml
More file actions
77 lines (77 loc) · 6.58 KB
/
Copy pathgauss_point_buoyancy.xml
File metadata and controls
77 lines (77 loc) · 6.58 KB
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
<?xml version='1.0' encoding='utf-8'?>
<testproblem>
<name>Test of buoyancy at gauss points.</name>
<owner userid="cwilson"/>
<problem_definition length="short" nprocs="1">
<command_line>fluidity -v3 -l p1dgp2_nodal_2d.flml; fluidity -v3 -l p1dgp2_gauss_2d.flml; fluidity -v3 -l p2p1_nodal_2d.flml; fluidity -v3 -l p2p1_gauss_2d.flml; fluidity -v3 -l p1dgp2_nodal_3d.flml; fluidity -v3 -l p1dgp2_gauss_3d.flml; fluidity -v3 -l p2p1_nodal_3d.flml; fluidity -v3 -l p2p1_gauss_3d.flml</command_line>
</problem_definition>
<variables>
<variable name="p1dgp2_nodal_2d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_nodal_2d.stat")
p1dgp2_nodal_2d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p1dgp2_gauss_2d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_gauss_2d.stat")
p1dgp2_gauss_2d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p1dgp2_nodal_2d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_nodal_2d.stat")
p1dgp2_nodal_2d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p1dgp2_gauss_2d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_gauss_2d.stat")
p1dgp2_gauss_2d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p2p1_nodal_2d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_nodal_2d.stat")
p2p1_nodal_2d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p2p1_gauss_2d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_gauss_2d.stat")
p2p1_gauss_2d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p2p1_nodal_2d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_nodal_2d.stat")
p2p1_nodal_2d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p2p1_gauss_2d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_gauss_2d.stat")
p2p1_gauss_2d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p1dgp2_nodal_3d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_nodal_3d.stat")
p1dgp2_nodal_3d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p1dgp2_gauss_3d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_gauss_3d.stat")
p1dgp2_gauss_3d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p1dgp2_nodal_3d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_nodal_3d.stat")
p1dgp2_nodal_3d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p1dgp2_gauss_3d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p1dgp2_gauss_3d.stat")
p1dgp2_gauss_3d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p2p1_nodal_3d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_nodal_3d.stat")
p2p1_nodal_3d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p2p1_gauss_3d_pressure_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_gauss_3d.stat")
p2p1_gauss_3d_pressure_error = stat["State"]["AbsoluteDifferencePressure"]["l2norm"][-1]</variable>
<variable name="p2p1_nodal_3d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_nodal_3d.stat")
p2p1_nodal_3d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
<variable name="p2p1_gauss_3d_velocity_error" language="python">from fluidity_tools import stat_parser
stat = stat_parser("p2p1_gauss_3d.stat")
p2p1_gauss_3d_velocity_error = stat["State"]["Velocity%magnitude"]["l2norm"][-1]</variable>
</variables>
<pass_tests>
<test name="p1dgp2_nodal_2d_pressure" language="python">assert abs(p1dgp2_nodal_2d_pressure_error - 1.5e-5) < 1.e-6</test>
<test name="p1dgp2_gauss_2d_pressure" language="python">assert p1dgp2_gauss_2d_pressure_error/p1dgp2_nodal_2d_pressure_error < 1.e-6</test>
<test name="p1dgp2_nodal_2d_velocity" language="python">assert abs(p1dgp2_nodal_2d_velocity_error - 1.5e-5) < 9.9e-5</test>
<test name="p1dgp2_gauss_2d_velocity" language="python">assert p1dgp2_gauss_2d_velocity_error/p1dgp2_nodal_2d_velocity_error < 1.e-6</test>
<test name="p2p1_nodal_2d_pressure" language="python">assert abs(p2p1_nodal_2d_pressure_error - 0.017) < 1.e-3</test>
<test name="p2p1_gauss_2d_pressure" language="python">assert p2p1_gauss_2d_pressure_error/p2p1_nodal_2d_pressure_error < 1.e-6</test>
<test name="p2p1_nodal_2d_velocity" language="python">assert abs(p2p1_nodal_2d_velocity_error - 0.035) < 1.e-3</test>
<test name="p2p1_gauss_2d_velocity" language="python">assert p2p1_gauss_2d_velocity_error/p2p1_nodal_2d_velocity_error < 1.e-6</test>
<test name="p1dgp2_nodal_3d_pressure" language="python">assert abs(p1dgp2_nodal_3d_pressure_error - 0.018) < 2.e-3</test>
<test name="p1dgp2_gauss_3d_pressure" language="python">assert p1dgp2_gauss_3d_pressure_error/p1dgp2_nodal_3d_pressure_error < 1.e-6</test>
<test name="p1dgp2_nodal_3d_velocity" language="python">assert abs(p1dgp2_nodal_3d_velocity_error - 0.15) < 1.e-2</test>
<test name="p1dgp2_gauss_3d_velocity" language="python">assert p1dgp2_gauss_3d_velocity_error/p1dgp2_nodal_3d_velocity_error < 1.e-6</test>
<test name="p2p1_nodal_3d_pressure" language="python">assert abs(p2p1_nodal_3d_pressure_error - 0.31) < 1.5e-2</test>
<test name="p2p1_gauss_3d_pressure" language="python">assert p2p1_gauss_3d_pressure_error/p2p1_nodal_3d_pressure_error < 1.e-6</test>
<test name="p2p1_nodal_3d_velocity" language="python">assert abs(p2p1_nodal_3d_velocity_error - 1.31) < 2.e-2</test>
<test name="p2p1_gauss_3d_velocity" language="python">assert p2p1_gauss_3d_velocity_error/p2p1_nodal_3d_velocity_error < 1.e-6</test>
</pass_tests>
<warn_tests/>
</testproblem>