← Back to branch summary

~fluidity-core/fluidity/shallow-water-dev

« back to all changes in this revision

Viewing changes to tests/cv_mass_matrix_check_p2/cv_mass_matrix_check_p2_1d.flml

  • Committer: colin.cotter at ac
  • Date: 2012-02-13 14:24:01 UTC
  • mfrom: (3565.1.362 fluidity)
  • Revision ID: colin.cotter@imperial.ac.uk-20120213142401-fz23s9jmvc8jjiy6
Merged from trunk.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/cv_mass_matrix_check_p2/cv_mass_matrix_check_p2_1d.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">cv_mass_matrix_check_p2_1d</string_value>
 
5
  </simulation_name>
 
6
  <problem_type>
 
7
    <string_value lines="1">fluids</string_value>
 
8
  </problem_type>
 
9
  <geometry>
 
10
    <dimension>
 
11
      <integer_value rank="0">1</integer_value>
 
12
    </dimension>
 
13
    <mesh name="CoordinateMesh">
 
14
      <from_file file_name="line">
 
15
        <format name="triangle"/>
 
16
        <stat>
 
17
          <include_in_stat/>
 
18
        </stat>
 
19
      </from_file>
 
20
    </mesh>
 
21
    <mesh name="DGP0">
 
22
      <from_mesh>
 
23
        <mesh name="CoordinateMesh"/>
 
24
        <mesh_shape>
 
25
          <polynomial_degree>
 
26
            <integer_value rank="0">0</integer_value>
 
27
          </polynomial_degree>
 
28
        </mesh_shape>
 
29
        <mesh_continuity>
 
30
          <string_value>discontinuous</string_value>
 
31
        </mesh_continuity>
 
32
        <stat>
 
33
          <exclude_from_stat/>
 
34
        </stat>
 
35
      </from_mesh>
 
36
    </mesh>
 
37
    <mesh name="CGP2">
 
38
      <from_mesh>
 
39
        <mesh name="CoordinateMesh"/>
 
40
        <mesh_shape>
 
41
          <polynomial_degree>
 
42
            <integer_value rank="0">2</integer_value>
 
43
          </polynomial_degree>
 
44
        </mesh_shape>
 
45
        <stat>
 
46
          <exclude_from_stat/>
 
47
        </stat>
 
48
      </from_mesh>
 
49
    </mesh>
 
50
    <mesh name="DGP2">
 
51
      <from_mesh>
 
52
        <mesh name="CoordinateMesh"/>
 
53
        <mesh_shape>
 
54
          <polynomial_degree>
 
55
            <integer_value rank="0">2</integer_value>
 
56
          </polynomial_degree>
 
57
        </mesh_shape>
 
58
        <mesh_continuity>
 
59
          <string_value>discontinuous</string_value>
 
60
        </mesh_continuity>
 
61
        <stat>
 
62
          <exclude_from_stat/>
 
63
        </stat>
 
64
      </from_mesh>
 
65
    </mesh>
 
66
    <quadrature>
 
67
      <degree>
 
68
        <integer_value rank="0">4</integer_value>
 
69
      </degree>
 
70
    </quadrature>
 
71
  </geometry>
 
72
  <io>
 
73
    <dump_format>
 
74
      <string_value>vtk</string_value>
 
75
    </dump_format>
 
76
    <dump_period_in_timesteps>
 
77
      <constant>
 
78
        <integer_value rank="0">0</integer_value>
 
79
      </constant>
 
80
    </dump_period_in_timesteps>
 
81
    <output_mesh name="DGP2"/>
 
82
    <stat/>
 
83
  </io>
 
84
  <timestepping>
 
85
    <current_time>
 
86
      <real_value rank="0">0.0</real_value>
 
87
    </current_time>
 
88
    <timestep>
 
89
      <real_value rank="0">1.0</real_value>
 
90
    </timestep>
 
91
    <finish_time>
 
92
      <real_value rank="0">1.0</real_value>
 
93
    </finish_time>
 
94
  </timestepping>
 
95
  <material_phase name="fluid">
 
96
    <vector_field name="Velocity" rank="1">
 
97
      <prescribed>
 
98
        <mesh name="CoordinateMesh"/>
 
99
        <value name="WholeMesh">
 
100
          <constant>
 
101
            <real_value shape="1" dim1="dim" rank="1">0.0</real_value>
 
102
          </constant>
 
103
        </value>
 
104
        <output/>
 
105
        <stat>
 
106
          <include_in_stat/>
 
107
        </stat>
 
108
        <detectors>
 
109
          <exclude_from_detectors/>
 
