1
=== removed directory 'tests/2d_periodic_p1p1_streamfunction'
2
=== removed file 'tests/2d_periodic_p1p1_streamfunction/2d_nonperiodic_p1p1_streamfunction.flml'
3
--- tests/2d_periodic_p1p1_streamfunction/2d_nonperiodic_p1p1_streamfunction.flml	2010-11-21 15:41:35 +0000
4
+++ tests/2d_periodic_p1p1_streamfunction/2d_nonperiodic_p1p1_streamfunction.flml	1970-01-01 00:00:00 +0000
5
@@ -1,159 +0,0 @@
6
1
<?xml version='1.0' encoding='utf-8'?>
7
2
<fluidity_options>
8
3
  <simulation_name>
9
4
    <string_value lines="1">2d_nonperiodic_p1p1_streamfunction</string_value>
10
5
  </simulation_name>
11
6
  <problem_type>
12
7
    <string_value lines="1">fluids</string_value>
13
8
  </problem_type>
14
9
  <geometry>
15
10
    <dimension>
16
11
      <integer_value rank="0">2</integer_value>
17
12
    </dimension>
18
13
    <mesh name="CoordinateMesh">
19
14
      <from_file file_name="2d_square">
20
15
        <format name="triangle"/>
21
16
        <stat>
22
17
          <include_in_stat/>
23
18
        </stat>
24
19
      </from_file>
25
20
    </mesh>
26
21
    <mesh name="VelocityMesh">
27
22
      <from_mesh>
28
23
        <mesh name="CoordinateMesh"/>
29
24
        <stat>
30
25
          <exclude_from_stat/>
31
26
        </stat>
32
27
      </from_mesh>
33
28
    </mesh>
34
29
    <mesh name="PressureMesh">
35
30
      <from_mesh>
36
31
        <mesh name="CoordinateMesh"/>
37
32
        <stat>
38
33
          <exclude_from_stat/>
39
34
        </stat>
40
35
      </from_mesh>
41
36
    </mesh>
42
37
    <quadrature>
43
38
      <degree>
44
39
        <integer_value rank="0">4</integer_value>
45
40
      </degree>
46
41
    </quadrature>
47
42
  </geometry>
48
43
  <io>
49
44
    <dump_format>
50
45
      <string_value>vtk</string_value>
51
46
    </dump_format>
52
47
    <dump_period>
53
48
      <constant>
54
49
        <real_value rank="0">1e27</real_value>
55
50
      </constant>
56
51
    </dump_period>
57
52
    <output_mesh name="CoordinateMesh"/>
58
53
    <stat/>
59
54
  </io>
60
55
  <timestepping>
61
56
    <current_time>
62
57
      <real_value rank="0">0.0</real_value>
63
58
    </current_time>
64
59
    <timestep>
65
60
      <real_value rank="0">1</real_value>
66
61
    </timestep>
67
62
    <finish_time>
68
63
      <real_value rank="0">20</real_value>
69
64
    </finish_time>
70
65
    <nonlinear_iterations>
71
66
      <integer_value rank="0">2</integer_value>
72
67
    </nonlinear_iterations>
73
68
  </timestepping>
74
69
  <physical_parameters/>
75
70
  <material_phase name="water">
76
71
    <equation_of_state>
77
72
      <fluids>
78
73
        <linear>
79
74
          <reference_density>
80
75
            <real_value rank="0">1.0</real_value>
81
76
          </reference_density>
82
77
        </linear>
83
78
      </fluids>
84
79
    </equation_of_state>
85
80
    <vector_field name="Velocity" rank="1">
86
81
      <prescribed>
87
82
        <mesh name="VelocityMesh"/>
88
83
        <value name="WholeMesh">
89
84
          <python>
90
85
            <string_value lines="20" type="python">def val(X,t):
91
86
  return [X[1], 0.0]</string_value>
92
87
          </python>
93
88
        </value>
94
89
        <output/>
95
90
        <stat>
96
91
          <include_in_stat/>
97
92
        </stat>
98
93
        <detectors>
99
94
          <exclude_from_detectors/>
100
95
        </detectors>
101
96
      </prescribed>
102
97
    </vector_field>
103
98
    <scalar_field depends="Velocity" name="MultiplyConnectedStreamFunction" rank="0">
104
99
      <prognostic>
105
100
        <mesh name="VelocityMesh"/>
106
101
        <solver>
107
102
          <iterative_method name="gmres">
108
103
            <restart>
109
104
              <integer_value rank="0">30</integer_value>
110
105
            </restart>
111
106
          </iterative_method>
112
107
          <preconditioner name="sor"/>
113
108
          <relative_error>
114
109
            <real_value rank="0">1e-7</real_value>
115
110
          </relative_error>
116
111
          <max_iterations>
117
112
            <integer_value rank="0">3000</integer_value>
118
113
          </max_iterations>
119
114
          <never_ignore_solver_failures/>
120
115
          <diagnostics>
121
116
            <monitors/>
122
117
          </diagnostics>
123
118
        </solver>
124
119
        <boundary_conditions name="bottom">
125
120
          <surface_ids>
126
121
            <integer_value shape="1" rank="1">7</integer_value>
127
122
          </surface_ids>
128
123
          <primary_boundary/>
129
124
          <type name="dirichlet">
130
125
            <internally_calculated/>
131
126
          </type>
132
127
        </boundary_conditions>
133
128
        <boundary_conditions name="top">
134
129
          <surface_ids>
135
130
            <integer_value shape="1" rank="1">9</integer_value>
136
131
          </surface_ids>
137
132
          <secondary_boundary>
138
133
            <primary_point>
139
134
              <real_value shape="2" dim1="dim" rank="1">0.5 0.0</real_value>
140
135
            </primary_point>
141
136
            <secondary_point>
142
137
              <real_value shape="2" dim1="dim" rank="1">0.5 1.0</real_value>
143
138
            </secondary_point>
144
139
          </secondary_boundary>
145
140
          <type name="dirichlet">
146
141
            <internally_calculated/>
147
142
          </type>
148
143
        </boundary_conditions>
149
144
        <exclude_from_checkpointing/>
150
145
        <output>
151
146
          <include_previous_time_step/>
152
147
        </output>
153
148
        <stat/>
154
149
        <detectors>
155
150
          <include_in_detectors/>
156
151
        </detectors>
157
152
        <steady_state>
158
153
          <include_in_steady_state/>
159
154
        </steady_state>
160
155
        <no_interpolation/>
161
156
      </prognostic>
162
157
    </scalar_field>
163
158
  </material_phase>
164
159
</fluidity_options>
165
160
0
166
=== removed file 'tests/2d_periodic_p1p1_streamfunction/2d_periodic_p1p1_streamfunction.flml'
167
--- tests/2d_periodic_p1p1_streamfunction/2d_periodic_p1p1_streamfunction.flml	2010-11-21 15:41:35 +0000
168
+++ tests/2d_periodic_p1p1_streamfunction/2d_periodic_p1p1_streamfunction.flml	1970-01-01 00:00:00 +0000
169
@@ -1,204 +0,0 @@
170
1
<?xml version='1.0' encoding='utf-8'?>
171
2
<fluidity_options>
172
3
  <simulation_name>
173
4
    <string_value lines="1">2d_periodic_p1p1_streamfunction</string_value>
174
5
  </simulation_name>
175
6
  <problem_type>
176
7
    <string_value lines="1">fluids</string_value>
177
8
  </problem_type>
178
9
  <geometry>
179
10
    <dimension>
180
11
      <integer_value rank="0">2</integer_value>
181
12
    </dimension>
182
13
    <mesh name="CoordinateMesh">
183
14
      <from_mesh>
184
15
        <mesh name="periodic_boundaries"/>
185
16
        <periodic_boundary_conditions name="east_west">
186
17
          <physical_boundary_ids>
187
18
            <integer_value shape="1" rank="1">10</integer_value>
188
19
          </physical_boundary_ids>
189
20
          <aliased_boundary_ids>
190
21
            <integer_value shape="1" rank="1">8</integer_value>
191
22
          </aliased_boundary_ids>
192
23
          <coordinate_map>
193
24
            <string_value lines="20" type="python">def val(x,t):
194
25
	return [x[0]-1,x[1]]</string_value>
195
26
          </coordinate_map>
196
27
          <inverse_coordinate_map>
197
28
            <string_value lines="20" type="python">def val(x,t):
198
29
	return [x[0]+1,x[1]]</string_value>
199
30
          </inverse_coordinate_map>
200
31
          <remove_periodicity/>
201
32
        </periodic_boundary_conditions>
202
33
        <stat>
203
34
          <exclude_from_stat/>
204
35
        </stat>
205
36
      </from_mesh>
206
37
    </mesh>
207
38
    <mesh name="VelocityMesh">
208
39
      <from_mesh>
209
40
        <mesh name="periodic_boundaries"/>
210
41
        <stat>
211
42
          <exclude_from_stat/>
212
43
        </stat>
213
44
      </from_mesh>
214
45
    </mesh>
215
46
    <mesh name="PressureMesh">
216
47
      <from_mesh>
217
48
        <mesh name="periodic_boundaries"/>
218
49
        <stat>
219
50
          <exclude_from_stat/>
220
51
        </stat>
221
52
      </from_mesh>
222
53
    </mesh>
223
54
    <mesh name="2D_mesh">
224
55
      <from_file file_name="2d_square">
225
56
        <format name="triangle"/>
226
57
        <stat>
227
58
          <include_in_stat/>
228
59
        </stat>
229
60
      </from_file>
230
61
    </mesh>
231
62
    <mesh name="periodic_boundaries">
232
63
      <from_mesh>
233
64
        <mesh name="2D_mesh"/>
234
65
        <periodic_boundary_conditions name="east_west">
235
66
          <physical_boundary_ids>
236
67
            <integer_value shape="1" rank="1">10</integer_value>
237
68
          </physical_boundary_ids>
238
69
          <aliased_boundary_ids>
239
70
            <integer_value shape="1" rank="1">8</integer_value>
240
71
          </aliased_boundary_ids>
241
72
          <coordinate_map>
242
73
            <string_value lines="20" type="python">def val(x,t):
243
74
	return [x[0]-1,x[1]]</string_value>
244
75
          </coordinate_map>
245
76
        </periodic_boundary_conditions>
246
77
        <stat>
247
78
          <exclude_from_stat/>
248
79
        </stat>
249
80
      </from_mesh>
250
81
    </mesh>
251
82
    <quadrature>
252
83
      <degree>
253
84
        <integer_value rank="0">4</integer_value>
254
85
      </degree>
255
86
    </quadrature>
256
87
  </geometry>
257
88
  <io>
258
89
    <dump_format>
259
90
      <string_value>vtk</string_value>
260
91
    </dump_format>
261
92
    <dump_period>
262
93
      <constant>
263
94
        <real_value rank="0">1e27</real_value>
264
95
      </constant>
265
96
    </dump_period>
266
97
    <output_mesh name="CoordinateMesh"/>
267
98
    <stat/>
268
99
  </io>
269
100
  <timestepping>
270
101
    <current_time>
271
102
      <real_value rank="0">0.0</real_value>
272
103
    </current_time>
273
104
    <timestep>
274
105
      <real_value rank="0">1</real_value>
275
106
    </timestep>
276
107
    <finish_time>
277
108
      <real_value rank="0">20</real_value>
278
109
    </finish_time>
279
110
    <nonlinear_iterations>
280
111
      <integer_value rank="0">2</integer_value>
281
112
    </nonlinear_iterations>
282
113
  </timestepping>
283
114
  <physical_parameters/>
284
115
  <material_phase name="water">
285
116
    <equation_of_state>
286
117
      <fluids>
287
118
        <linear>
288
119
          <reference_density>
289
120
            <real_value rank="0">1.0</real_value>
290
121
          </reference_density>
291
122
        </linear>
292
123
      </fluids>
293
124
    </equation_of_state>
294
125
    <vector_field name="Velocity" rank="1">
295
126
      <prescribed>
296
127
        <mesh name="VelocityMesh"/>
297
128
        <value name="WholeMesh">
298
129
          <python>
299
130
            <string_value lines="20" type="python">def val(X,t):
300
131
  return [X[1], 0.0]</string_value>
301
132
          </python>
302
133
        </value>
303
134
        <output/>
304
135
        <stat>
305
136
          <include_in_stat/>
306
137
        </stat>
307
138
        <detectors>
308
139
          <exclude_from_detectors/>
309
140
        </detectors>
310
141
      </prescribed>
311
142
    </vector_field>
312
143
    <scalar_field depends="Velocity" name="MultiplyConnectedStreamFunction" rank="0">
313
144
      <prognostic>
314
145
        <mesh name="VelocityMesh"/>
315
146
        <solver>
316
147
          <iterative_method name="gmres">
317
148
            <restart>
318
149
              <integer_value rank="0">30</integer_value>
319
150
            </restart>
320
151
          </iterative_method>
321
152
          <preconditioner name="sor"/>
322
153
          <relative_error>
323
154
            <real_value rank="0">1e-7</real_value>
324
155
          </relative_error>
325
156
          <max_iterations>
326
157
            <integer_value rank="0">3000</integer_value>
327
158
          </max_iterations>
328
159
          <never_ignore_solver_failures/>
329
160
          <diagnostics>
330
161
            <monitors/>
331
162
          </diagnostics>
332
163
        </solver>
333
164
        <boundary_conditions name="bottom">
334
165
          <surface_ids>
335
166
            <integer_value shape="1" rank="1">7</integer_value>
336
167
          </surface_ids>
337
168
          <primary_boundary/>
338
169
          <type name="dirichlet">
339
170
            <internally_calculated/>
340
171
          </type>
341
172
        </boundary_conditions>
342
173
        <boundary_conditions name="top">
343
174
          <surface_ids>
344
175
            <integer_value shape="1" rank="1">9</integer_value>
345
176
          </surface_ids>
346
177
          <secondary_boundary>
347
178
            <primary_point>
348
179
              <real_value shape="2" dim1="dim" rank="1">0.5 0.0</real_value>
349
180
            </primary_point>
350
181
            <secondary_point>
351
182
              <real_value shape="2" dim1="dim" rank="1">0.5 1.0</real_value>
352
183
            </secondary_point>
353
184
          </secondary_boundary>
354
185
          <type name="dirichlet">
355
186
            <internally_calculated/>
356
187
          </type>
357
188
        </boundary_conditions>
358
189
        <exclude_from_checkpointing/>
359
190
        <output>
360
191
          <include_previous_time_step/>
361
192
        </output>
362
193
        <stat/>
363
194
        <detectors>
364
195
          <include_in_detectors/>
365
196
        </detectors>
366
197
        <steady_state>
367
198
          <include_in_steady_state/>
368
199
        </steady_state>
369
200
        <no_interpolation/>
370
201
      </prognostic>
371
202
    </scalar_field>
372
203
  </material_phase>
373
204
</fluidity_options>
374
205
0
375
=== removed file 'tests/2d_periodic_p1p1_streamfunction/2d_periodic_p1p1_streamfunction.xml'
376
--- tests/2d_periodic_p1p1_streamfunction/2d_periodic_p1p1_streamfunction.xml	2011-10-21 10:05:25 +0000
377
+++ tests/2d_periodic_p1p1_streamfunction/2d_periodic_p1p1_streamfunction.xml	1970-01-01 00:00:00 +0000
378
@@ -1,106 +0,0 @@
379
1
<?xml version="1.0" encoding="UTF-8" ?>
380
2
381
3
<testproblem>
382
4
  <name>Periodic channel with streamfunction diagnostic check get correct boundary condition</name>
383
5
  <owner userid="amcvicar"/>
384
6
  <tags>flml zoltan parallel</tags>
385
7
  <problem_definition length="special" nprocs="1">
386
8
    <command_line> ../../bin/fluidity -v3 -l  2d_periodic_p1p1_streamfunction.flml &amp;&amp;
387
9
    ../../bin/fluidity -v3 -l  2d_nonperiodic_p1p1_streamfunction.flml &amp;&amp; 
388
10
    mpiexec -np 4  ../../bin/fluidity -v3 -l  2d_periodic_p1p1_streamfunction_periodised_flredecomp.flml </command_line>
389
11
  </problem_definition>
390
12
  <variables>
391
13
        <variable name="analytical_solution" language="python">
392
14
# compute analytical solution
393
15
import numpy
394
16
u_bottom = 0.0
395
17
u_top = 1.0
396
18
#divisions=12.0
397
19
#y=numpy.arange(0,(1+(u_top/divisions)),(u_top/divisions))
398
20
#velocity is a function of y
399
21
#x=[]
400
22
#for i in y:
401
23
#  x.append((i**2)/2)
402
24
analytical_solution=float((u_top**2.0)/2.0)
403
25
    </variable>
404
26
405
27
    <variable name="numerical_solution_periodic" language="python">
406
28
import vtktools
407
29
from fluidity_tools import stat_parser as stat
408
30
from numpy import array
409
31
# read in the .vtu
410
32
data=vtktools.vtu("2d_periodic_p1p1_streamfunction_1.vtu")
411
33
SF=data.GetField('MultiplyConnectedStreamFunction');
412
34
numerical_solution_periodic=[0,0,0]
413
35
numerical_solution_periodic[0]=max(abs(array(SF[:,0])));
414
36
numerical_solution_periodic[1]=min((array(SF[:,0])));
415
37
416
38
#final change in streamfunction
417
39
#numerical_solution_periodic[2]=stat("2d_periodic_p1p1_streamfunction.stat")["water"]["MultiplyConnectedStreamFunction"]["l2norm"][-1]
418
40
dt = stat("2d_periodic_p1p1_streamfunction.stat")["dt"]["value"][-1]
419
41
numerical_solution_periodic[2] = max(abs(data.GetScalarField("MultiplyConnectedStreamFunction")-data.GetScalarField("OldMultiplyConnectedStreamFunction"))/dt)
420
42
421
43
    </variable>
422
44
    <variable name="numerical_solution_nonperiodic" language="python">
423
45
import vtktools
424
46
from fluidity_tools import stat_parser as stat
425
47
from numpy import array
426
48
# read in the .vtu
427
49
data=vtktools.vtu("2d_nonperiodic_p1p1_streamfunction_1.vtu")
428
50
SF=data.GetField('MultiplyConnectedStreamFunction');
429
51
numerical_solution_nonperiodic=[0,0,0]
430
52
numerical_solution_nonperiodic[0]=max(abs(array(SF[:,0])));
431
53
numerical_solution_nonperiodic[1]=min((array(SF[:,0])));
432
54
433
55
#final change in streamfunction
434
56
#numerical_solution_nonperiodic[2]=stat("2d_nonperiodic_p1p1_streamfunction.stat")["water"]["MultiplyConnectedStreamFunction"]["l2norm"][-1]
435
57
dt = stat("2d_nonperiodic_p1p1_streamfunction.stat")["dt"]["value"][-1]
436
58
numerical_solution_nonperiodic[2] = max(abs(data.GetScalarField("MultiplyConnectedStreamFunction")-data.GetScalarField("OldMultiplyConnectedStreamFunction"))/dt)
437
59
438
60
    </variable>
439
61
440
62
    <variable name="numerical_solution_periodic_parallel" language="python">
441
63
import vtktools
442
64
from fluidity_tools import stat_parser as stat
443
65
from numpy import array
444
66
# read in the .vtu
445
67
data=vtktools.vtu("2d_periodic_p1p1_streamfunction_parallel_1.pvtu")
446
68
SF=data.GetField('MultiplyConnectedStreamFunction');
447
69
numerical_solution_periodic_parallel=[0,0,0]
448
70
numerical_solution_periodic_parallel[0]=max(abs(array(SF[:,0])));
449
71
numerical_solution_periodic_parallel[1]=min((array(SF[:,0])));
450
72
451
73
#final change in streamfunction
452
74
#numerical_solution_periodic_parallel[2]=stat("2d_periodic_p1p1_streamfunction_parallel.stat")["water"]["MultiplyConnectedStreamFunction"]["l2norm"][-1]
453
75
dt = stat("2d_periodic_p1p1_streamfunction_parallel.stat")["dt"]["value"][-1]
454
76
numerical_solution_periodic_parallel[2] = max(abs(data.GetScalarField("MultiplyConnectedStreamFunction")-data.GetScalarField("OldMultiplyConnectedStreamFunction"))/dt)
455
77
456
78
    </variable>
457
79
458
80
  </variables>
459
81
  <pass_tests>
460
82
    <test name="Compare max error of nonperiodic with analytical solution" language="python">
461
83
    assert abs(numerical_solution_periodic[0]-analytical_solution) &lt; 1e-6</test>
462
84
    <test name="Compare min error of periodic " language="python">
463
85
    assert abs(numerical_solution_periodic[1]) &lt; 1e-6 </test>
464
86
    <test name="Compare L2norm of periodic" language="python">
465
87
    assert abs(numerical_solution_periodic[2]) &lt; 1e-6 </test>
466
88
    
467
89
    <test name="Compare max error of nonperiodic with analytical solution" language="python">
468
90
    assert abs(numerical_solution_nonperiodic[0]-analytical_solution) &lt; 1e-6</test>
469
91
    <test name="Compare min error of nonperiodic " language="python">
470
92
    assert abs(numerical_solution_nonperiodic[1]) &lt; 1e-6 </test>
471
93
    <test name="Compare L2norm of nonperiodic" language="python">
472
94
    assert abs(numerical_solution_nonperiodic[2]) &lt; 1e-6 </test>
473
95
    
474
96
    <test name="Compare max error of periodic parallel with analytical solution" language="python">
475
97
    assert abs(numerical_solution_periodic_parallel[0]-analytical_solution) &lt; 1e-6</test>
476
98
    <test name="Compare min error of periodic parallel" language="python">
477
99
    assert abs(numerical_solution_periodic_parallel[1]) &lt; 1e-6 </test>
478
100
    <test name="Compare L2norm of periodic parallel" language="python">
479
101
    assert abs(numerical_solution_periodic_parallel[2]) &lt; 1e-6 </test>
480
102
    
481
103
  </pass_tests>
482
104
  <warn_tests>
483
105
  </warn_tests>
484
106
</testproblem>
485
107
0
486
=== removed file 'tests/2d_periodic_p1p1_streamfunction/Makefile'
487
--- tests/2d_periodic_p1p1_streamfunction/Makefile	2011-10-21 10:05:25 +0000
488
+++ tests/2d_periodic_p1p1_streamfunction/Makefile	1970-01-01 00:00:00 +0000
489
@@ -1,21 +0,0 @@
490
1
PROJECT = 2d_periodic_p1p1_streamfunction
491
2
492
3
input: clean
493
4
	gmsh -2 -o 2d_square.msh src/2d_square.geo
494
5
	../../bin/gmsh2triangle --2d 2d_square.msh
495
6
	
496
7
	
497
8
	#set up parallel run
498
9
	../../bin/periodise 2d_periodic_p1p1_streamfunction.flml
499
10
	mpiexec -np 4 ../../bin/flredecomp -i 1 -o 4 2d_periodic_p1p1_streamfunction_periodised 2d_periodic_p1p1_streamfunction_periodised_flredecomp
500
11
	spud-set 2d_periodic_p1p1_streamfunction_periodised_flredecomp.flml /simulation_name 2d_periodic_p1p1_streamfunction_parallel
501
12
502
13
	
503
14
clean:
504
15
	rm -f *.stat
505
16
	rm -f *vtu
506
17
	rm -f fluidity* matrixdump*
507
18
	rm -f *.ele *.face *.node *.edge
508
19
	rm -f *.msh
509
20
	rm -f *halo
510
21
	rm -f *periodise*
511
22
\ No newline at end of file
0
\ No newline at end of file
512
23
1
513
=== removed directory 'tests/2d_periodic_p1p1_streamfunction/src'
514
=== removed file 'tests/2d_periodic_p1p1_streamfunction/src/2d_square.geo'
515
--- tests/2d_periodic_p1p1_streamfunction/src/2d_square.geo	2010-06-13 13:10:46 +0000
516
+++ tests/2d_periodic_p1p1_streamfunction/src/2d_square.geo	1970-01-01 00:00:00 +0000
517
@@ -1,20 +0,0 @@
518
1
Point(1) = {0.0,0.0,0,0.08};
519
2
Point(2) = {1.0,0.0,0,0.08};
520
3
Point(3) = {1.0,1.0,0,0.08};
521
4
Point(4) = {0.0,1.0,0,0.08};
522
5
Line(1) = {1,2};
523
6
Line(2) = {2,3};
524
7
Line(3) = {3,4};
525
8
Line(4) = {4,1};
526
9
Line Loop(5) = {2,3,4,1};
527
10
Plane Surface(6) = {5};
528
11
// y const
529
12
Physical Line(7) = {1};
530
13
// x const
531
14
Physical Line(8) = {2};
532
15
// y const
533
16
Physical Line(9) = {3};
534
17
// x const
535
18
Physical Line(10) = {4};
536
19
// z const
537
20
Physical Surface(11) = {6};
538
21
0
539
=== removed directory 'tests/IGW'
540
=== removed file 'tests/IGW/IGW.xml'
541
--- tests/IGW/IGW.xml	2011-10-21 10:05:25 +0000
542
+++ tests/IGW/IGW.xml	1970-01-01 00:00:00 +0000
543
@@ -1,26 +0,0 @@
544
1
<?xml version="1.0" encoding="UTF-8" ?>
545
2
546
3
<testproblem>
547
4
  <name>IGW</name>
548
5
  <owner userid="colin"/>
549
6
  <tags>flml</tags>
550
7
  <problem_definition length="special" nprocs="1">
551
8
    <command_line>../../bin/IGW square1 IGW.log &amp;&gt; IGW.log</command_line>
552
9
  </problem_definition>
553
10
  <variables>
554
11
    <variable name="DeltaH" language="python">
555
12
import vtktools
556
13
start = vtktools.vtu("square1_0.vtu")
557
14
end = vtktools.vtu("square1_1.vtu")
558
15
DeltaH = max(start.GetScalarField("height") - end.GetScalarField("height"))
559
16
    </variable>
560
17
  </variables>
561
18
  <pass_tests>
562
19
    <!-- The height should not change -->
563
20
    <test name="Height difference (1.0e-12)" language="python">
564
21
assert DeltaH &lt; 1e-12
565
22
    </test>
566
23
  </pass_tests>
567
24
  <warn_tests>
568
25
  </warn_tests>
569
26
</testproblem>
570
27
0
571
=== removed file 'tests/IGW/Makefile'
572
--- tests/IGW/Makefile	2011-10-21 10:05:25 +0000
573
+++ tests/IGW/Makefile	1970-01-01 00:00:00 +0000
574
@@ -1,8 +0,0 @@
575
1
default: input
576
2
577
3
input: clean
578
4
	gmsh -2 square.geo > gmsh.log 2> gmsh_err.log
579
5
	../../bin/gmsh2triangle -2 square.msh
580
6
581
7
clean:
582
8
	rm -f *.ele *.edge *.face *.node *.vtu *.s *.stat *.log *.msh
583
9
0
584
=== removed file 'tests/IGW/basin.dat'
585
--- tests/IGW/basin.dat	2009-12-16 14:11:47 +0000
586
+++ tests/IGW/basin.dat	1970-01-01 00:00:00 +0000
587
@@ -1,3 +0,0 @@
588
1
&IGW_data tmax=10.0, tdump=1.0, dt=0.01, fr=1.0, ro=0.1, u_degree=1,
589
2
h_degree=2,u_cty=-1,steady_state=F, mesh="basin",
590
3
u_input="basin_u.py", h_input="basin_h.py", get_streamfunction=F/
591
4
\ No newline at end of file
0
\ No newline at end of file
592
5
1
593
=== removed file 'tests/IGW/basin.geo'
594
--- tests/IGW/basin.geo	2009-12-16 14:11:47 +0000
595
+++ tests/IGW/basin.geo	1970-01-01 00:00:00 +0000
596
@@ -1,9 +0,0 @@
597
1
Point(1) = {2e5,0,0,1e4};
598
2
Point(2) = {0,0,0,1e4};
599
3
Point(3) = {-2e5,0,0,1e4};
600
4
Circle(1) = {1,2,3};
601
5
Circle(2) = {3,2,1};
602
6
Physical Line(3) = {1,2};
603
7
Line Loop(4) = {1,2};
604
8
Plane Surface(5) = {4};
605
9
Physical Surface(6) = {5};
606
10
0
607
=== removed file 'tests/IGW/basin_h.py'
608
--- tests/IGW/basin_h.py	2009-12-16 14:11:47 +0000
609
+++ tests/IGW/basin_h.py	1970-01-01 00:00:00 +0000
610
@@ -1,8 +0,0 @@
611
1
def val(X,t):
612
2
   from math import exp
613
3
   r0 = 0.5
614
4
   r = ((X[0])**2 + (X[1]**2))**0.5
615
5
   if(r>0.0001):
616
6
      return exp((r-r0)/0.1)*X[0]/r
617
7
   else:
618
8
      return 0.0
619
9
0
620
=== removed file 'tests/IGW/basin_u.py'
621
--- tests/IGW/basin_u.py	2009-12-16 14:11:47 +0000
622
+++ tests/IGW/basin_u.py	1970-01-01 00:00:00 +0000
623
@@ -1,9 +0,0 @@
624
1
def val(X,t):
625
2
   from math import exp
626
3
   r0 = 0.5
627
4
   g = [0.0,0.0]
628
5
   r = ((X[0])**2 + (X[1])**2)**0.5
629
6
   if(r>0.0001):
630
7
      g[0] = -exp((r-r0)/0.1)*X[0]/r*X[1]/r
631
8
      g[1] = exp((r-r0)/0.1)*X[0]/r*X[0]/r
632
9
   return g
633
10
\ No newline at end of file
0
\ No newline at end of file
634
11
1
635
=== removed file 'tests/IGW/channel.1_u.py'
636
--- tests/IGW/channel.1_u.py	2009-12-16 14:11:47 +0000
637
+++ tests/IGW/channel.1_u.py	1970-01-01 00:00:00 +0000
638
@@ -1,6 +0,0 @@
639
1
def val(X,t):
640
2
   from math import exp
641
3
   g = [0.0,0.0]
642
4
   g[0] = 0.0
643
5
   g[1] = exp(-X[0]/0.1)*exp(-(X[1]-5.0)**2)
644
6
   return g
645
7
0
646
=== removed file 'tests/IGW/channel1.dat'
647
--- tests/IGW/channel1.dat	2009-12-16 14:11:47 +0000
648
+++ tests/IGW/channel1.dat	1970-01-01 00:00:00 +0000
649
@@ -1,1 +0,0 @@
650
1
&IGW_data tmax=10.0, tdump=1.0, dt=0.01, fr=1.0, ro=0.1 /
651
2
\ No newline at end of file
0
\ No newline at end of file
652
3
1
653
=== removed file 'tests/IGW/channel1_h.py'
654
--- tests/IGW/channel1_h.py	2009-12-16 14:11:47 +0000
655
+++ tests/IGW/channel1_h.py	1970-01-01 00:00:00 +0000
656
@@ -1,3 +0,0 @@
657
1
def val(X,t):
658
2
   from math import exp
659
3
   return -exp(-X[0]/0.1)*exp(-(X[1]-5.0)**2)
