fluidity/tests/explicit-hyperc-tets/explicit-hyperc-tets.xml at 5e8eef62f6c969c11cc54cfe4e7fc663f0c318c8 · FluidityProject/fluidity · GitHub

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
<?xml version="1.0" encoding="UTF-8" ?>

<testproblem>
  <name>Control volume, multi-material flow test</name>
  <owner userid="cwilson"/>
  <tags>flml solidity</tags>
    <problem_definition length="short" nprocs="1">
    <command_line>fluidity  channel_flow_input.flml || true &amp;&amp; fluidity  channel_flow_1_checkpoint.flml </command_line>
  </problem_definition>
  <variables>
    <variable name="solvers_converged" language="python">
import os
files = os.listdir("./")
solvers_converged = not "matrixdump" in files and not "matrixdump.info" in files
    </variable>
     <variable name="mxvfrac0" language="python">
import fluidity_tools
mxvfrac0 = fluidity_tools.stat_parser("channel_flow.stat")['Water']['Tracer']['max'][0]
    </variable>
     <variable name="mnvfrac0" language="python">
import fluidity_tools
mnvfrac0 = fluidity_tools.stat_parser("channel_flow.stat")['Water']['Tracer']['min'][0]
    </variable>
    <variable name="mxvfrac" language="python">
import fluidity_tools
mxvfrac = fluidity_tools.stat_parser("channel_flow.stat")['Water']['Tracer']['max'][-1]
    </variable>
     <variable name="mnvfrac" language="python">
import fluidity_tools
mnvfrac = fluidity_tools.stat_parser("channel_flow.stat")['Water']['Tracer']['min'][-1]
    </variable>
    <variable name="tracerbehindfront" language="python">
import vtktools
from numpy import array
coordinates=array([[0.005, 0.2, 0.005]])
vtufile=vtktools.vtu('channel_flow_3.vtu')
tracerbehindfront=vtktools.vtu.ProbeData(vtufile,coordinates,'Tracer')
    </variable>
     <variable name="tracerfrontfront" language="python">
import vtktools
from numpy import array
coordinates=array([[0.005, 0.3, 0.005]])
vtufile=vtktools.vtu('channel_flow_3.vtu')
tracerfrontfront=vtktools.vtu.ProbeData(vtufile,coordinates,'Tracer')
    </variable>
    <variable name="mxvelocity" language="python">
import fluidity_tools
mxvelocity = fluidity_tools.stat_parser("channel_flow.stat")['Water']['Velocity%magnitude']['max'][-1]
    </variable>
    <variable name="mxcfl" language="python">
import fluidity_tools
mxcfl = fluidity_tools.stat_parser("channel_flow.stat")['Water']['ControlVolumeCFLNumber']['max'][-1]
    </variable>
    <variable name="mncfl" language="python">
import fluidity_tools
mncfl = fluidity_tools.stat_parser("channel_flow.stat")['Water']['ControlVolumeCFLNumber']['min'][-1]
    </variable>
    <variable name="checkpointtracer" language="python">
import vtktools
vtu=vtktools.vtu('channel_flow_3.vtu')
vtu_checkpoint=vtktools.vtu('channel_flow_checkpoint_2.vtu')
vtu_diff=vtktools.VtuDiff(vtu,vtu_checkpoint)
checkpointtracer=max(vtu_diff.GetScalarRange("Tracer"))
    </variable>
  </variables>
    <pass_tests>
      <test name="Solvers converged" language="python">
        assert(solvers_converged)
      </test>
    <test name="check minimum volume fraction is 0.0" language="python">
   assert abs(mnvfrac) &lt; 1e-10
    </test>
    <test name="check maximum volume fraction is 1.0" language="python">
   assert abs(mxvfrac-1.0) &lt; 1e-10
    </test>
    <test name="check volume fraction in front of front is 0.0 (relaxed for simplex elements)" language="python">
   assert abs(tracerfrontfront) &lt; 1e-2
    </test>
    <test name="check volume fraction behind front is 1.0 (relaxed for simplex elements)" language="python">
   assert abs(tracerbehindfront-1.0) &lt; 1e-5
    </test>
    <test name="check initial minimum volume fraction is 0.0" language="python">
   assert abs(mnvfrac0) &lt; 1e-10
    </test>
    <test name="check initial maximum volume fraction is 1.0" language="python">
   assert abs(mxvfrac0-1.0) &lt; 1e-10
    </test>
    <test name="check maximum speed is 1.0" language="python">
	assert abs(mxvelocity-1.0) &lt; 1e-10
    </test>
    <test name="check maximum cv cfl" language="python">
      assert abs(mxcfl-1.6) &lt; 1e-10
    </test>
    <test name="check minimum cv cfl (lower bound)" language="python">
      assert mncfl-0.559 &gt; 1e-10
    </test>
    <test name="check difference is zero to checkpointed run" language="python">
      assert abs(checkpointtracer) &lt; 1e-10
    </test>
  </pass_tests>
  <warn_tests>
  </warn_tests>
</testproblem>