110
        </detectors>
 
111
      </prescribed>
 
112
    </vector_field>
 
113
    <scalar_field name="CVMassMatrix" rank="0">
 
114
      <diagnostic>
 
115
        <algorithm name="control_volume_mass_matrix" material_phase_support="single"/>
 
116
        <mesh name="CGP2"/>
 
117
        <output/>
 
118
        <stat/>
 
119
        <convergence>
 
120
          <include_in_convergence/>
 
121
        </convergence>
 
122
        <detectors>
 
123
          <include_in_detectors/>
 
124
        </detectors>
 
125
        <steady_state>
 
126
          <include_in_steady_state/>
 
127
        </steady_state>
 
128
      </diagnostic>
 
129
    </scalar_field>
 
130
    <scalar_field name="FELumpedMassMatrix" rank="0">
 
131
      <diagnostic>
 
132
        <algorithm name="finite_element_lumped_mass_matrix" material_phase_support="single"/>
 
133
        <mesh name="CGP2"/>
 
134
        <output/>
 
135
        <stat/>
 
136
        <convergence>
 
137
          <include_in_convergence/>
 
138
        </convergence>
 
139
        <detectors>
 
140
          <include_in_detectors/>
 
141
        </detectors>
 
142
        <steady_state>
 
143
          <include_in_steady_state/>
 
144
        </steady_state>
 
145
      </diagnostic>
 
146
    </scalar_field>
 
147
    <scalar_field name="ElementVolume" rank="0">
 
148
      <diagnostic>
 
149
        <algorithm name="scalar_python_diagnostic" material_phase_support="multiple">
 
150
          <string_value lines="20" type="python">coord = states["fluid"].vector_fields["Coordinate"]
 
151
 
 
152
for ele in range(field.element_count):
 
153
   
 
154
   # form element transform object
 
155
   tr_ele = Transform(ele, coord)
 
156
   
 
157
   # set volume via sum of detwei
 
158
   field.set(ele, sum(tr_ele.detwei))</string_value>
 
159
        </algorithm>
 
160
        <mesh name="DGP0"/>
 
161
        <output/>
 
162
        <stat/>
 
163
        <convergence>
 
164
          <include_in_convergence/>
 
165
        </convergence>
 
166
        <detectors>
 
167
          <include_in_detectors/>
 
168
        </detectors>
 
169
        <steady_state>
 
170
          <include_in_steady_state/>
 
171
        </steady_state>
 
172
      </diagnostic>
 
173
    </scalar_field>
 
174
    <scalar_field name="SumCVMassMatrix" rank="0">
 
175
      <diagnostic>
 
176
        <algorithm name="scalar_python_diagnostic" material_phase_support="multiple">
 
177
          <string_value lines="20" type="python">cvmass = states["fluid"].scalar_fields["CVMassMatrix"]
 
178
 
 
179
sum_cvmass = sum(cvmass.val)
 
180
 
 
181
for ele in range(field.element_count):
 
182
   
 
183
   field.set(ele, sum_cvmass)</string_value>
 
184
        </algorithm>
 
185
        <mesh name="DGP0"/>
 
186
        <output/>
 
187
        <stat/>
 
188
        <convergence>
 
189
          <include_in_convergence/>
 
190
        </convergence>
 
191
        <detectors>
 
192
          <include_in_detectors/>
 
193
        </detectors>
 
194
        <steady_state>
 
195
          <include_in_steady_state/>
 
196
        </steady_state>
 
197
      </diagnostic>
 
198
    </scalar_field>
 
199
    <scalar_field name="SumFELumpedMassMatrix" rank="0">
 
200
      <diagnostic>
 
201
        <algorithm name="scalar_python_diagnostic" material_phase_support="multiple">
 
202
          <string_value lines="20" type="python">felmass = states["fluid"].scalar_fields["FELumpedMassMatrix"]
 
203
 
 
204
sum_felmass = sum(felmass.val)
 
205
 
 
206
for ele in range(field.element_count):
 
207
   
 
208
   field.set(ele, sum_felmass)</string_value>
 
209
        </algorithm>
 
210
        <mesh name="DGP0"/>
 
211
        <output/>
 
212
        <stat/>
 
213
        <convergence>
 
214
          <include_in_convergence/>
 
215
        </convergence>
 
216
        <detectors>
 
217
          <include_in_detectors/>
 
218
        </detectors>
 
219
        <steady_state>
 
220
          <include_in_steady_state/>
 
221
        </steady_state>
 
222
      </diagnostic>
 
223
    </scalar_field>
 
224
  </material_phase>
 
225
</fluidity_options>