660
4
0
661
=== removed file 'tests/IGW/channel1_u.py'
662
--- tests/IGW/channel1_u.py	2009-12-16 14:11:47 +0000
663
+++ tests/IGW/channel1_u.py	1970-01-01 00:00:00 +0000
664
@@ -1,6 +0,0 @@
665
1
def val(X,t):
666
2
   from math import exp
667
3
   g = [0.0,0.0]
668
4
   g[0] = 0.0
669
5
   g[1] = exp(-X[0]/0.1)*exp(-(X[1]-5.0)**2)
670
6
   return g
671
7
\ No newline at end of file
0
\ No newline at end of file
672
8
1
673
=== removed file 'tests/IGW/gauss_h.py'
674
--- tests/IGW/gauss_h.py	2009-12-16 14:11:47 +0000
675
+++ tests/IGW/gauss_h.py	1970-01-01 00:00:00 +0000
676
@@ -1,11 +0,0 @@
677
1
def val(X,t):
678
2
   from math import exp, pi, cos
679
3
   x = X[0]
680
4
   y = X[1]
681
5
   rsq = (X[0]-0.5)**2 + (X[1]-0.5)**2
682
6
   r0sq = 0.3**2
683
7
   g = exp(-rsq/r0sq)
684
8
   if x>0.999 or x<0.001 or y>0.999 or y<0.001:
685
9
      g = 0.0
686
10
   
687
11
   return g
688
12
0
689
=== removed file 'tests/IGW/kelvin.dat'
690
--- tests/IGW/kelvin.dat	2009-12-16 14:11:47 +0000
691
+++ tests/IGW/kelvin.dat	1970-01-01 00:00:00 +0000
692
@@ -1,5 +0,0 @@
693
1
&IGW_data tmax=10.0, tdump=1.0, dt=0.001, fr=1.0, ro=0.1, u_degree=1,
694
2
h_degree=2,u_cty=-1,balanced_u=F, steady_state=F, mesh="kelvin_mod",
695
3
u_input="kelvin_u.py", h_input="kelvin_h.py", 
696
4
exact_h_input="kelvin_h.py", get_streamfunction=F,
697
5
projection_test=F,compute_error=F, exact_u_input="kelvin_u.py" /
698
6
\ No newline at end of file
0
\ No newline at end of file
699
7
1
700
=== removed file 'tests/IGW/kelvin.geo'
701
--- tests/IGW/kelvin.geo	2009-12-16 14:11:47 +0000
702
+++ tests/IGW/kelvin.geo	1970-01-01 00:00:00 +0000
703
@@ -1,12 +0,0 @@
704
1
Point(1) = {-10,0,0,0.01};
705
2
Point(2) = {10,0,0,0.01};
706
3
Point(3) = {-10.0,1.0,0,0.01};
707
4
Point(4) = {10.0,1.0,0,0.01};
708
5
Line(1) = {1,2};
709
6
Line(2) = {2,4};
710
7
Line(4) = {3,1};
711
8
Line(5) = {4,3};
712
9
Line Loop(8) = {1,2,4,5};
713
10
714
11
Plane Surface(9) = {8};
715
12
Physical Surface(10) = {9};
716
13
0
717
=== removed file 'tests/IGW/kelvin_h.py'
718
--- tests/IGW/kelvin_h.py	2009-12-16 14:11:47 +0000
719
+++ tests/IGW/kelvin_h.py	1970-01-01 00:00:00 +0000
720
@@ -1,5 +0,0 @@
721
1
def val(X,t):
722
2
   from math import exp
723
3
   X0 = 5.0 - t
724
4
   sig = 1.0
725
5
   return exp(-X[1]/0.1)*exp(-(X[0]-X0)**2/(sig*sig))
726
6
0
727
=== removed file 'tests/IGW/kelvin_mod.geo'
728
--- tests/IGW/kelvin_mod.geo	2009-12-16 14:11:47 +0000
729
+++ tests/IGW/kelvin_mod.geo	1970-01-01 00:00:00 +0000
730
@@ -1,23 +0,0 @@
731
1
Point(1) = {-15,0,0,2.0};
732
2
Point(2) = {15,0,0,2.0};
733
3
Point(3) = {-15.0,3.0,0,2.0};
734
4
Point(4) = {15.0,3.0,0,2.0};
735
5
Point(5) = {-8,0,0,0.15};
736
6
Point(6) = {8,0,0,0.15};
737
7
Point(7) = {-8,0.6,0,0.15};
738
8
Point(8) = {8,0.6,0,0.15};
739
9
Line(1) = {1,5};
740
10
Line(2) = {6,5};
741
11
Line(3) = {6,2};
742
12
Line(4) = {2,4};
743
13
Line(5) = {4,3};
744
14
Line(6) = {3,1};
745
15
Line(7) = {5,7};
746
16
Line(8) = {7,8};
747
17
Line(9) = {8,6};
748
18
Line Loop(10) = {5,6,1,7,8,9,3,4};
749
19
Line Loop(11) = {8,9,2,7};
750
20
Plane Surface(12) = {11};
751
21
Plane Surface(13) = {10};
752
22
Physical Line(14) = {5,6,1,2,3,4};
753
23
Physical Surface(15) = {13,12};
754
24
0
755
=== removed file 'tests/IGW/kelvin_u.py'
756
--- tests/IGW/kelvin_u.py	2009-12-16 14:11:47 +0000
757
+++ tests/IGW/kelvin_u.py	1970-01-01 00:00:00 +0000
758
@@ -1,7 +0,0 @@
759
1
def val(X,t):
760
2
   from math import exp
761
3
   X0 = 5.0 - t
762
4
   sig = 1.0
763
5
   g = [0.0,0.0]
764
6
   g[0] = exp(-X[1]/0.1)*exp(-(X[0]-X0)**2/(sig*sig))
765
7
   return g
766
8
0
767
=== removed file 'tests/IGW/rand_h.py'
768
--- tests/IGW/rand_h.py	2009-12-16 14:11:47 +0000
769
+++ tests/IGW/rand_h.py	1970-01-01 00:00:00 +0000
770
@@ -1,11 +0,0 @@
771
1
def val(X,t):
772
2
   from random import random
773
3
774
4
   x = X[0]
775
5
   y = X[1]
776
6
   
777
7
   g = random()
778
8
   if x>0.999 or x<0.001 or y>0.999 or y<0.001:
779
9
      g = 0.0
780
10
   
781
11
   return g
782
12
0
783
=== removed file 'tests/IGW/sin_h.py'
784
--- tests/IGW/sin_h.py	2009-12-16 14:11:47 +0000
785
+++ tests/IGW/sin_h.py	1970-01-01 00:00:00 +0000
786
@@ -1,3 +0,0 @@
787
1
def val(X,t):
788
2
   from math import sin, pi
789
3
   return sin(pi*X[0])*sin(pi*X[1])
790
4
0
791
=== removed file 'tests/IGW/square.geo'
792
--- tests/IGW/square.geo	2009-12-16 14:11:47 +0000
793
+++ tests/IGW/square.geo	1970-01-01 00:00:00 +0000
794
@@ -1,12 +0,0 @@
795
1
Point(1) = {0,0,0,0.1};
796
2
Point(2) = {1,0,0,0.1};
797
3
Point(3) = {1,1,0,0.1};
798
4
Point(4) = {0,1,0,0.1};
799
5
Line(1) = {1,2};
800
6
Line(2) = {2,3};
801
7
Line(3) = {3,4};
802
8
Line(4) = {4,1};
803
9
Physical Line(5) = {4,3,2,1};
804
10
Line Loop(6) = {3,4,1,2};
805
11
Plane Surface(7) = {6};
806
12
Physical Surface(8) = {7};
807
13
0
808
=== removed file 'tests/IGW/square1.dat'
809
--- tests/IGW/square1.dat	2009-12-16 14:11:47 +0000
810
+++ tests/IGW/square1.dat	1970-01-01 00:00:00 +0000
811
@@ -1,4 +0,0 @@
812
1
&IGW_data tmax=10.01, tdump=10.0, dt=0.01, fr=1.0, ro=0.1, u_degree=1,
813
2
h_degree=2,u_cty=-1,balanced_u=T,steady_state=T, mesh="square",
814
3
u_input="square1_u.py", h_input="rand_h.py", get_streamfunction=F,
815
4
projection_test=F, compute_error=F/
816
5
\ No newline at end of file
0
\ No newline at end of file
817
6
1
818
=== removed file 'tests/IGW/square1_h.py'
819
--- tests/IGW/square1_h.py	2009-12-16 14:11:47 +0000
820
+++ tests/IGW/square1_h.py	1970-01-01 00:00:00 +0000
821
@@ -1,6 +0,0 @@
822
1
def val(X,t):
823
2
   from random import random
824
3
   if X[0]<0.0001 or X[0]>0.9999 or X[1]<0.0001 or X[1]>0.9999:
825
4
      return 0.0
826
5
   else:
827
6
      return random()
828
7
0
829
=== removed file 'tests/IGW/square1_u.py'
830
--- tests/IGW/square1_u.py	2009-12-16 14:11:47 +0000
831
+++ tests/IGW/square1_u.py	1970-01-01 00:00:00 +0000
832
@@ -1,9 +0,0 @@
833
1
def val(X,t):
834
2
   from math import exp
835
3
   print "Hello, world"
836
4
   x = X[0]
837
5
   y = X[1]
838
6
   g = [None, None]
839
7
   g[0] =  0.555555555556*(y - 0.5)*exp(2.77777777778*(-(y - 0.5)**2 - (x - 0.5)**2))
840
8
   g[1] = -0.555555555556*(x - 0.5)*exp(2.77777777778*(-(y - 0.5)**2 - (x - 0.5)**2))
841
9
   return g
842
10
0
843
=== removed file 'tests/IGW/square_rnd_h.py'
844
--- tests/IGW/square_rnd_h.py	2009-12-16 14:11:47 +0000
845
+++ tests/IGW/square_rnd_h.py	1970-01-01 00:00:00 +0000
846
@@ -1,6 +0,0 @@
847
1
def val(X,t):
848
2
   from random import random
849
3
   if X[0]<0.0001 or X[0]>0.9999 or X[1]<0.0001 or X[1]>0.9999:
850
4
      return 0.0
851
5
   else:
852
6
      return random()
853
7
0
854
=== removed file 'tests/IGW/square_sin.dat'
855
--- tests/IGW/square_sin.dat	2009-12-16 14:11:47 +0000
856
+++ tests/IGW/square_sin.dat	1970-01-01 00:00:00 +0000
857
@@ -1,3 +0,0 @@
858
1
&IGW_data tmax=1.0, tdump=0.1, dt=0.01, fr=1.0, ro=0.1, u_degree=1,
859
2
h_degree=2,u_cty=-1,balanced_u=T,steady_state=T, mesh="square1",
860
3
u_input="square1_u.py", h_input="sin_h.py" /
861
4
\ No newline at end of file
0
\ No newline at end of file
862
5
1
863
=== removed file 'tests/IGW/square_weak.dat'
864
--- tests/IGW/square_weak.dat	2009-12-16 14:11:47 +0000
865
+++ tests/IGW/square_weak.dat	1970-01-01 00:00:00 +0000
866
@@ -1,4 +0,0 @@
867
1
&IGW_data tmax=0.1, tdump=0.01, dt=0.01, fr=1.0, ro=0.1, u_degree=1,
868
2
h_degree=2,u_cty=-1,balanced_u=T,steady_state=T, mesh="square.1",
869
3
u_input="square1_u.py", h_input="square_weak_h.py", get_streamfunction=F,
870
4
projection_test=F, compute_error=F/
871
5
0
872
=== removed file 'tests/IGW/squarep1p1.dat'
873
--- tests/IGW/squarep1p1.dat	2009-12-16 14:11:47 +0000
874
+++ tests/IGW/squarep1p1.dat	1970-01-01 00:00:00 +0000
875
@@ -1,3 +0,0 @@
876
1
&IGW_data tmax=0.02, tdump=0.01, dt=0.01, fr=1.0, ro=0.1, u_degree=1,
877
2
h_degree=2,u_cty=0,balanced_u=T,steady_state=T, mesh="square.1",
878
3
u_input="square1_u.py", h_input="gauss_h.py", get_streamfunction=T/
879
4
\ No newline at end of file
0
\ No newline at end of file
880
5
1
881
=== removed directory 'tests/Stokes_mms_p0p1cv'
882
=== removed file 'tests/Stokes_mms_p0p1cv/MMS_A.flml'
883
--- tests/Stokes_mms_p0p1cv/MMS_A.flml	2010-11-21 15:41:35 +0000
884
+++ tests/Stokes_mms_p0p1cv/MMS_A.flml	1970-01-01 00:00:00 +0000
885
@@ -1,478 +0,0 @@
886
1
<?xml version='1.0' encoding='utf-8'?>
887
2
<fluidity_options>
888
3
  <simulation_name>
889
4
    <string_value lines="1">MMS_A</string_value>
890
5
  </simulation_name>
891
6
  <problem_type>
892
7
    <string_value lines="1">fluids</string_value>
893
8
  </problem_type>
894
9
  <geometry>
895
10
    <dimension>
896
11
      <integer_value rank="0">2</integer_value>
897
12
    </dimension>
898
13
    <mesh name="CoordinateMesh">
899
14
      <from_mesh>
900
15
        <mesh name="InputMesh"/>
901
16
        <stat>
902
17
          <exclude_from_stat/>
903
18
        </stat>
904
19
      </from_mesh>
905
20
    </mesh>
906
21
    <mesh name="VelocityMesh">
907
22
      <from_mesh>
908
23
        <mesh name="InputMesh"/>
909
24
        <mesh_shape>
910
25
          <polynomial_degree>
911
26
            <integer_value rank="0">0</integer_value>
912
27
          </polynomial_degree>
913
28
        </mesh_shape>
914
29
        <mesh_continuity>
915
30
          <string_value>discontinuous</string_value>
916
31
        </mesh_continuity>
917
32
        <stat>
918
33
          <exclude_from_stat/>
919
34
        </stat>
920
35
      </from_mesh>
921
36
    </mesh>
922
37
    <mesh name="PressureMesh">
923
38
      <from_mesh>
924
39
        <mesh name="InputMesh"/>
925
40
        <stat>
926
41
          <exclude_from_stat/>
927
42
        </stat>
928
43
      </from_mesh>
929
44
    </mesh>
930
45
    <mesh name="InputMesh">
931
46
      <from_file file_name="MMS_A">
932
47
        <format name="triangle"/>
933
48
        <stat>
934
49
          <include_in_stat/>
935
50
        </stat>
936
51
      </from_file>
937
52
    </mesh>
938
53
    <quadrature>
939
54
      <degree>
940
55
        <integer_value rank="0">4</integer_value>
941
56
      </degree>
942
57
    </quadrature>
943
58
  </geometry>
944
59
  <io>
945
60
    <dump_format>
946
61
      <string_value>vtk</string_value>
947
62
    </dump_format>
948
63
    <dump_period>
949
64
      <constant>
950
65
        <real_value rank="0">15.0</real_value>
951
66
      </constant>
952
67
    </dump_period>
953
68
    <output_mesh name="VelocityMesh"/>
954
69
    <checkpointing>
955
70
      <checkpoint_period_in_dumps>
956
71
        <integer_value rank="0">10000</integer_value>
957
72
      </checkpoint_period_in_dumps>
958
73
      <checkpoint_at_end/>
959
74
    </checkpointing>
960
75
    <stat/>
961
76
  </io>
962
77
  <timestepping>
963
78
    <current_time>
964
79
      <real_value rank="0">0.0</real_value>
965
80
    </current_time>
966
81
    <timestep>
967
82
      <real_value rank="0">0.0072673325013204136</real_value>
968
83
      <comment>gives a max cfl number of 0.1</comment>
969
84
    </timestep>
970
85
    <finish_time>
971
86
      <real_value rank="0">10.0</real_value>
972
87
      <comment>10.0</comment>
973
88
    </finish_time>
974
89
    <steady_state>
975
90
      <tolerance>
976
91
        <real_value rank="0">1.E-8</real_value>
977
92
        <infinity_norm/>
978
93
      </tolerance>
979
94
    </steady_state>
980
95
  </timestepping>
981
96
  <physical_parameters/>
982
97
  <material_phase name="NS">
983
98
    <equation_of_state>
984
99
      <fluids>
985
100
        <linear>
986
101
          <reference_density>
987
102
            <real_value rank="0">1.0</real_value>
988
103
          </reference_density>
989
104
        </linear>
990
105
      </fluids>
991
106
    </equation_of_state>
992
107
    <scalar_field name="Pressure" rank="0">
993
108
      <prognostic>
994
109
        <mesh name="PressureMesh"/>
995
110
        <spatial_discretisation>
996
111
          <control_volumes/>
997
112
        </spatial_discretisation>
998
113
        <reference_node>
999
114
          <integer_value rank="0">1</integer_value>
1000
115
        </reference_node>
1001
116
        <scheme>
1002
117
          <poisson_pressure_solution>
1003
118
            <string_value lines="1">never</string_value>
1004
119
          </poisson_pressure_solution>
1005
120
          <use_projection_method>
1006
121
            <full_schur_complement>
1007
122
              <inner_matrix name="FullMomentumMatrix">
1008
123
                <solver>
1009
124
                  <iterative_method name="cg"/>
1010
125
                  <preconditioner name="sor"/>
1011
126
                  <relative_error>
1012
127
                    <real_value rank="0">1E-10</real_value>
1013
128
                  </relative_error>
1014
129
                  <max_iterations>
1015
130
                    <integer_value rank="0">10000</integer_value>
1016
131
                  </max_iterations>
1017
132
                  <start_from_zero/>
1018
133
                  <never_ignore_solver_failures/>
1019
134
                  <diagnostics>
1020
135
                    <monitors/>
1021
136
                  </diagnostics>
1022
137
                </solver>
1023
138
              </inner_matrix>
1024
139
              <preconditioner_matrix name="DiagonalSchurComplement"/>
1025
140
            </full_schur_complement>
1026
141
          </use_projection_method>
1027
142
        </scheme>
1028
143
        <solver>
1029
144
          <iterative_method name="fgmres"/>
1030
145
          <preconditioner name="mg"/>
1031
146
          <relative_error>
1032
147
            <real_value rank="0">1.0e-07</real_value>
1033
148
          </relative_error>
1034
149
          <max_iterations>
1035
150
            <integer_value rank="0">10000</integer_value>
1036
151
          </max_iterations>
1037
152
          <never_ignore_solver_failures/>
1038
153
          <diagnostics>
1039
154
            <monitors/>
1040
155
          </diagnostics>
1041
156
        </solver>
1042
157
        <output>
1043
158
          <include_previous_time_step/>
1044
159
        </output>
1045
160
        <stat/>
1046
161
        <convergence>
1047
162
          <include_in_convergence/>
1048
163
        </convergence>
1049
164
        <detectors>
1050
165
          <exclude_from_detectors/>
1051
166
        </detectors>
1052
167
        <steady_state>
1053
168
          <include_in_steady_state/>
1054
169
        </steady_state>
1055
170
        <consistent_interpolation/>
1056
171
      </prognostic>
1057
172
    </scalar_field>
1058
173
    <scalar_field name="Density" rank="0">
1059
174
      <diagnostic>
1060
175
        <algorithm name="Internal" material_phase_support="multiple"/>
1061
176
        <mesh name="VelocityMesh"/>
1062
177
        <output/>
1063
178
        <stat/>
1064
179
        <convergence>
1065
180
          <include_in_convergence/>
1066
181
        </convergence>
1067
182
        <detectors>
1068
183
          <include_in_detectors/>
1069
184
        </detectors>
1070
185
        <steady_state>
1071
186
          <include_in_steady_state/>
1072
187
        </steady_state>
1073
188
      </diagnostic>
1074
189
    </scalar_field>
1075
190
    <vector_field name="Velocity" rank="1">
1076
191
      <prognostic>
1077
192
        <mesh name="VelocityMesh"/>
1078
193
        <equation name="LinearMomentum"/>
1079
194
        <spatial_discretisation>
1080
195
          <discontinuous_galerkin>
1081
196
            <mass_terms>
1082
197
              <exclude_mass_terms/>
1083
198
            </mass_terms>
1084
199
            <viscosity_scheme>
1085
200
              <bassi_rebay/>
1086
201
            </viscosity_scheme>
1087
202
            <advection_scheme>
1088
203
              <none/>
1089
204
              <integrate_advection_by_parts>
1090
205
                <twice/>
1091
206
              </integrate_advection_by_parts>
1092
207
            </advection_scheme>
1093
208
          </discontinuous_galerkin>
1094
209
          <conservative_advection>
1095
210
            <real_value rank="0">0.0</real_value>
1096
211
          </conservative_advection>
1097
212
        </spatial_discretisation>
1098
213
        <temporal_discretisation>
1099
214
          <theta>
1100
215
            <real_value rank="0">1</real_value>
1101
216
          </theta>
1102
217
          <relaxation>
1103
218
            <real_value rank="0">1</real_value>
1104
219
          </relaxation>
1105
220
        </temporal_discretisation>
1106
221
        <solver>
1107
222
          <iterative_method name="cg"/>
1108
223
          <preconditioner name="eisenstat"/>
1109
224
          <relative_error>
1110
225
            <real_value rank="0">1.E-7</real_value>
1111
226
          </relative_error>
1112
227
          <max_iterations>
1113
228
            <integer_value rank="0">10000</integer_value>
1114
229
          </max_iterations>
1115
230
          <never_ignore_solver_failures/>
1116
231
          <diagnostics>
1117
232
            <monitors/>
1118
233
          </diagnostics>
1119
234
        </solver>
1120
235
        <initial_condition name="WholeMesh">
1121
236
          <constant>
1122
237
            <real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
1123
238
          </constant>
1124
239
        </initial_condition>
1125
240
        <boundary_conditions name="all_sides">
1126
241
          <surface_ids>
1127
242
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
1128
243
          </surface_ids>
1129
244
          <type name="dirichlet">
1130
245
            <apply_weakly/>
1131
246
            <align_bc_with_cartesian>
1132
247
              <x_component>
1133
248
                <python>
1134
249
                  <string_value lines="20" type="python">def val(X,t):
1135
250
   from math import sin, cos
1136
251
   # Shear rotation about origin.
1137
252
   u = sin(X[0])*cos(X[1])
1138
253
   return u</string_value>
1139
254
                </python>
1140
255
              </x_component>
1141
256
              <y_component>
1142
257
                <python>
1143
258
                  <string_value lines="20" type="python">def val(X,t):
1144
259
   from math import sin, cos
1145
260
   # Shear rotation about origin.
1146
261
   v = -1.0*cos(X[0])*sin(X[1])
1147
262
   return v</string_value>
1148
263
                </python>
1149
264
              </y_component>
1150
265
            </align_bc_with_cartesian>
1151
266
          </type>
1152
267
        </boundary_conditions>
1153
268
        <tensor_field name="Viscosity" rank="2">
1154
269
          <prescribed>
1155
270
            <value name="WholeMesh">
1156
271
              <isotropic>
1157
272
                <constant>
1158
273
                  <real_value rank="0">0.7</real_value>
1159
274
                </constant>
1160
275
              </isotropic>
1161
276
            </value>
1162
277
            <output/>
1163
278
          </prescribed>
1164
279
        </tensor_field>
1165
280
        <vector_field name="Source" rank="1">
1166
281
          <prescribed>
1167
282
            <value name="WholeMesh">
1168
283
              <python>
1169
284
                <string_value lines="20" type="python">def val(XX, t):
1170
285
   from math import sin,cos
1171
286
   nu = 0.7;
1172
287
   adv = 0.0;
1173
288
   beta = 0.0;
1174
289
   mass = 0.0;
1175
290
   rho = 1.0;
1176
291
   x = XX[0];
1177
292
   y = XX[1];
1178
293
   Su = adv*rho*(cos(x)*sin(x)*sin(y)**2 + cos(x)*sin(x)*cos(y)**2) + 2*nu*sin(x)*cos(y) - sin(x)*cos(y)
1179
294
   Sv = adv*rho*(sin(x)**2*cos(y)*sin(y) + cos(x)**2*cos(y)*sin(y)) - 2*nu*cos(x)*sin(y) - cos(x)*sin(y)
1180
295
   return (Su, Sv)</string_value>
1181
296
                <comment>u(x,y,rho,beta,nu,adv,mass,t) = sin(x)*cos(y)&#13;
1182
297
v(x,y,rho,beta,nu,adv,masst) = -cos(x)*sin(y)
1183
298
p(x,y,rho,beta,nu,adv,mass,t) = cos(x)*cos(y)&#13;
1184
299
1185
300
u_t = diff(u,t)
1186
301
v_t = diff(v,t)
1187
302
1188
303
u_x = diff(u,x)&#13;
1189
304
u_y = diff(u,y)
1190
305
v_x = diff(v,x)
1191
306
v_y = diff(v,y)
1192
307
1193
308
p_x = diff(p,x)
1194
309
p_y = diff(p,y)
1195
310
&#13;
1196
311
u_xx = diff(u_x,x)
1197
312
u_yy = diff(u_y,y)
1198
313
v_xx = diff(v_x,x)
1199
314
v_yy = diff(v_y,y)
1200
315
1201
316
mass*rho*u_t + adv*rho*(u*u_x+v*u_y + beta*(u*u_x + u*v_y)) + p_x - nu*u_xx - nu*u_yy
1202
317
mass*rho*v_t + adv*rho*(u*v_x+v*v_y + beta*(v*u_x + v*v_y)) + p_y - nu*v_xx - nu*v_yy</comment>
1203
318
              </python>
1204
319
            </value>
1205
320
            <output/>
1206
321
            <stat>
1207
322
              <include_in_stat/>
1208
323
            </stat>
1209
324
            <detectors>
1210
325
              <exclude_from_detectors/>
1211
326
            </detectors>
1212
327
          </prescribed>
1213
328
        </vector_field>
1214
329
        <output>
1215
330
          <include_previous_time_step/>
1216
331
        </output>
1217
332
        <stat>
1218
333
          <include_in_stat/>
1219
334
          <previous_time_step>
1220
335
            <exclude_from_stat/>
1221
336
          </previous_time_step>
1222
337
          <nonlinear_field>
1223
338
            <exclude_from_stat/>
1224
339
          </nonlinear_field>
1225
340
        </stat>
1226
341
        <convergence>
1227
342
          <include_in_convergence/>
1228
343
        </convergence>
1229
344
        <detectors>
1230
345
          <include_in_detectors/>
1231
346
        </detectors>
1232
347
        <steady_state>
1233
348
          <include_in_steady_state/>
1234
349
        </steady_state>
1235
350
        <consistent_interpolation/>
1236
351
      </prognostic>
1237
352
    </vector_field>
1238
353
    <scalar_field name="CFLNumber" rank="0">
1239
354
      <diagnostic>
1240
355
        <algorithm name="Internal" material_phase_support="multiple"/>
1241
356
        <mesh name="VelocityMesh"/>
1242
357
        <output/>
1243
358
        <stat/>
1244
359
        <convergence>
1245
360
          <include_in_convergence/>
1246
361
        </convergence>
1247
362
        <detectors>
1248
363
          <include_in_detectors/>
1249
364
        </detectors>
1250
365
        <steady_state>
1251
366
          <include_in_steady_state/>
1252
367
        </steady_state>
1253
368
      </diagnostic>
1254
369
    </scalar_field>
1255
370
    <scalar_field name="PressureAnalytical" rank="0">
1256
371
      <prescribed>
1257
372
        <mesh name="PressureMesh"/>
1258
373
        <value name="WholeMesh">
1259
374
          <python>
1260
375
            <string_value lines="20" type="python">def val(X,t):
1261
376
   from math import cos
1262
377
   p = cos(X[0])*cos(X[1])
1263
378
   return p</string_value>
1264
379
          </python>
1265
380
        </value>
1266
381
        <output/>
1267
382
        <stat/>
1268
383
        <detectors>
1269
384
          <exclude_from_detectors/>
1270
385
        </detectors>
1271
386
      </prescribed>
1272
387
    </scalar_field>
1273
388
    <scalar_field name="ScalarAbsoluteDifference" rank="0">
1274
389
      <diagnostic field_name_b="Pressure" field_name_a="PressureAnalytical">
1275
390
        <algorithm name="Internal" material_phase_support="multiple"/>
1276
391
        <mesh name="PressureMesh"/>
1277
392
        <output/>
1278
393
        <stat/>
1279
394
        <convergence>
1280
395
          <include_in_convergence/>
1281
396
        </convergence>
1282
397
        <detectors>
1283
398
          <include_in_detectors/>
1284
399
        </detectors>
1285
400
        <steady_state>
1286
401
          <include_in_steady_state/>
1287
402
        </steady_state>
1288
403
        <relative_to_average/>
1289
404
      </diagnostic>
1290
405
    </scalar_field>
1291
406
    <scalar_field name="FiniteElementDivergence" rank="0">
1292
407
      <diagnostic field_name="Velocity">
1293
408
        <algorithm legacy="true" name="Internal" material_phase_support="multiple"/>
1294
409
        <mesh name="PressureMesh"/>
1295
410
        <integrate_divergence_by_parts/>
1296
411
        <solver>
1297
412
          <iterative_method name="cg"/>
1298
413
          <preconditioner name="sor"/>
1299
414
          <relative_error>
1300
415
            <real_value rank="0">1.E-10</real_value>
1301
416
          </relative_error>
1302
417
          <max_iterations>
1303
418
            <integer_value rank="0">10000</integer_value>
1304
419
          </max_iterations>
1305
420
          <never_ignore_solver_failures/>
1306
421
          <diagnostics>
1307
422
            <monitors/>
1308
423
          </diagnostics>
1309
424
        </solver>
1310
425
        <output/>
1311
426
        <stat/>
1312
427
        <convergence>
1313
428
          <include_in_convergence/>
1314
429
        </convergence>
1315
430
        <detectors>
1316
431
          <include_in_detectors/>
1317
432
        </detectors>
1318
433
        <steady_state>
1319
434
          <include_in_steady_state/>
1320
435
        </steady_state>
1321
436
      </diagnostic>
1322
437
    </scalar_field>
1323
438
    <vector_field name="VelocityAnalytical" rank="1">
1324
439
      <prescribed>
1325
440
        <mesh name="VelocityMesh"/>
1326
441
        <value name="WholeMesh">
1327
442
          <python>
1328
443
            <string_value lines="20" type="python">def val(X,t):
1329
444
   from math import sin, cos
1330
445
   # Shear rotation about origin.
1331
446
   return (sin(X[0])*cos(X[1]), -1.0*cos(X[0])*sin(X[1]))</string_value>
1332
447
          </python>
1333
448
        </value>
1334
449
        <output/>
1335
450
        <stat>
1336
451
          <include_in_stat/>
1337
452
        </stat>
1338
453
        <detectors>
1339
454
          <exclude_from_detectors/>
1340
455
        </detectors>
1341
456
      </prescribed>
1342
457
    </vector_field>
1343
458
    <vector_field name="VectorAbsoluteDifference" rank="1">
1344
459
      <diagnostic field_name_b="Velocity" field_name_a="VelocityAnalytical">
1345
460
        <algorithm name="Internal" material_phase_support="multiple"/>
