← Back to branch summary

~fluidity-core/fluidity/bug_1183080

« back to all changes in this revision

Viewing changes to tests/Helmholtz-anisotropic-scalar-mms-p1p1cg-unstructured/MMS_B_unstructured.flml

  • Committer: Tim Greaves
  • Date: 2013-04-19 11:53:43 UTC
  • mfrom: (3463.2.94 turbulence-clean)
  • Revision ID: tim.greaves@imperial.ac.uk-20130419115343-qq0m21v905aachuy
This merge is a replay of a previous merge which was reviewed in detail by
Christian Jacobs. There have been no changes since this was approved, hence I'm
fast-tracking this through as the concerns that led to the previous uncommit
have been resolved.

This merge brings in general updates from Jonathan Bull to his turbulence codes
and updates the BFS examples and various related tests to reflect this.

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/Helmholtz-anisotropic-scalar-mms-p1p1cg-unstructured/MMS_B_unstructured.flml
 
1
<?xml version='1.0' encoding='utf-8'?>
 
2
<fluidity_options>
 
3
  <simulation_name>
 
4
    <string_value lines="1">MMS_B_unstructured</string_value>
 
5
  </simulation_name>
 
6
  <problem_type>
 
7
    <string_value lines="1">fluids</string_value>
 
8
  </problem_type>
 
9
  <geometry>
 
10
    <dimension>
 
11
      <integer_value rank="0">2</integer_value>
 
12
    </dimension>
 
13
    <mesh name="CoordinateMesh">
 
14
      <from_mesh>
 
15
        <mesh name="InputMesh"/>
 
16
        <stat>
 
17
          <exclude_from_stat/>
 
18
        </stat>
 
19
      </from_mesh>
 
20
    </mesh>
 
21
    <mesh name="VelocityMesh">
 
22
      <from_mesh>
 
23
        <mesh name="InputMesh"/>
 
24
        <mesh_shape>
 
25
          <polynomial_degree>
 
26
            <integer_value rank="0">1</integer_value>
 
27
          </polynomial_degree>
 
28
        </mesh_shape>
 
29
        <stat>
 
30
          <exclude_from_stat/>
 
31
        </stat>
 
32
      </from_mesh>
 
33
    </mesh>
 
34
    <mesh name="PressureMesh">
 
35
      <from_mesh>
 
36
        <mesh name="InputMesh"/>
 
37
        <stat>
 
38
          <exclude_from_stat/>
 
39
        </stat>
 
40
      </from_mesh>
 
41
    </mesh>
 
42
    <mesh name="InputMesh">
 
43
      <from_file file_name="MMS_B_unstructured">
 
44
        <format name="triangle"/>
 
45
        <stat>
 
46
          <include_in_stat/>
 
47
        </stat>
 
48
      </from_file>
 
49
    </mesh>
 
50
    <mesh name="ErrorMesh">
 
51
      <from_mesh>
 
52
        <mesh name="CoordinateMesh"/>
 
53
        <mesh_shape>
 
54
          <polynomial_degree>
 
55
            <integer_value rank="0">4</integer_value>
 
56
          </polynomial_degree>
 
57
        </mesh_shape>
 
58
        <mesh_continuity>
 
59
          <string_value>continuous</string_value>
 
60
        </mesh_continuity>
 
61
        <stat>
 
62
          <exclude_from_stat/>
 
63
        </stat>
 
64
      </from_mesh>
 
65
    </mesh>
 
66
    <quadrature>
 
67
      <degree>
 
68
        <integer_value rank="0">8</integer_value>
 
69
      </degree>
 
70
    </quadrature>
 
71
  </geometry>
 
72
  <io>
 
73
    <dump_format>
 
74
      <string_value>vtk</string_value>
 
75
    </dump_format>
 
76
    <dump_period>
 
77
      <constant>
 
78
        <real_value rank="0">1</real_value>
 
79
      </constant>
 
80
    </dump_period>
 
81
    <output_mesh name="VelocityMesh"/>
 
82
    <stat/>
 
83
  </io>
 
84
  <timestepping>
 
85
    <current_time>
 
86
      <real_value rank="0">0.0</real_value>
 
