1
<?xml version="1.0" encoding="UTF-8" ?>
2
<!DOCTYPE testproblem SYSTEM "regressiontest.dtd">
5
<name>1 material column under gravity using CG pressure with continuity tested with the CV dual</name>
6
<owner userid="btollit"/>
7
<tags>flml cg_press_test_cty_cv</tags>
8
<problem_definition length="short" nprocs="1">
10
fluidity gravity-col-press_cg_test_cty_cv_p0p1_2d.flml
11
fluidity gravity-col-press_cg_test_cty_cv_p1dgp2_2d.flml
12
fluidity gravity-col-press_cg_test_cty_cv_p2lumpedp1_2d.flml
13
fluidity gravity-col-press_cg_test_cty_cv_p2p1_2d.flml
17
<variable name="pressure_p0p1_2d" language="python">
19
pressure_p0p1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p0p1_2d.stat")['Material1']['Pressure']['max'][-1]
21
<variable name="minpressure_p0p1_2d" language="python">
23
minpressure_p0p1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p0p1_2d.stat")['Material1']['Pressure']['min'][-1]
25
<variable name="velocity_p0p1_2d" language="python">
27
velocity_p0p1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p0p1_2d.stat")['Material1']['Velocity%magnitude']['max'][-1]
29
<variable name="pressure_p1dgp2_2d" language="python">
31
pressure_p1dgp2_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p1dgp2_2d.stat")['Material1']['Pressure']['max'][-1]
33
<variable name="minpressure_p1dgp2_2d" language="python">
35
minpressure_p1dgp2_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p1dgp2_2d.stat")['Material1']['Pressure']['min'][-1]
37
<variable name="velocity_p1dgp2_2d" language="python">
39
velocity_p1dgp2_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p1dgp2_2d.stat")['Material1']['Velocity%magnitude']['max'][-1]
41
<variable name="pressure_p2lumpedp1_2d" language="python">
43
pressure_p2lumpedp1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p2lumpedp1_2d.stat")['Material1']['Pressure']['max'][-1]
45
<variable name="minpressure_p2lumpedp1_2d" language="python">
47
minpressure_p2lumpedp1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p2lumpedp1_2d.stat")['Material1']['Pressure']['min'][-1]
49
<variable name="velocity_p2lumpedp1_2d" language="python">
51
velocity_p2lumpedp1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p2lumpedp1_2d.stat")['Material1']['Velocity%magnitude']['max'][-1]
53
<variable name="pressure_p2p1_2d" language="python">
55
pressure_p2p1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p2p1_2d.stat")['Material1']['Pressure']['max'][-1]
57
<variable name="minpressure_p2p1_2d" language="python">
59
minpressure_p2p1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p2p1_2d.stat")['Material1']['Pressure']['min'][-1]
61
<variable name="velocity_p2p1_2d" language="python">
63
velocity_p2p1_2d = fluidity_tools.stat_parser("gravity-col-press_cg_test_cty_cv_p2p1_2d.stat")['Material1']['Velocity%magnitude']['max'][-1]
65
<variable name = "solvers_converged" language = "python">
67
files = os.listdir("./")
68
solvers_converged = not "matrixdump" in files and not "matrixdump.info" in files
72
<test name="check maximum pressure_p0p1_2d is 10.0 with tolerance 1.0e-10" language="python">
73
print 'max pressure error: ',abs(pressure_p0p1_2d-1e+1)
74
assert abs(pressure_p0p1_2d-1e+1) < 1.0e-10
76
<test name="check pressure_p0p1_2d range is 10.0 with tolerance 1.0e-10" language="python">
77
print 'pressure range error: ',abs(pressure_p0p1_2d-minpressure_p0p1_2d-1e+1)
78
assert abs(pressure_p0p1_2d-minpressure_p0p1_2d-1e+1) < 1.0e-10
80
<test name="check maximum speed_p0p1_2d is zero with tolerance 1.0e-10" language="python">
81
assert abs(velocity_p0p1_2d) < 1.0e-10
83
<test name="check maximum pressure_p1dgp2_2d is 10.0 with tolerance 1.0e-10" language="python">
84
print 'max pressure error: ',abs(pressure_p1dgp2_2d-1e+1)
85
assert abs(pressure_p1dgp2_2d-1e+1) < 1.0e-10
87
<test name="check pressure_p1dgp2_2d range is 10.0 with tolerance 1.0e-10" language="python">
88
print 'pressure range error: ',abs(pressure_p1dgp2_2d-minpressure_p1dgp2_2d-1e+1)
89
assert abs(pressure_p1dgp2_2d-minpressure_p1dgp2_2d-1e+1) < 1.0e-10
91
<test name="check maximum speed_p1dgp2_2d is zero with tolerance 1.0e-10" language="python">
92
assert abs(velocity_p1dgp2_2d) < 1.0e-10
94
<test name="check maximum pressure_p2lumpedp1_2d is 10.0 with tolerance 1.0e-10" language="python">
95
print 'max pressure error: ',abs(pressure_p2lumpedp1_2d-1e+1)
96
assert abs(pressure_p2lumpedp1_2d-1e+1) < 1.0e-10
98
<test name="check pressure_p2lumpedp1_2d range is 10.0 with tolerance 1.0e-10" language="python">
99
print 'pressure range error: ',abs(pressure_p2lumpedp1_2d-minpressure_p2lumpedp1_2d-1e+1)
100
assert abs(pressure_p2lumpedp1_2d-minpressure_p2lumpedp1_2d-1e+1) < 1.0e-10
102
<test name="check maximum speed_p2lumpedp1_2d is zero with tolerance 1.0e-10" language="python">
103
assert abs(velocity_p2lumpedp1_2d) < 1.0e-10
105
<test name="check maximum pressure_p2p1_2d is 10.0 with tolerance 1.0e-10" language="python">
106
print 'max pressure error: ',abs(pressure_p2p1_2d-1e+1)
107
assert abs(pressure_p2p1_2d-1e+1) < 1.0e-10
109
<test name="check pressure_p2p1_2d range is 10.0 with tolerance 1.0e-10" language="python">
110
print 'pressure range error: ',abs(pressure_p2p1_2d-minpressure_p2p1_2d-1e+1)
111
assert abs(pressure_p2p1_2d-minpressure_p2p1_2d-1e+1) < 1.0e-10
113
<test name="check maximum speed_p2p1_2d is zero with tolerance 1.0e-10" language="python">
114
assert abs(velocity_p2p1_2d) < 1.0e-10
116
<test name = "Solvers converged" language = "python">
117
assert(solvers_converged)