1346
461
        <mesh name="VelocityMesh"/>
1347
462
        <output/>
1348
463
        <stat>
1349
464
          <include_in_stat/>
1350
465
        </stat>
1351
466
        <convergence>
1352
467
          <include_in_convergence/>
1353
468
        </convergence>
1354
469
        <detectors>
1355
470
          <include_in_detectors/>
1356
471
        </detectors>
1357
472
        <steady_state>
1358
473
          <include_in_steady_state/>
1359
474
        </steady_state>
1360
475
      </diagnostic>
1361
476
    </vector_field>
1362
477
  </material_phase>
1363
478
</fluidity_options>
1364
479
0
1365
=== removed file 'tests/Stokes_mms_p0p1cv/MMS_B.flml'
1366
--- tests/Stokes_mms_p0p1cv/MMS_B.flml	2010-11-21 15:41:35 +0000
1367
+++ tests/Stokes_mms_p0p1cv/MMS_B.flml	1970-01-01 00:00:00 +0000
1368
@@ -1,478 +0,0 @@
1369
1
<?xml version='1.0' encoding='utf-8'?>
1370
2
<fluidity_options>
1371
3
  <simulation_name>
1372
4
    <string_value lines="1">MMS_B</string_value>
1373
5
  </simulation_name>
1374
6
  <problem_type>
1375
7
    <string_value lines="1">fluids</string_value>
1376
8
  </problem_type>
1377
9
  <geometry>
1378
10
    <dimension>
1379
11
      <integer_value rank="0">2</integer_value>
1380
12
    </dimension>
1381
13
    <mesh name="CoordinateMesh">
1382
14
      <from_mesh>
1383
15
        <mesh name="InputMesh"/>
1384
16
        <stat>
1385
17
          <exclude_from_stat/>
1386
18
        </stat>
1387
19
      </from_mesh>
1388
20
    </mesh>
1389
21
    <mesh name="VelocityMesh">
1390
22
      <from_mesh>
1391
23
        <mesh name="InputMesh"/>
1392
24
        <mesh_shape>
1393
25
          <polynomial_degree>
1394
26
            <integer_value rank="0">0</integer_value>
1395
27
          </polynomial_degree>
1396
28
        </mesh_shape>
1397
29
        <mesh_continuity>
1398
30
          <string_value>discontinuous</string_value>
1399
31
        </mesh_continuity>
1400
32
        <stat>
1401
33
          <exclude_from_stat/>
1402
34
        </stat>
1403
35
      </from_mesh>
1404
36
    </mesh>
1405
37
    <mesh name="PressureMesh">
1406
38
      <from_mesh>
1407
39
        <mesh name="InputMesh"/>
1408
40
        <stat>
1409
41
          <exclude_from_stat/>
1410
42
        </stat>
1411
43
      </from_mesh>
1412
44
    </mesh>
1413
45
    <mesh name="InputMesh">
1414
46
      <from_file file_name="MMS_B">
1415
47
        <format name="triangle"/>
1416
48
        <stat>
1417
49
          <include_in_stat/>
1418
50
        </stat>
1419
51
      </from_file>
1420
52
    </mesh>
1421
53
    <quadrature>
1422
54
      <degree>
1423
55
        <integer_value rank="0">4</integer_value>
1424
56
      </degree>
1425
57
    </quadrature>
1426
58
  </geometry>
1427
59
  <io>
1428
60
    <dump_format>
1429
61
      <string_value>vtk</string_value>
1430
62
    </dump_format>
1431
63
    <dump_period>
1432
64
      <constant>
1433
65
        <real_value rank="0">15.0</real_value>
1434
66
      </constant>
1435
67
    </dump_period>
1436
68
    <output_mesh name="VelocityMesh"/>
1437
69
    <checkpointing>
1438
70
      <checkpoint_period_in_dumps>
1439
71
        <integer_value rank="0">10000</integer_value>
1440
72
      </checkpoint_period_in_dumps>
1441
73
      <checkpoint_at_end/>
1442
74
    </checkpointing>
1443
75
    <stat/>
1444
76
  </io>
1445
77
  <timestepping>
1446
78
    <current_time>
1447
79
      <real_value rank="0">0.0</real_value>
1448
80
    </current_time>
1449
81
    <timestep>
1450
82
      <real_value rank="0">0.0040276830059454609</real_value>
1451
83
      <comment>gives a max cfl number of 0.1</comment>
1452
84
    </timestep>
1453
85
    <finish_time>
1454
86
      <real_value rank="0">10.0</real_value>
1455
87
      <comment>10.0</comment>
1456
88
    </finish_time>
1457
89
    <steady_state>
1458
90
      <tolerance>
1459
91
        <real_value rank="0">1.E-8</real_value>
1460
92
        <infinity_norm/>
1461
93
      </tolerance>
1462
94
    </steady_state>
1463
95
  </timestepping>
1464
96
  <physical_parameters/>
1465
97
  <material_phase name="NS">
1466
98
    <equation_of_state>
1467
99
      <fluids>
1468
100
        <linear>
1469
101
          <reference_density>
1470
102
            <real_value rank="0">1.0</real_value>
1471
103
          </reference_density>
1472
104
        </linear>
1473
105
      </fluids>
1474
106
    </equation_of_state>
1475
107
    <scalar_field name="Pressure" rank="0">
1476
108
      <prognostic>
1477
109
        <mesh name="PressureMesh"/>
1478
110
        <spatial_discretisation>
1479
111
          <control_volumes/>
1480
112
        </spatial_discretisation>
1481
113
        <reference_node>
1482
114
          <integer_value rank="0">1</integer_value>
1483
115
        </reference_node>
1484
116
        <scheme>
1485
117
          <poisson_pressure_solution>
1486
118
            <string_value lines="1">never</string_value>
1487
119
          </poisson_pressure_solution>
1488
120
          <use_projection_method>
1489
121
            <full_schur_complement>
1490
122
              <inner_matrix name="FullMomentumMatrix">
1491
123
                <solver>
1492
124
                  <iterative_method name="cg"/>
1493
125
                  <preconditioner name="sor"/>
1494
126
                  <relative_error>
1495
127
                    <real_value rank="0">1E-10</real_value>
1496
128
                  </relative_error>
1497
129
                  <max_iterations>
1498
130
                    <integer_value rank="0">10000</integer_value>
1499
131
                  </max_iterations>
1500
132
                  <start_from_zero/>
1501
133
                  <never_ignore_solver_failures/>
1502
134
                  <diagnostics>
1503
135
                    <monitors/>
1504
136
                  </diagnostics>
1505
137
                </solver>
1506
138
              </inner_matrix>
1507
139
              <preconditioner_matrix name="DiagonalSchurComplement"/>
1508
140
            </full_schur_complement>
1509
141
          </use_projection_method>
1510
142
        </scheme>
1511
143
        <solver>
1512
144
          <iterative_method name="fgmres"/>
1513
145
          <preconditioner name="mg"/>
1514
146
          <relative_error>
1515
147
            <real_value rank="0">1.0e-07</real_value>
1516
148
          </relative_error>
1517
149
          <max_iterations>
1518
150
            <integer_value rank="0">10000</integer_value>
1519
151
          </max_iterations>
1520
152
          <never_ignore_solver_failures/>
1521
153
          <diagnostics>
1522
154
            <monitors/>
1523
155
          </diagnostics>
1524
156
        </solver>
1525
157
        <output>
1526
158
          <include_previous_time_step/>
1527
159
        </output>
1528
160
        <stat/>
1529
161
        <convergence>
1530
162
          <include_in_convergence/>
1531
163
        </convergence>
1532
164
        <detectors>
1533
165
          <exclude_from_detectors/>
1534
166
        </detectors>
1535
167
        <steady_state>
1536
168
          <include_in_steady_state/>
1537
169
        </steady_state>
1538
170
        <consistent_interpolation/>
1539
171
      </prognostic>
1540
172
    </scalar_field>
1541
173
    <scalar_field name="Density" rank="0">
1542
174
      <diagnostic>
1543
175
        <algorithm name="Internal" material_phase_support="multiple"/>
1544
176
        <mesh name="VelocityMesh"/>
1545
177
        <output/>
1546
178
        <stat/>
1547
179
        <convergence>
1548
180
          <include_in_convergence/>
1549
181
        </convergence>
1550
182
        <detectors>
1551
183
          <include_in_detectors/>
1552
184
        </detectors>
1553
185
        <steady_state>
1554
186
          <include_in_steady_state/>
1555
187
        </steady_state>
1556
188
      </diagnostic>
1557
189
    </scalar_field>
1558
190
    <vector_field name="Velocity" rank="1">
1559
191
      <prognostic>
1560
192
        <mesh name="VelocityMesh"/>
1561
193
        <equation name="LinearMomentum"/>
1562
194
        <spatial_discretisation>
1563
195
          <discontinuous_galerkin>
1564
196
            <mass_terms>
1565
197
              <exclude_mass_terms/>
1566
198
            </mass_terms>
1567
199
            <viscosity_scheme>
1568
200
              <bassi_rebay/>
1569
201
            </viscosity_scheme>
1570
202
            <advection_scheme>
1571
203
              <none/>
1572
204
              <integrate_advection_by_parts>
1573
205
                <twice/>
1574
206
              </integrate_advection_by_parts>
1575
207
            </advection_scheme>
1576
208
          </discontinuous_galerkin>
1577
209
          <conservative_advection>
1578
210
            <real_value rank="0">0.0</real_value>
1579
211
          </conservative_advection>
1580
212
        </spatial_discretisation>
1581
213
        <temporal_discretisation>
1582
214
          <theta>
1583
215
            <real_value rank="0">1</real_value>
1584
216
          </theta>
1585
217
          <relaxation>
1586
218
            <real_value rank="0">1</real_value>
1587
219
          </relaxation>
1588
220
        </temporal_discretisation>
1589
221
        <solver>
1590
222
          <iterative_method name="cg"/>
1591
223
          <preconditioner name="eisenstat"/>
1592
224
          <relative_error>
1593
225
            <real_value rank="0">1.E-07</real_value>
1594
226
          </relative_error>
1595
227
          <max_iterations>
1596
228
            <integer_value rank="0">10000</integer_value>
1597
229
          </max_iterations>
1598
230
          <never_ignore_solver_failures/>
1599
231
          <diagnostics>
1600
232
            <monitors/>
1601
233
          </diagnostics>
1602
234
        </solver>
1603
235
        <initial_condition name="WholeMesh">
1604
236
          <constant>
1605
237
            <real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
1606
238
          </constant>
1607
239
        </initial_condition>
1608
240
        <boundary_conditions name="all_sides">
1609
241
          <surface_ids>
1610
242
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
1611
243
          </surface_ids>
1612
244
          <type name="dirichlet">
1613
245
            <apply_weakly/>
1614
246
            <align_bc_with_cartesian>
1615
247
              <x_component>
1616
248
                <python>
1617
249
                  <string_value lines="20" type="python">def val(X,t):
1618
250
   from math import sin, cos
1619
251
   # Shear rotation about origin.
1620
252
   u = sin(X[0])*cos(X[1])
1621
253
   return u</string_value>
1622
254
                </python>
1623
255
              </x_component>
1624
256
              <y_component>
1625
257
                <python>
1626
258
                  <string_value lines="20" type="python">def val(X,t):
1627
259
   from math import sin, cos
1628
260
   # Shear rotation about origin.
1629
261
   v = -1.0*cos(X[0])*sin(X[1])
1630
262
   return v</string_value>
1631
263
                </python>
1632
264
              </y_component>
1633
265
            </align_bc_with_cartesian>
1634
266
          </type>
1635
267
        </boundary_conditions>
1636
268
        <tensor_field name="Viscosity" rank="2">
1637
269
          <prescribed>
1638
270
            <value name="WholeMesh">
1639
271
              <isotropic>
1640
272
                <constant>
1641
273
                  <real_value rank="0">0.7</real_value>
1642
274
                </constant>
1643
275
              </isotropic>
1644
276
            </value>
1645
277
            <output/>
1646
278
          </prescribed>
1647
279
        </tensor_field>
1648
280
        <vector_field name="Source" rank="1">
1649
281
          <prescribed>
1650
282
            <value name="WholeMesh">
1651
283
              <python>
1652
284
                <string_value lines="20" type="python">def val(XX, t):
1653
285
   from math import sin,cos
1654
286
   nu = 0.7;
1655
287
   adv = 0.0;
1656
288
   beta = 0.0;
1657
289
   mass = 0.0;
1658
290
   rho = 1.0;
1659
291
   x = XX[0];
1660
292
   y = XX[1];
1661
293
   Su = adv*rho*(cos(x)*sin(x)*sin(y)**2 + cos(x)*sin(x)*cos(y)**2) + 2*nu*sin(x)*cos(y) - sin(x)*cos(y)
1662
294
   Sv = adv*rho*(sin(x)**2*cos(y)*sin(y) + cos(x)**2*cos(y)*sin(y)) - 2*nu*cos(x)*sin(y) - cos(x)*sin(y)
1663
295
   return (Su, Sv)</string_value>
1664
296
                <comment>u(x,y,rho,beta,nu,adv,mass,t) = sin(x)*cos(y)&#13;
1665
297
v(x,y,rho,beta,nu,adv,masst) = -cos(x)*sin(y)
1666
298
p(x,y,rho,beta,nu,adv,mass,t) = cos(x)*cos(y)&#13;
1667
299
1668
300
u_t = diff(u,t)
1669
301
v_t = diff(v,t)
1670
302
1671
303
u_x = diff(u,x)&#13;
1672
304
u_y = diff(u,y)
1673
305
v_x = diff(v,x)
1674
306
v_y = diff(v,y)
1675
307
1676
308
p_x = diff(p,x)
1677
309
p_y = diff(p,y)
1678
310
&#13;
1679
311
u_xx = diff(u_x,x)
1680
312
u_yy = diff(u_y,y)
1681
313
v_xx = diff(v_x,x)
1682
314
v_yy = diff(v_y,y)
1683
315
1684
316
mass*rho*u_t + adv*rho*(u*u_x+v*u_y + beta*(u*u_x + u*v_y)) + p_x - nu*u_xx - nu*u_yy
1685
317
mass*rho*v_t + adv*rho*(u*v_x+v*v_y + beta*(v*u_x + v*v_y)) + p_y - nu*v_xx - nu*v_yy</comment>
1686
318
              </python>
1687
319
            </value>
1688
320
            <output/>
1689
321
            <stat>
1690
322
              <include_in_stat/>
1691
323
            </stat>
1692
324
            <detectors>
1693
325
              <exclude_from_detectors/>
1694
326
            </detectors>
1695
327
          </prescribed>
1696
328
        </vector_field>
1697
329
        <output>
1698
330
          <include_previous_time_step/>
1699
331
        </output>
1700
332
        <stat>
1701
333
          <include_in_stat/>
1702
334
          <previous_time_step>
1703
335
            <exclude_from_stat/>
1704
336
          </previous_time_step>
1705
337
          <nonlinear_field>
1706
338
            <exclude_from_stat/>
1707
339
          </nonlinear_field>
1708
340
        </stat>
1709
341
        <convergence>
1710
342
          <include_in_convergence/>
1711
343
        </convergence>
1712
344
        <detectors>
1713
345
          <include_in_detectors/>
1714
346
        </detectors>
1715
347
        <steady_state>
1716
348
          <include_in_steady_state/>
1717
349
        </steady_state>
1718
350
        <consistent_interpolation/>
1719
351
      </prognostic>
1720
352
    </vector_field>
1721
353
    <scalar_field name="CFLNumber" rank="0">
1722
354
      <diagnostic>
1723
355
        <algorithm name="Internal" material_phase_support="multiple"/>
1724
356
        <mesh name="VelocityMesh"/>
1725
357
        <output/>
1726
358
        <stat/>
1727
359
        <convergence>
1728
360
          <include_in_convergence/>
1729
361
        </convergence>
1730
362
        <detectors>
1731
363
          <include_in_detectors/>
1732
364
        </detectors>
1733
365
        <steady_state>
1734
366
          <include_in_steady_state/>
1735
367
        </steady_state>
1736
368
      </diagnostic>
1737
369
    </scalar_field>
1738
370
    <scalar_field name="PressureAnalytical" rank="0">
1739
371
      <prescribed>
1740
372
        <mesh name="PressureMesh"/>
1741
373
        <value name="WholeMesh">
1742
374
          <python>
1743
375
            <string_value lines="20" type="python">def val(X,t):
1744
376
   from math import cos
1745
377
   p = cos(X[0])*cos(X[1])
1746
378
   return p</string_value>
1747
379
          </python>
1748
380
        </value>
1749
381
        <output/>
1750
382
        <stat/>
1751
383
        <detectors>
1752
384
          <exclude_from_detectors/>
1753
385
        </detectors>
1754
386
      </prescribed>
1755
387
    </scalar_field>
1756
388
    <scalar_field name="ScalarAbsoluteDifference" rank="0">
1757
389
      <diagnostic field_name_b="Pressure" field_name_a="PressureAnalytical">
1758
390
        <algorithm name="Internal" material_phase_support="multiple"/>
1759
391
        <mesh name="PressureMesh"/>
1760
392
        <output/>
1761
393
        <stat/>
1762
394
        <convergence>
1763
395
          <include_in_convergence/>
1764
396
        </convergence>
1765
397
        <detectors>
1766
398
          <include_in_detectors/>
1767
399
        </detectors>
1768
400
        <steady_state>
1769
401
          <include_in_steady_state/>
1770
402
        </steady_state>
1771
403
        <relative_to_average/>
1772
404
      </diagnostic>
1773
405
    </scalar_field>
1774
406
    <scalar_field name="FiniteElementDivergence" rank="0">
1775
407
      <diagnostic field_name="Velocity">
1776
408
        <algorithm legacy="true" name="Internal" material_phase_support="multiple"/>
1777
409
        <mesh name="PressureMesh"/>
1778
410
        <integrate_divergence_by_parts/>
1779
411
        <solver>
1780
412
          <iterative_method name="cg"/>
1781
413
          <preconditioner name="sor"/>
1782
414
          <relative_error>
1783
415
            <real_value rank="0">1.E-10</real_value>
1784
416
          </relative_error>
1785
417
          <max_iterations>
1786
418
            <integer_value rank="0">10000</integer_value>
1787
419
          </max_iterations>
1788
420
          <never_ignore_solver_failures/>
1789
421
          <diagnostics>
1790
422
            <monitors/>
1791
423
          </diagnostics>
1792
424
        </solver>
1793
425
        <output/>
1794
426
        <stat/>
1795
427
        <convergence>
1796
428
          <include_in_convergence/>
1797
429
        </convergence>
1798
430
        <detectors>
1799
431
          <include_in_detectors/>
1800
432
        </detectors>
1801
433
        <steady_state>
1802
434
          <include_in_steady_state/>
1803
435
        </steady_state>
1804
436
      </diagnostic>
1805
437
    </scalar_field>
1806
438
    <vector_field name="VelocityAnalytical" rank="1">
1807
439
      <prescribed>
1808
440
        <mesh name="VelocityMesh"/>
1809
441
        <value name="WholeMesh">
1810
442
          <python>
1811
443
            <string_value lines="20" type="python">def val(X,t):
1812
444
   from math import sin, cos
1813
445
   # Shear rotation about origin.
1814
446
   return (sin(X[0])*cos(X[1]), -1.0*cos(X[0])*sin(X[1]))</string_value>
1815
447
          </python>
1816
448
        </value>
1817
449
        <output/>
1818
450
        <stat>
1819
451
          <include_in_stat/>
1820
452
        </stat>
1821
453
        <detectors>
1822
454
          <exclude_from_detectors/>
1823
455
        </detectors>
1824
456
      </prescribed>
1825
457
    </vector_field>
1826
458
    <vector_field name="VectorAbsoluteDifference" rank="1">
1827
459
      <diagnostic field_name_b="Velocity" field_name_a="VelocityAnalytical">
1828
460
        <algorithm name="Internal" material_phase_support="multiple"/>
1829
461
        <mesh name="VelocityMesh"/>
1830
462
        <output/>
1831
463
        <stat>
1832
464
          <include_in_stat/>
1833
465
        </stat>
1834
466
        <convergence>
1835
467
          <include_in_convergence/>
1836
468
        </convergence>
1837
469
        <detectors>
1838
470
          <include_in_detectors/>
1839
471
        </detectors>
1840
472
        <steady_state>
1841
473
          <include_in_steady_state/>
1842
474
        </steady_state>
1843
475
      </diagnostic>
1844
476
    </vector_field>
1845
477
  </material_phase>
1846
478
</fluidity_options>
1847
479
0
1848
=== removed file 'tests/Stokes_mms_p0p1cv/MMS_C.flml'
1849
--- tests/Stokes_mms_p0p1cv/MMS_C.flml	2010-11-21 15:41:35 +0000
1850
+++ tests/Stokes_mms_p0p1cv/MMS_C.flml	1970-01-01 00:00:00 +0000
1851
@@ -1,478 +0,0 @@
1852
1
<?xml version='1.0' encoding='utf-8'?>
1853
2
<fluidity_options>
1854
3
  <simulation_name>
1855
4
    <string_value lines="1">MMS_C</string_value>
1856
5
  </simulation_name>
1857
6
  <problem_type>
1858
7
    <string_value lines="1">fluids</string_value>
1859
8
  </problem_type>
1860
9
  <geometry>
1861
10
    <dimension>
1862
11
      <integer_value rank="0">2</integer_value>
1863
12
    </dimension>
1864
13
    <mesh name="CoordinateMesh">
1865
14
      <from_mesh>
1866
15
        <mesh name="InputMesh"/>
1867
16
        <stat>
1868
17
          <exclude_from_stat/>
1869
18
        </stat>
1870
19
      </from_mesh>
1871
20
    </mesh>
1872
21
    <mesh name="VelocityMesh">
1873
22
      <from_mesh>
1874
23
        <mesh name="InputMesh"/>
1875
24
        <mesh_shape>
1876
25
          <polynomial_degree>
1877
26
            <integer_value rank="0">0</integer_value>
1878
27
          </polynomial_degree>
1879
28
        </mesh_shape>
1880
29
        <mesh_continuity>
1881
30
          <string_value>discontinuous</string_value>
1882
31
        </mesh_continuity>
1883
32
        <stat>
1884
33
          <exclude_from_stat/>
1885
34
        </stat>
1886
35
      </from_mesh>
1887
36
    </mesh>
1888
37
    <mesh name="PressureMesh">
1889
38
      <from_mesh>
1890
39
        <mesh name="InputMesh"/>
1891
40
        <stat>
1892
41
          <exclude_from_stat/>
1893
42
        </stat>
1894
43
      </from_mesh>
1895
44
    </mesh>
1896
45
    <mesh name="InputMesh">
1897
46
      <from_file file_name="MMS_C">
1898
47
        <format name="triangle"/>
1899
48
        <stat>
1900
49
          <include_in_stat/>
1901
50
        </stat>
1902
51
      </from_file>
1903
52
    </mesh>
1904
53
    <quadrature>
1905
54
      <degree>
1906
55
        <integer_value rank="0">4</integer_value>
1907
56
      </degree>
1908
57
    </quadrature>
1909
58
  </geometry>
1910
59
  <io>
1911
60
    <dump_format>
1912
61
      <string_value>vtk</string_value>
1913
62
    </dump_format>
1914
63
    <dump_period>
1915
64
      <constant>
1916
65
        <real_value rank="0">15.0</real_value>
1917
66
      </constant>
1918
67
    </dump_period>
1919
68
    <output_mesh name="VelocityMesh"/>
1920
69
    <checkpointing>
1921
70
      <checkpoint_period_in_dumps>
1922
71
        <integer_value rank="0">10000</integer_value>
1923
72
      </checkpoint_period_in_dumps>
1924
73
      <checkpoint_at_end/>
1925
74
    </checkpointing>
1926
75
    <stat/>
1927
76
  </io>
1928
77
  <timestepping>
1929
78
    <current_time>
1930
79
      <real_value rank="0">0.0</real_value>
1931
80
    </current_time>
1932
81
    <timestep>
1933
82
      <real_value rank="0">0.0020011146041869813</real_value>
1934
83
      <comment>gives a max cfl number of 0.1</comment>
1935
84
    </timestep>
1936
85
    <finish_time>
1937
86
      <real_value rank="0">10.0</real_value>
1938
87
      <comment>10.0</comment>
1939
88
    </finish_time>
1940
89
    <steady_state>
1941
90
      <tolerance>
1942
91
        <real_value rank="0">1.E-8</real_value>
1943
92
        <infinity_norm/>
1944
93
      </tolerance>
1945
94
    </steady_state>
1946
95
  </timestepping>
1947
96
  <physical_parameters/>
1948
97
  <material_phase name="NS">
1949
98
    <equation_of_state>
1950
99
      <fluids>
1951
100
        <linear>
1952
101
          <reference_density>
1953
102
            <real_value rank="0">1.0</real_value>
1954
103
          </reference_density>
1955
104
        </linear>
1956
105
      </fluids>
1957
106
    </equation_of_state>
1958
107
    <scalar_field name="Pressure" rank="0">
1959
108
      <prognostic>
1960
109
        <mesh name="PressureMesh"/>
1961
110
        <spatial_discretisation>
1962
111
          <control_volumes/>
1963
112
        </spatial_discretisation>
1964
113
        <reference_node>
1965
114
          <integer_value rank="0">1</integer_value>
1966
115
        </reference_node>
1967
116
        <scheme>
1968
117
          <poisson_pressure_solution>
1969
118
            <string_value lines="1">never</string_value>
1970
119
          </poisson_pressure_solution>
1971
120
          <use_projection_method>
1972
121
            <full_schur_complement>
1973
122
              <inner_matrix name="FullMomentumMatrix">
1974
123
                <solver>
1975
124
                  <iterative_method name="cg"/>
1976
125
                  <preconditioner name="sor"/>
1977
126
                  <relative_error>
1978
127
                    <real_value rank="0">1E-10</real_value>
1979
128
                  </relative_error>
1980
129
                  <max_iterations>
1981
130
                    <integer_value rank="0">10000</integer_value>
1982
131
                  </max_iterations>
1983
132
                  <start_from_zero/>
1984
133
                  <never_ignore_solver_failures/>
1985
134
                  <diagnostics>
1986
135
                    <monitors/>
1987
136
                  </diagnostics>
1988
137
                </solver>
1989
138
              </inner_matrix>
1990
139
              <preconditioner_matrix name="DiagonalSchurComplement"/>
1991
140
            </full_schur_complement>
1992
141
          </use_projection_method>
1993
142
        </scheme>
1994
143
        <solver>
1995
144
          <iterative_method name="fgmres"/>
1996
145
          <preconditioner name="mg"/>
1997
146
          <relative_error>
1998
147
            <real_value rank="0">1.0e-07</real_value>
1999
148
          </relative_error>
2000
149
          <max_iterations>
2001
150
            <integer_value rank="0">10000</integer_value>
2002
151
          </max_iterations>
2003
152
          <never_ignore_solver_failures/>
2004
153
          <diagnostics>
2005
154
            <monitors/>
2006
155
          </diagnostics>
2007
156
        </solver>
2008
157
        <output>
2009
158
          <include_previous_time_step/>
2010
159
        </output>
2011
160
        <stat/>
2012
161
        <convergence>
2013
162
          <include_in_convergence/>
2014
163
        </convergence>
2015
164
        <detectors>
2016
165
          <exclude_from_detectors/>
2017
166
        </detectors>
2018
167
        <steady_state>
2019
168
          <include_in_steady_state/>
2020
169
        </steady_state>
2021
170
        <consistent_interpolation/>
2022
171
      </prognostic>
2023
172
    </scalar_field>
2024
173
    <scalar_field name="Density" rank="0">
2025
174
      <diagnostic>
2026
175
        <algorithm name="Internal" material_phase_support="multiple"/>
2027
176
        <mesh name="VelocityMesh"/>
2028
177
        <output/>
2029
178
        <stat/>
2030
179
        <convergence>
2031
180
          <include_in_convergence/>
2032
181
        </convergence>
2033
182
        <detectors>
2034
183
          <include_in_detectors/>
2035
184
        </detectors>
2036
185
        <steady_state>
2037
186
          <include_in_steady_state/>
2038
187
        </steady_state>
2039
188
      </diagnostic>
2040
189
    </scalar_field>
2041
190
    <vector_field name="Velocity" rank="1">
2042
191
      <prognostic>
2043
192
        <mesh name="VelocityMesh"/>
2044
193
        <equation name="LinearMomentum"/>
2045
194
        <spatial_discretisation>
2046
195
          <discontinuous_galerkin>
2047
196
            <mass_terms>
2048
197
              <exclude_mass_terms/>
2049
198
            </mass_terms>
2050
199
            <viscosity_scheme>
2051
200
              <bassi_rebay/>
2052
201
            </viscosity_scheme>
2053
202
            <advection_scheme>
2054
203
              <none/>
2055
204
              <integrate_advection_by_parts>
2056
205
                <twice/>
2057
206
              </integrate_advection_by_parts>
2058
207
            </advection_scheme>
2059
208
          </discontinuous_galerkin>
2060
209
          <conservative_advection>
2061
210
            <real_value rank="0">0.0</real_value>
2062
211
          </conservative_advection>
2063
212
        </spatial_discretisation>
2064
213
        <temporal_discretisation>
2065
214
          <theta>
2066
215
            <real_value rank="0">1</real_value>
2067
216
          </theta>
2068
217
          <relaxation>
2069
218
            <real_value rank="0">1</real_value>
2070
219
          </relaxation>
2071
220
        </temporal_discretisation>
2072
221
        <solver>
2073
222
          <iterative_method name="cg"/>
2074
223
          <preconditioner name="eisenstat"/>
2075
224
          <relative_error>
2076
225
            <real_value rank="0">1.E-07</real_value>
2077
226
          </relative_error>
2078
227
          <max_iterations>
2079
228
            <integer_value rank="0">10000</integer_value>
2080
229
          </max_iterations>
2081
230
          <never_ignore_solver_failures/>
2082
231
          <diagnostics>
2083
232
            <monitors/>
2084
233
          </diagnostics>
2085
234
        </solver>
2086
235
        <initial_condition name="WholeMesh">
2087
236
          <constant>
2088
237
            <real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
2089
238
          </constant>
2090
239
        </initial_condition>
2091
240
        <boundary_conditions name="all_sides">
2092
241
          <surface_ids>
2093
242
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
2094
243
          </surface_ids>
2095
244
          <type name="dirichlet">
2096
245
            <apply_weakly/>
2097
246
            <align_bc_with_cartesian>
2098
247
              <x_component>
2099
248
                <python>
2100
249
                  <string_value lines="20" type="python">def val(X,t):
2101
250
   from math import sin, cos
2102
251
   # Shear rotation about origin.
2103
252
   u = sin(X[0])*cos(X[1])
2104
253
   return u</string_value>
2105
254
                </python>
2106
255
              </x_component>
2107
256
              <y_component>
2108
257
                <python>
2109
258
                  <string_value lines="20" type="python">def val(X,t):
2110
259
   from math import sin, cos
2111
260
   # Shear rotation about origin.
2112
261
   v = -1.0*cos(X[0])*sin(X[1])