87
    </current_time>
 
88
    <timestep>
 
89
      <real_value rank="0">0.1</real_value>
 
90
    </timestep>
 
91
    <finish_time>
 
92
      <real_value rank="0">1</real_value>
 
93
    </finish_time>
 
94
    <nonlinear_iterations>
 
95
      <integer_value rank="0">2</integer_value>
 
96
    </nonlinear_iterations>
 
97
  </timestepping>
 
98
  <material_phase name="NS">
 
99
    <vector_field name="Velocity" rank="1">
 
100
      <prescribed>
 
101
        <mesh name="VelocityMesh"/>
 
102
        <value name="WholeMesh">
 
103
          <python>
 
104
            <string_value lines="20" type="code" language="python">def val(X,t):
 
105
    import helmholtz_tools as g
 
106
    return g.velocity(X)</string_value>
 
107
          </python>
 
108
        </value>
 
109
        <output/>
 
110
        <stat>
 
111
          <include_in_stat/>
 
112
        </stat>
 
113
        <detectors>
 
114
          <exclude_from_detectors/>
 
115
        </detectors>
 
116
      </prescribed>
 
117
    </vector_field>
 
118
    <scalar_field name="Tracer" rank="0">
 
119
      <prognostic>
 
120
        <mesh name="VelocityMesh"/>
 
121
        <equation name="AdvectionDiffusion"/>
 
122
        <spatial_discretisation>
 
123
          <continuous_galerkin>
 
124
            <stabilisation>
 
125
              <no_stabilisation/>
 
126
            </stabilisation>
 
127
            <advection_terms/>
 
128
            <mass_terms/>
 
129
          </continuous_galerkin>
 
130
          <conservative_advection>
 
131
            <real_value rank="0">0.</real_value>
 
132
          </conservative_advection>
 
133
        </spatial_discretisation>
 
134
        <temporal_discretisation>
 
135
          <theta>
 
136
            <real_value rank="0">0.5</real_value>
 
137
          </theta>
 
138
        </temporal_discretisation>
 
139
        <solver>
 
140
          <iterative_method name="gmres">
 
141
            <restart>
 
142
              <integer_value rank="0">30</integer_value>
 
143
            </restart>
 
144
          </iterative_method>
 
145
          <preconditioner name="sor"/>
 
146
          <relative_error>
 
147
            <real_value rank="0">1.e-7</real_value>
 
148
          </relative_error>
 
149
          <max_iterations>
 
150
            <integer_value rank="0">1000</integer_value>
 
151
          </max_iterations>
 
152
          <never_ignore_solver_failures/>
 
153
          <diagnostics>
 
154
            <monitors/>
 
155
          </diagnostics>
 
156
        </solver>
 
157
        <initial_condition name="WholeMesh">
 
158
          <python>
 
159
            <string_value lines="20" type="code" language="python">def val(X,t):
 
160
    import helmholtz_tools as g
 
161
    return g.gaussian(X)</string_value>
 
162
          </python>
 
163
        </initial_condition>
 
164
        <boundary_conditions name="sides">
 
165
          <surface_ids>
 
166
            <integer_value shape="4" rank="1">7 8 9 10</integer_value>
 
167
          </surface_ids>
 
168
          <type name="dirichlet">
 
169
            <python>
 
170
              <string_value lines="20" type="code" language="python">def val(X,t):
 
171
    import helmholtz_tools as g
 
172
    return g.gaussian(X)</string_value>
 
173
            </python>
 
174
          </type>
 
175
        </boundary_conditions>
 
176
        <output/>
 
177
        <stat/>
 
178
        <convergence>
 
179
          <include_in_convergence/>
 
180
        </convergence>
 
181
        <detectors>
 
182
          <include_in_detectors/>
 
183
        </detectors>
 
184
        <steady_state>
 
185
          <include_in_steady_state/>
 
186
        </steady_state>
 
187
        <consistent_interpolation/>
 
188
      </prognostic>
 
189
    </scalar_field>
 
190
    <scalar_field name="TracerProjection" rank="0">
 
191
      <diagnostic>
 
192
        <algorithm source_field_type="scalar" material_phase_support="single" name="scalar_galerkin_projection" source_field_name="Tracer">
 
