From 402eba50626849ddd3562a9b4fb4591dd9380505 Mon Sep 17 00:00:00 2001
From: Kilian Weishaupt <kilian.weishaupt@iws.uni-stuttgart.de>
Date: Thu, 30 Nov 2017 11:02:47 +0100
Subject: [PATCH] [1p2c][test] Make mpfa tests pass
* Use reference file without velocity
---
.../1pnc/implicit/1p2cniconductionproblem.hh | 1 +
.../1pnc/implicit/CMakeLists.txt | 29 +++--
.../1p2cniccmpfaconduction-reference.vtu | 104 ++++++++--------
test/references/1p2ctestmpfa-reference.vtu | 113 ++++++++++++++++++
4 files changed, 186 insertions(+), 61 deletions(-)
create mode 100644 test/references/1p2ctestmpfa-reference.vtu
diff --git a/test/porousmediumflow/1pnc/implicit/1p2cniconductionproblem.hh b/test/porousmediumflow/1pnc/implicit/1p2cniconductionproblem.hh
index a8fdaa7b6e..7fcdc95235 100644
--- a/test/porousmediumflow/1pnc/implicit/1p2cniconductionproblem.hh
+++ b/test/porousmediumflow/1pnc/implicit/1p2cniconductionproblem.hh
@@ -45,6 +45,7 @@ namespace Properties
{
NEW_TYPE_TAG(OnePTwoCNIConductionProblem, INHERITS_FROM(OnePNCNI));
NEW_TYPE_TAG(OnePTwoCNIConductionCCTpfaProblem, INHERITS_FROM(CCTpfaModel, OnePTwoCNIConductionProblem));
+NEW_TYPE_TAG(OnePTwoCNIConductionCCMpfaProblem, INHERITS_FROM(CCMpfaModel, OnePTwoCNIConductionProblem));
NEW_TYPE_TAG(OnePTwoCNIConductionBoxProblem, INHERITS_FROM(BoxModel, OnePTwoCNIConductionProblem));
// Set the grid type
diff --git a/test/porousmediumflow/1pnc/implicit/CMakeLists.txt b/test/porousmediumflow/1pnc/implicit/CMakeLists.txt
index 6631bb494d..7c57dcd523 100644
--- a/test/porousmediumflow/1pnc/implicit/CMakeLists.txt
+++ b/test/porousmediumflow/1pnc/implicit/CMakeLists.txt
@@ -20,17 +20,27 @@ dune_add_test(NAME test_1p2c_tpfa
--command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2c_tpfa test_1p2c_fv.input -Problem.Name 1pnctest_tpfa"
--zeroThreshold {"velocity_w \(m/s\)_1":1e-15})
+#Note: The mpfa model does not support velocity output yet. We therefore deactivate it for the test and use a copy of the tpfa reference
+# solution where the velocities have been removed.
dune_add_test(NAME test_1p2c_mpfa
SOURCES test_1p2c_fv.cc
COMPILE_DEFINITIONS TYPETAG=OnePTwoCTestCCMpfaProblem
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
CMD_ARGS --script fuzzy
- --files ${CMAKE_SOURCE_DIR}/test/references/1p2ctestcc-reference.vtu
+ --files ${CMAKE_SOURCE_DIR}/test/references/1p2ctestmpfa-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/1pnctest_mpfa-00009.vtu
- --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2c_mpfa test_1p2c_fv.input -Problem.Name 1pnctest_mpfa"
- --zeroThreshold {"velocity_w \(m/s\)_1":1e-15})
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2c_mpfa test_1p2c_fv.input -Problem.Name 1pnctest_mpfa -Vtk.AddVelocity false")
# non-isothermal tests
+dune_add_test(NAME test_1p2cni_conduction_box
+ SOURCES test_1p2cni_conduction_fv.cc
+ COMPILE_DEFINITIONS TYPETAG=OnePTwoCNIConductionBoxProblem
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/1p2cniboxconduction-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/1p2cni_conductiontest_box-00006.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2cni_conduction_box test_1p2cni_conduction_fv.input -Problem.Name 1p2cni_conductiontest_box"
+ --zeroThreshold {"velocity_w \(m/s\)":1e-9})
dune_add_test(NAME test_1p2cni_conduction_tpfa
SOURCES test_1p2cni_conduction_fv.cc
@@ -42,12 +52,13 @@ dune_add_test(NAME test_1p2cni_conduction_tpfa
--command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2cni_conduction_tpfa test_1p2cni_conduction_fv.input -Problem.Name 1p2cni_conductiontest_tpfa"
--zeroThreshold {"velocity_w \(m/s\)_0":1e-9})
-dune_add_test(NAME test_1p2cni_conduction_box
+#Note: The mpfa model does not support velocity output yet. We therefore deactivate it for the test and use a copy of the tpfa reference
+# solution where the velocities have been removed.
+dune_add_test(NAME test_1p2cni_conduction_mpfa
SOURCES test_1p2cni_conduction_fv.cc
- COMPILE_DEFINITIONS TYPETAG=OnePTwoCNIConductionBoxProblem
+ COMPILE_DEFINITIONS TYPETAG=OnePTwoCNIConductionCCMpfaProblem
COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
CMD_ARGS --script fuzzy
- --files ${CMAKE_SOURCE_DIR}/test/references/1p2cniboxconduction-reference.vtu
- ${CMAKE_CURRENT_BINARY_DIR}/1p2cni_conductiontest_box-00006.vtu
- --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2cni_conduction_box test_1p2cni_conduction_fv.input -Problem.Name 1p2cni_conductiontest_box"
- --zeroThreshold {"velocity_w \(m/s\)":1e-9})
+ --files ${CMAKE_SOURCE_DIR}/test/references/1p2cniccmpfaconduction-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/1p2cni_conductiontest_mpfa-00006.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p2cni_conduction_mpfa test_1p2cni_conduction_fv.input -Problem.Name 1p2cni_conductiontest_mpfa -Vtk.AddVelocity false")
diff --git a/test/references/1p2cniccmpfaconduction-reference.vtu b/test/references/1p2cniccmpfaconduction-reference.vtu
index 144b4f4118..6e7bc268f3 100644
--- a/test/references/1p2cniccmpfaconduction-reference.vtu
+++ b/test/references/1p2cniccmpfaconduction-reference.vtu
@@ -2,8 +2,8 @@
<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
<UnstructuredGrid>
<Piece NumberOfCells="200" NumberOfPoints="402">
- <CellData Scalars="P">
- <DataArray type="Float32" Name="P" NumberOfComponents="1" format="ascii">
+ <CellData Scalars="pressure">
+ <DataArray type="Float32" Name="pressure" NumberOfComponents="1" format="ascii">
100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
@@ -23,23 +23,23 @@
100000 100000 100000 100000 100000 100000 100000 100000
</DataArray>
<DataArray type="Float32" Name="delp" NumberOfComponents="1" format="ascii">
- 0.000125888 0.000379777 0.000633389 0.000884312 0.00113124 0.00137247 0.00160609 0.00183034 0.00204366 0.00224472 0.00243243 0.00260818
- 0.00277185 0.00292308 0.00306169 0.00318771 0.00330128 0.00340272 0.00349243 0.00357092 0.00363874 0.00369653 0.00374492 0.00378458
- 0.00381637 0.00384174 0.00386168 0.00387667 0.00388718 0.00389365 0.00389649 0.00389607 0.00389275 0.00388684 0.00387865 0.00386842
- 0.00385641 0.00384282 0.00382785 0.00381167 0.00379443 0.00377626 0.00375729 0.00373761 0.00371733 0.00369652 0.00367526 0.0036536
- 0.0036316 0.00360931 0.00358676 0.003564 0.00354104 0.00351793 0.00349468 0.00347132 0.00344785 0.0034243 0.00340067 0.00337698
- 0.00335323 0.00332944 0.00330561 0.00328175 0.00325786 0.00323395 0.00321001 0.00318606 0.00316209 0.00313811 0.00311412 0.00309012
- 0.00306611 0.0030421 0.00301808 0.00299405 0.00297003 0.00294599 0.00292196 0.00289792 0.00287388 0.00284984 0.0028258 0.00280176
- 0.00277772 0.00275367 0.00272962 0.00270558 0.00268153 0.00265748 0.00263344 0.00260939 0.00258534 0.00256129 0.00253724 0.00251319
- 0.00248914 0.0024651 0.00244105 0.002417 0.00239295 0.0023689 0.00234485 0.0023208 0.00229675 0.0022727 0.00224865 0.0022246
- 0.00220055 0.0021765 0.00215245 0.0021284 0.00210435 0.0020803 0.00205625 0.0020322 0.00200815 0.0019841 0.00196005 0.001936
- 0.00191195 0.0018879 0.00186385 0.00183981 0.00181576 0.00179171 0.00176766 0.00174361 0.00171956 0.00169551 0.00167146 0.00164741
- 0.00162336 0.00159931 0.00157526 0.00155121 0.00152716 0.00150311 0.00147906 0.00145501 0.00143096 0.00140691 0.00138286 0.00135881
- 0.00133476 0.00131071 0.00128666 0.00126261 0.00123856 0.00121451 0.00119046 0.00116641 0.00114236 0.00111831 0.00109426 0.00107021
- 0.00104616 0.00102211 0.000998064 0.000974015 0.000949965 0.000925915 0.000901865 0.000877816 0.000853766 0.000829716 0.000805666 0.000781617
- 0.000757567 0.000733517 0.000709467 0.000685418 0.000661368 0.000637318 0.000613268 0.000589219 0.000565169 0.000541119 0.000517069 0.00049302
- 0.00046897 0.00044492 0.00042087 0.000396821 0.000372771 0.000348721 0.000324672 0.000300622 0.000276572 0.000252522 0.000228473 0.000204423
- 0.000180373 0.000156323 0.000132274 0.000108224 8.41741e-05 6.01244e-05 3.60746e-05 1.20249e-05
+ 0.000125533 0.000378711 0.000631602 0.000881784 0.00112799 0.00136848 0.00160133 0.00182478 0.00203727 0.00223747 0.00242431 0.00259927
+ 0.00276214 0.00291258 0.00305041 0.00317565 0.00328848 0.00338919 0.00347821 0.00355604 0.00362324 0.00368045 0.0037283 0.00376747
+ 0.00379885 0.00382391 0.00384356 0.0038583 0.0038686 0.00387488 0.00387755 0.003877 0.00387357 0.00386758 0.00385932 0.00384905
+ 0.00383701 0.00382341 0.00380845 0.00379228 0.00377507 0.00375694 0.00373802 0.00371841 0.0036982 0.00367746 0.00365628 0.0036347
+ 0.0036128 0.0035906 0.00356815 0.00354549 0.00352265 0.00349964 0.0034765 0.00345325 0.0034299 0.00340646 0.00338295 0.00335938
+ 0.00333575 0.00331208 0.00328837 0.00326463 0.00324086 0.00321707 0.00319326 0.00316943 0.00314559 0.00312173 0.00309786 0.00307399
+ 0.00305011 0.00302622 0.00300232 0.00297842 0.00295452 0.00293061 0.0029067 0.00288279 0.00285888 0.00283496 0.00281105 0.00278713
+ 0.00276321 0.00273929 0.00271537 0.00269145 0.00266753 0.00264361 0.00261968 0.00259576 0.00257184 0.00254792 0.00252399 0.00250007
+ 0.00247615 0.00245222 0.0024283 0.00240438 0.00238045 0.00235653 0.0023326 0.00230868 0.00228476 0.00226083 0.00223691 0.00221298
+ 0.00218906 0.00216514 0.00214121 0.00211729 0.00209336 0.00206944 0.00204551 0.00202159 0.00199767 0.00197374 0.00194982 0.00192589
+ 0.00190197 0.00187805 0.00185412 0.0018302 0.00180627 0.00178235 0.00175843 0.0017345 0.00171058 0.00168665 0.00166273 0.0016388
+ 0.00161488 0.00159096 0.00156703 0.00154311 0.00151918 0.00149526 0.00147134 0.00144741 0.00142349 0.00139956 0.00137564 0.00135171
+ 0.00132779 0.00130387 0.00127994 0.00125602 0.00123209 0.00120817 0.00118425 0.00116032 0.0011364 0.00111247 0.00108855 0.00106462
+ 0.0010407 0.00101678 0.000992852 0.000968928 0.000945004 0.00092108 0.000897156 0.000873232 0.000849307 0.000825383 0.000801459 0.000777535
+ 0.000753611 0.000729687 0.000705762 0.000681838 0.000657914 0.00063399 0.000610066 0.000586142 0.000562218 0.000538293 0.000514369 0.000490445
+ 0.000466521 0.000442597 0.000418673 0.000394748 0.000370824 0.0003469 0.000322976 0.000299052 0.000275128 0.000251204 0.000227279 0.000203355
+ 0.000179431 0.000155507 0.000131583 0.000107659 8.37345e-05 5.98104e-05 3.58862e-05 1.19621e-05
</DataArray>
<DataArray type="Float32" Name="x_H2O" NumberOfComponents="1" format="ascii">
0.99999 0.99999 0.99999 0.99999 0.99999 0.99999 0.99999 0.99999 0.99999 0.99999 0.99999 0.99999
@@ -118,31 +118,31 @@
1.55498e-05 1.55498e-05 1.55498e-05 1.55498e-05 1.55498e-05 1.55498e-05 1.55498e-05 1.55498e-05
</DataArray>
<DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
- 996.629 996.789 996.947 997.102 997.237 997.367 997.492 997.612 997.727 997.836 997.939 998.023
- 998.099 998.171 998.237 998.297 998.353 998.404 998.451 998.493 998.531 998.566 998.597 998.624
- 998.649 998.669 998.685 998.699 998.711 998.722 998.731 998.74 998.747 998.753 998.759 998.764
- 998.768 998.772 998.775 998.778 998.78 998.782 998.784 998.786 998.787 998.788 998.789 998.79
- 998.791 998.791 998.792 998.792 998.793 998.793 998.793 998.793 998.794 998.794 998.794 998.794
- 998.794 998.794 998.794 998.794 998.794 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
- 998.795 998.795 998.795 998.795 998.795 998.795 998.795 998.795
+ 996.624 996.784 996.943 997.099 997.234 997.364 997.489 997.609 997.724 997.833 997.936 998.02
+ 998.097 998.168 998.234 998.294 998.35 998.401 998.448 998.49 998.528 998.562 998.593 998.62
+ 998.645 998.665 998.68 998.694 998.706 998.717 998.726 998.735 998.742 998.748 998.754 998.759
+ 998.763 998.767 998.77 998.773 998.775 998.777 998.779 998.78 998.781 998.783 998.784 998.784
+ 998.785 998.786 998.786 998.787 998.787 998.787 998.788 998.788 998.788 998.788 998.788 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
+ 998.789 998.789 998.789 998.789 998.789 998.789 998.789 998.789
</DataArray>
<DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii">
- 0.000860464 0.000871799 0.000883027 0.00089408 0.000906248 0.000918226 0.000929801 0.00094092 0.000951541 0.000961628 0.000971156 0.000981222
- 0.000990825 0.000999761 0.00100804 0.00101567 0.00102267 0.00102907 0.00103491 0.0010402 0.00104499 0.0010493 0.00105317 0.00105664
- 0.00105973 0.00106267 0.0010655 0.00106801 0.00107022 0.00107217 0.00107389 0.0010754 0.00107672 0.00107787 0.00107888 0.00107976
- 0.00108052 0.00108118 0.00108176 0.00108226 0.00108269 0.00108306 0.00108338 0.00108366 0.0010839 0.0010841 0.00108428 0.00108443
- 0.00108456 0.00108467 0.00108477 0.00108485 0.00108492 0.00108498 0.00108503 0.00108507 0.00108511 0.00108514 0.00108516 0.00108519
- 0.00108521 0.00108522 0.00108524 0.00108525 0.00108526 0.00108527 0.00108527 0.00108528 0.00108529 0.00108529 0.00108529 0.0010853
+ 0.000860486 0.000871863 0.000883131 0.000894219 0.000906442 0.000918451 0.000930051 0.00094119 0.000951826 0.000961923 0.000971457 0.000981568
+ 0.000991169 0.0010001 0.00100836 0.00101598 0.00102297 0.00102936 0.00103518 0.00104046 0.00104522 0.00104952 0.00105338 0.00105683
+ 0.00105991 0.00106286 0.00106567 0.00106816 0.00107036 0.0010723 0.001074 0.0010755 0.0010768 0.00107795 0.00107895 0.00107982
+ 0.00108058 0.00108123 0.0010818 0.0010823 0.00108272 0.00108309 0.00108341 0.00108368 0.00108392 0.00108412 0.00108429 0.00108444
+ 0.00108457 0.00108468 0.00108477 0.00108485 0.00108492 0.00108498 0.00108503 0.00108507 0.00108511 0.00108514 0.00108517 0.00108519
+ 0.00108521 0.00108522 0.00108524 0.00108525 0.00108526 0.00108527 0.00108528 0.00108528 0.00108529 0.00108529 0.00108529 0.0010853
0.0010853 0.0010853 0.0010853 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531
0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531
0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531 0.00108531
@@ -175,11 +175,11 @@
0 0 0 0 0 0 0 0
</DataArray>
<DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
- 299.704 299.114 298.529 297.954 297.391 296.844 296.315 295.807 295.322 294.861 294.426 294.017
- 293.634 293.279 292.949 292.645 292.366 292.111 291.879 291.668 291.477 291.306 291.151 291.013
- 290.89 290.78 290.683 290.596 290.52 290.453 290.394 290.342 290.296 290.257 290.222 290.192
- 290.165 290.142 290.123 290.105 290.09 290.078 290.067 290.057 290.049 290.042 290.036 290.03
- 290.026 290.022 290.019 290.016 290.014 290.012 290.01 290.008 290.007 290.006 290.005 290.004
+ 299.703 299.111 298.524 297.947 297.382 296.834 296.304 295.795 295.309 294.848 294.412 294.003
+ 293.621 293.265 292.936 292.633 292.354 292.1 291.868 291.658 291.468 291.297 291.143 291.006
+ 290.883 290.774 290.677 290.591 290.515 290.449 290.39 290.338 290.293 290.254 290.219 290.189
+ 290.163 290.141 290.121 290.104 290.089 290.077 290.066 290.056 290.048 290.041 290.035 290.03
+ 290.026 290.022 290.019 290.016 290.013 290.011 290.01 290.008 290.007 290.006 290.005 290.004
290.004 290.003 290.003 290.002 290.002 290.002 290.001 290.001 290.001 290.001 290.001 290.001
290 290 290 290 290 290 290 290 290 290 290 290
290 290 290 290 290 290 290 290 290 290 290 290
@@ -194,11 +194,11 @@
290 290 290 290 290 290 290 290
</DataArray>
<DataArray type="Float32" Name="temperatureExact" NumberOfComponents="1" format="ascii">
- 299.723 299.171 298.623 298.082 297.549 297.028 296.521 296.029 295.555 295.099 294.664 294.25
- 293.859 293.49 293.145 292.823 292.524 292.248 291.994 291.762 291.55 291.359 291.186 291.031
- 290.892 290.769 290.66 290.564 290.48 290.407 290.344 290.289 290.242 290.201 290.167 290.138
- 290.113 290.093 290.076 290.061 290.05 290.04 290.032 290.025 290.02 290.016 290.013 290.01
- 290.008 290.006 290.005 290.004 290.003 290.002 290.002 290.001 290.001 290.001 290 290
+ 299.723 299.17 298.621 298.079 297.545 297.023 296.515 296.023 295.548 295.092 294.656 294.243
+ 293.851 293.482 293.137 292.815 292.516 292.24 291.986 291.754 291.543 291.352 291.18 291.025
+ 290.887 290.764 290.656 290.56 290.477 290.404 290.341 290.286 290.239 290.199 290.165 290.136
+ 290.112 290.092 290.075 290.061 290.049 290.039 290.031 290.025 290.02 290.016 290.012 290.01
+ 290.008 290.006 290.004 290.003 290.003 290.002 290.002 290.001 290.001 290.001 290 290
290 290 290 290 290 290 290 290 290 290 290 290
290 290 290 290 290 290 290 290 290 290 290 290
290 290 290 290 290 290 290 290 290 290 290 290
diff --git a/test/references/1p2ctestmpfa-reference.vtu b/test/references/1p2ctestmpfa-reference.vtu
new file mode 100644
index 0000000000..284fc0208e
--- /dev/null
+++ b/test/references/1p2ctestmpfa-reference.vtu
@@ -0,0 +1,113 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="40" NumberOfPoints="63">
+ <CellData Scalars="pressure">
+ <DataArray type="Float32" Name="pressure" NumberOfComponents="1" format="ascii">
+ 197500 192500 187500 182501 177501 172501 167501 162501 157501 152501 147501 142501
+ 137501 132501 127501 122501 117501 112500 107500 102500 197500 192500 187500 182501
+ 177501 172501 167501 162501 157501 152501 147501 142501 137501 132501 127501 122501
+ 117501 112500 107500 102500
+ </DataArray>
+ <DataArray type="Float32" Name="delp" NumberOfComponents="1" format="ascii">
+ 97500.1 92500.3 87500.4 82500.6 77500.7 72500.8 67500.8 62500.9 57500.9 52501 47501 42500.9
+ 37500.9 32500.8 27500.8 22500.7 17500.6 12500.4 7500.27 2500.1 97500.1 92500.3 87500.4 82500.6
+ 77500.7 72500.8 67500.8 62500.9 57500.9 52501 47501 42500.9 37500.9 32500.8 27500.8 22500.7
+ 17500.6 12500.4 7500.27 2500.1
+ </DataArray>
+ <DataArray type="Float32" Name="x_H2O" NumberOfComponents="1" format="ascii">
+ 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998
+ 0.99998 0.99998 0.99998 0.999981 0.999981 0.999981 0.999981 0.999981 0.99998 0.99998 0.99998 0.99998
+ 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.99998 0.999981
+ 0.999981 0.999981 0.999981 0.999981
+ </DataArray>
+ <DataArray type="Float32" Name="x_N2" NumberOfComponents="1" format="ascii">
+ 2e-05 1.99999e-05 1.99995e-05 1.99985e-05 1.99962e-05 1.99919e-05 1.99844e-05 1.99724e-05 1.99544e-05 1.99287e-05 1.98936e-05 1.98473e-05
+ 1.97882e-05 1.97146e-05 1.96251e-05 1.95186e-05 1.9394e-05 1.92505e-05 1.90878e-05 1.89055e-05 2e-05 1.99999e-05 1.99995e-05 1.99985e-05
+ 1.99962e-05 1.99919e-05 1.99844e-05 1.99724e-05 1.99544e-05 1.99287e-05 1.98936e-05 1.98473e-05 1.97882e-05 1.97146e-05 1.96251e-05 1.95186e-05
+ 1.9394e-05 1.92505e-05 1.90878e-05 1.89055e-05
+ </DataArray>
+ <DataArray type="Float32" Name="X_H2O" NumberOfComponents="1" format="ascii">
+ 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969
+ 0.999969 0.999969 0.999969 0.99997 0.99997 0.99997 0.99997 0.999971 0.999969 0.999969 0.999969 0.999969
+ 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.999969 0.99997
+ 0.99997 0.99997 0.99997 0.999971
+ </DataArray>
+ <DataArray type="Float32" Name="X_N2" NumberOfComponents="1" format="ascii">
+ 3.10994e-05 3.10992e-05 3.10986e-05 3.1097e-05 3.10935e-05 3.10868e-05 3.10752e-05 3.10565e-05 3.10285e-05 3.09886e-05 3.0934e-05 3.0862e-05
+ 3.077e-05 3.06556e-05 3.05165e-05 3.03508e-05 3.01571e-05 2.9934e-05 2.9681e-05 2.93976e-05 3.10994e-05 3.10992e-05 3.10986e-05 3.1097e-05
+ 3.10935e-05 3.10868e-05 3.10752e-05 3.10565e-05 3.10285e-05 3.09886e-05 3.0934e-05 3.0862e-05 3.077e-05 3.06556e-05 3.05165e-05 3.03508e-05
+ 3.01571e-05 2.9934e-05 2.9681e-05 2.93976e-05
+ </DataArray>
+ <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+ 998.235 998.233 998.231 998.229 998.226 998.224 998.222 998.219 998.217 998.215 998.213 998.21
+ 998.208 998.206 998.203 998.201 998.199 998.197 998.194 998.192 998.235 998.233 998.231 998.229
+ 998.226 998.224 998.222 998.219 998.217 998.215 998.213 998.21 998.208 998.206 998.203 998.201
+ 998.199 998.197 998.194 998.192
+ </DataArray>
+ <DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii">
+ 0.00100296 0.00100296 0.00100296 0.00100296 0.00100297 0.00100297 0.00100297 0.00100297 0.00100297 0.00100297 0.00100297 0.00100298
+ 0.00100298 0.00100298 0.00100298 0.00100298 0.00100298 0.00100299 0.00100299 0.00100299 0.00100296 0.00100296 0.00100296 0.00100296
+ 0.00100297 0.00100297 0.00100297 0.00100297 0.00100297 0.00100297 0.00100297 0.00100298 0.00100298 0.00100298 0.00100298 0.00100298
+ 0.00100298 0.00100299 0.00100299 0.00100299
+ </DataArray>
+ <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 0.05 0 0 0 0.5 0 0.05 0.5 0
+ 0.1 0 0 0.1 0.5 0 0.15 0 0 0.15 0.5 0
+ 0.2 0 0 0.2 0.5 0 0.25 0 0 0.25 0.5 0
+ 0.3 0 0 0.3 0.5 0 0.35 0 0 0.35 0.5 0
+ 0.4 0 0 0.4 0.5 0 0.45 0 0 0.45 0.5 0
+ 0.5 0 0 0.5 0.5 0 0.55 0 0 0.55 0.5 0
+ 0.6 0 0 0.6 0.5 0 0.65 0 0 0.65 0.5 0
+ 0.7 0 0 0.7 0.5 0 0.75 0 0 0.75 0.5 0
+ 0.8 0 0 0.8 0.5 0 0.85 0 0 0.85 0.5 0
+ 0.9 0 0 0.9 0.5 0 0.95 0 0 0.95 0.5 0
+ 1 0 0 1 0.5 0 0 1 0 0.05 1 0
+ 0.1 1 0 0.15 1 0 0.2 1 0 0.25 1 0
+ 0.3 1 0 0.35 1 0 0.4 1 0 0.45 1 0
+ 0.5 1 0 0.55 1 0 0.6 1 0 0.65 1 0
+ 0.7 1 0 0.75 1 0 0.8 1 0 0.85 1 0
+ 0.9 1 0 0.95 1 0 1 1 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 24 26 27 25 26 28 29 27 28 30 31 29
+ 30 32 33 31 32 34 35 33 34 36 37 35
+ 36 38 39 37 38 40 41 39 2 3 43 42
+ 3 5 44 43 5 7 45 44 7 9 46 45
+ 9 11 47 46 11 13 48 47 13 15 49 48
+ 15 17 50 49 17 19 51 50 19 21 52 51
+ 21 23 53 52 23 25 54 53 25 27 55 54
+ 27 29 56 55 29 31 57 56 31 33 58 57
+ 33 35 59 58 35 37 60 59 37 39 61 60
+ 39 41 62 61
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
--
GitLab