2113
262
   return v</string_value>
2114
263
                </python>
2115
264
              </y_component>
2116
265
            </align_bc_with_cartesian>
2117
266
          </type>
2118
267
        </boundary_conditions>
2119
268
        <tensor_field name="Viscosity" rank="2">
2120
269
          <prescribed>
2121
270
            <value name="WholeMesh">
2122
271
              <isotropic>
2123
272
                <constant>
2124
273
                  <real_value rank="0">0.7</real_value>
2125
274
                </constant>
2126
275
              </isotropic>
2127
276
            </value>
2128
277
            <output/>
2129
278
          </prescribed>
2130
279
        </tensor_field>
2131
280
        <vector_field name="Source" rank="1">
2132
281
          <prescribed>
2133
282
            <value name="WholeMesh">
2134
283
              <python>
2135
284
                <string_value lines="20" type="python">def val(XX, t):
2136
285
   from math import sin,cos
2137
286
   nu = 0.7;
2138
287
   adv = 0.0;
2139
288
   beta = 0.0;
2140
289
   mass = 0.0;
2141
290
   rho = 1.0;
2142
291
   x = XX[0];
2143
292
   y = XX[1];
2144
293
   Su = adv*rho*(cos(x)*sin(x)*sin(y)**2 + cos(x)*sin(x)*cos(y)**2) + 2*nu*sin(x)*cos(y) - sin(x)*cos(y)
2145
294
   Sv = adv*rho*(sin(x)**2*cos(y)*sin(y) + cos(x)**2*cos(y)*sin(y)) - 2*nu*cos(x)*sin(y) - cos(x)*sin(y)
2146
295
   return (Su, Sv)</string_value>
2147
296
                <comment>u(x,y,rho,beta,nu,adv,mass,t) = sin(x)*cos(y)&#13;
2148
297
v(x,y,rho,beta,nu,adv,masst) = -cos(x)*sin(y)
2149
298
p(x,y,rho,beta,nu,adv,mass,t) = cos(x)*cos(y)&#13;
2150
299
2151
300
u_t = diff(u,t)
2152
301
v_t = diff(v,t)
2153
302
2154
303
u_x = diff(u,x)&#13;
2155
304
u_y = diff(u,y)
2156
305
v_x = diff(v,x)
2157
306
v_y = diff(v,y)
2158
307
2159
308
p_x = diff(p,x)
2160
309
p_y = diff(p,y)
2161
310
&#13;
2162
311
u_xx = diff(u_x,x)
2163
312
u_yy = diff(u_y,y)
2164
313
v_xx = diff(v_x,x)
2165
314
v_yy = diff(v_y,y)
2166
315
2167
316
mass*rho*u_t + adv*rho*(u*u_x+v*u_y + beta*(u*u_x + u*v_y)) + p_x - nu*u_xx - nu*u_yy
2168
317
mass*rho*v_t + adv*rho*(u*v_x+v*v_y + beta*(v*u_x + v*v_y)) + p_y - nu*v_xx - nu*v_yy</comment>
2169
318
              </python>
2170
319
            </value>
2171
320
            <output/>
2172
321
            <stat>
2173
322
              <include_in_stat/>
2174
323
            </stat>
2175
324
            <detectors>
2176
325
              <exclude_from_detectors/>
2177
326
            </detectors>
2178
327
          </prescribed>
2179
328
        </vector_field>
2180
329
        <output>
2181
330
          <include_previous_time_step/>
2182
331
        </output>
2183
332
        <stat>
2184
333
          <include_in_stat/>
2185
334
          <previous_time_step>
2186
335
            <exclude_from_stat/>
2187
336
          </previous_time_step>
2188
337
          <nonlinear_field>
2189
338
            <exclude_from_stat/>
2190
339
          </nonlinear_field>
2191
340
        </stat>
2192
341
        <convergence>
2193
342
          <include_in_convergence/>
2194
343
        </convergence>
2195
344
        <detectors>
2196
345
          <include_in_detectors/>
2197
346
        </detectors>
2198
347
        <steady_state>
2199
348
          <include_in_steady_state/>
2200
349
        </steady_state>
2201
350
        <consistent_interpolation/>
2202
351
      </prognostic>
2203
352
    </vector_field>
2204
353
    <scalar_field name="CFLNumber" rank="0">
2205
354
      <diagnostic>
2206
355
        <algorithm name="Internal" material_phase_support="multiple"/>
2207
356
        <mesh name="VelocityMesh"/>
2208
357
        <output/>
2209
358
        <stat/>
2210
359
        <convergence>
2211
360
          <include_in_convergence/>
2212
361
        </convergence>
2213
362
        <detectors>
2214
363
          <include_in_detectors/>
2215
364
        </detectors>
2216
365
        <steady_state>
2217
366
          <include_in_steady_state/>
2218
367
        </steady_state>
2219
368
      </diagnostic>
2220
369
    </scalar_field>
2221
370
    <scalar_field name="PressureAnalytical" rank="0">
2222
371
      <prescribed>
2223
372
        <mesh name="PressureMesh"/>
2224
373
        <value name="WholeMesh">
2225
374
          <python>
2226
375
            <string_value lines="20" type="python">def val(X,t):
2227
376
   from math import cos
2228
377
   p = cos(X[0])*cos(X[1])
2229
378
   return p</string_value>
2230
379
          </python>
2231
380
        </value>
2232
381
        <output/>
2233
382
        <stat/>
2234
383
        <detectors>
2235
384
          <exclude_from_detectors/>
2236
385
        </detectors>
2237
386
      </prescribed>
2238
387
    </scalar_field>
2239
388
    <scalar_field name="ScalarAbsoluteDifference" rank="0">
2240
389
      <diagnostic field_name_b="Pressure" field_name_a="PressureAnalytical">
2241
390
        <algorithm name="Internal" material_phase_support="multiple"/>
2242
391
        <mesh name="PressureMesh"/>
2243
392
        <output/>
2244
393
        <stat/>
2245
394
        <convergence>
2246
395
          <include_in_convergence/>
2247
396
        </convergence>
2248
397
        <detectors>
2249
398
          <include_in_detectors/>
2250
399
        </detectors>
2251
400
        <steady_state>
2252
401
          <include_in_steady_state/>
2253
402
        </steady_state>
2254
403
        <relative_to_average/>
2255
404
      </diagnostic>
2256
405
    </scalar_field>
2257
406
    <scalar_field name="FiniteElementDivergence" rank="0">
2258
407
      <diagnostic field_name="Velocity">
2259
408
        <algorithm legacy="true" name="Internal" material_phase_support="multiple"/>
2260
409
        <mesh name="PressureMesh"/>
2261
410
        <integrate_divergence_by_parts/>
2262
411
        <solver>
2263
412
          <iterative_method name="cg"/>
2264
413
          <preconditioner name="sor"/>
2265
414
          <relative_error>
2266
415
            <real_value rank="0">1.E-10</real_value>
2267
416
          </relative_error>
2268
417
          <max_iterations>
2269
418
            <integer_value rank="0">10000</integer_value>
2270
419
          </max_iterations>
2271
420
          <never_ignore_solver_failures/>
2272
421
          <diagnostics>
2273
422
            <monitors/>
2274
423
          </diagnostics>
2275
424
        </solver>
2276
425
        <output/>
2277
426
        <stat/>
2278
427
        <convergence>
2279
428
          <include_in_convergence/>
2280
429
        </convergence>
2281
430
        <detectors>
2282
431
          <include_in_detectors/>
2283
432
        </detectors>
2284
433
        <steady_state>
2285
434
          <include_in_steady_state/>
2286
435
        </steady_state>
2287
436
      </diagnostic>
2288
437
    </scalar_field>
2289
438
    <vector_field name="VelocityAnalytical" rank="1">
2290
439
      <prescribed>
2291
440
        <mesh name="VelocityMesh"/>
2292
441
        <value name="WholeMesh">
2293
442
          <python>
2294
443
            <string_value lines="20" type="python">def val(X,t):
2295
444
   from math import sin, cos
2296
445
   # Shear rotation about origin.
2297
446
   return (sin(X[0])*cos(X[1]), -1.0*cos(X[0])*sin(X[1]))</string_value>
2298
447
          </python>
2299
448
        </value>
2300
449
        <output/>
2301
450
        <stat>
2302
451
          <include_in_stat/>
2303
452
        </stat>
2304
453
        <detectors>
2305
454
          <exclude_from_detectors/>
2306
455
        </detectors>
2307
456
      </prescribed>
2308
457
    </vector_field>
2309
458
    <vector_field name="VectorAbsoluteDifference" rank="1">
2310
459
      <diagnostic field_name_b="Velocity" field_name_a="VelocityAnalytical">
2311
460
        <algorithm name="Internal" material_phase_support="multiple"/>
2312
461
        <mesh name="VelocityMesh"/>
2313
462
        <output/>
2314
463
        <stat>
2315
464
          <include_in_stat/>
2316
465
        </stat>
2317
466
        <convergence>
2318
467
          <include_in_convergence/>
2319
468
        </convergence>
2320
469
        <detectors>
2321
470
          <include_in_detectors/>
2322
471
        </detectors>
2323
472
        <steady_state>
2324
473
          <include_in_steady_state/>
2325
474
        </steady_state>
2326
475
      </diagnostic>
2327
476
    </vector_field>
2328
477
  </material_phase>
2329
478
</fluidity_options>
2330
479
0
2331
=== removed file 'tests/Stokes_mms_p0p1cv/MMS_D.flml'
2332
--- tests/Stokes_mms_p0p1cv/MMS_D.flml	2010-11-21 15:41:35 +0000
2333
+++ tests/Stokes_mms_p0p1cv/MMS_D.flml	1970-01-01 00:00:00 +0000
2334
@@ -1,478 +0,0 @@
2335
1
<?xml version='1.0' encoding='utf-8'?>
2336
2
<fluidity_options>
2337
3
  <simulation_name>
2338
4
    <string_value lines="1">MMS_D</string_value>
2339
5
  </simulation_name>
2340
6
  <problem_type>
2341
7
    <string_value lines="1">fluids</string_value>
2342
8
  </problem_type>
2343
9
  <geometry>
2344
10
    <dimension>
2345
11
      <integer_value rank="0">2</integer_value>
2346
12
    </dimension>
2347
13
    <mesh name="CoordinateMesh">
2348
14
      <from_mesh>
2349
15
        <mesh name="InputMesh"/>
2350
16
        <stat>
2351
17
          <exclude_from_stat/>
2352
18
        </stat>
2353
19
      </from_mesh>
2354
20
    </mesh>
2355
21
    <mesh name="VelocityMesh">
2356
22
      <from_mesh>
2357
23
        <mesh name="InputMesh"/>
2358
24
        <mesh_shape>
2359
25
          <polynomial_degree>
2360
26
            <integer_value rank="0">0</integer_value>
2361
27
          </polynomial_degree>
2362
28
        </mesh_shape>
2363
29
        <mesh_continuity>
2364
30
          <string_value>discontinuous</string_value>
2365
31
        </mesh_continuity>
2366
32
        <stat>
2367
33
          <exclude_from_stat/>
2368
34
        </stat>
2369
35
      </from_mesh>
2370
36
    </mesh>
2371
37
    <mesh name="PressureMesh">
2372
38
      <from_mesh>
2373
39
        <mesh name="InputMesh"/>
2374
40
        <stat>
2375
41
          <exclude_from_stat/>
2376
42
        </stat>
2377
43
      </from_mesh>
2378
44
    </mesh>
2379
45
    <mesh name="InputMesh">
2380
46
      <from_file file_name="MMS_D">
2381
47
        <format name="triangle"/>
2382
48
        <stat>
2383
49
          <include_in_stat/>
2384
50
        </stat>
2385
51
      </from_file>
2386
52
    </mesh>
2387
53
    <quadrature>
2388
54
      <degree>
2389
55
        <integer_value rank="0">4</integer_value>
2390
56
      </degree>
2391
57
    </quadrature>
2392
58
  </geometry>
2393
59
  <io>
2394
60
    <dump_format>
2395
61
      <string_value>vtk</string_value>
2396
62
    </dump_format>
2397
63
    <dump_period>
2398
64
      <constant>
2399
65
        <real_value rank="0">15.0</real_value>
2400
66
      </constant>
2401
67
    </dump_period>
2402
68
    <output_mesh name="VelocityMesh"/>
2403
69
    <checkpointing>
2404
70
      <checkpoint_period_in_dumps>
2405
71
        <integer_value rank="0">10000</integer_value>
2406
72
      </checkpoint_period_in_dumps>
2407
73
      <checkpoint_at_end/>
2408
74
    </checkpointing>
2409
75
    <stat/>
2410
76
  </io>
2411
77
  <timestepping>
2412
78
    <current_time>
2413
79
      <real_value rank="0">0.0</real_value>
2414
80
    </current_time>
2415
81
    <timestep>
2416
82
      <real_value rank="0">0.00073714354843638517</real_value>
2417
83
      <comment>gives a max cfl number of 0.1</comment>
2418
84
    </timestep>
2419
85
    <finish_time>
2420
86
      <real_value rank="0">10.0</real_value>
2421
87
      <comment>10.0</comment>
2422
88
    </finish_time>
2423
89
    <steady_state>
2424
90
      <tolerance>
2425
91
        <real_value rank="0">1.E-8</real_value>
2426
92
        <infinity_norm/>
2427
93
      </tolerance>
2428
94
    </steady_state>
2429
95
  </timestepping>
2430
96
  <physical_parameters/>
2431
97
  <material_phase name="NS">
2432
98
    <equation_of_state>
2433
99
      <fluids>
2434
100
        <linear>
2435
101
          <reference_density>
2436
102
            <real_value rank="0">1.0</real_value>
2437
103
          </reference_density>
2438
104
        </linear>
2439
105
      </fluids>
2440
106
    </equation_of_state>
2441
107
    <scalar_field name="Pressure" rank="0">
2442
108
      <prognostic>
2443
109
        <mesh name="PressureMesh"/>
2444
110
        <spatial_discretisation>
2445
111
          <control_volumes/>
2446
112
        </spatial_discretisation>
2447
113
        <reference_node>
2448
114
          <integer_value rank="0">1</integer_value>
2449
115
        </reference_node>
2450
116
        <scheme>
2451
117
          <poisson_pressure_solution>
2452
118
            <string_value lines="1">never</string_value>
2453
119
          </poisson_pressure_solution>
2454
120
          <use_projection_method>
2455
121
            <full_schur_complement>
2456
122
              <inner_matrix name="FullMomentumMatrix">
2457
123
                <solver>
2458
124
                  <iterative_method name="cg"/>
2459
125
                  <preconditioner name="sor"/>
2460
126
                  <relative_error>
2461
127
                    <real_value rank="0">1E-10</real_value>
2462
128
                  </relative_error>
2463
129
                  <max_iterations>
2464
130
                    <integer_value rank="0">10000</integer_value>
2465
131
                  </max_iterations>
2466
132
                  <start_from_zero/>
2467
133
                  <never_ignore_solver_failures/>
2468
134
                  <diagnostics>
2469
135
                    <monitors/>
2470
136
                  </diagnostics>
2471
137
                </solver>
2472
138
              </inner_matrix>
2473
139
              <preconditioner_matrix name="DiagonalSchurComplement"/>
2474
140
            </full_schur_complement>
2475
141
          </use_projection_method>
2476
142
        </scheme>
2477
143
        <solver>
2478
144
          <iterative_method name="fgmres"/>
2479
145
          <preconditioner name="mg"/>
2480
146
          <relative_error>
2481
147
            <real_value rank="0">1.0e-07</real_value>
2482
148
          </relative_error>
2483
149
          <max_iterations>
2484
150
            <integer_value rank="0">10000</integer_value>
2485
151
          </max_iterations>
2486
152
          <never_ignore_solver_failures/>
2487
153
          <diagnostics>
2488
154
            <monitors/>
2489
155
          </diagnostics>
2490
156
        </solver>
2491
157
        <output>
2492
158
          <include_previous_time_step/>
2493
159
        </output>
2494
160
        <stat/>
2495
161
        <convergence>
2496
162
          <include_in_convergence/>
2497
163
        </convergence>
2498
164
        <detectors>
2499
165
          <exclude_from_detectors/>
2500
166
        </detectors>
2501
167
        <steady_state>
2502
168
          <include_in_steady_state/>
2503
169
        </steady_state>
2504
170
        <consistent_interpolation/>
2505
171
      </prognostic>
2506
172
    </scalar_field>
2507
173
    <scalar_field name="Density" rank="0">
2508
174
      <diagnostic>
2509
175
        <algorithm name="Internal" material_phase_support="multiple"/>
2510
176
        <mesh name="VelocityMesh"/>
2511
177
        <output/>
2512
178
        <stat/>
2513
179
        <convergence>
2514
180
          <include_in_convergence/>
2515
181
        </convergence>
2516
182
        <detectors>
2517
183
          <include_in_detectors/>
2518
184
        </detectors>
2519
185
        <steady_state>
2520
186
          <include_in_steady_state/>
2521
187
        </steady_state>
2522
188
      </diagnostic>
2523
189
    </scalar_field>
2524
190
    <vector_field name="Velocity" rank="1">
2525
191
      <prognostic>
2526
192
        <mesh name="VelocityMesh"/>
2527
193
        <equation name="LinearMomentum"/>
2528
194
        <spatial_discretisation>
2529
195
          <discontinuous_galerkin>
2530
196
            <mass_terms>
2531
197
              <exclude_mass_terms/>
2532
198
            </mass_terms>
2533
199
            <viscosity_scheme>
2534
200
              <bassi_rebay/>
2535
201
            </viscosity_scheme>
2536
202
            <advection_scheme>
2537
203
              <none/>
2538
204
              <integrate_advection_by_parts>
2539
205
                <twice/>
2540
206
              </integrate_advection_by_parts>
2541
207
            </advection_scheme>
2542
208
          </discontinuous_galerkin>
2543
209
          <conservative_advection>
2544
210
            <real_value rank="0">0.0</real_value>
2545
211
          </conservative_advection>
2546
212
        </spatial_discretisation>
2547
213
        <temporal_discretisation>
2548
214
          <theta>
2549
215
            <real_value rank="0">1</real_value>
2550
216
          </theta>
2551
217
          <relaxation>
2552
218
            <real_value rank="0">1</real_value>
2553
219
          </relaxation>
2554
220
        </temporal_discretisation>
2555
221
        <solver>
2556
222
          <iterative_method name="cg"/>
2557
223
          <preconditioner name="eisenstat"/>
2558
224
          <relative_error>
2559
225
            <real_value rank="0">1.E-07</real_value>
2560
226
          </relative_error>
2561
227
          <max_iterations>
2562
228
            <integer_value rank="0">10000</integer_value>
2563
229
          </max_iterations>
2564
230
          <never_ignore_solver_failures/>
2565
231
          <diagnostics>
2566
232
            <monitors/>
2567
233
          </diagnostics>
2568
234
        </solver>
2569
235
        <initial_condition name="WholeMesh">
2570
236
          <constant>
2571
237
            <real_value shape="2" dim1="dim" rank="1">0.0 0.0</real_value>
2572
238
          </constant>
2573
239
        </initial_condition>
2574
240
        <boundary_conditions name="all_sides">
2575
241
          <surface_ids>
2576
242
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
2577
243
          </surface_ids>
2578
244
          <type name="dirichlet">
2579
245
            <apply_weakly/>
2580
246
            <align_bc_with_cartesian>
2581
247
              <x_component>
2582
248
                <python>
2583
249
                  <string_value lines="20" type="python">def val(X,t):
2584
250
   from math import sin, cos
2585
251
   # Shear rotation about origin.
2586
252
   u = sin(X[0])*cos(X[1])
2587
253
   return u</string_value>
2588
254
                </python>
2589
255
              </x_component>
2590
256
              <y_component>
2591
257
                <python>
2592
258
                  <string_value lines="20" type="python">def val(X,t):
2593
259
   from math import sin, cos
2594
260
   # Shear rotation about origin.
2595
261
   v = -1.0*cos(X[0])*sin(X[1])
2596
262
   return v</string_value>
2597
263
                </python>
2598
264
              </y_component>
2599
265
            </align_bc_with_cartesian>
2600
266
          </type>
2601
267
        </boundary_conditions>
2602
268
        <tensor_field name="Viscosity" rank="2">
2603
269
          <prescribed>
2604
270
            <value name="WholeMesh">
2605
271
              <isotropic>
2606
272
                <constant>
2607
273
                  <real_value rank="0">0.7</real_value>
2608
274
                </constant>
2609
275
              </isotropic>
2610
276
            </value>
2611
277
            <output/>
2612
278
          </prescribed>
2613
279
        </tensor_field>
2614
280
        <vector_field name="Source" rank="1">
2615
281
          <prescribed>
2616
282
            <value name="WholeMesh">
2617
283
              <python>
2618
284
                <string_value lines="20" type="python">def val(XX, t):
2619
285
   from math import sin,cos
2620
286
   nu = 0.7;
2621
287
   adv = 0.0;
2622
288
   beta = 0.0;
2623
289
   mass = 0.0;
2624
290
   rho = 1.0;
2625
291
   x = XX[0];
2626
292
   y = XX[1];
2627
293
   Su = adv*rho*(cos(x)*sin(x)*sin(y)**2 + cos(x)*sin(x)*cos(y)**2) + 2*nu*sin(x)*cos(y) - sin(x)*cos(y)
2628
294
   Sv = adv*rho*(sin(x)**2*cos(y)*sin(y) + cos(x)**2*cos(y)*sin(y)) - 2*nu*cos(x)*sin(y) - cos(x)*sin(y)
2629
295
   return (Su, Sv)</string_value>
2630
296
                <comment>u(x,y,rho,beta,nu,adv,mass,t) = sin(x)*cos(y)&#13;
2631
297
v(x,y,rho,beta,nu,adv,masst) = -cos(x)*sin(y)
2632
298
p(x,y,rho,beta,nu,adv,mass,t) = cos(x)*cos(y)&#13;
2633
299
2634
300
u_t = diff(u,t)
2635
301
v_t = diff(v,t)
2636
302
2637
303
u_x = diff(u,x)&#13;
2638
304
u_y = diff(u,y)
2639
305
v_x = diff(v,x)
2640
306
v_y = diff(v,y)
2641
307
2642
308
p_x = diff(p,x)
2643
309
p_y = diff(p,y)
2644
310
&#13;
2645
311
u_xx = diff(u_x,x)
2646
312
u_yy = diff(u_y,y)
2647
313
v_xx = diff(v_x,x)
2648
314
v_yy = diff(v_y,y)
2649
315
2650
316
mass*rho*u_t + adv*rho*(u*u_x+v*u_y + beta*(u*u_x + u*v_y)) + p_x - nu*u_xx - nu*u_yy
2651
317
mass*rho*v_t + adv*rho*(u*v_x+v*v_y + beta*(v*u_x + v*v_y)) + p_y - nu*v_xx - nu*v_yy</comment>
2652
318
              </python>
2653
319
            </value>
2654
320
            <output/>
2655
321
            <stat>
2656
322
              <include_in_stat/>
2657
323
            </stat>
2658
324
            <detectors>
2659
325
              <exclude_from_detectors/>
2660
326
            </detectors>
2661
327
          </prescribed>
2662
328
        </vector_field>
2663
329
        <output>
2664
330
          <include_previous_time_step/>
2665
331
        </output>
2666
332
        <stat>
2667
333
          <include_in_stat/>
2668
334
          <previous_time_step>
2669
335
            <exclude_from_stat/>
2670
336
          </previous_time_step>
2671
337
          <nonlinear_field>
2672
338
            <exclude_from_stat/>
2673
339
          </nonlinear_field>
2674
340
        </stat>
2675
341
        <convergence>
2676
342
          <include_in_convergence/>
2677
343
        </convergence>
2678
344
        <detectors>
2679
345
          <include_in_detectors/>
2680
346
        </detectors>
2681
347
        <steady_state>
2682
348
          <include_in_steady_state/>
2683
349
        </steady_state>
2684
350
        <consistent_interpolation/>
2685
351
      </prognostic>
2686
352
    </vector_field>
2687
353
    <scalar_field name="CFLNumber" rank="0">
2688
354
      <diagnostic>
2689
355
        <algorithm name="Internal" material_phase_support="multiple"/>
2690
356
        <mesh name="VelocityMesh"/>
2691
357
        <output/>
2692
358
        <stat/>
2693
359
        <convergence>
2694
360
          <include_in_convergence/>
2695
361
        </convergence>
2696
362
        <detectors>
2697
363
          <include_in_detectors/>
2698
364
        </detectors>
2699
365
        <steady_state>
2700
366
          <include_in_steady_state/>
2701
367
        </steady_state>
2702
368
      </diagnostic>
2703
369
    </scalar_field>
2704
370
    <scalar_field name="PressureAnalytical" rank="0">
2705
371
      <prescribed>
2706
372
        <mesh name="PressureMesh"/>
2707
373
        <value name="WholeMesh">
2708
374
          <python>
2709
375
            <string_value lines="20" type="python">def val(X,t):
2710
376
   from math import cos
2711
377
   p = cos(X[0])*cos(X[1])
2712
378
   return p</string_value>
2713
379
          </python>
2714
380
        </value>
2715
381
        <output/>
2716
382
        <stat/>
2717
383
        <detectors>
2718
384
          <exclude_from_detectors/>
2719
385
        </detectors>
2720
386
      </prescribed>
2721
387
    </scalar_field>
2722
388
    <scalar_field name="ScalarAbsoluteDifference" rank="0">
2723
389
      <diagnostic field_name_b="Pressure" field_name_a="PressureAnalytical">
2724
390
        <algorithm name="Internal" material_phase_support="multiple"/>
2725
391
        <mesh name="PressureMesh"/>
2726
392
        <output/>
2727
393
        <stat/>
2728
394
        <convergence>
2729
395
          <include_in_convergence/>
2730
396
        </convergence>
2731
397
        <detectors>
2732
398
          <include_in_detectors/>
2733
399
        </detectors>
2734
400
        <steady_state>
2735
401
          <include_in_steady_state/>
2736
402
        </steady_state>
2737
403
        <relative_to_average/>
2738
404
      </diagnostic>
2739
405
    </scalar_field>
2740
406
    <scalar_field name="FiniteElementDivergence" rank="0">
2741
407
      <diagnostic field_name="Velocity">
2742
408
        <algorithm legacy="true" name="Internal" material_phase_support="multiple"/>
2743
409
        <mesh name="PressureMesh"/>
2744
410
        <integrate_divergence_by_parts/>
2745
411
        <solver>
2746
412
          <iterative_method name="cg"/>
2747
413
          <preconditioner name="sor"/>
2748
414
          <relative_error>
2749
415
            <real_value rank="0">1.E-10</real_value>
2750
416
          </relative_error>
2751
417
          <max_iterations>
2752
418
            <integer_value rank="0">10000</integer_value>
2753
419
          </max_iterations>
2754
420
          <never_ignore_solver_failures/>
2755
421
          <diagnostics>
2756
422
            <monitors/>
2757
423
          </diagnostics>
2758
424
        </solver>
2759
425
        <output/>
2760
426
        <stat/>
2761
427
        <convergence>
2762
428
          <include_in_convergence/>
2763
429
        </convergence>
2764
430
        <detectors>
2765
431
          <include_in_detectors/>
2766
432
        </detectors>
2767
433
        <steady_state>
2768
434
          <include_in_steady_state/>
2769
435
        </steady_state>
2770
436
      </diagnostic>
2771
437
    </scalar_field>
2772
438
    <vector_field name="VelocityAnalytical" rank="1">
2773
439
      <prescribed>
2774
440
        <mesh name="VelocityMesh"/>
2775
441
        <value name="WholeMesh">
2776
442
          <python>
2777
443
            <string_value lines="20" type="python">def val(X,t):
2778
444
   from math import sin, cos
2779
445
   # Shear rotation about origin.
2780
446
   return (sin(X[0])*cos(X[1]), -1.0*cos(X[0])*sin(X[1]))</string_value>
2781
447
          </python>
2782
448
        </value>
2783
449
        <output/>
2784
450
        <stat>
2785
451
          <include_in_stat/>
2786
452
        </stat>
2787
453
        <detectors>
2788
454
          <exclude_from_detectors/>
2789
455
        </detectors>
2790
456
      </prescribed>
2791
457
    </vector_field>
2792
458
    <vector_field name="VectorAbsoluteDifference" rank="1">
2793
459
      <diagnostic field_name_b="Velocity" field_name_a="VelocityAnalytical">
2794
460
        <algorithm name="Internal" material_phase_support="multiple"/>
2795
461
        <mesh name="VelocityMesh"/>
2796
462
        <output/>
2797
463
        <stat>
2798
464
          <include_in_stat/>
2799
465
        </stat>
2800
466
        <convergence>
2801
467
          <include_in_convergence/>
2802
468
        </convergence>
2803
469
        <detectors>
2804
470
          <include_in_detectors/>
2805
471
        </detectors>
2806
472
        <steady_state>
2807
473
          <include_in_steady_state/>
2808
474
        </steady_state>
2809
475
      </diagnostic>
2810
476
    </vector_field>
2811
477
  </material_phase>
2812
478
</fluidity_options>
2813
479
0
2814
=== removed file 'tests/Stokes_mms_p0p1cv/Makefile'
2815
--- tests/Stokes_mms_p0p1cv/Makefile	2011-10-21 10:05:25 +0000
2816
+++ tests/Stokes_mms_p0p1cv/Makefile	1970-01-01 00:00:00 +0000
2817
@@ -1,8 +0,0 @@
2818
1
input:	clean
2819
2
	../../bin/gmsh2triangle --2d src/MMS_A.msh
2820
3
	../../bin/gmsh2triangle --2d src/MMS_B.msh
2821
4
	../../bin/gmsh2triangle --2d src/MMS_C.msh
2822
5
	../../bin/gmsh2triangle --2d src/MMS_D.msh
2823
6
2824
7
clean:
2825
8
	rm -f *.vtu *.stat *.s *.log *.node *.ele *.edge *checkpoint*
2826
9
0
2827
=== removed file 'tests/Stokes_mms_p0p1cv/Stokes_mms_p0p1cv.xml'
2828
--- tests/Stokes_mms_p0p1cv/Stokes_mms_p0p1cv.xml	2011-10-21 10:05:25 +0000
2829
+++ tests/Stokes_mms_p0p1cv/Stokes_mms_p0p1cv.xml	1970-01-01 00:00:00 +0000
2830
@@ -1,301 +0,0 @@
2831
1
<?xml version="1.0" encoding="UTF-8" ?>
2832
2
<testproblem>
2833
3
  <name>p0p1cv stokes spatial convergence test</name>
2834
4
  <owner userid="rhodrid"/>
2835
5
  <tags>flml</tags>
2836
6
  <problem_definition length="special" nprocs="1">
2837
7
    <command_line>fluidity -v2 -l MMS_A.flml &amp;&amp; fluidity MMS_B.flml &amp;&amp; fluidity MMS_C.flml &amp;&amp; fluidity MMS_D.flml </command_line>
2838
8
  </problem_definition>
2839
9
  <variables>
2840
10
    <variable name="ab_convergence_vel" language="python">