193
          <solver>
 
194
            <iterative_method name="cg"/>
 
195
            <preconditioner name="sor"/>
 
196
            <relative_error>
 
197
              <real_value rank="0">1.e-10</real_value>
 
198
            </relative_error>
 
199
            <max_iterations>
 
200
              <integer_value rank="0">1000</integer_value>
 
201
            </max_iterations>
 
202
            <never_ignore_solver_failures/>
 
203
            <diagnostics>
 
204
              <monitors/>
 
205
            </diagnostics>
 
206
          </solver>
 
207
        </algorithm>
 
208
        <mesh name="ErrorMesh"/>
 
209
        <output/>
 
210
        <stat/>
 
211
        <convergence>
 
212
          <include_in_convergence/>
 
213
        </convergence>
 
214
        <detectors>
 
215
          <include_in_detectors/>
 
216
        </detectors>
 
217
        <steady_state>
 
218
          <include_in_steady_state/>
 
219
        </steady_state>
 
220
      </diagnostic>
 
221
    </scalar_field>
 
222
    <scalar_field name="IsoFilteredTracer" rank="0">
 
223
      <diagnostic>
 
224
        <algorithm source_field_type="component" material_phase_support="single" name="helmholtz_smoothed_scalar" source_field_name="Tracer">
 
225
          <smoothing_scale_factor>
 
226
            <real_value rank="0">2.0</real_value>
 
227
          </smoothing_scale_factor>
 
228
          <solver>
 
229
            <iterative_method name="cg"/>
 
230
            <preconditioner name="sor"/>
 
231
            <relative_error>
 
232
              <real_value rank="0">1.e-7</real_value>
 
233
            </relative_error>
 
234
            <max_iterations>
 
235
              <integer_value rank="0">1000</integer_value>
 
236
            </max_iterations>
 
237
            <never_ignore_solver_failures/>
 
238
            <diagnostics>
 
239
              <monitors/>
 
240
            </diagnostics>
 
241
          </solver>
 
242
        </algorithm>
 
243
        <mesh name="VelocityMesh"/>
 
244
        <output/>
 
245
        <stat/>
 
246
        <convergence>
 
247
          <include_in_convergence/>
 
248
        </convergence>
 
249
        <detectors>
 
250
          <include_in_detectors/>
 
251
        </detectors>
 
252
        <steady_state>
 
253
          <include_in_steady_state/>
 
254
        </steady_state>
 
255
      </diagnostic>
 
256
    </scalar_field>
 
257
    <scalar_field name="IsoFilteredTracerProjection" rank="0">
 
258
      <diagnostic>
 
259
        <algorithm source_field_type="scalar" material_phase_support="single" name="scalar_galerkin_projection" source_field_name="IsoFilteredTracer">
 
260
          <solver>
 
261
            <iterative_method name="cg"/>
 
262
            <preconditioner name="sor"/>
 
263
            <relative_error>
 
264
              <real_value rank="0">1.e-10</real_value>
 
265
            </relative_error>
 
266
            <max_iterations>
 
267
              <integer_value rank="0">1000</integer_value>
 
268
            </max_iterations>
 
269
            <never_ignore_solver_failures/>
 
270
            <diagnostics>
 
271
              <monitors/>
 
272
            </diagnostics>
 
273
          </solver>
 
274
        </algorithm>
 
275
        <mesh name="ErrorMesh"/>
 
276
        <output/>
 
277
        <stat/>
 
278
        <convergence>
 
279
          <include_in_convergence/>
 
280
        </convergence>
 
281
        <detectors>
 
282
          <include_in_detectors/>
 
283
        </detectors>
 
284
        <steady_state>
 
285
          <include_in_steady_state/>
 
286
        </steady_state>
 
287
      </diagnostic>
 
288
    </scalar_field>
 
289
    <scalar_field name="IsoFilterDifference" rank="0">
 
290
      <diagnostic>
 
291
        <algorithm source_field_2_type="scalar" name="scalar_difference" source_field_1_name="TracerProjection" source_field_2_name="IsoFilteredTracerProjection" material_phase_support="single" source_field_1_type="scalar"/>
 
