From bf1a4aed05dc2d53ebe9696da9cd3c5d4170af67 Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?Dennis=20Gl=C3=A4ser?= <dennis.glaeser@iws.uni-stuttgart.de>
Date: Fri, 23 Nov 2018 12:15:17 +0100
Subject: [PATCH] [test][1p] test mpfa scv velocity reconstruction

---
 .../1p/implicit/incompressible/CMakeLists.txt |  4 +-
 .../1p/implicit/incompressible/main.cc        | 32 +++++++++
 .../test_1p_mpfa_scvvelocities-reference.vtp  | 67 +++++++++++++++++++
 3 files changed, 102 insertions(+), 1 deletion(-)
 create mode 100644 test/references/test_1p_mpfa_scvvelocities-reference.vtp

diff --git a/test/porousmediumflow/1p/implicit/incompressible/CMakeLists.txt b/test/porousmediumflow/1p/implicit/incompressible/CMakeLists.txt
index 30dae4f00b..10477bb979 100644
--- a/test/porousmediumflow/1p/implicit/incompressible/CMakeLists.txt
+++ b/test/porousmediumflow/1p/implicit/incompressible/CMakeLists.txt
@@ -41,7 +41,9 @@ dune_add_test(NAME test_1p_incompressible_mpfa
               CMD_ARGS  --script fuzzy
                         --files ${CMAKE_SOURCE_DIR}/test/references/test_1p_cc-reference.vtu
                                 ${CMAKE_CURRENT_BINARY_DIR}/test_1p_incompressible_mpfa-00001.vtu
-                        --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p_incompressible_mpfa_anadiff params.input -Problem.Name test_1p_incompressible_mpfa")
+                                ${CMAKE_SOURCE_DIR}/test/references/test_1p_mpfa_scvvelocities-reference.vtp
+                                ${CMAKE_CURRENT_BINARY_DIR}/mpfa_scv_velocities.vtp
+                        --command "${CMAKE_CURRENT_BINARY_DIR}/test_1p_incompressible_mpfa_anadiff params.input -Problem.Name test_1p_incompressible_mpfa -IO.WriteMpfaVelocities true")
 
 # using box and analytical Jacobian
 dune_add_test(NAME test_1p_incompressible_box
diff --git a/test/porousmediumflow/1p/implicit/incompressible/main.cc b/test/porousmediumflow/1p/implicit/incompressible/main.cc
index bb67afdb45..a3bee4023d 100644
--- a/test/porousmediumflow/1p/implicit/incompressible/main.cc
+++ b/test/porousmediumflow/1p/implicit/incompressible/main.cc
@@ -45,11 +45,39 @@
 
 #include <dumux/io/vtkoutputmodule.hh>
 #include <dumux/io/grid/gridmanager.hh>
+#include <dumux/io/pointcloudvtkwriter.hh>
+
+#include <dumux/discretization/method.hh>
+#include <dumux/discretization/cellcentered/mpfa/scvgradients.hh>
 
 #include <dumux/assembly/fvassembler.hh>
 
 #include "problem.hh"
 
+//! function to write out the scv-wise velocities (overload for mpfa)
+template<class FVGridGeometry, class GridVariables, class Sol,
+         std::enable_if_t<FVGridGeometry::discMethod == Dumux::DiscretizationMethod::ccmpfa, int> = 0>
+void writeMpfaVelocities(const FVGridGeometry& fvGridGeometry,
+                         const GridVariables& gridVariables,
+                         const Sol& x)
+{
+    using Scalar = typename GridVariables::Scalar;
+    using GlobalPos = typename FVGridGeometry::SubControlVolume::GlobalPosition;
+
+    const auto velocities = Dumux::CCMpfaScvGradients::computeVelocities(fvGridGeometry, gridVariables, x, /*phaseIdx*/0);
+    Dumux::PointCloudVtkWriter<Scalar, GlobalPos> writer(velocities.first);
+    writer.addPointData(velocities.second, "velocity (m/s)");
+    writer.write("mpfa_scv_velocities");
+}
+
+//! function to write out the scv-wise velocities (overload for NOT mpfa)
+template<class FVGridGeometry, class GridVariables, class Sol,
+         std::enable_if_t<FVGridGeometry::discMethod != Dumux::DiscretizationMethod::ccmpfa, int> = 0>
+void writeMpfaVelocities(const FVGridGeometry& fvGridGeometry,
+                         const GridVariables& gridVariables,
+                         const Sol& x)
+{}
+
 int main(int argc, char** argv) try
 {
     using namespace Dumux;
@@ -191,6 +219,10 @@ int main(int argc, char** argv) try
                   << comm.size() << " processes.\n"
                   << "The cumulative CPU time was " << timer.elapsed()*comm.size() << " seconds.\n";
 
+    // For the mpfa test, write out the gradients in the scv centers
+    if (getParam<bool>("IO.WriteMpfaVelocities", false))
+        writeMpfaVelocities(*fvGridGeometry, *gridVariables, x);
+
     if (mpiHelper.rank() == 0)
         Parameters::print();
 
diff --git a/test/references/test_1p_mpfa_scvvelocities-reference.vtp b/test/references/test_1p_mpfa_scvvelocities-reference.vtp
new file mode 100644
index 0000000000..da52ede32f
--- /dev/null
+++ b/test/references/test_1p_mpfa_scvvelocities-reference.vtp
@@ -0,0 +1,67 @@
+<?xml version="1.0"?>
+<VTKFile type="PolyData" version="0.1" byte_order="LittleEndian">
+<PolyData>
+<Piece NumberOfLines="0" NumberOfPoints="400">
+<Points>
+<DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+0.025 0.025 0 0.025 0.075 0 0.025 0.125 0 0.075 0.025 0 0.125 0.025 0 0.075 0.075 0 0.125 0.075 0 0.075 0.125 0 0.125 0.125 0 0.175 0.025 0 0.225 0.025 0 0.175 0.075 0 0.225 0.075 0 0.175 0.125 0 0.225 0.125 0 0.275 0.025 0 
+0.325 0.025 0 0.275 0.075 0 0.325 0.075 0 0.275 0.125 0 0.325 0.125 0 0.375 0.025 0 0.425 0.025 0 0.375 0.075 0 0.425 0.075 0 0.375 0.125 0 0.425 0.125 0 0.475 0.025 0 0.525 0.025 0 0.475 0.075 0 0.525 0.075 0 0.475 0.125 0 
+0.525 0.125 0 0.575 0.025 0 0.625 0.025 0 0.575 0.075 0 0.625 0.075 0 0.575 0.125 0 0.625 0.125 0 0.675 0.025 0 0.725 0.025 0 0.675 0.075 0 0.725 0.075 0 0.675 0.125 0 0.725 0.125 0 0.775 0.025 0 0.825 0.025 0 0.775 0.075 0 
+0.825 0.075 0 0.775 0.125 0 0.825 0.125 0 0.875 0.025 0 0.925 0.025 0 0.875 0.075 0 0.925 0.075 0 0.875 0.125 0 0.925 0.125 0 0.975 0.025 0 0.975 0.075 0 0.975 0.125 0 0.025 0.175 0 0.025 0.225 0 0.075 0.175 0 0.125 0.175 0 
+0.075 0.225 0 0.125 0.225 0 0.175 0.175 0 0.225 0.175 0 0.175 0.225 0 0.225 0.225 0 0.275 0.175 0 0.325 0.175 0 0.275 0.225 0 0.325 0.225 0 0.375 0.175 0 0.425 0.175 0 0.375 0.225 0 0.425 0.225 0 0.475 0.175 0 0.525 0.175 0 
+0.475 0.225 0 0.525 0.225 0 0.575 0.175 0 0.625 0.175 0 0.575 0.225 0 0.625 0.225 0 0.675 0.175 0 0.725 0.175 0 0.675 0.225 0 0.725 0.225 0 0.775 0.175 0 0.825 0.175 0 0.775 0.225 0 0.825 0.225 0 0.875 0.175 0 0.925 0.175 0 
+0.875 0.225 0 0.925 0.225 0 0.975 0.175 0 0.975 0.225 0 0.025 0.275 0 0.025 0.325 0 0.075 0.275 0 0.125 0.275 0 0.075 0.325 0 0.125 0.325 0 0.175 0.275 0 0.225 0.275 0 0.175 0.325 0 0.225 0.325 0 0.275 0.275 0 0.325 0.275 0 
+0.275 0.325 0 0.325 0.325 0 0.375 0.275 0 0.425 0.275 0 0.375 0.325 0 0.425 0.325 0 0.475 0.275 0 0.525 0.275 0 0.475 0.325 0 0.525 0.325 0 0.575 0.275 0 0.625 0.275 0 0.575 0.325 0 0.625 0.325 0 0.675 0.275 0 0.725 0.275 0 
+0.675 0.325 0 0.725 0.325 0 0.775 0.275 0 0.825 0.275 0 0.775 0.325 0 0.825 0.325 0 0.875 0.275 0 0.925 0.275 0 0.875 0.325 0 0.925 0.325 0 0.975 0.275 0 0.975 0.325 0 0.025 0.375 0 0.025 0.425 0 0.075 0.375 0 0.125 0.375 0 
+0.075 0.425 0 0.125 0.425 0 0.175 0.375 0 0.225 0.375 0 0.175 0.425 0 0.225 0.425 0 0.275 0.375 0 0.325 0.375 0 0.275 0.425 0 0.325 0.425 0 0.375 0.375 0 0.425 0.375 0 0.375 0.425 0 0.425 0.425 0 0.475 0.375 0 0.525 0.375 0 
+0.475 0.425 0 0.525 0.425 0 0.575 0.375 0 0.625 0.375 0 0.575 0.425 0 0.625 0.425 0 0.675 0.375 0 0.725 0.375 0 0.675 0.425 0 0.725 0.425 0 0.775 0.375 0 0.825 0.375 0 0.775 0.425 0 0.825 0.425 0 0.875 0.375 0 0.925 0.375 0 
+0.875 0.425 0 0.925 0.425 0 0.975 0.375 0 0.975 0.425 0 0.025 0.475 0 0.025 0.525 0 0.075 0.475 0 0.125 0.475 0 0.075 0.525 0 0.125 0.525 0 0.175 0.475 0 0.225 0.475 0 0.175 0.525 0 0.225 0.525 0 0.275 0.475 0 0.325 0.475 0 
+0.275 0.525 0 0.325 0.525 0 0.375 0.475 0 0.425 0.475 0 0.375 0.525 0 0.425 0.525 0 0.475 0.475 0 0.525 0.475 0 0.475 0.525 0 0.525 0.525 0 0.575 0.475 0 0.625 0.475 0 0.575 0.525 0 0.625 0.525 0 0.675 0.475 0 0.725 0.475 0 
+0.675 0.525 0 0.725 0.525 0 0.775 0.475 0 0.825 0.475 0 0.775 0.525 0 0.825 0.525 0 0.875 0.475 0 0.925 0.475 0 0.875 0.525 0 0.925 0.525 0 0.975 0.475 0 0.975 0.525 0 0.025 0.575 0 0.025 0.625 0 0.075 0.575 0 0.125 0.575 0 
+0.075 0.625 0 0.125 0.625 0 0.175 0.575 0 0.225 0.575 0 0.175 0.625 0 0.225 0.625 0 0.275 0.575 0 0.325 0.575 0 0.275 0.625 0 0.325 0.625 0 0.375 0.575 0 0.425 0.575 0 0.375 0.625 0 0.425 0.625 0 0.475 0.575 0 0.525 0.575 0 
+0.475 0.625 0 0.525 0.625 0 0.575 0.575 0 0.625 0.575 0 0.575 0.625 0 0.625 0.625 0 0.675 0.575 0 0.725 0.575 0 0.675 0.625 0 0.725 0.625 0 0.775 0.575 0 0.825 0.575 0 0.775 0.625 0 0.825 0.625 0 0.875 0.575 0 0.925 0.575 0 
+0.875 0.625 0 0.925 0.625 0 0.975 0.575 0 0.975 0.625 0 0.025 0.675 0 0.025 0.725 0 0.075 0.675 0 0.125 0.675 0 0.075 0.725 0 0.125 0.725 0 0.175 0.675 0 0.225 0.675 0 0.175 0.725 0 0.225 0.725 0 0.275 0.675 0 0.325 0.675 0 
+0.275 0.725 0 0.325 0.725 0 0.375 0.675 0 0.425 0.675 0 0.375 0.725 0 0.425 0.725 0 0.475 0.675 0 0.525 0.675 0 0.475 0.725 0 0.525 0.725 0 0.575 0.675 0 0.625 0.675 0 0.575 0.725 0 0.625 0.725 0 0.675 0.675 0 0.725 0.675 0 
+0.675 0.725 0 0.725 0.725 0 0.775 0.675 0 0.825 0.675 0 0.775 0.725 0 0.825 0.725 0 0.875 0.675 0 0.925 0.675 0 0.875 0.725 0 0.925 0.725 0 0.975 0.675 0 0.975 0.725 0 0.025 0.775 0 0.025 0.825 0 0.075 0.775 0 0.125 0.775 0 
+0.075 0.825 0 0.125 0.825 0 0.175 0.775 0 0.225 0.775 0 0.175 0.825 0 0.225 0.825 0 0.275 0.775 0 0.325 0.775 0 0.275 0.825 0 0.325 0.825 0 0.375 0.775 0 0.425 0.775 0 0.375 0.825 0 0.425 0.825 0 0.475 0.775 0 0.525 0.775 0 
+0.475 0.825 0 0.525 0.825 0 0.575 0.775 0 0.625 0.775 0 0.575 0.825 0 0.625 0.825 0 0.675 0.775 0 0.725 0.775 0 0.675 0.825 0 0.725 0.825 0 0.775 0.775 0 0.825 0.775 0 0.775 0.825 0 0.825 0.825 0 0.875 0.775 0 0.925 0.775 0 
+0.875 0.825 0 0.925 0.825 0 0.975 0.775 0 0.975 0.825 0 0.025 0.875 0 0.025 0.925 0 0.075 0.875 0 0.125 0.875 0 0.075 0.925 0 0.125 0.925 0 0.175 0.875 0 0.225 0.875 0 0.175 0.925 0 0.225 0.925 0 0.275 0.875 0 0.325 0.875 0 
+0.275 0.925 0 0.325 0.925 0 0.375 0.875 0 0.425 0.875 0 0.375 0.925 0 0.425 0.925 0 0.475 0.875 0 0.525 0.875 0 0.475 0.925 0 0.525 0.925 0 0.575 0.875 0 0.625 0.875 0 0.575 0.925 0 0.625 0.925 0 0.675 0.875 0 0.725 0.875 0 
+0.675 0.925 0 0.725 0.925 0 0.775 0.875 0 0.825 0.875 0 0.775 0.925 0 0.825 0.925 0 0.875 0.875 0 0.925 0.875 0 0.875 0.925 0 0.925 0.925 0 0.975 0.875 0 0.975 0.925 0 0.025 0.975 0 0.075 0.975 0 0.125 0.975 0 0.175 0.975 0 
+0.225 0.975 0 0.275 0.975 0 0.325 0.975 0 0.375 0.975 0 0.425 0.975 0 0.475 0.975 0 0.525 0.975 0 0.575 0.975 0 0.625 0.975 0 0.675 0.975 0 0.725 0.975 0 0.775 0.975 0 0.825 0.975 0 0.875 0.975 0 0.925 0.975 0 0.975 0.975 0 
+
+</DataArray>
+</Points>
+<PointData Vectors="gradP">
+<DataArray type="Float32" Name="gradP" NumberOfComponents="3" format="ascii">
+0 -0.00817984 0 0 -0.00873783 0 0 -0.00873783 0 0.000557991 -0.00817984 0 0.000557991 -0.00706385 0 0.000557991 -0.00873783 0 0.000557991 -0.0077827 0 0.00151312 -0.00873783 0 0.00151312 -0.0077827 0 0.00127684 -0.00706385 0 0.00127684 -0.00451017 0 0.00127684 -0.0077827 0 0.00127684 -0.00400127 0 0.00505827 -0.0077827 0 0.00505827 -0.00400127 0 0.000767938 -0.00451017 0 
+0.000767938 -0.0029743 0 0.000767938 -0.00400127 0 0.000767938 -0.00253139 0 0.00223782 -0.00400127 0 0.00223782 -0.00253139 0 0.00032503 -0.0029743 0 0.00032503 -0.00232424 0 0.00032503 -0.00253139 0 0.00032503 -0.00199921 0 0.000857211 -0.00253139 0 0.000857211 -0.00199921 0 -4.65661e-16 -0.00232424 0 -4.65661e-16 -0.00232424 0 -4.65661e-16 -0.00199921 0 -4.65661e-16 -0.00199921 0 4.07454e-16 -0.00199921 0 
+4.07454e-16 -0.00199921 0 -0.00032503 -0.00232424 0 -0.00032503 -0.0029743 0 -0.00032503 -0.00199921 0 -0.00032503 -0.00253139 0 -0.000857211 -0.00199921 0 -0.000857211 -0.00253139 0 -0.000767938 -0.0029743 0 -0.000767938 -0.00451017 0 -0.000767938 -0.00253139 0 -0.000767938 -0.00400127 0 -0.00223782 -0.00253139 0 -0.00223782 -0.00400127 0 -0.00127684 -0.00451017 0 -0.00127684 -0.00706385 0 -0.00127684 -0.00400127 0 
+-0.00127684 -0.0077827 0 -0.00505827 -0.00400127 0 -0.00505827 -0.0077827 0 -0.000557991 -0.00706385 0 -0.000557991 -0.00817984 0 -0.000557991 -0.0077827 0 -0.000557991 -0.00873783 0 -0.00151312 -0.0077827 0 -0.00151312 -0.00873783 0 0 -0.00817984 0 0 -0.00873783 0 0 -0.00873783 0 0 -0.0102509 0 0 -0.0102509 0 0.00151312 -0.0102509 0 0.00151312 -0.0113279 0 
+0.000436197 -0.0102509 0 0.000436197 -0.0113279 0 0.00505827 -0.0113279 0 0.00505827 -0.000219242 0 0.000105236 -0.0113279 0 0.000105236 -0.000219242 0 0.00223782 -0.000219242 0 0.00223782 -0.000189209 0 3.75448e-05 -0.000219242 0 3.75448e-05 -0.000189209 0 0.000857211 -0.000189209 0 0.000857211 -0.000180423 0 1.3009e-05 -0.000189209 0 1.3009e-05 -0.000180423 0 4.90109e-14 -0.000180423 0 -4.81959e-14 -0.000180423 0 
+6.35046e-16 -0.000180423 0 -3.35276e-16 -0.000180423 0 -0.000857211 -0.000180423 0 -0.000857211 -0.000189209 0 -1.3009e-05 -0.000180423 0 -1.3009e-05 -0.000189209 0 -0.00223782 -0.000189209 0 -0.00223782 -0.000219242 0 -3.75448e-05 -0.000189209 0 -3.75448e-05 -0.000219242 0 -0.00505827 -0.000219242 0 -0.00505827 -0.0113279 0 -0.000105236 -0.000219242 0 -0.000105236 -0.0113279 0 -0.00151312 -0.0113279 0 -0.00151312 -0.0102509 0 
+-0.000436197 -0.0113279 0 -0.000436197 -0.0102509 0 0 -0.0102509 0 0 -0.0102509 0 0 -0.0106871 0 -5.82077e-17 -0.0106871 0 0.000436197 -0.0106871 0 0.000436197 -0.0109969 0 0.000126437 -0.0106871 0 0.000126437 -0.0109969 0 0.000105236 -0.0109969 0 0.000105236 -0.000141936 0 4.19355e-05 -0.0109969 0 4.19355e-05 -0.000141936 0 3.75448e-05 -0.000141936 0 3.75448e-05 -0.000155058 0 
+2.44225e-05 -0.000141936 0 2.44225e-05 -0.000155058 0 1.3009e-05 -0.000155058 0 1.3009e-05 -0.000157799 0 1.02684e-05 -0.000155058 0 1.02684e-05 -0.000157799 0 1.49594e-16 -0.000157799 0 1.50176e-16 -0.000157799 0 9.02219e-17 -0.000157799 0 8.96398e-17 -0.000157799 0 -1.3009e-05 -0.000157799 0 -1.3009e-05 -0.000155058 0 -1.02684e-05 -0.000157799 0 -1.02684e-05 -0.000155058 0 -3.75448e-05 -0.000155058 0 -3.75448e-05 -0.000141936 0 
+-2.44225e-05 -0.000155058 0 -2.44225e-05 -0.000141936 0 -0.000105236 -0.000141936 0 -0.000105236 -0.0109969 0 -4.19355e-05 -0.000141936 0 -4.19355e-05 -0.0109969 0 -0.000436197 -0.0109969 0 -0.000436197 -0.0106871 0 -0.000126437 -0.0109969 0 -0.000126437 -0.0106871 0 0 -0.0106871 0 0 -0.0106871 0 0 -0.0108136 0 0 -0.0108136 0 0.000126437 -0.0108136 0 0.000126437 -0.0109124 0 
+2.76157e-05 -0.0108136 0 2.76157e-05 -0.0109124 0 4.19355e-05 -0.0109124 0 4.19355e-05 -0.000124422 0 1.16415e-05 -0.0109124 0 1.16415e-05 -0.000124422 0 2.44225e-05 -0.000124422 0 2.44225e-05 -0.000140904 0 7.94117e-06 -0.000124422 0 7.94117e-06 -0.000140904 0 1.02684e-05 -0.000140904 0 1.02684e-05 -0.00014753 0 3.64191e-06 -0.000140904 0 3.64191e-06 -0.00014753 0 9.02219e-17 -0.00014753 0 8.96398e-17 -0.00014753 0 
+1.49012e-16 -0.00014753 0 1.4843e-16 -0.00014753 0 -1.02684e-05 -0.00014753 0 -1.02684e-05 -0.000140904 0 -3.64191e-06 -0.00014753 0 -3.64191e-06 -0.000140904 0 -2.44225e-05 -0.000140904 0 -2.44225e-05 -0.000124422 0 -7.94117e-06 -0.000140904 0 -7.94117e-06 -0.000124422 0 -4.19355e-05 -0.000124422 0 -4.19355e-05 -0.0109124 0 -1.16415e-05 -0.000124422 0 -1.16415e-05 -0.0109124 0 -0.000126437 -0.0109124 0 -0.000126437 -0.0108136 0 
+-2.76157e-05 -0.0109124 0 -2.76157e-05 -0.0108136 0 0 -0.0108136 0 0 -0.0108136 0 -5.82077e-17 -0.0108412 0 0 -0.0108412 0 2.76157e-05 -0.0108412 0 2.76157e-05 -0.0108964 0 -2.76157e-05 -0.0108412 0 -2.76157e-05 -0.0108964 0 1.16415e-05 -0.0108964 0 1.16415e-05 -0.000120722 0 -1.16415e-05 -0.0108964 0 -1.16415e-05 -0.000120722 0 7.94117e-06 -0.000120722 0 7.94117e-06 -0.000136604 0 
+-7.94117e-06 -0.000120722 0 -7.94117e-06 -0.000136604 0 3.64191e-06 -0.000136604 0 3.64191e-06 -0.000143888 0 -3.64191e-06 -0.000136604 0 -3.64191e-06 -0.000143888 0 1.49012e-16 -0.000143888 0 1.4843e-16 -0.000143888 0 2.37487e-16 -0.000143888 0 2.38651e-16 -0.000143888 0 -3.64191e-06 -0.000143888 0 -3.64191e-06 -0.000136604 0 3.64191e-06 -0.000143888 0 3.64191e-06 -0.000136604 0 -7.94117e-06 -0.000136604 0 -7.94117e-06 -0.000120722 0 
+7.94117e-06 -0.000136604 0 7.94117e-06 -0.000120722 0 -1.16415e-05 -0.000120722 0 -1.16415e-05 -0.0108964 0 1.16415e-05 -0.000120722 0 1.16415e-05 -0.0108964 0 -2.76157e-05 -0.0108964 0 -2.76157e-05 -0.0108412 0 2.76157e-05 -0.0108964 0 2.76157e-05 -0.0108412 0 -5.82077e-17 -0.0108412 0 0 -0.0108412 0 0 -0.0108136 0 0 -0.0108136 0 -2.76157e-05 -0.0108136 0 -2.76157e-05 -0.0109124 0 
+-0.000126437 -0.0108136 0 -0.000126437 -0.0109124 0 -1.16415e-05 -0.0109124 0 -1.16415e-05 -0.000124422 0 -4.19355e-05 -0.0109124 0 -4.19355e-05 -0.000124422 0 -7.94117e-06 -0.000124422 0 -7.94117e-06 -0.000140904 0 -2.44225e-05 -0.000124422 0 -2.44225e-05 -0.000140904 0 -3.64191e-06 -0.000140904 0 -3.64191e-06 -0.00014753 0 -1.02684e-05 -0.000140904 0 -1.02684e-05 -0.00014753 0 2.37487e-16 -0.00014753 0 2.38651e-16 -0.00014753 0 
+1.74623e-16 -0.00014753 0 1.74914e-16 -0.00014753 0 3.64191e-06 -0.00014753 0 3.64191e-06 -0.000140904 0 1.02684e-05 -0.00014753 0 1.02684e-05 -0.000140904 0 7.94117e-06 -0.000140904 0 7.94117e-06 -0.000124422 0 2.44225e-05 -0.000140904 0 2.44225e-05 -0.000124422 0 1.16415e-05 -0.000124422 0 1.16415e-05 -0.0109124 0 4.19355e-05 -0.000124422 0 4.19355e-05 -0.0109124 0 2.76157e-05 -0.0109124 0 2.76157e-05 -0.0108136 0 
+0.000126437 -0.0109124 0 0.000126437 -0.0108136 0 0 -0.0108136 0 0 -0.0108136 0 0 -0.0106871 0 0 -0.0106871 0 -0.000126437 -0.0106871 0 -0.000126437 -0.0109969 0 -0.000436197 -0.0106871 0 -0.000436197 -0.0109969 0 -4.19355e-05 -0.0109969 0 -4.19355e-05 -0.000141936 0 -0.000105236 -0.0109969 0 -0.000105236 -0.000141936 0 -2.44225e-05 -0.000141936 0 -2.44225e-05 -0.000155058 0 
+-3.75448e-05 -0.000141936 0 -3.75448e-05 -0.000155058 0 -1.02684e-05 -0.000155058 0 -1.02684e-05 -0.000157799 0 -1.3009e-05 -0.000155058 0 -1.3009e-05 -0.000157799 0 1.75205e-16 -0.000157799 0 1.74332e-16 -0.000157799 0 8.11997e-17 -0.000157799 0 8.09086e-17 -0.000157799 0 1.02684e-05 -0.000157799 0 1.02684e-05 -0.000155058 0 1.3009e-05 -0.000157799 0 1.3009e-05 -0.000155058 0 2.44225e-05 -0.000155058 0 2.44225e-05 -0.000141936 0 
+3.75448e-05 -0.000155058 0 3.75448e-05 -0.000141936 0 4.19355e-05 -0.000141936 0 4.19355e-05 -0.0109969 0 0.000105236 -0.000141936 0 0.000105236 -0.0109969 0 0.000126437 -0.0109969 0 0.000126437 -0.0106871 0 0.000436197 -0.0109969 0 0.000436197 -0.0106871 0 -5.82077e-17 -0.0106871 0 -2.91038e-17 -0.0106871 0 0 -0.0102509 0 0 -0.0102509 0 -0.000436197 -0.0102509 0 -0.000436197 -0.0113279 0 
+-0.00151312 -0.0102509 0 -0.00151312 -0.0113279 0 -0.000105236 -0.0113279 0 -0.000105236 -0.000219242 0 -0.00505827 -0.0113279 0 -0.00505827 -0.000219242 0 -3.75448e-05 -0.000219242 0 -3.75448e-05 -0.000189209 0 -0.00223782 -0.000219242 0 -0.00223782 -0.000189209 0 -1.3009e-05 -0.000189209 0 -1.3009e-05 -0.000180423 0 -0.000857211 -0.000189209 0 -0.000857211 -0.000180423 0 -4.04834e-16 -0.000180423 0 5.66943e-16 -0.000180423 0 
+-4.83706e-14 -0.000180423 0 4.87489e-14 -0.000180423 0 1.3009e-05 -0.000180423 0 1.3009e-05 -0.000189209 0 0.000857211 -0.000180423 0 0.000857211 -0.000189209 0 3.75448e-05 -0.000189209 0 3.75448e-05 -0.000219242 0 0.00223782 -0.000189209 0 0.00223782 -0.000219242 0 0.000105236 -0.000219242 0 0.000105236 -0.0113279 0 0.00505827 -0.000219242 0 0.00505827 -0.0113279 0 0.000436197 -0.0113279 0 0.000436197 -0.0102509 0 
+0.00151312 -0.0113279 0 0.00151312 -0.0102509 0 0 -0.0102509 0 0 -0.0102509 0 0 -0.00873783 0 0 -0.00873783 0 -0.00151312 -0.00873783 0 -0.00151312 -0.0077827 0 -0.000557991 -0.00873783 0 -0.000557991 -0.0077827 0 -0.00505827 -0.0077827 0 -0.00505827 -0.00400127 0 -0.00127684 -0.0077827 0 -0.00127684 -0.00400127 0 -0.00223782 -0.00400127 0 -0.00223782 -0.00253139 0 
+-0.000767938 -0.00400127 0 -0.000767938 -0.00253139 0 -0.000857211 -0.00253139 0 -0.000857211 -0.00199921 0 -0.00032503 -0.00253139 0 -0.00032503 -0.00199921 0 2.03727e-16 -0.00199921 0 1.74623e-16 -0.00199921 0 1.45519e-16 -0.00199921 0 1.16415e-16 -0.00199921 0 0.000857211 -0.00199921 0 0.000857211 -0.00253139 0 0.00032503 -0.00199921 0 0.00032503 -0.00253139 0 0.00223782 -0.00253139 0 0.00223782 -0.00400127 0 
+0.000767938 -0.00253139 0 0.000767938 -0.00400127 0 0.00505827 -0.00400127 0 0.00505827 -0.0077827 0 0.00127684 -0.00400127 0 0.00127684 -0.0077827 0 0.00151312 -0.0077827 0 0.00151312 -0.00873783 0 0.000557991 -0.0077827 0 0.000557991 -0.00873783 0 -2.91038e-17 -0.00873783 0 -2.91038e-17 -0.00873783 0 0 -0.00817984 0 -0.000557991 -0.00817984 0 -0.000557991 -0.00706385 0 -0.00127684 -0.00706385 0 
+-0.00127684 -0.00451017 0 -0.000767938 -0.00451017 0 -0.000767938 -0.0029743 0 -0.00032503 -0.0029743 0 -0.00032503 -0.00232424 0 1.16415e-16 -0.00232424 0 1.45519e-16 -0.00232424 0 0.00032503 -0.00232424 0 0.00032503 -0.0029743 0 0.000767938 -0.0029743 0 0.000767938 -0.00451017 0 0.00127684 -0.00451017 0 0.00127684 -0.00706385 0 0.000557991 -0.00706385 0 0.000557991 -0.00817984 0 0 -0.00817984 0 
+
+</DataArray>
+</PointData>
+</Piece>
+</PolyData>
+</VTKFile>
\ No newline at end of file
-- 
GitLab