2841
11
from fluidity_tools import stat_parser as stat
2842
12
from math import log
2843
13
2844
14
a_error_x = stat("MMS_A.stat")["NS"]["VectorAbsoluteDifference%1"]["l2norm"][-1]
2845
15
b_error_x = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%1"]["l2norm"][-1]
2846
16
a_error_y = stat("MMS_A.stat")["NS"]["VectorAbsoluteDifference%2"]["l2norm"][-1]
2847
17
b_error_y = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%2"]["l2norm"][-1]
2848
18
2849
19
a_error_inf_x = stat("MMS_A.stat")["NS"]["VectorAbsoluteDifference%1"]["max"][-1]
2850
20
b_error_inf_x = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%1"]["max"][-1]
2851
21
a_error_inf_y = stat("MMS_A.stat")["NS"]["VectorAbsoluteDifference%2"]["max"][-1]
2852
22
b_error_inf_y = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%2"]["max"][-1]
2853
23
2854
24
print a_error_x
2855
25
print b_error_x
2856
26
print a_error_y
2857
27
print b_error_y
2858
28
2859
29
print a_error_inf_x
2860
30
print b_error_inf_x
2861
31
print a_error_inf_y
2862
32
print b_error_inf_y
2863
33
2864
34
ab_ratio_x = a_error_x / b_error_x
2865
35
ab_ratio_y = a_error_y / b_error_y
2866
36
ab_ratio_inf_x = a_error_inf_x / b_error_inf_x
2867
37
ab_ratio_inf_y = a_error_inf_y / b_error_inf_y
2868
38
ab_convergence_vel = [[log(ab_ratio_x, 2), log(ab_ratio_inf_x, 2)], [log(ab_ratio_y, 2), log(ab_ratio_inf_y, 2)]]
2869
39
    </variable>
2870
40
    <variable name="ab_convergence_p" language="python">
2871
41
from fluidity_tools import stat_parser as stat
2872
42
from math import log
2873
43
2874
44
a_error = stat("MMS_A.stat")["NS"]["ScalarAbsoluteDifference"]["l2norm"][-1]
2875
45
b_error = stat("MMS_B.stat")["NS"]["ScalarAbsoluteDifference"]["l2norm"][-1]
2876
46
2877
47
a_error_inf = stat("MMS_A.stat")["NS"]["ScalarAbsoluteDifference"]["max"][-1]
2878
48
b_error_inf = stat("MMS_B.stat")["NS"]["ScalarAbsoluteDifference"]["max"][-1]
2879
49
2880
50
print a_error
2881
51
print b_error
2882
52
2883
53
print a_error_inf
2884
54
print b_error_inf
2885
55
2886
56
ab_ratio = a_error / b_error
2887
57
ab_ratio_inf = a_error_inf / b_error_inf
2888
58
ab_convergence_p = [log(ab_ratio, 2), log(ab_ratio_inf, 2)]
2889
59
    </variable>
2890
60
    <variable name="bc_convergence_vel" language="python">
2891
61
from fluidity_tools import stat_parser as stat
2892
62
from math import log
2893
63
2894
64
b_error_x = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%1"]["l2norm"][-1]
2895
65
c_error_x = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%1"]["l2norm"][-1]
2896
66
b_error_y = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%2"]["l2norm"][-1]
2897
67
c_error_y = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%2"]["l2norm"][-1]
2898
68
2899
69
b_error_inf_x = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%1"]["max"][-1]
2900
70
c_error_inf_x = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%1"]["max"][-1]
2901
71
b_error_inf_y = stat("MMS_B.stat")["NS"]["VectorAbsoluteDifference%2"]["max"][-1]
2902
72
c_error_inf_y = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%2"]["max"][-1]
2903
73
2904
74
print b_error_x
2905
75
print c_error_x
2906
76
print b_error_y
2907
77
print c_error_y
2908
78
2909
79
print b_error_inf_x
2910
80
print c_error_inf_x
2911
81
print b_error_inf_y
2912
82
print c_error_inf_y
2913
83
2914
84
bc_ratio_x = b_error_x / c_error_x
2915
85
bc_ratio_y = b_error_y / c_error_y
2916
86
bc_ratio_inf_x = b_error_inf_x / c_error_inf_x
2917
87
bc_ratio_inf_y = b_error_inf_y / c_error_inf_y
2918
88
bc_convergence_vel = [[log(bc_ratio_x, 2), log(bc_ratio_inf_x, 2)], [log(bc_ratio_y, 2), log(bc_ratio_inf_y, 2)]]
2919
89
    </variable>
2920
90
    <variable name="bc_convergence_p" language="python">
2921
91
from fluidity_tools import stat_parser as stat
2922
92
from math import log
2923
93
2924
94
b_error = stat("MMS_B.stat")["NS"]["ScalarAbsoluteDifference"]["l2norm"][-1]
2925
95
c_error = stat("MMS_C.stat")["NS"]["ScalarAbsoluteDifference"]["l2norm"][-1]
2926
96
2927
97
b_error_inf = stat("MMS_B.stat")["NS"]["ScalarAbsoluteDifference"]["max"][-1]
2928
98
c_error_inf = stat("MMS_C.stat")["NS"]["ScalarAbsoluteDifference"]["max"][-1]
2929
99
2930
100
print b_error
2931
101
print c_error
2932
102
2933
103
print b_error_inf
2934
104
print c_error_inf
2935
105
2936
106
bc_ratio = b_error / c_error
2937
107
bc_ratio_inf = b_error_inf / c_error_inf
2938
108
bc_convergence_p = [log(bc_ratio, 2), log(bc_ratio_inf, 2)]
2939
109
    </variable>
2940
110
    <variable name="cd_convergence_vel" language="python">
2941
111
from fluidity_tools import stat_parser as stat
2942
112
from math import log
2943
113
2944
114
c_error_x = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%1"]["l2norm"][-1]
2945
115
d_error_x = stat("MMS_D.stat")["NS"]["VectorAbsoluteDifference%1"]["l2norm"][-1]
2946
116
c_error_y = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%2"]["l2norm"][-1]
2947
117
d_error_y = stat("MMS_D.stat")["NS"]["VectorAbsoluteDifference%2"]["l2norm"][-1]
2948
118
2949
119
c_error_inf_x = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%1"]["max"][-1]
2950
120
d_error_inf_x = stat("MMS_D.stat")["NS"]["VectorAbsoluteDifference%1"]["max"][-1]
2951
121
c_error_inf_y = stat("MMS_C.stat")["NS"]["VectorAbsoluteDifference%2"]["max"][-1]
2952
122
d_error_inf_y = stat("MMS_D.stat")["NS"]["VectorAbsoluteDifference%2"]["max"][-1]
2953
123
2954
124
print c_error_x
2955
125
print d_error_x
2956
126
print c_error_y
2957
127
print d_error_y
2958
128
2959
129
print c_error_inf_x
2960
130
print d_error_inf_x
2961
131
print c_error_inf_y
2962
132
print d_error_inf_y
2963
133
2964
134
cd_ratio_x = c_error_x / d_error_x
2965
135
cd_ratio_y = c_error_y / d_error_y
2966
136
cd_ratio_inf_x = c_error_inf_x / d_error_inf_x
2967
137
cd_ratio_inf_y = c_error_inf_y / d_error_inf_y
2968
138
cd_convergence_vel = [[log(cd_ratio_x, 2), log(cd_ratio_inf_x, 2)], [log(cd_ratio_y, 2), log(cd_ratio_inf_y, 2)]]
2969
139
    </variable>
2970
140
    <variable name="cd_convergence_p" language="python">
2971
141
from fluidity_tools import stat_parser as stat
2972
142
from math import log
2973
143
2974
144
c_error = stat("MMS_C.stat")["NS"]["ScalarAbsoluteDifference"]["l2norm"][-1]
2975
145
d_error = stat("MMS_D.stat")["NS"]["ScalarAbsoluteDifference"]["l2norm"][-1]
2976
146
2977
147
c_error_inf = stat("MMS_C.stat")["NS"]["ScalarAbsoluteDifference"]["max"][-1]
2978
148
d_error_inf = stat("MMS_D.stat")["NS"]["ScalarAbsoluteDifference"]["max"][-1]
2979
149
2980
150
print c_error
2981
151
print d_error
2982
152
2983
153
print c_error_inf
2984
154
print d_error_inf
2985
155
2986
156
cd_ratio = c_error / d_error
2987
157
cd_ratio_inf = c_error_inf / d_error_inf
2988
158
cd_convergence_p = [log(cd_ratio, 2), log(cd_ratio_inf, 2)]
2989
159
    </variable>
2990
160
   <variable name="a_finish_time" language="python">
2991
161
from fluidity_tools import stat_parser as stat
2992
162
a_finish_time = stat("MMS_A.stat")["ElapsedTime"]["value"][-1]
2993
163
    </variable>
2994
164
    <variable name="b_finish_time" language="python">
2995
165
from fluidity_tools import stat_parser as stat
2996
166
b_finish_time = stat("MMS_B.stat")["ElapsedTime"]["value"][-1]
2997
167
    </variable>
2998
168
    <variable name="c_finish_time" language="python">
2999
169
from fluidity_tools import stat_parser as stat
3000
170
c_finish_time = stat("MMS_C.stat")["ElapsedTime"]["value"][-1]
3001
171
    </variable>
3002
172
    <variable name="d_finish_time" language="python">
3003
173
from fluidity_tools import stat_parser as stat
3004
174
d_finish_time = stat("MMS_D.stat")["ElapsedTime"]["value"][-1]
3005
175
    </variable>
3006
176
    <variable name="a_final_change_vel" language="python">
3007
177
import vtktools
3008
178
from fluidity_tools import stat_parser as stat
3009
179
vtu = vtktools.vtu("MMS_A_1.vtu")
3010
180
dt = stat("MMS_A.stat")["dt"]["value"][-1]
3011
181
a_final_change_vel = max(max(abs(vtu.GetVectorField("Velocity")[:,0]-vtu.GetVectorField("OldVelocity")[:,0])/dt), max(abs(vtu.GetVectorField("Velocity")[:,1]-vtu.GetVectorField("OldVelocity")[:,1])/dt))
3012
182
    </variable>
3013
183
    <variable name="b_final_change_vel" language="python">
3014
184
import vtktools
3015
185
from fluidity_tools import stat_parser as stat
3016
186
vtu = vtktools.vtu("MMS_B_1.vtu")
3017
187
dt = stat("MMS_B.stat")["dt"]["value"][-1]
3018
188
b_final_change_vel = max(max(abs(vtu.GetVectorField("Velocity")[:,0]-vtu.GetVectorField("OldVelocity")[:,0])/dt), max(abs(vtu.GetVectorField("Velocity")[:,1]-vtu.GetVectorField("OldVelocity")[:,1])/dt))
3019
189
    </variable>
3020
190
    <variable name="c_final_change_vel" language="python">
3021
191
import vtktools
3022
192
from fluidity_tools import stat_parser as stat
3023
193
vtu = vtktools.vtu("MMS_C_1.vtu")
3024
194
dt = stat("MMS_C.stat")["dt"]["value"][-1]
3025
195
c_final_change_vel = max(max(abs(vtu.GetVectorField("Velocity")[:,0]-vtu.GetVectorField("OldVelocity")[:,0])/dt), max(abs(vtu.GetVectorField("Velocity")[:,1]-vtu.GetVectorField("OldVelocity")[:,1])/dt))
3026
196
    </variable>
3027
197
    <variable name="d_final_change_vel" language="python">
3028
198
import vtktools
3029
199
from fluidity_tools import stat_parser as stat
3030
200
vtu = vtktools.vtu("MMS_D_1.vtu")
3031
201
dt = stat("MMS_D.stat")["dt"]["value"][-1]
3032
202
d_final_change_vel = max(max(abs(vtu.GetVectorField("Velocity")[:,0]-vtu.GetVectorField("OldVelocity")[:,0])/dt), max(abs(vtu.GetVectorField("Velocity")[:,1]-vtu.GetVectorField("OldVelocity")[:,1])/dt))
3033
203
    </variable>
3034
204
    <variable name="a_final_change_p" language="python">
3035
205
import vtktools
3036
206
from fluidity_tools import stat_parser as stat
3037
207
vtu = vtktools.vtu("MMS_A_1.vtu")
3038
208
dt = stat("MMS_A.stat")["dt"]["value"][-1]
3039
209
a_final_change_p = max(abs(vtu.GetScalarField("Pressure")-vtu.GetScalarField("OldPressure"))/dt)
3040
210
    </variable>
3041
211
    <variable name="b_final_change_p" language="python">
3042
212
import vtktools
3043
213
from fluidity_tools import stat_parser as stat
3044
214
vtu = vtktools.vtu("MMS_B_1.vtu")
3045
215
dt = stat("MMS_B.stat")["dt"]["value"][-1]
3046
216
b_final_change_p = max(abs(vtu.GetScalarField("Pressure")-vtu.GetScalarField("OldPressure"))/dt)
3047
217
    </variable>
3048
218
    <variable name="c_final_change_p" language="python">
3049
219
import vtktools
3050
220
from fluidity_tools import stat_parser as stat
3051
221
vtu = vtktools.vtu("MMS_C_1.vtu")
3052
222
dt = stat("MMS_C.stat")["dt"]["value"][-1]
3053
223
c_final_change_p = max(abs(vtu.GetScalarField("Pressure")-vtu.GetScalarField("OldPressure"))/dt)
3054
224
    </variable>
3055
225
    <variable name="d_final_change_p" language="python">
3056
226
import vtktools
3057
227
from fluidity_tools import stat_parser as stat
3058
228
vtu = vtktools.vtu("MMS_D_1.vtu")
3059
229
dt = stat("MMS_D.stat")["dt"]["value"][-1]
3060
230
d_final_change_p = max(abs(vtu.GetScalarField("Pressure")-vtu.GetScalarField("OldPressure"))/dt)
3061
231
    </variable>
3062
232
  </variables>
3063
233
  <pass_tests>
3064
234
    <test name="ab_convergence_vel_x: L2 order > 0.9" language="python">
3065
235
assert(ab_convergence_vel[0][0] &gt; 0.9)
3066
236
    </test>
3067
237
    <test name="ab_convergence_vel_y: L2 order > 0.9" language="python">
3068
238
assert(ab_convergence_vel[1][0] &gt; 0.9)
3069
239
    </test>
3070
240
    <test name="ab_convergence_p: L2 order > 0.9" language="python">
3071
241
assert(ab_convergence_p[0] &gt; 0.9)
3072
242
    </test>
3073
243
    <test name="bc_convergence_vel_x: L2 order > 0.9" language="python">
3074
244
assert(bc_convergence_vel[0][0] &gt; 0.9)
3075
245
    </test>
3076
246
    <test name="bc_convergence_vel_y: L2 order > 0.9" language="python">
3077
247
assert(bc_convergence_vel[1][0] &gt; 0.9)
3078
248
    </test>
3079
249
    <test name="bc_convergence_p: L2 order > 0.9" language="python">
3080
250
assert(bc_convergence_p[0] &gt; 0.9)
3081
251
    </test>
3082
252
    <test name="cd_convergence_vel_x: L2 order > 0.9" language="python">
3083
253
assert(cd_convergence_vel[0][0] &gt; 0.9)
3084
254
    </test>
3085
255
    <test name="cd_convergence_vel_y: L2 order > 0.9" language="python">
3086
256
assert(cd_convergence_vel[1][0] &gt; 0.9)
3087
257
    </test>
3088
258
    <test name="cd_convergence_p: L2 order > 0.9" language="python">
3089
259
assert(cd_convergence_p[0] &gt; 0.9)
3090
260
    </test>
3091
261
    <test name="checking A finished in less than 0.03" language="python">
3092
262
assert(a_finish_time -0.03 &lt; 1.E-10)
3093
263
    </test>
3094
264
    <test name="checking B finished in less than 0.03" language="python">
3095
265
assert(b_finish_time -0.03 &lt; 1.E-10)
3096
266
    </test>
3097
267
    <test name="checking C finished in less than 0.03" language="python">
3098
268
assert(c_finish_time -0.03 &lt; 1.E-10)
3099
269
    </test>
3100
270
    <test name="checking D finished in less than 0.03" language="python">
3101
271
assert(d_finish_time -0.03 &lt; 1.E-10)
3102
272
    </test>
3103
273
    <test name="checking steady state was reached in A vel" language="python">
3104
274
assert(a_final_change_vel &lt; 1.E-3)
3105
275
    </test>
3106
276
    <test name="checking steady state was reached in B vel" language="python">
3107
277
assert(b_final_change_vel &lt; 1.E-3)
3108
278
    </test>
3109
279
    <test name="checking steady state was reached in C vel" language="python">
3110
280
assert(c_final_change_vel &lt; 1.E-3)
3111
281
    </test>
3112
282
    <test name="checking steady state was reached in D vel" language="python">
3113
283
assert(d_final_change_vel &lt; 1.E-3)
3114
284
    </test>
3115
285
    <test name="checking steady state was reached in A p" language="python">
3116
286
assert(a_final_change_p &lt; 1.E-3)
3117
287
    </test>
3118
288
    <test name="checking steady state was reached in B p" language="python">
3119
289
assert(b_final_change_p &lt; 1.E-3)
3120
290
    </test>
3121
291
    <test name="checking steady state was reached in C p" language="python">
3122
292
assert(c_final_change_p &lt; 1.E-3)
3123
293
    </test>
3124
294
    <test name="checking steady state was reached in D p" language="python">
3125
295
assert(d_final_change_p &lt; 1.E-3)
3126
296
    </test>
3127
297
  </pass_tests>
3128
298
  <warn_tests>
3129
299
  </warn_tests>