292
        <mesh name="ErrorMesh"/>
 
293
        <output/>
 
294
        <stat/>
 
295
        <convergence>
 
296
          <include_in_convergence/>
 
297
        </convergence>
 
298
        <detectors>
 
299
          <include_in_detectors/>
 
300
        </detectors>
 
301
        <steady_state>
 
302
          <include_in_steady_state/>
 
303
        </steady_state>
 
304
      </diagnostic>
 
305
    </scalar_field>
 
306
    <scalar_field name="AnisoFilteredTracer" rank="0">
 
307
      <diagnostic>
 
308
        <algorithm source_field_type="scalar" material_phase_support="single" name="helmholtz_anisotropic_smoothed_scalar" source_field_name="Tracer">
 
309
          <smoothing_scale_factor>
 
310
            <real_value rank="0">2.0</real_value>
 
311
          </smoothing_scale_factor>
 
312
          <solver>
 
313
            <iterative_method name="cg"/>
 
314
            <preconditioner name="sor"/>
 
315
            <relative_error>
 
316
              <real_value rank="0">1.e-7</real_value>
 
317
            </relative_error>
 
318
            <max_iterations>
 
319
              <integer_value rank="0">1000</integer_value>
 
320
            </max_iterations>
 
321
            <never_ignore_solver_failures/>
 
322
            <diagnostics>
 
323
              <monitors/>
 
324
            </diagnostics>
 
325
          </solver>
 
326
        </algorithm>
 
327
        <mesh name="VelocityMesh"/>
 
328
        <output/>
 
329
        <stat/>
 
330
        <convergence>
 
331
          <include_in_convergence/>
 
332
        </convergence>
 
333
        <detectors>
 
334
          <include_in_detectors/>
 
335
        </detectors>
 
336
        <steady_state>
 
337
          <include_in_steady_state/>
 
338
        </steady_state>
 
339
      </diagnostic>
 
340
    </scalar_field>
 
341
    <scalar_field name="AnisoFilteredTracerProjection" rank="0">
 
342
      <diagnostic>
 
343
        <algorithm source_field_type="scalar" material_phase_support="single" name="scalar_galerkin_projection" source_field_name="AnisoFilteredTracer">
 
344
          <solver>
 
345
            <iterative_method name="cg"/>
 
346
            <preconditioner name="sor"/>
 
347
            <relative_error>
 
348
              <real_value rank="0">1.e-10</real_value>
 
349
            </relative_error>
 
350
            <max_iterations>
 
351
              <integer_value rank="0">1000</integer_value>
 
352
            </max_iterations>
 
353
            <never_ignore_solver_failures/>
 
354
            <diagnostics>
 
355
              <monitors/>
 
356
            </diagnostics>
 
357
          </solver>
 
358
        </algorithm>
 
359
        <mesh name="ErrorMesh"/>
 
360
        <output/>
 
361
        <stat/>
 
362
        <convergence>
 
363
          <include_in_convergence/>
 
364
        </convergence>
 
365
        <detectors>
 
366
          <include_in_detectors/>
 
367
        </detectors>
 
368
        <steady_state>
 
369
          <include_in_steady_state/>
 
370
        </steady_state>
 
371
      </diagnostic>
 
372
    </scalar_field>
 
373
    <scalar_field name="AnisoFilterDifference" rank="0">
 
374
      <diagnostic>
 
375
        <algorithm source_field_2_type="scalar" name="scalar_difference" source_field_1_name="TracerProjection" source_field_2_name="AnisoFilteredTracerProjection" material_phase_support="single" source_field_1_type="scalar"/>
 
376
        <mesh name="ErrorMesh"/>
 
377
        <output/>
 
378
        <stat/>
 
379
        <convergence>
 
380
          <include_in_convergence/>
 
381
        </convergence>
 
382
        <detectors>
 
383
          <include_in_detectors/>
 
384
        </detectors>
 
385
        <steady_state>
 
386
          <include_in_steady_state/>
 
387
        </steady_state>
 
388
      </diagnostic>
 
389
    </scalar_field>
 
390
    <scalar_field name="MLFilteredTracer" rank="0">
 
