1
<?xml version="1.0" encoding="UTF-8" ?>
2
<!DOCTYPE testproblem SYSTEM "regressiontest.dtd">
5
<name>heat_2d_dg</name>
6
<owner userid="wence"/>
7
<tags>flml parallel</tags>
8
<problem_definition length="short" nprocs="1">
9
<command_line>mpiexec -np 16 ../../bin/fluidity -v2 -l heat.flml </command_line>
10
<!-- Two dimensional diffusion only problem. -->
13
<variable name="tracer1_integral0" language="python">
14
from fluidity_tools import stat_parser
15
s = stat_parser("heat.stat")
16
tracer1_integral0 = s["Water"]["Velocity%1"]["integral"][0]
18
<variable name="tracer2_integral0" language="python">
19
from fluidity_tools import stat_parser
20
s = stat_parser("heat.stat")
21
tracer2_integral0 = s["Water"]["Velocity%2"]["integral"][0]
23
<variable name="tracer1_integral" language="python">
24
from fluidity_tools import stat_parser
25
s = stat_parser("heat.stat")
26
tracer1_integral = s["Water"]["Velocity%1"]["integral"][-1]
28
<variable name="tracer2_integral" language="python">
29
from fluidity_tools import stat_parser
30
s = stat_parser("heat.stat")
31
tracer2_integral = s["Water"]["Velocity%2"]["integral"][-1]
33
<variable name="tracer1_error0" language="python">
34
from fluidity_tools import stat_parser
35
s = stat_parser("heat.stat")
36
tracer1_error0 = s["Water"]["VectorAbsoluteDifference%1"]["max"][0]
38
<variable name="tracer1_error" language="python">
39
from fluidity_tools import stat_parser
40
s = stat_parser("heat.stat")
41
tracer1_error = s["Water"]["VectorAbsoluteDifference%1"]["max"][-1]
43
<variable name="tracer2_error0" language="python">
44
from fluidity_tools import stat_parser
45
s = stat_parser("heat.stat")
46
tracer2_error0 = s["Water"]["VectorAbsoluteDifference%2"]["max"][0]
48
<variable name="tracer2_error" language="python">
49
from fluidity_tools import stat_parser
50
s = stat_parser("heat.stat")
51
tracer2_error = s["Water"]["VectorAbsoluteDifference%2"]["max"][-1]
55
<test name="tracer1 integral is 1.0" language="python">
56
assert abs(tracer1_integral - 1.0) < 1e-6
58
<test name="tracer2 integral is 1.0" language="python">
59
assert abs(tracer2_integral - 1.0) < 1e-6
61
<test name="initial tracer1 integral is 1.0" language="python">
62
assert abs(tracer1_integral0 - 1.0) < 1e-3
64
<test name="initial tracer2 integral is 1.0" language="python">
65
assert abs(tracer2_integral0 - 1.0) < 1e-3
67
<test name="initial tracer1 error is small relative to serial" language="python">
68
assert abs(tracer1_error0 - 0.00877193063655) < 1e-6
70
<test name="tracer1 error is small relative to serial" language="python">
71
assert abs(tracer1_error - 0.00877192982458) < 5e-9
73
<test name="initial tracer2 error is small relative to serial" language="python">
74
assert abs(tracer2_error0 - 0.00877193063655) < 1e-6
76
<test name="tracer2 error is small relative to serial" language="python">
77
assert abs(tracer2_error - 0.00877192982458) < 5e-9