3130
300
</testproblem>
3131
301
3132
302
0
3133
=== removed directory 'tests/Stokes_mms_p0p1cv/src'
3134
=== removed file 'tests/Stokes_mms_p0p1cv/src/MMS_A.geo'
3135
--- tests/Stokes_mms_p0p1cv/src/MMS_A.geo	2010-05-07 20:23:24 +0000
3136
+++ tests/Stokes_mms_p0p1cv/src/MMS_A.geo	1970-01-01 00:00:00 +0000
3137
@@ -1,15 +0,0 @@
3138
1
Point(1) = {0.0,0.0,0,0.32};
3139
2
Point(2) = {3.1415926535897931,0.0,0,0.32};
3140
3
Point(3) = {3.1415926535897931,3.1415926535897931,0,0.32};
3141
4
Point(4) = {0.0,3.1415926535897931,0,0.32};
3142
5
Line(1) = {4,3};
3143
6
Line(2) = {3,2};
3144
7
Line(3) = {2,1};
3145
8
Line(4) = {1,4};
3146
9
Line Loop(5) = {1,2,3,4};
3147
10
Plane Surface(6) = {5};
3148
11
Physical Line(7) = {3};
3149
12
Physical Line(8) = {2};
3150
13
Physical Line(9) = {1};
3151
14
Physical Line(10) = {4};
3152
15
Physical Surface(12) = {6};
3153
16
0
3154
=== removed file 'tests/Stokes_mms_p0p1cv/src/MMS_A.msh'
3155
--- tests/Stokes_mms_p0p1cv/src/MMS_A.msh	2010-05-07 20:23:24 +0000
3156
+++ tests/Stokes_mms_p0p1cv/src/MMS_A.msh	1970-01-01 00:00:00 +0000
3157
@@ -1,373 +0,0 @@
3158
1
$MeshFormat
3159
2
2 0 8
3160
3
$EndMeshFormat
3161
4
$Nodes
3162
5
122
3163
6
1 0 0 0
3164
7
2 3.141592653589793 0 0
3165
8
3 3.141592653589793 3.141592653589793 0
3166
9
4 0 3.141592653589793 0
3167
10
5 0.349065850398648 3.141592653589793 0
3168
11
6 0.698131700797296 3.141592653589793 0
3169
12
7 1.04719755119647 3.141592653589793 0
3170
13
8 1.39626340159582 3.141592653589793 0
3171
14
9 1.745329251995478 3.141592653589793 0
3172
15
10 2.094395102395445 3.141592653589793 0
3173
16
11 2.443460952794795 3.141592653589793 0
3174
17
12 2.792526803192294 3.141592653589793 0
3175
18
13 3.141592653589793 2.79252680319031 0
3176
19
14 3.141592653589793 2.443460952790828 0
3177
20
15 3.141592653589793 2.094395102391345 0
3178
21
16 3.141592653589793 1.745329251991862 0
3179
22
17 3.141592653589793 1.39626340159238 0
3180
23
18 3.141592653589793 1.047197551192897 0
3181
24
19 3.141592653589793 0.6981317007940313 0
3182
25
20 3.141592653589793 0.3490658503970154 0
3183
26
21 2.79252680319031 0 0
3184
27
22 2.443460952790828 0 0
3185
28
23 2.094395102391345 0 0
3186
29
24 1.745329251991862 0 0
3187
30
25 1.39626340159238 0 0
3188
31
26 1.047197551192897 0 0
3189
32
27 0.6981317007940313 0 0
3190
33
28 0.3490658503970154 0 0
3191
34
29 0 0.349065850398648 0
3192
35
30 0 0.698131700797296 0
3193
36
31 0 1.04719755119647 0
3194
37
32 0 1.39626340159582 0
3195
38
33 0 1.745329251995478 0
3196
39
34 0 2.094395102395445 0
3197
40
35 0 2.443460952794795 0
3198
41
36 0 2.792526803192294 0
3199
42
37 2.245331016203158 2.237477093448061 0
3200
43
38 1.939212850091832 1.488667447826441 0
3201
44
39 0.4805987052845769 1.481941171402636 0
3202
45
40 1.375039289697466 2.531093011873753 0
3203
46
41 1.748454674281334 2.006877314733054 0
3204
47
42 1.8413234790713 0.3414603770972079 0
3205
48
43 0.2647906118507061 0.8690486946178108 0
3206
49
44 1.171580759472431 1.202053915551927 0
3207
50
45 2.010303042677967 2.681011374302387 0
3208
51
46 2.782496476113133 2.551712022195696 0
3209
52
47 0.7850460661416393 2.837550150787113 0
3210
53
48 0.9124421070782256 2.132069388734857 0
3211
54
49 0.2226014046652034 1.958442683592825 0
3212
55
50 0.241860467851656 2.553231261758386 0
3213
56
51 2.493755826600072 1.399197545233538 0
3214
57
52 1.684307131411933 1.652051739326319 0
3215
58
53 2.123029470445046 1.995647508820058 0
3216
59
54 0.2350335459538366 1.491064457657457 0
3217
60
55 0.9804807296768365 1.51485526305302 0
3218
61
56 0.5363395983201773 1.993552794837632 0
3219
62
57 1.435777128333467 2.876001160134476 0
3220
63
58 1.346020536528309 2.147871846007586 0
3221
64
59 1.947386564510164 2.369844994536343 0
3222
65
60 0.8436424848152182 2.513204947469001 0
3223
66
61 1.458433528545092 1.815345563296753 0
3224
67
62 2.527862670756731 1.082068937537555 0
3225
68
63 2.024545858803479 1.172807122982141 0
3226
69
64 2.537203282389552 2.241741256935653 0
3227
70
65 2.845190229270507 0.9032782271684311 0
3228
71
66 2.196647011217093 0.2907308185295358 0
3229
72
67 2.26945679302129 0.9079902985898242 0
3230
73
68 0.3812530977240584 0.6134794414300975 0
3231
74
69 2.360921103780199 1.720989337040876 0
3232
75
70 0.3607654427034569 2.885456200345249 0
3233
76
71 0.2761773190375531 2.274367702446057 0
3234
77
72 1.375582385360377 1.506097837016471 0
3235
78
73 0.3203185168137341 1.203516306009632 0
3236
79
74 0.3494632020933188 1.71876393612117 0
3237
80
75 0.6738549898012032 1.297248441378809 0
3238
81
76 1.684217955851304 2.70643444402808 0
3239
82
77 1.599991567676058 2.377063147082594 0
3240
83
78 1.126113714829085 2.43915168813852 0
3241
84
79 0.7728374833067017 0.6267909252667454 0
3242
85
80 1.689420976530337 1.29367521559253 0
3243
86
81 2.822497879170552 1.272041449704428 0
3244
87
82 2.281684470019294 1.140700436727976 0
3245
88
83 2.214421965622866 1.410336338841596 0
3246
89
84 2.369606313367625 2.018880092251258 0
3247
90
85 2.535355324434343 0.7496541746228972 0
3248
91
86 1.110636296117287 0.5848766808981929 0
3249
92
87 1.730037176915837 0.7344628669538847 0
3250
93
88 1.18185425377986 0.2541827504744055 0
3251
94
89 1.505061768143136 0.2283365892598895 0
3252
95
90 2.696800710726369 1.603790382737259 0
3253
96
91 2.058874293757006 1.684998687961457 0
3254
97
92 2.629210384431764 1.92457545984789 0
3255
98
93 2.339377584302134 2.540772651868505 0
3256
99
94 0.5162073646051875 2.568841404134913 0
3257
100
95 1.152327727137135 2.79629955199542 0
3258
101
96 0.5923536438084981 0.9208652715981994 0
3259
102
97 2.855494253163295 2.181711512090013 0
3260
103
98 0.6716279419198647 1.657495117542256 0
3261
104
99 1.160478616613587 1.835174145137185 0
3262
105
100 1.453686683434912 0.5069572601828125 0
3263
106
101 0.5810966210891484 2.297418966154744 0
3264
107
102 2.809070899053095 0.218901494312215 0
3265
108
103 2.487657761817074 0.4163827666359571 0
3266
109
104 2.223176270471563 2.864617743636775 0
3267
110
105 1.913747817514841 2.918272375122196 0
3268
111
106 1.793906512578577 1.003757118136154 0
3269
112
107 0.9804986137608793 0.9008312905398148 0
3270
113
108 0.8938343370954145 0.3066847625821323 0
3271
114
109 2.829593115791567 0.5600474173912698 0
3272
115
110 2.898634635889831 2.889041028524417 0
3273
116
111 2.581080305526595 2.829905597609742 0
3274
117
112 0.8407956278336383 1.849857759732165 0
3275
118
113 2.905913288695151 1.835562116980959 0
3276
119
114 1.561116438886176 1.039730365788334 0
3277
120
115 0.53725735778748 0.2993180933067856 0
3278
121
116 2.004220755858769 0.8385154074300597 0
3279
122
117 2.93579076722317 1.546198917550701 0
3280
123
118 1.422474588967255 1.215876736816864 0
3281
124
119 2.139327863413657 0.5942903098409081 0
3282
125
120 0.8765942977129461 1.137678729790205 0
3283
126
121 1.323189517933578 0.9135733591372802 0
3284
127
122 0.199019453500453 0.2759548735034221 0
3285
128
$EndNodes
3286
129
$Elements
3287
130
242
3288
131
1 1 3 9 1 0 4 5
3289
132
2 1 3 9 1 0 5 6
3290
133
3 1 3 9 1 0 6 7
3291
134
4 1 3 9 1 0 7 8
3292
135
5 1 3 9 1 0 8 9
3293
136
6 1 3 9 1 0 9 10
3294
137
7 1 3 9 1 0 10 11
3295
138
8 1 3 9 1 0 11 12
3296
139
9 1 3 9 1 0 12 3
3297
140
10 1 3 8 2 0 3 13
3298
141
11 1 3 8 2 0 13 14
3299
142
12 1 3 8 2 0 14 15
3300
143
13 1 3 8 2 0 15 16
3301
144
14 1 3 8 2 0 16 17
3302
145
15 1 3 8 2 0 17 18
3303
146
16 1 3 8 2 0 18 19
3304
147
17 1 3 8 2 0 19 20
3305
148
18 1 3 8 2 0 20 2
3306
149
19 1 3 7 3 0 2 21
3307
150
20 1 3 7 3 0 21 22
3308
151
21 1 3 7 3 0 22 23
3309
152
22 1 3 7 3 0 23 24
3310
153
23 1 3 7 3 0 24 25
3311
154
24 1 3 7 3 0 25 26
3312
155
25 1 3 7 3 0 26 27
3313
156
26 1 3 7 3 0 27 28
3314
157
27 1 3 7 3 0 28 1
3315
158
28 1 3 10 4 0 1 29
3316
159
29 1 3 10 4 0 29 30
3317
160
30 1 3 10 4 0 30 31
3318
161
31 1 3 10 4 0 31 32
3319
162
32 1 3 10 4 0 32 33
3320
163
33 1 3 10 4 0 33 34
3321
164
34 1 3 10 4 0 34 35
3322
165
35 1 3 10 4 0 35 36
3323
166
36 1 3 10 4 0 36 4
3324
167
37 2 3 12 6 0 24 42 23
3325
168
38 2 3 12 6 0 43 30 31
3326
169
39 2 3 12 6 0 47 6 7
3327
170
40 2 3 12 6 0 34 49 33
3328
171
41 2 3 12 6 0 41 61 58
3329
172
42 2 3 12 6 0 61 41 52
3330
173
43 2 3 12 6 0 70 36 4
3331
174
44 2 3 12 6 0 36 50 35
3332
175
45 2 3 12 6 0 36 70 50
3333
176
46 2 3 12 6 0 32 33 54
3334
177
47 2 3 12 6 0 4 5 70
3335
178
48 2 3 12 6 0 71 49 34
3336
179
49 2 3 12 6 0 73 43 31
3337
180
50 2 3 12 6 0 74 54 33
3338
181
51 2 3 12 6 0 49 74 33
3339
182
52 2 3 12 6 0 38 80 52
3340
183
53 2 3 12 6 0 62 81 65
3341
184
54 2 3 12 6 0 82 62 67
3342
185
55 2 3 12 6 0 84 53 37
3343
186
56 2 3 12 6 0 64 84 37
3344
187
57 2 3 12 6 0 37 93 64
3345
188
58 2 3 12 6 0 94 50 70
3346
189
59 2 3 12 6 0 95 7 8
3347
190
60 2 3 12 6 0 57 95 8
3348
191
61 2 3 12 6 0 96 68 43
3349
192
62 2 3 12 6 0 97 14 15
3350
193
63 2 3 12 6 0 8 9 57
3351
194
64 2 3 12 6 0 76 57 9
3352
195
65 2 3 12 6 0 5 6 70
3353
196
66 2 3 12 6 0 49 71 56
3354
197
67 2 3 12 6 0 72 44 55
3355
198
68 2 3 12 6 0 72 61 52
3356
199
69 2 3 12 6 0 80 72 52
3357
200
70 2 3 12 6 0 39 73 54
3358
201
71 2 3 12 6 0 39 75 73
3359
202
72 2 3 12 6 0 74 39 54
3360
203
73 2 3 12 6 0 74 49 56
3361
204
74 2 3 12 6 0 76 40 57
3362
205
75 2 3 12 6 0 77 40 76
3363
206
76 2 3 12 6 0 76 45 59
3364
207
77 2 3 12 6 0 59 77 76
3365
208
78 2 3 12 6 0 40 77 58
3366
209
79 2 3 12 6 0 41 77 59
3367
210
80 2 3 12 6 0 77 41 58
3368
211
81 2 3 12 6 0 78 40 58
3369
212
82 2 3 12 6 0 78 48 60
3370
213
83 2 3 12 6 0 48 78 58
3371
214
84 2 3 12 6 0 68 96 79
3372
215
85 2 3 12 6 0 80 38 63
3373
216
86 2 3 12 6 0 51 81 62
3374
217
87 2 3 12 6 0 51 90 81
3375
218
88 2 3 12 6 0 82 51 62
3376
219
89 2 3 12 6 0 83 51 82
3377
220
90 2 3 12 6 0 63 83 82
3378
221
91 2 3 12 6 0 82 67 63
3379
222
92 2 3 12 6 0 38 83 63
3380
223
93 2 3 12 6 0 38 91 83
3381
224
94 2 3 12 6 0 51 83 69
3382
225
95 2 3 12 6 0 91 69 83
3383
226
96 2 3 12 6 0 53 84 69
3384
227
97 2 3 12 6 0 92 69 84
3385
228
98 2 3 12 6 0 23 42 66
3386
229
99 2 3 12 6 0 17 18 81
3387
230
100 2 3 12 6 0 18 65 81
3388
231
101 2 3 12 6 0 62 85 67
3389
232
102 2 3 12 6 0 85 62 65
3390
233
103 2 3 12 6 0 25 26 88
3391
234
104 2 3 12 6 0 89 25 88
3392
235
105 2 3 12 6 0 90 51 69
3393
236
106 2 3 12 6 0 69 92 90
3394
237
107 2 3 12 6 0 91 38 52
3395
238
108 2 3 12 6 0 91 53 69
3396
239
109 2 3 12 6 0 91 41 53
3397
240
110 2 3 12 6 0 41 91 52
3398
241
111 2 3 12 6 0 94 47 60
3399
242
112 2 3 12 6 0 94 71 50
3400
243
113 2 3 12 6 0 31 32 73
3401
244
114 2 3 12 6 0 32 54 73
3402
245
115 2 3 12 6 0 24 25 89
3403
246
116 2 3 12 6 0 42 24 89
3404
247
117 2 3 12 6 0 95 47 7
3405
248
118 2 3 12 6 0 40 95 57
3406
249
119 2 3 12 6 0 40 78 95
3407
250
120 2 3 12 6 0 18 19 65
3408
251
121 2 3 12 6 0 46 97 64
3409
252
122 2 3 12 6 0 97 46 14
3410
253
123 2 3 12 6 0 34 35 71
3411
254
124 2 3 12 6 0 35 50 71
3412
255
125 2 3 12 6 0 99 61 72
3413
256
126 2 3 12 6 0 55 99 72
3414
257
127 2 3 12 6 0 86 100 88
3415
258
128 2 3 12 6 0 100 89 88
3416
259
129 2 3 12 6 0 94 101 71
3417
260
130 2 3 12 6 0 101 56 71
3418
261
131 2 3 12 6 0 2 102 20
3419
262
132 2 3 12 6 0 102 2 21
3420
263
133 2 3 12 6 0 39 98 75
3421
264
134 2 3 12 6 0 98 55 75
3422
265
135 2 3 12 6 0 93 46 64
3423
266
136 2 3 12 6 0 45 104 93
3424
267
137 2 3 12 6 0 99 48 58
3425
268
138 2 3 12 6 0 61 99 58
3426
269
139 2 3 12 6 0 48 101 60
3427
270
140 2 3 12 6 0 101 48 56
3428
271
141 2 3 12 6 0 101 94 60
3429
272
142 2 3 12 6 0 104 10 11
3430
273
143 2 3 12 6 0 9 105 76
3431
274
144 2 3 12 6 0 45 105 104
3432
275
145 2 3 12 6 0 105 45 76
3433
276
146 2 3 12 6 0 27 108 26
3434
277
147 2 3 12 6 0 108 88 26
3435
278
148 2 3 12 6 0 20 109 19
3436
279
149 2 3 12 6 0 110 12 3
3437
280
150 2 3 12 6 0 13 110 3
3438
281
151 2 3 12 6 0 93 111 46
3439
282
152 2 3 12 6 0 92 84 64
3440
283
153 2 3 12 6 0 97 92 64
3441
284
154 2 3 12 6 0 98 39 74
3442
285
155 2 3 12 6 0 56 98 74
3443
286
156 2 3 12 6 0 63 106 80
3444
287
157 2 3 12 6 0 79 107 86
3445
288
158 2 3 12 6 0 88 108 86
3446
289
159 2 3 12 6 0 108 79 86
3447
290
160 2 3 12 6 0 9 10 105
3448
291
161 2 3 12 6 0 10 104 105
3449
292
162 2 3 12 6 0 48 112 56
3450
293
163 2 3 12 6 0 112 98 56
3451
294
164 2 3 12 6 0 108 115 79
3452
295
165 2 3 12 6 0 115 68 79
3453
296
166 2 3 12 6 0 43 73 96
3454
297
167 2 3 12 6 0 75 96 73
3455
298
168 2 3 12 6 0 107 79 96
3456
299
169 2 3 12 6 0 100 87 42
3457
300
170 2 3 12 6 0 42 89 100
3458
301
171 2 3 12 6 0 109 85 65
3459
302
172 2 3 12 6 0 19 109 65
3460
303
173 2 3 12 6 0 12 111 11
3461
304
174 2 3 12 6 0 111 104 11
3462
305
175 2 3 12 6 0 104 111 93
3463
306
176 2 3 12 6 0 112 48 99
3464
307
177 2 3 12 6 0 55 112 99
3465
308
178 2 3 12 6 0 112 55 98
3466
309
179 2 3 12 6 0 113 97 15
3467
310
180 2 3 12 6 0 97 113 92
3468
311
181 2 3 12 6 0 16 113 15
3469
312
182 2 3 12 6 0 113 90 92
3470
313
183 2 3 12 6 0 67 116 63
3471
314
184 2 3 12 6 0 116 106 63
3472
315
185 2 3 12 6 0 117 17 81
3473
316
186 2 3 12 6 0 117 90 113
3474
317
187 2 3 12 6 0 90 117 81
3475
318
188 2 3 12 6 0 85 109 103
3476
319
189 2 3 12 6 0 110 111 12
3477
320
190 2 3 12 6 0 111 110 46
3478
321
191 2 3 12 6 0 68 30 43
3479
322
192 2 3 12 6 0 106 114 80
3480
323
193 2 3 12 6 0 87 114 106
3481
324
194 2 3 12 6 0 116 87 106
3482
325
195 2 3 12 6 0 16 17 117
3483
326
196 2 3 12 6 0 113 16 117
3484
327
197 2 3 12 6 0 118 72 80
3485
328
198 2 3 12 6 0 114 118 80
3486
329
199 2 3 12 6 0 72 118 44
3487
330
200 2 3 12 6 0 21 22 102
3488
331
201 2 3 12 6 0 109 20 102
3489
332
202 2 3 12 6 0 109 102 103
3490
333
203 2 3 12 6 0 46 13 14
3491
334
204 2 3 12 6 0 13 46 110
3492
335
205 2 3 12 6 0 115 27 28
3493
336
206 2 3 12 6 0 115 108 27
3494
337
207 2 3 12 6 0 119 103 66
3495
338
208 2 3 12 6 0 103 119 85
3496
339
209 2 3 12 6 0 67 119 116
3497
340
210 2 3 12 6 0 119 67 85
3498
341
211 2 3 12 6 0 119 42 87
3499
342
212 2 3 12 6 0 42 119 66
3500
343
213 2 3 12 6 0 116 119 87
3501
344
214 2 3 12 6 0 75 120 96
3502
345
215 2 3 12 6 0 120 107 96
3503
346
216 2 3 12 6 0 44 121 107
3504
347
217 2 3 12 6 0 103 22 66
3505
348
218 2 3 12 6 0 22 103 102
3506
349
219 2 3 12 6 0 23 66 22
3507
350
220 2 3 12 6 0 44 120 55
3508
351
221 2 3 12 6 0 120 44 107
3509
352
222 2 3 12 6 0 120 75 55
3510
353
223 2 3 12 6 0 100 121 87
3511
354
224 2 3 12 6 0 86 121 100
3512
355
225 2 3 12 6 0 121 114 87
3513
356
226 2 3 12 6 0 107 121 86
3514
357
227 2 3 12 6 0 118 121 44
3515
358
228 2 3 12 6 0 118 114 121
3516
359
229 2 3 12 6 0 60 95 78
3517
360
230 2 3 12 6 0 95 60 47
3518
361
231 2 3 12 6 0 59 93 37
3519
362
232 2 3 12 6 0 93 59 45
3520
363
233 2 3 12 6 0 47 70 6
3521
364
234 2 3 12 6 0 70 47 94
3522
365
235 2 3 12 6 0 53 59 37
3523
366
236 2 3 12 6 0 41 59 53
3524
367
237 2 3 12 6 0 122 28 1
3525
368
238 2 3 12 6 0 29 122 1
3526
369
239 2 3 12 6 0 122 29 30
3527
370
240 2 3 12 6 0 122 115 28
3528
371
241 2 3 12 6 0 122 68 115
3529
372
242 2 3 12 6 0 68 122 30
3530
373
$EndElements
3531
374
0
3532
=== removed file 'tests/Stokes_mms_p0p1cv/src/MMS_B.geo'
3533
--- tests/Stokes_mms_p0p1cv/src/MMS_B.geo	2010-05-07 20:23:24 +0000
3534
+++ tests/Stokes_mms_p0p1cv/src/MMS_B.geo	1970-01-01 00:00:00 +0000
3535
@@ -1,15 +0,0 @@
3536
1
Point(1) = {0.0,0.0,0,0.16};
3537
2
Point(2) = {3.1415926535897931,0.0,0,0.16};
3538
3
Point(3) = {3.1415926535897931,3.1415926535897931,0,0.16};
3539
4
Point(4) = {0.0,3.1415926535897931,0,0.16};
3540
5
Line(1) = {4,3};
3541
6
Line(2) = {3,2};
3542
7
Line(3) = {2,1};
3543
8
Line(4) = {1,4};
3544
9
Line Loop(5) = {1,2,3,4};
3545
10
Plane Surface(6) = {5};
3546
11
Physical Line(7) = {3};
3547
12
Physical Line(8) = {2};
3548
13
Physical Line(9) = {1};
3549
14
Physical Line(10) = {4};
3550
15
Physical Surface(12) = {6};
3551
16
0
3552
=== removed file 'tests/Stokes_mms_p0p1cv/src/MMS_B.msh'
3553
--- tests/Stokes_mms_p0p1cv/src/MMS_B.msh	2010-05-07 20:23:24 +0000
3554
+++ tests/Stokes_mms_p0p1cv/src/MMS_B.msh	1970-01-01 00:00:00 +0000
3555
@@ -1,1291 +0,0 @@
3556
1
$MeshFormat
3557
2
2 0 8
3558
3
$EndMeshFormat
3559
4
$Nodes
3560
5
428
3561
6
1 0 0 0
3562
7
2 3.141592653589793 0 0
3563
8
3 3.141592653589793 3.141592653589793 0
3564
9
4 0 3.141592653589793 0
3565
10
5 0.1653469817677806 3.141592653589793 0
3566
11
6 0.3306939635355612 3.141592653589793 0
3567
12
7 0.4960409453033418 3.141592653589793 0
3568
13
8 0.6613879270711225 3.141592653589793 0
3569
14
9 0.8267349088389864 3.141592653589793 0
3570
15
10 0.9920818906070992 3.141592653589793 0
3571
16
11 1.157428872375212 3.141592653589793 0
3572
17
12 1.322775854143325 3.141592653589793 0
3573
18
13 1.488122835911438 3.141592653589793 0
3574
19
14 1.653469817679697 3.141592653589793 0
3575
20
15 1.818816799448103 3.141592653589793 0
3576
21
16 1.984163781216508 3.141592653589793 0
3577
22
17 2.149510762984913 3.141592653589793 0
3578
23
18 2.314857744753318 3.141592653589793 0
3579
24
19 2.480204726520847 3.141592653589793 0
3580
25
20 2.645551708288084 3.141592653589793 0
3581
26
21 2.81089869005532 3.141592653589793 0
3582
27
22 2.976245671822556 3.141592653589793 0
3583
28
23 3.141592653589793 2.976245671821617 0
3584
29
24 3.141592653589793 2.810898690053441 0
3585
30
25 3.141592653589793 2.645551708285265 0
3586
31
26 3.141592653589793 2.480204726517089 0
3587
32
27 3.141592653589793 2.314857744748913 0
3588
33
28 3.141592653589793 2.149510762980737 0
3589
34
29 3.141592653589793 1.984163781212561 0
3590
35
30 3.141592653589793 1.818816799444385 0
3591
36
31 3.141592653589793 1.653469817676209 0
3592
37
32 3.141592653589793 1.488122835908033 0
3593
38
33 3.141592653589793 1.322775854139857 0
3594
39
34 3.141592653589793 1.157428872371681 0
3595
40
35 3.141592653589793 0.9920818906035049 0
3596
41
36 3.141592653589793 0.8267349088353289 0
3597
42
37 3.141592653589793 0.6613879270680294 0
3598
43
38 3.141592653589793 0.4960409453010217 0
3599
44
39 3.141592653589793 0.3306939635340145 0
3600
45
40 3.141592653589793 0.1653469817670068 0
3601
46
41 2.976245671821617 0 0
3602
47
42 2.810898690053441 0 0
3603
48
43 2.645551708285265 0 0
3604
49
44 2.480204726517089 0 0
3605
50
45 2.314857744748913 0 0
3606
51
46 2.149510762980737 0 0
3607
52
47 1.984163781212561 0 0
3608
53
48 1.818816799444385 0 0
3609
54
49 1.653469817676209 0 0
3610
55
50 1.488122835908033 0 0
3611
56
51 1.322775854139857 0 0
3612
57
52 1.157428872371681 0 0
3613
58
53 0.9920818906035049 0 0
3614
59
54 0.8267349088353289 0 0
3615
60
55 0.6613879270680294 0 0
3616
61
56 0.4960409453010217 0 0
3617
62
57 0.3306939635340145 0 0
3618
63
58 0.1653469817670068 0 0
3619
64
59 0 0.1653469817677806 0
3620
65
60 0 0.3306939635355612 0
3621
66
61 0 0.4960409453033418 0
3622
67
62 0 0.6613879270711225 0
3623
68
63 0 0.8267349088389864 0
3624
69
64 0 0.9920818906070992 0
3625
70
65 0 1.157428872375212 0
3626
71
66 0 1.322775854143325 0
3627
72
67 0 1.488122835911438 0
3628
73
68 0 1.653469817679697 0
3629
74
69 0 1.818816799448103 0
3630
75
70 0 1.984163781216508 0
3631
76
71 0 2.149510762984913 0
3632
77
72 0 2.314857744753318 0
3633
78
73 0 2.480204726520847 0
3634
79
74 0 2.645551708288084 0
3635
80
75 0 2.81089869005532 0
3636
81
76 0 2.976245671822556 0
3637
82
77 2.710154345613917 0.7124035191059989 0
3638
83
78 1.335705447144501 1.537291746312253 0
3639
84
79 0.7792697275134636 0.99155755697381 0
3640
85
80 0.473361672149427 0.5418702562722421 0
3641
86
81 2.823087642901273 0.3079710975474592 0
3642
87
82 2.767078285855207 2.593662908315215 0
3643
88
83 1.658230909855657 0.6993179675070205 0
3644
89
84 2.477046825221157 0.8983989183155836 0
3645
90
85 2.286546527540184 0.5291155451054952 0
3646
91
86 2.763441830817289 0.1540079993047146 0
3647
92
87 0.9310433468591642 0.3631462188368522 0
3648
93
88 2.880590855368471 0.1914635091263185 0
3649
94
89 2.044444915452922 2.598831894750799 0
3650
95
90 2.453220504963881 2.694332345363849 0
3651
96
91 2.942829198104365 2.488434356447135 0
3652
97
92 2.540561744397069 2.974795387608982 0
3653
98
93 0.1671901984448033 2.205538451778842 0
3654
99
94 0.9240428927373574 2.055219669428184 0
3655
100
95 1.564809982517981 1.831524450922345 0
3656
101
96 2.206327931818251 2.138406934106811 0
3657
102
97 0.9811357034587043 2.441202898590433 0
3658
103
98 1.632942267552848 2.193404877477654 0
3659
104
99 1.299583009832575 1.039221487677734 0
3660
105
100 0.3035032222759849 1.695110777719873 0
3661
106
101 1.96947052174502 1.144479511012929 0
3662
107
102 1.58269508073549 1.26184070387476 0
3663
108
103 2.425109744138061 1.180221257733585 0
3664
109
104 2.723716816817247 0.979342777164218 0
3665
110
105 2.095691154277549 0.7791856729143436 0
3666
111
106 2.353379159360377 0.8459424898903043 0
3667
112
107 2.515199651387261 0.7244646397736486 0
3668
113
108 2.553648424475276 1.050966575365818 0
3669
114
109 1.11071832585682 0.9840325084045711 0
3670
115
110 0.9289740643403908 0.1585447704921896 0
3671
116
111 2.294460657733453 2.476400761853999 0
3672
117
112 1.951287116773943 2.283306621468514 0
3673
118
113 2.109760924474123 1.791468361544839 0
3674
119
114 2.688912087705139 1.892009783199424 0
3675
120
115 2.401227236877782 2.577292805750611 0
3676
121
116 2.640159959712249 2.433090295279496 0
3677
122
117 0.1090382934827166 1.900967598307153 0
3678
123
118 0.4678553322226026 2.405625195819564 0
3679
124
119 1.089494724561111 1.752787835249928 0
3680
125
120 0.7176582363183167 2.304965125626813 0
3681
126
121 0.1694144837145111 2.746174807704512 0
3682
127
122 2.888436227529274 0.9062277134744589 0
3683
128
123 2.243924545238765 0.7838910206692115 0
3684
129
124 0.1669969730218628 1.421715088176076 0
3685
130
125 0.5914734535192498 1.021337589027003 0
3686
131
126 2.66972695866243 0.2910078877299449 0
3687
132
127 0.1652780543481153 1.085153939649221 0
3688
133
128 2.884253770498661 0.7470890018964499 0
3689
134
129 3.002027911171429 1.718940992262709 0
3690
135
130 2.15761823550225 2.918575049244855 0
3691
136
131 2.181887167942433 2.734039762339985 0
3692
137
132 2.35433176432498 2.980987417517035 0
3693
138
133 2.661778219387621 2.859399783359583 0
3694
139
134 2.695772180708321 2.270698610323994 0
3695
140
135 2.938604072962539 2.932107199689162 0
3696
141
136 2.969844772018041 2.63124359395824 0
3697
142
137 3.036448711113867 3.032732175523536 0
3698
143
138 3.022534085109895 2.76676177911272 0
3699
144
139 1.686607241758467 2.790087479645788 0
3700
145
140 1.533433507278696 2.441656265951141 0
3701
146
141 1.919822703006044 1.727223384093773 0
3702
147
142 1.813923339990162 1.976987356695722 0
3703
148
143 1.455707771832521 2.070047537570578 0
3704
149
144 1.995978352592519 2.091349223262582 0
3705
150
145 2.108298862994522 2.423650675329858 0
3706
151
146 2.284547605920695 1.869199921200383 0
3707
152
147 1.345235098133746 2.66114913664867 0
3708
153
148 1.604268139703333 2.612495845366983 0
3709
154
149 0.7480073345023147 2.496672244404683 0
3710
155
150 0.9157863229348679 2.255055440610529 0
3711
156
151 1.44274325799834 0.2630952735574613 0
3712
157
152 1.607676446880413 0.5385375169385637 0
3713
158
153 0.9346302377518008 1.871397278260167 0
3714
159
154 0.3745910538022397 2.261605443231008 0
3715
160
155 0.7456713452935297 2.124121853483894 0
3716
161
156 2.115804446734388 1.094312619518252 0
3717
162
157 1.66238584006893 1.553670206395838 0
3718
163
158 2.253275892353718 1.163631830974447 0
3719
164
159 2.736922052829414 1.120090998796195 0
3720
165
160 2.60352163203479 1.216270600568891 0
3721
166
161 2.454186324495676 0.536611096451201 0
3722
167
162 1.285133235198161 1.223330086912121 0
3723
168
163 2.315756432187607 2.671310410452097 0
3724
169
164 2.812657142602099 2.396292962413941 0
3725
170
165 0.136852048977916 2.034883063678378 0
3726
171
166 0.5418281626293912 2.308493077677789 0
3727
172
167 0.443536649234231 2.545627860579855 0
3728
173
168 1.100311822943816 1.946022059238785 0
3729
174
169 0.8539321934869285 2.40197529637043 0
3730
175
170 1.186237453043325 2.288334197967351 0
3731
176
171 0.7652856616157222 1.943984178166221 0
3732
177
172 0.6098194245086899 2.985574538487605 0
3733
178
173 2.539456843526181 2.031240610128277 0
3734
179
174 1.833344553071007 0.898979730527855 0
3735
180
175 1.769719959305706 1.194304879895215 0
3736
181
176 1.518261546581936 0.408405483080395 0
3737
182
177 1.408367353727112 1.353815562665538 0
3738
183
178 1.932180696962042 0.7459977838391768 0
3739
184
179 2.3420997702715 0.6759969493163754 0
3740
185
180 2.211904838264805 0.9456172545347595 0
3741
186
181 2.182619609175662 0.6338366984510819 0
3742
187
182 0.6660988608060098 0.8519320040293099 0
3743
188
183 2.824632978175203 0.5699197156284113 0
3744
189
184 2.615721458511351 0.5489570794216285 0
3745
190
185 2.196430927761531 2.591538455810981 0
3746
191
186 2.0240788720558 2.772637017881345 0
3747
192
187 2.465862732904224 2.420217185392864 0
3748
193
188 2.568220559737175 2.600182804832682 0
3749
194
189 0.3687128050060238 2.783987860321025 0
3750
195
190 2.988992031784484 0.1172880830341935 0
3751
196
191 0.3059174942391116 0.1417643898823577 0
3752
197
192 0.126996997968816 0.329995198176525 0
3753
198
193 3.006955753660114 2.306404028916538 0
3754
199
194 0.1425553326982193 2.924554374540718 0
3755
200
195 0.2907902169237725 2.895253588119128 0
3756
201
196 2.153054405587897 0.465094686595338 0
3757
202
197 0.3155916798373833 0.5015827428176055 0
3758
203
198 0.179903986983802 1.250989382466576 0
3759
204
199 0.5315144135057466 0.7232508287032177 0
3760
205
200 0.5079504532573573 0.8928950868888056 0
3761
206
201 0.3449217198866198 0.8610334168910246 0
3762
207
202 2.579334723265651 0.1766533029758535 0
3763
208
203 2.414110454944511 2.820266676899157 0
3764
209
204 2.69212459811774 2.724412480751462 0
3765
210
205 2.711428055510238 2.099743093885627 0
3766
211
206 1.805763900278314 2.146832042461267 0
3767
212
207 1.353056848686425 2.455241293224084 0
3768
213
208 1.613843243420686 2.333732165484082 0
3769
214
209 1.84944905926157 2.781396682077941 0
3770
215
210 1.96323967183751 1.89679428432906 0
3771
216
211 2.040658004318893 1.602334017109482 0
3772
217
212 1.843169709479856 1.534817886649233 0
3773
218
213 2.161738227762281 1.484011418071472 0
3774
219
214 1.646199642992884 2.037888001505255 0
3775
220
215 1.729033014210835 1.730430744173156 0
3776
221
216 1.456899565782772 2.279117622191337 0
3777
222
217 2.340388289307997 2.006693281832478 0
3778
223
218 2.084450381363953 2.244585644332628 0
3779
224
219 2.135363929179332 1.960323208439774 0
3780
225
220 1.588628191525907 2.964236847596425 0
3781
226
221 1.470144871603232 2.575409651244652 0
3782
227
222 1.268170660571357 2.004953917311955 0
3783
228
223 1.254028559084924 0.2715361565992326 0
3784
229
224 0.1470838008660966 1.764796140457005 0
3785
230
225 0.1593462005375125 1.586491973739174 0
3786
231
226 0.4352774855837946 1.792906898802827 0
3787
232
227 0.3507362945887816 1.359618152516893 0
3788
233
228 0.3395649867638134 1.517021777660599 0
3789
234
229 0.3066316055909449 2.068973059576606 0
3790
235
230 0.2737261496659326 1.865906700757079 0
3791
236
231 0.5613504165317504 1.170292178863893 0
3792
237
232 0.4563596603613842 2.097353188592915 0
3793
238
233 2.007728511021513 0.9549956123321318 0
3794
239
234 0.7267482545634016 2.834948776194259 0
3795
240
235 0.3253391498805185 2.624366646627285 0
3796
241
236 2.358627197546403 1.010584251933394 0
3797
242
237 2.336460474883364 1.472537901961013 0
3798
243
238 2.479946954196146 1.342373798035041 0
3799
244
239 2.407160364816117 1.639950922364161 0
3800
245
240 0.7774898959576743 2.997773012818021 0
3801
246
241 0.429755363134886 2.968536834507562 0
3802
247
242 0.841700969612011 0.5402764443878689 0
3803
248
243 1.074292031581404 2.990127302842237 0
3804
249
244 0.1649770749113588 2.37846011276108 0
3805
250
245 0.4294893176197124 1.025176751379064 0
3806
251
246 2.769865645042701 0.8425557550515876 0
3807
252
247 1.19722712608645 1.365215187482435 0
3808
253
248 1.41888751620493 1.178776675558878 0
3809
254
249 1.01600305704664 1.61892336628441 0
3810
255
250 1.023299239274799 1.099376148391514 0
3811
256
251 0.7094616098701299 1.146355733472533 0
3812
257
252 0.8375570011288848 1.30017994094187 0
3813
258
253 1.40046851928461 0.8850515394496341 0
3814
259
254 0.7574604210817141 1.75807980481268 0
3815
260
255 0.5369151619932577 1.450604254046114 0
3816
261
256 1.084158188555139 0.2569001986084773 0
3817
262
257 0.5776611781110469 2.531623442938306 0
3818
263
258 1.102251000331462 2.14688264034314 0
3819
264
259 2.56613476204501 2.178438629672458 0
3820
265
260 2.551233083018216 1.754393000862182 0
3821
266
261 2.718645215727205 1.713386235550407 0
3822
267
262 1.828158880816279 1.343316234877889 0
3823
268
263 1.499290546069746 1.549275088469224 0
3824
269
264 1.520749295326505 1.008432435162344 0
3825
270
265 1.549142593648994 0.8185190690796893 0
3826
271
266 1.708851074954128 1.030162227796758 0
3827
272
267 2.451621018712137 2.26069077209076 0
3828
273
268 0.1444235686033272 0.2104049248412986 0
3829
274
269 0.2820022552685038 0.3230138617005892 0
3830
275
270 0.1839590994501898 0.4510570537117202 0
3831
276
271 0.3425521250128601 0.6690147704806004 0
3832
277
272 2.287412426038745 2.819484384956406 0
3833
278
273 2.724785911251578 2.989319848466447 0
3834
279
274 2.824932499814135 2.789116392863983 0
3835
280
275 2.565990396946163 2.764732077967902 0
3836
281
276 3.009808231531253 1.886577471600333 0
3837
282
277 3.056007070619947 2.897572887613173 0
3838
283
278 2.917628401610219 1.188558014457213 0
3839
284
279 2.122722288949147 0.302574109213219 0
3840
285
280 2.763165731824999 1.268594204164478 0
3841
286
281 1.958000548259156 2.961262599328205 0
3842
287
282 0.1553356301265957 2.548246547728349 0
3843
288
283 1.777702181350695 2.330116968837266 0
3844
289
284 1.680292733716306 2.472177953840805 0
3845
290
285 1.851499212068124 2.52467123334242 0
3846
291
286 1.20814977253769 2.578169806855605 0
3847
292
287 1.392352206339733 1.876937404819547 0
3848
293
288 0.499890186166176 1.292292622650896 0
3849
294
289 0.7542572053621431 1.451695067728285 0
3850
295
290 0.6546138884760674 1.305862776640875 0
3851
296
291 0.5548962248034992 2.193629241092517 0
3852
297
292 0.6027882673202325 2.017601331865698 0
3853
298
293 0.4910627065190025 2.662288043412755 0
3854
299
294 2.273005319549937 1.747779133720325 0
3855
300
295 0.5530824535235275 2.815412226353714 0
3856
301
296 0.7620380704874339 0.2589209365668249 0
3857
302
297 0.7364730363668714 0.6900940035646084 0
3858
303
298 0.9246022946897132 3.001030182296128 0
3859
304
299 0.3200299973484256 2.455228633808153 0
3860
305
300 0.9290423920253242 1.726431709340246 0
3861
306
301 1.150718936248786 1.568604424763654 0
3862
307
302 0.8914629311843184 1.131937106054937 0
3863
308
303 1.154664763582443 1.170138034437028 0
3864
309
304 0.854665300346682 1.597275241603588 0
3865
310
305 1.410721445098315 0.5629256465494171 0
3866
311
306 0.6104694414399311 2.411542658017817 0
3867
312
307 2.466514684951866 1.889882086334573 0
3868
313
308 1.93806535792229 1.283000199568838 0
3869
314
309 1.987488134914787 1.429230321501737 0
3870
315
310 1.568284230829065 1.416113341296378 0
3871
316
311 1.855891735188512 1.057798544484445 0
3872
317
312 1.619632038141906 1.14129891089492 0
3873
318
313 2.565965918470776 2.296260696734317 0
3874
319
314 2.321016000401902 2.288575618181463 0
3875
320
315 2.663087866447886 1.411312838732698 0
3876
321
316 2.532811941169788 0.3710875147649104 0
3877
322
317 2.895052860809531 3.048501750156999 0
3878
323
318 2.856541782263733 1.791395735065057 0
3879
324
319 2.617749083627754 0.8845501063299834 0
3880
325
320 2.697219297628306 0.4193152077994686 0
3881
326
321 0.2556832434807254 3.021145788653913 0
3882
327
322 1.296685312196859 2.331677341888401 0
3883
328
323 1.06237936441321 2.31677290980887 0
3884
329
324 1.500729908018818 2.77335675278342 0
3885
330
325 1.391049094471233 0.7417524850937569 0
3886
331
326 2.232065711744161 1.356021105282206 0
3887
332
327 1.674107597316582 0.8822445806230195 0
3888
333
328 1.54675580061779 0.6830568032054767 0
3889
334
329 0.1757174282502098 0.7559230239156856 0
3890
335
330 0.15932529930038 0.9227960540648468 0
3891
336
331 2.862262164284501 1.038203126666982 0
3892
337
332 2.056605044693401 0.135847050210787 0
3893
338
333 1.16067571736549 0.1156698764675235 0
3894
339
334 1.292857304692313 2.181632517433974 0
3895
340
335 1.57977957601495 1.692703976734909 0
3896
341
336 0.8534573824871737 2.687797426763937 0
3897
342
337 1.002271508751479 2.825548258824476 0
3898
343
338 0.3604143317251066 1.173606927815346 0
3899
344
339 0.5992769246286508 1.860392519688679 0
3900
345
340 0.4493439884107928 1.950759679505188 0
3901
346
341 2.203956420728705 1.644102800232776 0
3902
347
342 0.6067093771829247 0.1585828075188501 0
3903
348
343 0.772738374189719 0.1182309707316591 0
3904
349
344 0.9790361305650853 1.454463681016988 0
3905
350
345 1.018552862012585 1.259424765306716 0
3906
351
346 1.342362599820749 0.4050057887797586 0
3907
352
347 2.395551665363852 1.787227490965682 0
3908
353
348 1.708996837581472 1.373603827058711 0
3909
354
349 2.383805587455049 2.142898858986513 0
3910
355
350 0.8472755822110353 0.8317564635518965 0
3911
356
351 3.002774439871097 2.082169293500771 0
3912
357
352 2.853641435177462 2.192865411227553 0
3913
358
353 2.614364464430576 1.59544094010444 0
3914
359
354 2.887648503285594 0.4213835509232221 0
3915
360
355 2.238690472772636 0.1694610244221905 0
3916
361
356 1.34488922768415 2.842534551691272 0
3917
362
357 2.203849954645166 2.321222635808788 0
3918
363
358 2.326131458112737 0.3696924785243472 0
3919
364
359 0.681666871544073 1.60125037665294 0
3920
365
360 1.110042668600366 0.6632433477829441 0
3921
366
361 2.995988225231864 1.543180837001202 0
3922
367
362 1.357893380725764 0.1296425574983743 0
3923
368
363 1.890369584957931 2.6503947001601 0
3924
369
364 1.752366164911217 2.623376405020105 0
3925
370
365 1.958274883400448 2.448185981660364 0
3926
371
366 2.029390110009729 0.5919371397059507 0
3927
372
367 1.246647941414314 2.986884284685245 0
3928
373
368 1.150594164092547 2.436245361905836 0
3929
374
369 2.345167919667988 1.315368777794578 0
3930
375
370 2.096965906954571 1.311378914074521 0
3931
376
371 0.2903550052120112 1.015410723538438 0
3932
377
372 3.002292328422895 1.021277073295557 0
3933
378
373 1.047673625418839 2.622397830087023 0
3934
379
374 1.176985830813522 2.791365555264978 0
3935
380
375 0.8728091635283713 2.877755689978193 0
3936
381
376 0.4816092231819856 1.632254565845142 0
3937
382
377 0.4292873395726138 0.3956349384765626 0
3938
383
378 0.7513490314779386 0.4234232354093667 0
3939
384
379 0.6341042373478326 0.568270478040689 0
3940
385
380 1.22234407148227 0.8300171394864007 0
3941
386
381 2.857623085556361 1.988197178704509 0
3942
387
382 2.507158237500938 1.472652862906713 0
3943
388
383 2.409527621986015 0.1859400598560078 0
3944
389
384 1.764698773955381 2.965686815609694 0
3945
390
385 0.8997417750110508 2.54443163563051 0
3946
391
386 1.416793031411793 2.979766293462839 0
3947
392
387 1.525513981616164 0.1449723232592194 0
3948
393
388 0.6733841155748304 2.649380725568614 0
3949
394
389 1.031423079732368 0.5247476212108881 0
3950
395
390 0.9314933951417547 0.6653814442277431 0
3951
396
391 1.266624125517718 0.6623633819621698 0
3952
397
392 2.817746465499031 1.601336967556995 0
3953
398
393 1.614530770409514 0.2710620393770625 0
3954
399
394 1.71636688628676 0.1396713013247922 0
3955
400
395 1.282708898714251 1.7174725730963 0
3956
401
396 1.447212987524484 1.711750206278115 0
3957
402
397 0.5767722250337984 1.731728970816225 0
3958
403
398 0.5913272474677719 0.3732866975702252 0
3959
404
399 2.989585509747945 0.3005544744224926 0
3960
405
400 1.783426395629271 0.7390310606086292 0
3961
406
401 0.1939614271025184 0.6024078510995812 0
3962
407
402 0.9372148501050993 0.9811426585705684 0
3963
408
403 2.864348317610181 1.399970656621234 0
3964
409
404 1.754522340738722 0.4587208708514305 0
3965
410
405 1.150300942140853 0.3994409683546766 0
3966
411
406 1.023248075413499 0.817594015684914 0
3967
412
407 1.242779111659225 1.859539070769294 0
3968
413
408 0.4486014544794452 0.2439224749890189 0
3969
414
409 3.015503673251204 0.4990226676685885 0
3970
415
410 3.019697961286836 0.8619059836157521 0
3971
416
411 3.030260467539848 1.34333347067228 0
3972
417
412 1.996468216947338 0.4138691345250425 0
3973
418
413 1.778567667136346 0.2808503494953919 0
3974
419
414 1.868562216861688 0.5618806000914466 0
3975
420
415 1.663393210490735 0.3902015376645016 0
3976
421
416 1.218439071129093 0.5384544506519755 0
3977
422
417 3.011265819676178 0.6931410673899771 0
3978
423
418 1.942644646849937 0.2488856793463607 0
3979
424
419 1.880570585646719 0.369815449419561 0
3980
425
420 1.537037228489016 1.950578716018899 0
3981
426
421 2.221266302057239 3.045599411745064 0
3982
427
422 1.829599455239953 1.83583038649709 0
3983
428
423 0.1027025713541935 3.046990539861628 0
3984
429
424 1.684157697788862 1.899067369585522 0
3985
430
425 1.729317902377592 0.5907287518243273 0
3986
431
426 0.4493613937975903 0.1108083172592962 0
3987
432
427 2.086150493609806 3.040684625350802 0
3988
433
428 1.879303317923693 0.1245686623238263 0
3989
434
$EndNodes
3990
435
$Elements
3991
436
854
3992
437
1 1 3 9 1 0 4 5
3993
438
2 1 3 9 1 0 5 6
3994
439
3 1 3 9 1 0 6 7
3995
440
4 1 3 9 1 0 7 8
3996
441
5 1 3 9 1 0 8 9
3997
442
6 1 3 9 1 0 9 10
3998
443
7 1 3 9 1 0 10 11
3999
444
8 1 3 9 1 0 11 12
4000
445
9 1 3 9 1 0 12 13
4001
446
10 1 3 9 1 0 13 14
4002
447
11 1 3 9 1 0 14 15
4003
448
12 1 3 9 1 0 15 16
4004
449
13 1 3 9 1 0 16 17
4005
450
14 1 3 9 1 0 17 18
4006
451
15 1 3 9 1 0 18 19
4007
452
16 1 3 9 1 0 19 20
4008
453
17 1 3 9 1 0 20 21
4009
454
18 1 3 9 1 0 21 22
4010
455
19 1 3 9 1 0 22 3
4011
456
20 1 3 8 2 0 3 23
4012
457
21 1 3 8 2 0 23 24
4013
458
22 1 3 8 2 0 24 25
4014
459
23 1 3 8 2 0 25 26
4015
460
24 1 3 8 2 0 26 27
4016
461
25 1 3 8 2 0 27 28
4017
462
26 1 3 8 2 0 28 29
4018
463
27 1 3 8 2 0 29 30
4019
464
28 1 3 8 2 0 30 31
4020
465
29 1 3 8 2 0 31 32
4021
466
30 1 3 8 2 0 32 33
4022
467
31 1 3 8 2 0 33 34
4023
468
32 1 3 8 2 0 34 35
4024
469
33 1 3 8 2 0 35 36
4025
470
34 1 3 8 2 0 36 37
4026
471
35 1 3 8 2 0 37 38
4027
472
36 1 3 8 2 0 38 39
4028
473
37 1 3 8 2 0 39 40
4029
474
38 1 3 8 2 0 40 2
4030
475
39 1 3 7 3 0 2 41
4031
476
40 1 3 7 3 0 41 42
4032
477
41 1 3 7 3 0 42 43
4033
478
42 1 3 7 3 0 43 44
4034
479
43 1 3 7 3 0 44 45
4035
480
44 1 3 7 3 0 45 46
4036
481
45 1 3 7 3 0 46 47
4037
482
46 1 3 7 3 0 47 48
4038
483
47 1 3 7 3 0 48 49
4039
484
48 1 3 7 3 0 49 50
4040
485
49 1 3 7 3 0 50 51
4041
486
50 1 3 7 3 0 51 52
4042
487
51 1 3 7 3 0 52 53
4043
488
52 1 3 7 3 0 53 54
4044
489
53 1 3 7 3 0 54 55
4045
490
54 1 3 7 3 0 55 56
4046
491
55 1 3 7 3 0 56 57
4047
492
56 1 3 7 3 0 57 58
4048
493
57 1 3 7 3 0 58 1
4049
494
58 1 3 10 4 0 1 59
4050
495
59 1 3 10 4 0 59 60
4051
496
60 1 3 10 4 0 60 61
4052
497
61 1 3 10 4 0 61 62
4053
498
62 1 3 10 4 0 62 63
4054
499
63 1 3 10 4 0 63 64
4055
500
64 1 3 10 4 0 64 65
4056
501
65 1 3 10 4 0 65 66
4057
502
66 1 3 10 4 0 66 67
4058
503
67 1 3 10 4 0 67 68
4059
504
68 1 3 10 4 0 68 69
4060
505
69 1 3 10 4 0 69 70
4061
506
70 1 3 10 4 0 70 71
4062
507
71 1 3 10 4 0 71 72
4063
508
72 1 3 10 4 0 72 73
4064
509
73 1 3 10 4 0 73 74
4065
510
74 1 3 10 4 0 74 75
4066
511
75 1 3 10 4 0 75 76
4067
512
76 1 3 10 4 0 76 4
4068
513
77 2 3 12 6 0 81 88 86
4069
514
78 2 3 12 6 0 92 19 20
4070
515
79 2 3 12 6 0 137 22 3
4071
516
80 2 3 12 6 0 23 137 3
4072
517
81 2 3 12 6 0 53 54 110
4073
518
82 2 3 12 6 0 126 81 86
4074
519
83 2 3 12 6 0 31 129 30
4075
520
84 2 3 12 6 0 18 19 132
4076
521
85 2 3 12 6 0 19 92 132
4077
522
86 2 3 12 6 0 136 91 82
4078
523
87 2 3 12 6 0 69 70 117
4079
524
88 2 3 12 6 0 25 138 24
4080
525
89 2 3 12 6 0 138 25 136
4081
526
90 2 3 12 6 0 116 164 134
4082
527
91 2 3 12 6 0 165 117 70
4083
528
92 2 3 12 6 0 118 166 154
4084
529
93 2 3 12 6 0 168 119 153
4085
530
94 2 3 12 6 0 120 169 150
4086
531
95 2 3 12 6 0 169 120 149
4087
532
96 2 3 12 6 0 169 97 150
4088
533
97 2 3 12 6 0 171 94 153
4089
534
98 2 3 12 6 0 94 171 155
4090
535
99 2 3 12 6 0 90 188 115
4091
536
100 2 3 12 6 0 26 193 91
4092
537
101 2 3 12 6 0 195 121 194
4093
538
102 2 3 12 6 0 200 125 182
4094
539
103 2 3 12 6 0 183 184 77
4095
540
104 2 3 12 6 0 183 77 128
4096
541
105 2 3 12 6 0 190 41 42
4097
542
106 2 3 12 6 0 42 86 88
4098
543
107 2 3 12 6 0 88 190 42
4099
544
108 2 3 12 6 0 159 104 108
4100
545
109 2 3 12 6 0 108 160 159
4101
546
110 2 3 12 6 0 103 160 108
4102
547
111 2 3 12 6 0 85 179 161
4103
548
112 2 3 12 6 0 64 65 127
4104
549
113 2 3 12 6 0 163 90 115
4105
550
114 2 3 12 6 0 111 163 115
4106
551
115 2 3 12 6 0 111 185 163
4107
552
116 2 3 12 6 0 185 131 163
4108
553
117 2 3 12 6 0 205 114 173
4109
554
118 2 3 12 6 0 85 181 179
4110
555
119 2 3 12 6 0 106 179 123
4111
556
120 2 3 12 6 0 181 123 179
4112
557
121 2 3 12 6 0 107 161 179
4113
558
122 2 3 12 6 0 105 180 123
4114
559
123 2 3 12 6 0 180 106 123
4115
560
124 2 3 12 6 0 181 105 123
4116
561
125 2 3 12 6 0 125 79 182
4117
562
126 2 3 12 6 0 77 184 107
4118
563
127 2 3 12 6 0 185 89 131
4119
564
128 2 3 12 6 0 89 185 145
4120
565
129 2 3 12 6 0 185 111 145
4121
566
130 2 3 12 6 0 89 186 131
4122
567
131 2 3 12 6 0 131 186 130
4123
568
132 2 3 12 6 0 187 111 115
4124
569
133 2 3 12 6 0 116 188 82
4125
570
134 2 3 12 6 0 116 187 188
4126
571
135 2 3 12 6 0 187 115 188
4127
572
136 2 3 12 6 0 2 190 40
4128
573
137 2 3 12 6 0 190 2 41
4129
574
138 2 3 12 6 0 191 57 58
4130
575
139 2 3 12 6 0 60 61 192
4131
576
140 2 3 12 6 0 27 193 26
4132
577
141 2 3 12 6 0 75 76 194
4133
578
142 2 3 12 6 0 121 75 194
4134
579
143 2 3 12 6 0 195 189 121
4135
580
144 2 3 12 6 0 198 127 65
4136
581
145 2 3 12 6 0 66 198 65
4137
582
146 2 3 12 6 0 42 43 86
4138
583
147 2 3 12 6 0 202 86 43
4139
584
148 2 3 12 6 0 202 126 86
4140
585
149 2 3 12 6 0 92 203 132
4141
586
150 2 3 12 6 0 163 203 90
4142
587
151 2 3 12 6 0 188 204 82
4143
588
152 2 3 12 6 0 70 71 165
4144
589
153 2 3 12 6 0 74 75 121
4145
590
154 2 3 12 6 0 221 140 207
4146
591
155 2 3 12 6 0 147 221 207
4147
592
156 2 3 12 6 0 117 224 69
4148
593
157 2 3 12 6 0 225 124 67
4149
594
158 2 3 12 6 0 93 229 165
4150
595
159 2 3 12 6 0 233 156 180
4151
596
160 2 3 12 6 0 105 233 180
4152
597
161 2 3 12 6 0 158 236 180
4153
598
162 2 3 12 6 0 236 106 180
4154
599
163 2 3 12 6 0 122 246 104
4155
600
164 2 3 12 6 0 162 247 177
4156
601
165 2 3 12 6 0 247 78 177
4157
602
166 2 3 12 6 0 248 162 177
4158
603
167 2 3 12 6 0 162 248 99
4159
604
168 2 3 12 6 0 125 251 79
4160
605
169 2 3 12 6 0 173 259 205
4161
606
170 2 3 12 6 0 259 134 205
4162
607
171 2 3 12 6 0 192 269 268
4163
608
172 2 3 12 6 0 269 191 268
4164
609
173 2 3 12 6 0 131 272 163
4165
610
174 2 3 12 6 0 272 203 163
4166
611
175 2 3 12 6 0 273 92 20
4167
612
176 2 3 12 6 0 92 273 133
4168
613
177 2 3 12 6 0 90 275 188
4169
614
178 2 3 12 6 0 275 204 188
4170
615
179 2 3 12 6 0 44 202 43
4171
616
180 2 3 12 6 0 158 180 156
4172
617
181 2 3 12 6 0 184 161 107
4173
618
182 2 3 12 6 0 85 196 181
4174
619
183 2 3 12 6 0 199 200 182
4175
620
184 2 3 12 6 0 164 82 91
4176
621
185 2 3 12 6 0 164 116 82
4177
622
186 2 3 12 6 0 7 8 172
4178
623
187 2 3 12 6 0 172 241 7
4179
624
188 2 3 12 6 0 210 113 141
4180
625
189 2 3 12 6 0 219 113 210
4181
626
190 2 3 12 6 0 144 219 210
4182
627
191 2 3 12 6 0 113 211 141
4183
628
192 2 3 12 6 0 208 98 216
4184
629
193 2 3 12 6 0 207 140 216
4185
630
194 2 3 12 6 0 140 208 216
4186
631
195 2 3 12 6 0 218 96 144
4187
632
196 2 3 12 6 0 218 112 145
4188
633
197 2 3 12 6 0 96 219 144
4189
634
198 2 3 12 6 0 113 219 146
4190
635
199 2 3 12 6 0 140 221 148
4191
636
200 2 3 12 6 0 230 100 224
4192
637
201 2 3 12 6 0 100 228 225
4193
638
202 2 3 12 6 0 100 230 226
4194
639
203 2 3 12 6 0 198 124 227
4195
640
204 2 3 12 6 0 229 93 154
4196
641
205 2 3 12 6 0 154 232 229
4197
642
206 2 3 12 6 0 251 125 231
4198
643
207 2 3 12 6 0 233 101 156
4199
644
208 2 3 12 6 0 172 240 234
4200
645
209 2 3 12 6 0 189 235 121
4201
646
210 2 3 12 6 0 103 236 158
4202
647
211 2 3 12 6 0 236 103 108
4203
648
212 2 3 12 6 0 236 84 106
4204
649
213 2 3 12 6 0 84 236 108
4205
650
214 2 3 12 6 0 103 238 160
4206
651
215 2 3 12 6 0 240 8 9
4207
652
216 2 3 12 6 0 8 240 172
4208
653
217 2 3 12 6 0 244 93 72
4209
654
218 2 3 12 6 0 125 200 245
4210
655
219 2 3 12 6 0 231 125 245
4211
656
220 2 3 12 6 0 200 201 245
4212
657
221 2 3 12 6 0 77 246 128
4213
658
222 2 3 12 6 0 246 122 128
4214
659
223 2 3 12 6 0 99 264 253
4215
660
224 2 3 12 6 0 87 256 110
4216
661
225 2 3 12 6 0 167 257 118
4217
662
226 2 3 12 6 0 258 94 150
4218
663
227 2 3 12 6 0 114 261 260
4219
664
228 2 3 12 6 0 30 276 29
4220
665
229 2 3 12 6 0 268 191 58
4221
666
230 2 3 12 6 0 192 270 269
4222
667
231 2 3 12 6 0 270 197 269
4223
668
232 2 3 12 6 0 59 60 268
4224
669
233 2 3 12 6 0 60 192 268
4225
670
234 2 3 12 6 0 67 68 225
4226
671
235 2 3 12 6 0 61 270 192
4227
672
236 2 3 12 6 0 66 67 124
4228
673
237 2 3 12 6 0 198 66 124
4229
674
238 2 3 12 6 0 197 271 80
4230
675
239 2 3 12 6 0 271 199 80
4231
676
240 2 3 12 6 0 272 130 132
4232
677
241 2 3 12 6 0 130 272 131
4233
678
242 2 3 12 6 0 203 272 132
4234
679
243 2 3 12 6 0 204 133 274
4235
680
244 2 3 12 6 0 274 136 82
4236
681
245 2 3 12 6 0 204 274 82
4237
682
246 2 3 12 6 0 203 275 90
4238
683
247 2 3 12 6 0 275 203 92
4239
684
248 2 3 12 6 0 275 133 204
4240
685
249 2 3 12 6 0 133 275 92
4241
686
250 2 3 12 6 0 129 276 30
4242
687
251 2 3 12 6 0 24 277 23
4243
688
252 2 3 12 6 0 74 282 73
4244
689
253 2 3 12 6 0 283 98 208
4245
690
254 2 3 12 6 0 140 284 208
4246
691
255 2 3 12 6 0 166 291 154
4247
692
256 2 3 12 6 0 291 232 154
4248
693
257 2 3 12 6 0 235 293 167
4249
694
258 2 3 12 6 0 298 9 10
4250
695
259 2 3 12 6 0 243 298 10
4251
696
260 2 3 12 6 0 299 118 154
4252
697
261 2 3 12 6 0 300 153 119
4253
698
262 2 3 12 6 0 249 300 119
4254
699
263 2 3 12 6 0 249 304 300
4255
700
264 2 3 12 6 0 304 254 300
4256
701
265 2 3 12 6 0 306 166 118
4257
702
266 2 3 12 6 0 257 306 118
4258
703
267 2 3 12 6 0 308 101 175
4259
704
268 2 3 12 6 0 262 308 175
4260
705
269 2 3 12 6 0 310 102 177
4261
706
270 2 3 12 6 0 101 311 175
4262
707
271 2 3 12 6 0 311 266 175
4263
708
272 2 3 12 6 0 312 102 175
4264
709
273 2 3 12 6 0 266 312 175
4265
710
274 2 3 12 6 0 313 134 259
4266
711
275 2 3 12 6 0 267 313 259
4267
712
276 2 3 12 6 0 317 21 22
4268
713
277 2 3 12 6 0 138 277 24
4269
714
278 2 3 12 6 0 298 240 9
4270
715
279 2 3 12 6 0 137 317 22
4271
716
280 2 3 12 6 0 137 135 317
4272
717
281 2 3 12 6 0 277 137 23
4273
718
282 2 3 12 6 0 135 137 277
4274
719
283 2 3 12 6 0 241 6 7
4275
720
284 2 3 12 6 0 135 138 274
4276
721
285 2 3 12 6 0 135 277 138
4277
722
286 2 3 12 6 0 136 274 138
4278
723
287 2 3 12 6 0 195 241 189
4279
724
288 2 3 12 6 0 293 235 189
4280
725
289 2 3 12 6 0 241 295 189
4281
726
290 2 3 12 6 0 71 93 165
4282
727
291 2 3 12 6 0 93 71 72
4283
728
292 2 3 12 6 0 124 225 228
4284
729
293 2 3 12 6 0 227 124 228
4285
730
294 2 3 12 6 0 248 264 99
4286
731
295 2 3 12 6 0 177 102 248
4287
732
296 2 3 12 6 0 288 227 255
4288
733
297 2 3 12 6 0 255 227 228
4289
734
298 2 3 12 6 0 264 265 253
4290
735
299 2 3 12 6 0 273 20 21
4291
736
300 2 3 12 6 0 21 317 273
4292
737
301 2 3 12 6 0 274 133 273
4293
738
302 2 3 12 6 0 273 135 274
4294
739
303 2 3 12 6 0 317 135 273
4295
740
304 2 3 12 6 0 159 280 278
4296
741
305 2 3 12 6 0 160 280 159
4297
742
306 2 3 12 6 0 160 315 280
4298
743
307 2 3 12 6 0 121 282 74
4299
744
308 2 3 12 6 0 282 121 235
4300
745
309 2 3 12 6 0 206 98 283
4301
746
310 2 3 12 6 0 112 206 283
4302
747
311 2 3 12 6 0 284 140 148
4303
748
312 2 3 12 6 0 222 143 287
4304
749
313 2 3 12 6 0 252 290 289
4305
750
314 2 3 12 6 0 291 120 155
4306
751
315 2 3 12 6 0 120 291 166
4307
752
316 2 3 12 6 0 155 292 291
4308
753
317 2 3 12 6 0 292 232 291
4309
754
318 2 3 12 6 0 171 292 155
4310
755
319 2 3 12 6 0 294 113 146
4311
756
320 2 3 12 6 0 172 295 241
4312
757
321 2 3 12 6 0 234 295 172
4313
758
322 2 3 12 6 0 296 87 110
4314
759
323 2 3 12 6 0 297 199 182
4315
760
324 2 3 12 6 0 167 118 299
4316
761
325 2 3 12 6 0 299 235 167
4317
762
326 2 3 12 6 0 282 235 299
4318
763
327 2 3 12 6 0 299 244 282
4319
764
328 2 3 12 6 0 301 249 119
4320
765
329 2 3 12 6 0 109 303 99
4321
766
330 2 3 12 6 0 303 109 250
4322
767
331 2 3 12 6 0 303 162 99
4323
768
332 2 3 12 6 0 176 305 152
4324
769
333 2 3 12 6 0 120 306 149
4325
770
334 2 3 12 6 0 306 120 166
4326
771
335 2 3 12 6 0 306 257 149
4327
772
336 2 3 12 6 0 173 307 217
4328
773
337 2 3 12 6 0 217 307 146
4329
774
338 2 3 12 6 0 309 212 211
4330
775
339 2 3 12 6 0 262 212 309
4331
776
340 2 3 12 6 0 157 310 263
4332
777
341 2 3 12 6 0 311 101 233
4333
778
342 2 3 12 6 0 174 311 233
4334
779
343 2 3 12 6 0 311 174 266
4335
780
344 2 3 12 6 0 102 312 248
4336
781
345 2 3 12 6 0 264 312 266
4337
782
346 2 3 12 6 0 312 264 248
4338
783
347 2 3 12 6 0 134 313 116
4339
784
348 2 3 12 6 0 313 187 116
4340
785
349 2 3 12 6 0 187 313 267
4341
786
350 2 3 12 6 0 267 314 187
4342
787
351 2 3 12 6 0 314 111 187
4343
788
352 2 3 12 6 0 315 160 238
4344
789
353 2 3 12 6 0 316 126 202
4345
790
354 2 3 12 6 0 316 161 184
4346
791
355 2 3 12 6 0 318 261 114
4347
792
356 2 3 12 6 0 129 318 276
4348
793
357 2 3 12 6 0 319 84 108
4349
794
358 2 3 12 6 0 104 319 108
4350
795
359 2 3 12 6 0 320 316 184
4351
796
360 2 3 12 6 0 316 320 126
4352
797
361 2 3 12 6 0 321 6 241
4353
798
362 2 3 12 6 0 195 321 241
4354
799
363 2 3 12 6 0 323 170 258
4355
800
364 2 3 12 6 0 150 323 258
4356
801
365 2 3 12 6 0 148 324 139
4357
802
366 2 3 12 6 0 325 253 265
4358
803
367 2 3 12 6 0 327 265 264
4359
804
368 2 3 12 6 0 266 327 264
4360
805
369 2 3 12 6 0 152 328 83
4361
806
370 2 3 12 6 0 328 265 83
4362
807
371 2 3 12 6 0 329 62 63
4363
808
372 2 3 12 6 0 64 330 63
4364
809
373 2 3 12 6 0 331 104 159
4365
810
374 2 3 12 6 0 278 331 159
4366
811
375 2 3 12 6 0 334 216 143
4367
812
376 2 3 12 6 0 215 335 95
4368
813
377 2 3 12 6 0 335 215 157
4369
814
378 2 3 12 6 0 245 338 231
4370
815
379 2 3 12 6 0 338 288 231
4371
816
380 2 3 12 6 0 254 339 171
4372
817
381 2 3 12 6 0 339 292 171
4373
818
382 2 3 12 6 0 341 211 113
4374
819
383 2 3 12 6 0 294 341 113
4375
820
384 2 3 12 6 0 54 343 110
4376
821
385 2 3 12 6 0 343 296 110
4377
822
386 2 3 12 6 0 347 294 146
4378
823
387 2 3 12 6 0 307 347 146
4379
824
388 2 3 12 6 0 96 349 217
4380
825
389 2 3 12 6 0 349 96 314
4381
826
390 2 3 12 6 0 353 260 261
4382
827
391 2 3 12 6 0 247 301 78
4383
828
392 2 3 12 6 0 162 303 247
4384
829
393 2 3 12 6 0 303 345 247
4385
830
394 2 3 12 6 0 302 79 251
4386
831
395 2 3 12 6 0 290 252 251
4387
832
396 2 3 12 6 0 302 251 252
4388
833
397 2 3 12 6 0 290 251 231
4389
834
398 2 3 12 6 0 281 130 186
4390
835
399 2 3 12 6 0 209 281 186
4391
836
400 2 3 12 6 0 52 333 51
4392
837
401 2 3 12 6 0 286 147 207
4393
838
402 2 3 12 6 0 290 288 255
4394
839
403 2 3 12 6 0 255 289 290
4395
840
404 2 3 12 6 0 290 231 288
4396
841
405 2 3 12 6 0 257 167 293
4397
842
406 2 3 12 6 0 295 293 189
4398
843
407 2 3 12 6 0 73 244 72
4399
844
408 2 3 12 6 0 244 73 282
4400
845
409 2 3 12 6 0 243 10 11
4401
846
410 2 3 12 6 0 262 309 308
4402
847
411 2 3 12 6 0 246 319 104
4403
848
412 2 3 12 6 0 319 107 84
4404
849
413 2 3 12 6 0 319 77 107
4405
850
414 2 3 12 6 0 246 77 319
4406
851
415 2 3 12 6 0 6 321 5
4407
852
416 2 3 12 6 0 194 321 195
4408
853
417 2 3 12 6 0 216 322 207
4409
854
418 2 3 12 6 0 216 334 322
4410
855
419 2 3 12 6 0 97 323 150
4411
856
420 2 3 12 6 0 326 213 237
4412
857
421 2 3 12 6 0 174 327 266
4413
858
422 2 3 12 6 0 327 83 265
4414
859
423 2 3 12 6 0 127 330 64
4415
860
424 2 3 12 6 0 331 122 104
4416
861
425 2 3 12 6 0 52 53 333
4417
862
426 2 3 12 6 0 333 256 223
4418
863
427 2 3 12 6 0 298 243 337
4419
864
428 2 3 12 6 0 338 227 288
4420
865
429 2 3 12 6 0 227 338 198
4421
866
430 2 3 12 6 0 341 213 211
4422
867
431 2 3 12 6 0 213 341 237
4423
868
432 2 3 12 6 0 239 341 294
4424
869
433 2 3 12 6 0 341 239 237
4425
870
434 2 3 12 6 0 296 343 342
4426
871
435 2 3 12 6 0 344 249 301
4427
872
436 2 3 12 6 0 249 344 304
4428
873
437 2 3 12 6 0 151 346 176
4429
874
438 2 3 12 6 0 346 151 223
4430
875
439 2 3 12 6 0 346 305 176
4431
876
440 2 3 12 6 0 239 347 260
4432
877
441 2 3 12 6 0 347 239 294
4433
878
442 2 3 12 6 0 347 307 260
4434
879
443 2 3 12 6 0 102 348 175
4435
880
444 2 3 12 6 0 348 102 310
4436
881
445 2 3 12 6 0 348 262 175
4437
882
446 2 3 12 6 0 348 157 212
4438
883
447 2 3 12 6 0 157 348 310
4439
884
448 2 3 12 6 0 262 348 212
4440
885
449 2 3 12 6 0 79 350 182
4441
886
450 2 3 12 6 0 350 297 182
4442
887
451 2 3 12 6 0 351 29 276
4443
888
452 2 3 12 6 0 351 28 29
4444
889
453 2 3 12 6 0 134 352 205
4445
890
454 2 3 12 6 0 352 134 164
4446
891
455 2 3 12 6 0 353 239 260
4447
892
456 2 3 12 6 0 46 355 45
4448
893
457 2 3 12 6 0 332 279 355
4449
894
458 2 3 12 6 0 357 145 111
4450
895
459 2 3 12 6 0 314 357 111
4451
896
460 2 3 12 6 0 359 304 289
4452
897
461 2 3 12 6 0 304 359 254
4453
898
462 2 3 12 6 0 364 285 284
4454
899
463 2 3 12 6 0 148 364 284
4455
900
464 2 3 12 6 0 145 365 89
4456
901
465 2 3 12 6 0 365 285 89
4457
902
466 2 3 12 6 0 367 11 12
4458
903
467 2 3 12 6 0 238 369 237
4459
904
468 2 3 12 6 0 369 326 237
4460
905
469 2 3 12 6 0 330 371 201
4461
906
470 2 3 12 6 0 371 245 201
4462
907
471 2 3 12 6 0 375 240 298
4463
908
472 2 3 12 6 0 337 375 298
4464
909
473 2 3 12 6 0 228 376 255
4465
910
474 2 3 12 6 0 199 379 80
4466
911
475 2 3 12 6 0 379 199 297
4467
912
476 2 3 12 6 0 238 382 315
4468
913
477 2 3 12 6 0 382 353 315
4469
914
478 2 3 12 6 0 202 383 316
4470
915
479 2 3 12 6 0 212 141 211
4471
916
480 2 3 12 6 0 13 14 220
4472
917
481 2 3 12 6 0 284 283 208
4473
918
482 2 3 12 6 0 284 285 283
4474
919
483 2 3 12 6 0 211 213 309
4475
920
484 2 3 12 6 0 370 309 213
4476
921
485 2 3 12 6 0 326 370 213
4477
922
486 2 3 12 6 0 81 126 320
4478
923
487 2 3 12 6 0 354 81 320
4479
924
488 2 3 12 6 0 183 320 184
4480
925
489 2 3 12 6 0 354 320 183
4481
926
490 2 3 12 6 0 305 328 152
4482
927
491 2 3 12 6 0 328 325 265
4483
928
492 2 3 12 6 0 305 325 328
4484
929
493 2 3 12 6 0 329 201 271
4485
930
494 2 3 12 6 0 329 330 201
4486
931
495 2 3 12 6 0 330 329 63
4487
932
496 2 3 12 6 0 51 362 50
4488
933
497 2 3 12 6 0 143 222 334
4489
934
498 2 3 12 6 0 170 334 258
4490
935
499 2 3 12 6 0 170 322 334
4491
936
500 2 3 12 6 0 335 157 263
4492
937
501 2 3 12 6 0 226 340 339
4493
938
502 2 3 12 6 0 340 292 339
4494
939
503 2 3 12 6 0 232 340 229
4495
940
504 2 3 12 6 0 340 232 292
4496
941
505 2 3 12 6 0 173 349 259
4497
942
506 2 3 12 6 0 349 173 217
4498
943
507 2 3 12 6 0 349 267 259
4499
944
508 2 3 12 6 0 314 267 349
4500
945
509 2 3 12 6 0 345 302 252
4501
946
510 2 3 12 6 0 345 250 302
4502
947
511 2 3 12 6 0 250 345 303
4503
948
512 2 3 12 6 0 357 218 145
4504
949
513 2 3 12 6 0 218 357 96
4505
950
514 2 3 12 6 0 357 314 96
4506
951
515 2 3 12 6 0 358 85 161
4507
952
516 2 3 12 6 0 316 358 161
4508
953
517 2 3 12 6 0 383 358 316
4509
954
518 2 3 12 6 0 31 361 129
4510
955
519 2 3 12 6 0 361 31 32
4511
956
520 2 3 12 6 0 362 51 333
4512
957
521 2 3 12 6 0 223 362 333
4513
958
522 2 3 12 6 0 362 223 151
4514
959
523 2 3 12 6 0 186 89 363
4515
960
524 2 3 12 6 0 89 285 363
4516
961
525 2 3 12 6 0 363 209 186
4517
962
526 2 3 12 6 0 139 364 148
4518
963
527 2 3 12 6 0 364 139 209
4519
964
528 2 3 12 6 0 209 363 364
4520
965
529 2 3 12 6 0 363 285 364
4521
966
530 2 3 12 6 0 112 365 145
4522
967
531 2 3 12 6 0 365 112 283
4523
968
532 2 3 12 6 0 285 365 283
4524
969
533 2 3 12 6 0 196 366 181
4525
970
534 2 3 12 6 0 11 367 243
4526
971
535 2 3 12 6 0 368 170 323
4527
972
536 2 3 12 6 0 368 322 170
4528
973
537 2 3 12 6 0 368 286 207
4529
974
538 2 3 12 6 0 322 368 207
4530
975
539 2 3 12 6 0 97 368 323
4531
976
540 2 3 12 6 0 369 103 158
4532
977
541 2 3 12 6 0 103 369 238
4533
978
542 2 3 12 6 0 326 369 158
4534
979
543 2 3 12 6 0 370 308 309
4535
980
544 2 3 12 6 0 101 308 370
4536
981
545 2 3 12 6 0 326 158 370
4537
982
546 2 3 12 6 0 156 370 158
4538
983
547 2 3 12 6 0 101 370 156
4539
984
548 2 3 12 6 0 127 371 330
4540
985
549 2 3 12 6 0 245 371 338
4541
986
550 2 3 12 6 0 372 122 331
4542
987
551 2 3 12 6 0 35 372 34
4543
988
552 2 3 12 6 0 278 372 331
4544
989
553 2 3 12 6 0 372 278 34
4545
990
554 2 3 12 6 0 368 97 373
4546
991
555 2 3 12 6 0 336 337 373
4547
992
556 2 3 12 6 0 286 368 373
4548
993
557 2 3 12 6 0 375 234 240
4549
994
558 2 3 12 6 0 234 375 336
4550
995
559 2 3 12 6 0 375 337 336
4551
996
560 2 3 12 6 0 100 376 228
4552
997
561 2 3 12 6 0 100 226 376
4553
998
562 2 3 12 6 0 377 269 197
4554
999
563 2 3 12 6 0 80 377 197
4555
1000
564 2 3 12 6 0 54 55 343
4556
1001
565 2 3 12 6 0 55 342 343
4557
1002
566 2 3 12 6 0 296 378 87
4558
1003
567 2 3 12 6 0 99 380 109
4559
1004
568 2 3 12 6 0 99 253 380
4560
1005
569 2 3 12 6 0 380 253 325
4561
1006
570 2 3 12 6 0 318 381 276
4562
1007
571 2 3 12 6 0 381 318 114
4563
1008
572 2 3 12 6 0 381 351 276
4564
1009
573 2 3 12 6 0 205 381 114
4565
1010
574 2 3 12 6 0 27 28 193
4566
1011
575 2 3 12 6 0 382 237 239
4567
1012
576 2 3 12 6 0 237 382 238
4568
1013
577 2 3 12 6 0 353 382 239
4569
1014
578 2 3 12 6 0 44 383 202
4570
1015
579 2 3 12 6 0 46 47 332
4571
1016
580 2 3 12 6 0 355 46 332
4572
1017
581 2 3 12 6 0 384 139 220
4573
1018
582 2 3 12 6 0 385 336 373
4574
1019
583 2 3 12 6 0 336 385 149
4575
1020
584 2 3 12 6 0 97 385 373
4576
1021
585 2 3 12 6 0 386 13 220
4577
1022
586 2 3 12 6 0 387 50 362
4578
1023
587 2 3 12 6 0 151 387 362
4579
1024
588 2 3 12 6 0 388 149 257
4580
1025
589 2 3 12 6 0 149 388 336
4581
1026
590 2 3 12 6 0 390 360 389
4582
1027
591 2 3 12 6 0 242 390 389
4583
1028
592 2 3 12 6 0 391 360 380
4584
1029
593 2 3 12 6 0 325 391 380
4585
1030
594 2 3 12 6 0 392 261 318
4586
1031
595 2 3 12 6 0 393 151 176
4587
1032
596 2 3 12 6 0 151 393 387
4588
1033
597 2 3 12 6 0 301 395 78
4589
1034
598 2 3 12 6 0 95 396 287
4590
1035
599 2 3 12 6 0 396 95 335
4591
1036
600 2 3 12 6 0 339 397 226
4592
1037
601 2 3 12 6 0 397 376 226
4593
1038
602 2 3 12 6 0 398 377 80
4594
1039
603 2 3 12 6 0 169 385 97
4595
1040
604 2 3 12 6 0 169 149 385
4596
1041
605 2 3 12 6 0 117 165 230
4597
1042
606 2 3 12 6 0 224 117 230
4598
1043
607 2 3 12 6 0 229 230 165
4599
1044
608 2 3 12 6 0 229 340 230
4600
1045
609 2 3 12 6 0 340 226 230
4601
1046
610 2 3 12 6 0 344 345 252
4602
1047
611 2 3 12 6 0 289 255 359
4603
1048
612 2 3 12 6 0 376 359 255
4604
1049
613 2 3 12 6 0 374 337 243
4605
1050
614 2 3 12 6 0 367 374 243
4606
1051
615 2 3 12 6 0 286 374 147
4607
1052
616 2 3 12 6 0 286 373 374
4608
1053
617 2 3 12 6 0 373 337 374
4609
1054
618 2 3 12 6 0 139 384 209
4610
1055
619 2 3 12 6 0 384 281 209
4611
1056
620 2 3 12 6 0 13 386 12
4612
1057
621 2 3 12 6 0 350 390 297
4613
1058
622 2 3 12 6 0 390 242 297
4614
1059
623 2 3 12 6 0 305 391 325
4615
1060
624 2 3 12 6 0 394 48 49
4616
1061
625 2 3 12 6 0 119 395 301
4617
1062
626 2 3 12 6 0 254 397 339
4618
1063
627 2 3 12 6 0 397 254 359
4619
1064
628 2 3 12 6 0 376 397 359
4620
1065
629 2 3 12 6 0 398 296 342
4621
1066
630 2 3 12 6 0 378 296 398
4622
1067
631 2 3 12 6 0 399 39 40
4623
1068
632 2 3 12 6 0 399 88 81
4624
1069
633 2 3 12 6 0 354 399 81
4625
1070
634 2 3 12 6 0 190 399 40
4626
1071
635 2 3 12 6 0 88 399 190
4627
1072
636 2 3 12 6 0 174 400 327
4628
1073
637 2 3 12 6 0 400 174 178
4629
1074
638 2 3 12 6 0 400 83 327
4630
1075
639 2 3 12 6 0 401 271 197
4631
1076
640 2 3 12 6 0 271 401 329
4632
1077
641 2 3 12 6 0 402 302 250
4633
1078
642 2 3 12 6 0 406 390 350
4634
1079
643 2 3 12 6 0 407 222 287
4635
1080
644 2 3 12 6 0 395 407 287
4636
1081
645 2 3 12 6 0 408 398 342
4637
1082
646 2 3 12 6 0 354 409 399
4638
1083
647 2 3 12 6 0 68 69 224
4639
1084
648 2 3 12 6 0 68 224 225
4640
1085
649 2 3 12 6 0 100 225 224
4641
1086
650 2 3 12 6 0 378 242 87
4642
1087
651 2 3 12 6 0 389 87 242
4643
1088
652 2 3 12 6 0 374 367 356
4644
1089
653 2 3 12 6 0 374 356 147
4645
1090
654 2 3 12 6 0 242 379 297
4646
1091
655 2 3 12 6 0 242 378 379
4647
1092
656 2 3 12 6 0 379 398 80
4648
1093
657 2 3 12 6 0 379 378 398
4649
1094
658 2 3 12 6 0 388 295 234
4650
1095
659 2 3 12 6 0 295 388 293
4651
1096
660 2 3 12 6 0 388 234 336
4652
1097
661 2 3 12 6 0 388 257 293
4653
1098
662 2 3 12 6 0 393 394 387
4654
1099
663 2 3 12 6 0 396 395 287
4655
1100
664 2 3 12 6 0 396 335 263
4656
1101
665 2 3 12 6 0 383 45 355
4657
1102
666 2 3 12 6 0 45 383 44
4658
1103
667 2 3 12 6 0 402 79 302
4659
1104
668 2 3 12 6 0 79 402 350
4660
1105
669 2 3 12 6 0 403 280 315
4661
1106
670 2 3 12 6 0 223 405 346
4662
1107
671 2 3 12 6 0 256 405 223
4663
1108
672 2 3 12 6 0 406 360 390
4664
1109
673 2 3 12 6 0 406 380 360
4665
1110
674 2 3 12 6 0 380 406 109
4666
1111
675 2 3 12 6 0 49 50 387
4667
1112
676 2 3 12 6 0 394 49 387
4668
1113
677 2 3 12 6 0 168 407 119
4669
1114
678 2 3 12 6 0 407 168 222
4670
1115
679 2 3 12 6 0 407 395 119
4671
1116
680 2 3 12 6 0 377 408 269
4672
1117
681 2 3 12 6 0 408 377 398
4673
1118
682 2 3 12 6 0 408 191 269
4674
1119
683 2 3 12 6 0 38 409 37
4675
1120
684 2 3 12 6 0 409 354 183
4676
1121
685 2 3 12 6 0 410 122 372
4677
1122
686 2 3 12 6 0 411 33 34
4678
1123
687 2 3 12 6 0 278 411 34
4679
1124
688 2 3 12 6 0 342 55 56
4680
1125
689 2 3 12 6 0 220 139 324
4681
1126
690 2 3 12 6 0 324 221 147
4682
1127
691 2 3 12 6 0 356 324 147
4683
1128
692 2 3 12 6 0 221 324 148
4684
1129
693 2 3 12 6 0 392 353 261
4685
1130
694 2 3 12 6 0 403 392 361
4686
1131
695 2 3 12 6 0 401 61 62
4687
1132
696 2 3 12 6 0 61 401 270
4688
1133
697 2 3 12 6 0 62 329 401
4689
1134
698 2 3 12 6 0 270 401 197
4690
1135
699 2 3 12 6 0 250 109 402
4691
1136
700 2 3 12 6 0 406 402 109
4692
1137
701 2 3 12 6 0 402 406 350
4693
1138
702 2 3 12 6 0 38 39 409
4694
1139
703 2 3 12 6 0 39 399 409
4695
1140
704 2 3 12 6 0 14 15 384
4696
1141
705 2 3 12 6 0 220 14 384
4697
1142
706 2 3 12 6 0 410 128 122
4698
1143
707 2 3 12 6 0 35 36 410
4699
1144
708 2 3 12 6 0 372 35 410
4700
1145
709 2 3 12 6 0 411 32 33
4701
1146
710 2 3 12 6 0 412 366 196
4702
1147
711 2 3 12 6 0 414 366 412
4703
1148
712 2 3 12 6 0 279 412 196
4704
1149
713 2 3 12 6 0 413 394 393
4705
1150
714 2 3 12 6 0 393 415 413
4706
1151
715 2 3 12 6 0 414 178 366
4707
1152
716 2 3 12 6 0 400 178 414
4708
1153
717 2 3 12 6 0 415 176 152
4709
1154
718 2 3 12 6 0 176 415 393
4710
1155
719 2 3 12 6 0 404 415 152
4711
1156
720 2 3 12 6 0 416 305 346
4712
1157
721 2 3 12 6 0 305 416 391
4713
1158
722 2 3 12 6 0 405 416 346
4714
1159
723 2 3 12 6 0 360 416 389
4715
1160
724 2 3 12 6 0 416 360 391
4716
1161
725 2 3 12 6 0 416 405 389
4717
1162
726 2 3 12 6 0 366 105 181
4718
1163
727 2 3 12 6 0 105 366 178
4719
1164
728 2 3 12 6 0 98 143 216
4720
1165
729 2 3 12 6 0 143 98 214
4721
1166
730 2 3 12 6 0 254 153 300
4722
1167
731 2 3 12 6 0 153 254 171
4723
1168
732 2 3 12 6 0 214 206 142
4724
1169
733 2 3 12 6 0 206 214 98
4725
1170
734 2 3 12 6 0 307 114 260
4726
1171
735 2 3 12 6 0 114 307 173
4727
1172
736 2 3 12 6 0 338 127 198
4728
1173
737 2 3 12 6 0 127 338 371
4729
1174
738 2 3 12 6 0 355 358 383
4730
1175
739 2 3 12 6 0 358 355 279
4731
1176
740 2 3 12 6 0 381 352 351
4732
1177
741 2 3 12 6 0 352 381 205
4733
1178
742 2 3 12 6 0 168 94 258
4734
1179
743 2 3 12 6 0 94 168 153
4735
1180
744 2 3 12 6 0 222 258 334
4736
1181
745 2 3 12 6 0 258 222 168
4737
1182
746 2 3 12 6 0 196 358 279
4738
1183
747 2 3 12 6 0 358 196 85
4739
1184
748 2 3 12 6 0 177 263 310
4740
1185
749 2 3 12 6 0 263 177 78
4741
1186
750 2 3 12 6 0 78 396 263
4742
1187
751 2 3 12 6 0 396 78 395
4743
1188
752 2 3 12 6 0 405 87 389
4744
1189
753 2 3 12 6 0 87 405 256
4745
1190
754 2 3 12 6 0 155 150 94
4746
1191
755 2 3 12 6 0 155 120 150
4747
1192
756 2 3 12 6 0 217 219 96
4748
1193
757 2 3 12 6 0 217 146 219
4749
1194
758 2 3 12 6 0 178 233 105
4750
1195
759 2 3 12 6 0 178 174 233
4751
1196
760 2 3 12 6 0 107 106 84
4752
1197
761 2 3 12 6 0 107 179 106
4753
1198
762 2 3 12 6 0 112 144 206
4754
1199
763 2 3 12 6 0 218 144 112
4755
1200
764 2 3 12 6 0 193 164 91
4756
1201
765 2 3 12 6 0 352 164 193
4757
1202
766 2 3 12 6 0 386 367 12
4758
1203
767 2 3 12 6 0 386 356 367
4759
1204
768 2 3 12 6 0 417 128 410
4760
1205
769 2 3 12 6 0 143 420 287
4761
1206
770 2 3 12 6 0 420 143 214
4762
1207
771 2 3 12 6 0 420 95 287
4763
1208
772 2 3 12 6 0 421 18 132
4764
1209
773 2 3 12 6 0 130 421 132
4765
1210
774 2 3 12 6 0 18 421 17
4766
1211
775 2 3 12 6 0 210 422 142
4767
1212
776 2 3 12 6 0 215 422 141
4768
1213
777 2 3 12 6 0 422 210 141
4769
1214
778 2 3 12 6 0 423 76 4
4770
1215
779 2 3 12 6 0 76 423 194
4771
1216
780 2 3 12 6 0 5 423 4
4772
1217
781 2 3 12 6 0 95 424 215
4773
1218
782 2 3 12 6 0 142 424 214
4774
1219
783 2 3 12 6 0 321 423 5
4775
1220
784 2 3 12 6 0 321 194 423
4776
1221
785 2 3 12 6 0 422 424 142
4777
1222
786 2 3 12 6 0 422 215 424
4778
1223
787 2 3 12 6 0 420 424 95
4779
1224
788 2 3 12 6 0 420 214 424
4780
1225
789 2 3 12 6 0 318 129 392
4781
1226
790 2 3 12 6 0 129 361 392
4782
1227
791 2 3 12 6 0 128 417 183
4783
1228
792 2 3 12 6 0 409 417 37
4784
1229
793 2 3 12 6 0 417 409 183
4785
1230
794 2 3 12 6 0 36 37 417
4786
1231
795 2 3 12 6 0 410 36 417
4787
1232
796 2 3 12 6 0 332 418 279
4788
1233
797 2 3 12 6 0 418 412 279
4789
1234
798 2 3 12 6 0 419 412 418
4790
1235
799 2 3 12 6 0 413 419 418
4791
1236
800 2 3 12 6 0 419 404 414
4792
1237
801 2 3 12 6 0 412 419 414
4793
1238
802 2 3 12 6 0 404 413 415
4794
1239
803 2 3 12 6 0 413 404 419
4795
1240
804 2 3 12 6 0 142 144 210
4796
1241
805 2 3 12 6 0 144 142 206
4797
1242
806 2 3 12 6 0 83 425 152
4798
1243
807 2 3 12 6 0 404 425 414
4799
1244
808 2 3 12 6 0 425 404 152
4800
1245
809 2 3 12 6 0 56 426 342
4801
1246
810 2 3 12 6 0 408 426 191
4802
1247
811 2 3 12 6 0 426 408 342
4803
1248
812 2 3 12 6 0 427 16 17
4804
1249
813 2 3 12 6 0 427 130 281
4805
1250
814 2 3 12 6 0 16 427 281
4806
1251
815 2 3 12 6 0 421 427 17
4807
1252
816 2 3 12 6 0 421 130 427
4808
1253
817 2 3 12 6 0 400 425 83
4809
1254
818 2 3 12 6 0 400 414 425
4810
1255
819 2 3 12 6 0 56 57 426
4811
1256
820 2 3 12 6 0 57 191 426
4812
1257
821 2 3 12 6 0 48 428 47
4813
1258
822 2 3 12 6 0 428 48 394
4814
1259
823 2 3 12 6 0 413 428 394
4815
1260
824 2 3 12 6 0 413 418 428
4816
1261
825 2 3 12 6 0 428 332 47
4817
1262
826 2 3 12 6 0 418 332 428
4818
1263
827 2 3 12 6 0 344 247 345
4819
1264
828 2 3 12 6 0 247 344 301
4820
1265
829 2 3 12 6 0 386 324 356
4821
1266
830 2 3 12 6 0 324 386 220
4822
1267
831 2 3 12 6 0 315 392 403
4823
1268
832 2 3 12 6 0 392 315 353
4824
1269
833 2 3 12 6 0 411 361 32
4825
1270
834 2 3 12 6 0 361 411 403
4826
1271
835 2 3 12 6 0 244 154 93
4827
1272
836 2 3 12 6 0 154 244 299
4828
1273
837 2 3 12 6 0 351 193 28
4829
1274
838 2 3 12 6 0 193 351 352
4830
1275
839 2 3 12 6 0 110 333 53
4831
1276
840 2 3 12 6 0 333 110 256
4832
1277
841 2 3 12 6 0 15 281 384
4833
1278
842 2 3 12 6 0 16 281 15
4834
1279
843 2 3 12 6 0 136 26 91
4835
1280
844 2 3 12 6 0 26 136 25
4836
1281
845 2 3 12 6 0 215 212 157
4837
1282
846 2 3 12 6 0 212 215 141
4838
1283
847 2 3 12 6 0 201 199 271
4839
1284
848 2 3 12 6 0 199 201 200
4840
1285
849 2 3 12 6 0 344 289 304
4841
1286
850 2 3 12 6 0 289 344 252
4842
1287
851 2 3 12 6 0 278 403 411
4843
1288
852 2 3 12 6 0 403 278 280
4844
1289
853 2 3 12 6 0 59 58 1
4845
1290
854 2 3 12 6 0 58 59 268
4846
1291
$EndElements
4847
1292
0
4848
=== removed file 'tests/Stokes_mms_p0p1cv/src/MMS_C.geo'
4849
--- tests/Stokes_mms_p0p1cv/src/MMS_C.geo	2010-05-07 20:23:24 +0000
4850
+++ tests/Stokes_mms_p0p1cv/src/MMS_C.geo	1970-01-01 00:00:00 +0000
4851
@@ -1,15 +0,0 @@
4852
1
Point(1) = {0.0,0.0,0,0.08};
4853
2
Point(2) = {3.1415926535897931,0.0,0,0.08};
4854
3
Point(3) = {3.1415926535897931,3.1415926535897931,0,0.08};
4855
4
Point(4) = {0.0,3.1415926535897931,0,0.08};
4856
5
Line(1) = {4,3};
4857
6
Line(2) = {3,2};
4858
7
Line(3) = {2,1};
4859
8
Line(4) = {1,4};
4860
9
Line Loop(5) = {1,2,3,4};
4861
10
Plane Surface(6) = {5};
4862
11
Physical Line(7) = {3};
4863
12
Physical Line(8) = {2};
4864
13
Physical Line(9) = {1};
4865
14
Physical Line(10) = {4};
4866
15
Physical Surface(12) = {6};
4867
16
0
4868
=== removed file 'tests/Stokes_mms_p0p1cv/src/MMS_C.msh'
4869
--- tests/Stokes_mms_p0p1cv/src/MMS_C.msh	2010-05-07 20:23:24 +0000
4870
+++ tests/Stokes_mms_p0p1cv/src/MMS_C.msh	1970-01-01 00:00:00 +0000
4871
@@ -1,4828 +0,0 @@
4872
1
$MeshFormat
4873
2
2 0 8
4874
3
$EndMeshFormat
4875
4
$Nodes
4876
5
1607
4877
6
1 0 0 0
4878
7
2 3.141592653589793 0 0
4879
8
3 3.141592653589793 3.141592653589793 0
4880
9
4 0 3.141592653589793 0
4881
10
5 0.08055365778430337 3.141592653589793 0
4882
11
6 0.1611073155686067 3.141592653589793 0
4883
12
7 0.2416609733529101 3.141592653589793 0
4884
13
8 0.3222146311372135 3.141592653589793 0
4885
14
9 0.4027682889215169 3.141592653589793 0
4886
15
10 0.4833219467058202 3.141592653589793 0
4887
16
11 0.5638756044901236 3.141592653589793 0
4888
17
12 0.644429262274427 3.141592653589793 0
4889
18
13 0.7249829200587304 3.141592653589793 0
4890
19
14 0.8055365778430743 3.141592653589793 0
4891
20
15 0.8860902356275395 3.141592653589793 0
4892
21
16 0.9666438934120049 3.141592653589793 0
4893
22
17 1.04719755119647 3.141592653589793 0
4894
23
18 1.127751208980936 3.141592653589793 0
4895
24
19 1.208304866765401 3.141592653589793 0
4896
25
20 1.288858524549866 3.141592653589793 0
4897
26
21 1.369412182334331 3.141592653589793 0
4898
27
22 1.449965840118797 3.141592653589793 0
4899
28
23 1.530519497903262 3.141592653589793 0
4900
29
24 1.611073155687799 3.141592653589793 0
4901
30
25 1.691626813472406 3.141592653589793 0
4902
31
26 1.772180471257014 3.141592653589793 0
4903
32
27 1.852734129041621 3.141592653589793 0
4904
33
28 1.933287786826229 3.141592653589793 0
4905
34
29 2.013841444610837 3.141592653589793 0
4906
35
30 2.094395102395445 3.141592653589793 0
4907
36
31 2.174948760180052 3.141592653589793 0
4908
37
32 2.25550241796466 3.141592653589793 0
4909
38
33 2.336056075749267 3.141592653589793 0
4910
39
34 2.416609733533448 3.141592653589793 0
4911
40
35 2.497163391317486 3.141592653589793 0
4912
41
36 2.577717049101525 3.141592653589793 0
4913
42
37 2.658270706885563 3.141592653589793 0
4914
43
38 2.738824364669601 3.141592653589793 0
4915
44
39 2.81937802245364 3.141592653589793 0
4916
45
40 2.899931680237678 3.141592653589793 0
4917
46
41 2.980485338021716 3.141592653589793 0
4918
47
42 3.061038995805755 3.141592653589793 0
4919
48
43 3.141592653589793 3.061038995805297 0
4920
49
44 3.141592653589793 2.980485338020801 0
4921
50
45 3.141592653589793 2.899931680236305 0
4922
51
46 3.141592653589793 2.819378022451809 0
4923
52
47 3.141592653589793 2.738824364667313 0
4924
53
48 3.141592653589793 2.658270706882817 0
4925
54
49 3.141592653589793 2.577717049098321 0
4926
55
50 3.141592653589793 2.497163391313825 0
4927
56
51 3.141592653589793 2.416609733529329 0
4928
57
52 3.141592653589793 2.336056075744833 0
4929
58
53 3.141592653589793 2.255502417960337 0
4930
59
54 3.141592653589793 2.174948760175841 0
4931
60
55 3.141592653589793 2.094395102391345 0
4932
61
56 3.141592653589793 2.013841444606849 0
4933
62
57 3.141592653589793 1.933287786822353 0
4934
63
58 3.141592653589793 1.852734129037857 0
4935
64
59 3.141592653589793 1.772180471253361 0
4936
65
60 3.141592653589793 1.691626813468865 0
4937
66
61 3.141592653589793 1.611073155684369 0
4938
67
62 3.141592653589793 1.530519497899874 0
4939
68
63 3.141592653589793 1.449965840115377 0
4940
69
64 3.141592653589793 1.369412182330881 0
4941
70
65 3.141592653589793 1.288858524546385 0
4942
71
66 3.141592653589793 1.208304866761889 0
4943
72
67 3.141592653589793 1.127751208977393 0
4944
73
68 3.141592653589793 1.047197551192897 0
4945
74
69 3.141592653589793 0.966643893408401 0
4946
75
70 3.141592653589793 0.8860902356239051 0
4947
76
71 3.141592653589793 0.8055365778394092 0
4948
77
72 3.141592653589793 0.7249829200553402 0
4949
78
73 3.141592653589793 0.6444292622714136 0
4950
79
74 3.141592653589793 0.5638756044874866 0
4951
80
75 3.141592653589793 0.4833219467035597 0
4952
81
76 3.141592653589793 0.4027682889196336 0
4953
82
77 3.141592653589793 0.3222146311357066 0
4954
83
78 3.141592653589793 0.2416609733517801 0
4955
84
79 3.141592653589793 0.1611073155678531 0
4956
85
80 3.141592653589793 0.08055365778392698 0
4957
86
81 3.061038995805297 0 0
4958
87
82 2.980485338020801 0 0
4959
88
83 2.899931680236305 0 0
4960
89
84 2.819378022451809 0 0
4961
90
85 2.738824364667313 0 0
4962
91
86 2.658270706882817 0 0
4963
92
87 2.577717049098321 0 0
4964
93
88 2.497163391313825 0 0
4965
94
89 2.416609733529329 0 0
4966
95
90 2.336056075744833 0 0
4967
96
91 2.255502417960337 0 0
4968
97
92 2.174948760175841 0 0
4969
98
93 2.094395102391345 0 0
4970
99
94 2.013841444606849 0 0
4971
100
95 1.933287786822353 0 0
4972
101
96 1.852734129037857 0 0
4973
102
97 1.772180471253361 0 0
4974
103
98 1.691626813468865 0 0
4975
104
99 1.611073155684369 0 0
4976
105
100 1.530519497899874 0 0
4977
106
101 1.449965840115377 0 0
4978
107
102 1.369412182330881 0 0
4979
108
103 1.288858524546385 0 0
4980
109
104 1.208304866761889 0 0
4981
110
105 1.127751208977393 0 0
4982
111
106 1.047197551192897 0 0
4983
112
107 0.966643893408401 0 0
4984
113
108 0.8860902356239051 0 0
4985
114
109 0.8055365778394092 0 0
4986
115
110 0.7249829200553402 0 0
4987
116
111 0.6444292622714136 0 0
4988
117
112 0.5638756044874866 0 0
4989
118
113 0.4833219467035597 0 0
4990
119
114 0.4027682889196336 0 0
4991
120
115 0.3222146311357066 0 0
4992
121
116 0.2416609733517801 0 0
4993
122
117 0.1611073155678531 0 0
4994
123
118 0.08055365778392698 0 0
4995
124
119 0 0.08055365778430337 0
4996
125
120 0 0.1611073155686067 0
4997
126
121 0 0.2416609733529101 0
4998
127
122 0 0.3222146311372135 0
4999
128
123 0 0.4027682889215169 0
5000
129
124 0 0.4833219467058202 0
Reviewer	Review Type	Date Requested	Status
Jon Hill		2011-11-28	Approve on 2011-11-28
Review via email: mp+83605@code.launchpad.net