391
      <diagnostic>
 
392
        <algorithm source_field_type="scalar" material_phase_support="single" name="lumped_mass_smoothed_scalar" source_field_name="Tracer"/>
 
393
        <mesh name="VelocityMesh"/>
 
394
        <output/>
 
395
        <stat/>
 
396
        <convergence>
 
397
          <include_in_convergence/>
 
398
        </convergence>
 
399
        <detectors>
 
400
          <include_in_detectors/>
 
401
        </detectors>
 
402
        <steady_state>
 
403
          <include_in_steady_state/>
 
404
        </steady_state>
 
405
      </diagnostic>
 
406
    </scalar_field>
 
407
    <scalar_field name="MLFilteredTracerProjection" rank="0">
 
408
      <diagnostic>
 
409
        <algorithm source_field_type="scalar" material_phase_support="single" name="scalar_galerkin_projection" source_field_name="MLFilteredTracer">
 
410
          <solver>
 
411
            <iterative_method name="cg"/>
 
412
            <preconditioner name="sor"/>
 
413
            <relative_error>
 
414
              <real_value rank="0">1.e-10</real_value>
 
415
            </relative_error>
 
416
            <max_iterations>
 
417
              <integer_value rank="0">1000</integer_value>
 
418
            </max_iterations>
 
419
            <never_ignore_solver_failures/>
 
420
            <diagnostics>
 
421
              <monitors/>
 
422
            </diagnostics>
 
423
          </solver>
 
424
        </algorithm>
 
425
        <mesh name="ErrorMesh"/>
 
426
        <output/>
 
427
        <stat/>
 
428
        <convergence>
 
429
          <include_in_convergence/>
 
430
        </convergence>
 
431
        <detectors>
 
432
          <include_in_detectors/>
 
433
        </detectors>
 
434
        <steady_state>
 
435
          <include_in_steady_state/>
 
436
        </steady_state>
 
437
      </diagnostic>
 
438
    </scalar_field>
 
439
    <scalar_field name="MLFilterDifference" rank="0">
 
440
      <diagnostic>
 
441
        <algorithm source_field_2_type="scalar" name="scalar_difference" source_field_1_name="TracerProjection" source_field_2_name="MLFilteredTracerProjection" material_phase_support="single" source_field_1_type="scalar"/>
 
442
        <mesh name="ErrorMesh"/>
 
443
        <output/>
 
444
        <stat/>
 
445
        <convergence>
 
446
          <include_in_convergence/>
 
447
        </convergence>
 
448
        <detectors>
 
449
          <include_in_detectors/>
 
450
        </detectors>
 
451
        <steady_state>
 
452
          <include_in_steady_state/>
 
453
        </steady_state>
 
454
      </diagnostic>
 
455
    </scalar_field>
 
456
    <vector_field name="GradTracer" rank="1">
 
457
      <diagnostic>
 
458
        <algorithm source_field_type="scalar" material_phase_support="single" name="grad" source_field_name="Tracer"/>
 
459
        <mesh name="VelocityMesh"/>
 
460
        <output/>
 
461
        <stat>
 
462
          <include_in_stat/>
 
463
        </stat>
 
464
        <convergence>
 
465
          <include_in_convergence/>
 
466
        </convergence>
 
467
        <detectors>
 
468
          <include_in_detectors/>
 
469
        </detectors>
 
470
        <steady_state>
 
471
          <include_in_steady_state/>
 
472
        </steady_state>
 
473
      </diagnostic>
 
474
    </vector_field>
 
475
    <vector_field name="IsoFilteredGradTracer" rank="1">
 
476
      <diagnostic>
 
477
        <algorithm source_field_type="vector" material_phase_support="single" name="helmholtz_smoothed_vector" source_field_name="GradTracer">
 
478
          <smoothing_scale_factor>
 
479
            <real_value rank="0">2.0</real_value>
 
480
          </smoothing_scale_factor>
 
481
          <solver>
 
482
            <iterative_method name="cg"/>
 
483
            <preconditioner name="sor"/>
 
484
            <relative_error>
 
485
              <real_value rank="0">1.e-7</real_value>
 
486
            </relative_error>
 
487
            <max_iterations>
 
488
              <integer_value rank="0">1000</integer_value>
 
489
            </max_iterations>
 
490
            <never_ignore_solver_failures/>
 
491
            <diagnostics>
 
492
              <monitors/>
 
493
            </diagnostics>
 
494
          </solver>
 
495
        </algorithm>
 
496
        <mesh name="VelocityMesh"/>
 
497
        <output/>
 
498
        <stat>
 
499
          <include_in_stat/>
 
500
        </stat>
 
501
        <convergence>
 
502
          <include_in_convergence/>
 
503
        </convergence>
 
504
        <detectors>
 
505
          <include_in_detectors/>
 
506
        </detectors>
 
507
        <steady_state>
 
508
          <include_in_steady_state/>
 
509
        </steady_state>
 
510
      </diagnostic>
 
511
    </vector_field>
 
512
    <vector_field name="GradIsoFilteredTracer" rank="1">
 
513
      <diagnostic>
 
514
        <algorithm source_field_type="scalar" material_phase_support="single" name="grad" source_field_name="IsoFilteredTracer"/>
 
515
        <mesh name="VelocityMesh"/>
 
516
        <output/>
 
517
        <stat>
 
518
          <include_in_stat/>
 
519
        </stat>
 
520
        <convergence>
 
521
          <include_in_convergence/>
 
522
        </convergence>
 
523
        <detectors>
 
524
          <include_in_detectors/>
 
525
        </detectors>
 
526
        <steady_state>
 
527
          <include_in_steady_state/>
 
528
        </steady_state>
 
529
      </diagnostic>
 
530
    </vector_field>
 
531
    <vector_field name="IsoCommError" rank="1">
 
532
      <diagnostic>
 
533
        <algorithm source_field_2_type="vector" name="vector_difference" source_field_1_name="IsoFilteredGradTracer" source_field_2_name="GradIsoFilteredTracer" material_phase_support="single" source_field_1_type="vector">
 
534
          <absolute_difference/>
 
535
        </algorithm>
 
536
        <mesh name="ErrorMesh"/>
 
537
        <output/>
 
538
        <stat>
 
539
          <include_in_stat/>
 
540
        </stat>
 
541
        <convergence>
 
542
          <include_in_convergence/>
 
543
        </convergence>
 
544
        <detectors>
 
545
          <include_in_detectors/>
 
546
        </detectors>
 
547
        <steady_state>
 
548
          <include_in_steady_state/>
 
549
        </steady_state>
 
550
      </diagnostic>
 
551
    </vector_field>
 
552
    <vector_field name="AnisoFilteredGradTracer" rank="1">
 
553
      <diagnostic>
 
554
        <algorithm source_field_type="vector" material_phase_support="single" name="helmholtz_anisotropic_smoothed_vector" source_field_name="GradTracer">
 
555
          <smoothing_scale_factor>
 
556
            <real_value rank="0">2.0</real_value>
 
557
          </smoothing_scale_factor>
 
558
          <solver>
 
559
            <iterative_method name="cg"/>
 
560
            <preconditioner name="sor"/>
 
561
            <relative_error>
 
562
              <real_value rank="0">1.e-7</real_value>
 
563
            </relative_error>
 
564
            <max_iterations>
 
565
              <integer_value rank="0">1000</integer_value>
 
566
            </max_iterations>
 
567
            <never_ignore_solver_failures/>
 
568
            <diagnostics>
 
569
              <monitors/>
 
570
            </diagnostics>
 
571
          </solver>
 
572
        </algorithm>
 
573
        <mesh name="VelocityMesh"/>
 
574
        <output/>
 
575
        <stat>
 
576
          <include_in_stat/>
 
577
        </stat>
 
578
        <convergence>
 
579
          <include_in_convergence/>
 
580
        </convergence>
 
581
        <detectors>
 
582
          <include_in_detectors/>
 
583
        </detectors>
 
584
        <steady_state>
 
585
          <include_in_steady_state/>
 
586
        </steady_state>
 
587
      </diagnostic>
 
588
    </vector_field>
 
589
    <vector_field name="GradAnisoFilteredTracer" rank="1">
 
590
      <diagnostic>
 
591
        <algorithm source_field_type="scalar" material_phase_support="single" name="grad" source_field_name="AnisoFilteredTracer"/>
 
592
        <mesh name="VelocityMesh"/>
 
593
        <output/>
 
594
        <stat>
 
595
          <include_in_stat/>
 
596
        </stat>
 
597
        <convergence>
 
598
          <include_in_convergence/>
 
599
        </convergence>
 
600
        <detectors>
 
601
          <include_in_detectors/>
 
602
        </detectors>
 
603
        <steady_state>
 
604
          <include_in_steady_state/>
 
605
        </steady_state>
 
606
      </diagnostic>
 
607
    </vector_field>
 
608
    <vector_field name="AnisoCommError" rank="1">
 
609
      <diagnostic>
 
610
        <algorithm source_field_2_type="vector" name="vector_difference" source_field_1_name="AnisoFilteredGradTracer" source_field_2_name="GradAnisoFilteredTracer" material_phase_support="single" source_field_1_type="vector">
 
611
          <absolute_difference/>
 
612
        </algorithm>
 
613
        <mesh name="ErrorMesh"/>
 
614
        <output/>
 
615
        <stat>
 
616
          <include_in_stat/>
 
617
        </stat>
 
618
        <convergence>
 
619
          <include_in_convergence/>
 
620
        </convergence>
 
621
        <detectors>
 
622
          <include_in_detectors/>
 
623
        </detectors>
 
624
        <steady_state>
 
625
          <include_in_steady_state/>
 
626
        </steady_state>
 
627
      </diagnostic>
 
628
    </vector_field>
 
629
    <vector_field name="MLFilteredGradTracer" rank="1">
 
630
      <diagnostic>
 
631
        <algorithm source_field_type="vector" material_phase_support="single" name="lumped_mass_smoothed_vector" source_field_name="GradTracer"/>
 
632
        <mesh name="VelocityMesh"/>
 
633
        <output/>
 
634
        <stat>
 
635
          <include_in_stat/>
 
636
        </stat>
 
637
        <convergence>
 
638
          <include_in_convergence/>
 
639
        </convergence>
 
640
        <detectors>
 
641
          <include_in_detectors/>
 
642
        </detectors>
 
643
        <steady_state>
 
644
          <include_in_steady_state/>
 
645
        </steady_state>
 
646
      </diagnostic>
 
647
    </vector_field>
 
648
    <vector_field name="GradMLFilteredTracer" rank="1">
 
649
      <diagnostic>
 
650
        <algorithm source_field_type="scalar" material_phase_support="single" name="grad" source_field_name="MLFilteredTracer"/>
 
651
        <mesh name="VelocityMesh"/>
 
652
        <output/>
 
653
        <stat>
 
654
          <include_in_stat/>
 
655
        </stat>
 
656
        <convergence>
 
657
          <include_in_convergence/>
 
658
        </convergence>
 
659
        <detectors>
 
660
          <include_in_detectors/>
 
661
        </detectors>
 
662
        <steady_state>
 
663
          <include_in_steady_state/>
 
664
        </steady_state>
 
665
      </diagnostic>
 
666
    </vector_field>
 
667
    <vector_field name="MLCommError" rank="1">
 
668
      <diagnostic>
 
669
        <algorithm source_field_2_type="vector" name="vector_difference" source_field_1_name="MLFilteredGradTracer" source_field_2_name="GradMLFilteredTracer" material_phase_support="single" source_field_1_type="vector">
 
670
          <absolute_difference/>
 
671
        </algorithm>
 
672
        <mesh name="ErrorMesh"/>
 
673
        <output/>
 
674
        <stat>
 
675
          <include_in_stat/>
 
676
        </stat>
 
677
        <convergence>
 
678
          <include_in_convergence/>
 
679
        </convergence>
 
680
        <detectors>
 
681
          <include_in_detectors/>
 
682
        </detectors>
 
683
        <steady_state>
 
684
          <include_in_steady_state/>
 
685
        </steady_state>
 
686
      </diagnostic>
 
687
    </vector_field>
 
688
  </material_phase>
 
689
</fluidity_options>