Merge branch 'fix/multidomainNormalMomentum' into 'master'

Fix/multidomain normal momentum

See merge request !252

dumux/multidomain/2cnistokes2p2cni/stokesncnicouplinglocalresidual.hh

@@ -147,10 +147,8 @@ public:
                     // if p.n comes from the pm
                     if (bcTypes.isCouplingNeumann(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx))
                     {
-                        DimVector pressureCorrection(intersection.centerUnitOuterNormal());
-                        pressureCorrection *= volVars.pressure();
-                        this->residual_[scvIdx][momentumYIdx] += pressureCorrection[momentumYIdx]
-                                                                  * boundaryVars.face().area;
+                        this->residual_[scvIdx][momentumYIdx] -= volVars.pressure()
+                                                                 * boundaryVars.face().normal[momentumYIdx];
                         Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]);
                     }
 
@@ -158,16 +156,13 @@ public:
                     for (int compIdx = 0; compIdx < numComponents; compIdx++)
                     {
                         int eqIdx =  conti0EqIdx + compIdx;
-                        if (eqIdx != massBalanceIdx)
+                        if ((eqIdx != massBalanceIdx) && bcTypes.isCouplingDirichlet(eqIdx))
                         {
-                            if (bcTypes.isCouplingDirichlet(eqIdx))
-                            {
-                                if(useMoles)
-                                    this->residual_[scvIdx][eqIdx] = volVars.moleFraction(compIdx);
-                                else
-                                    this->residual_[scvIdx][eqIdx] = volVars.massFraction(compIdx);
-                                Valgrind::CheckDefined(this->residual_[scvIdx][compIdx]);
-                            }
+                            if(useMoles)
+                                this->residual_[scvIdx][eqIdx] = volVars.moleFraction(compIdx);
+                            else
+                                this->residual_[scvIdx][eqIdx] = volVars.massFraction(compIdx);
+                            Valgrind::CheckDefined(this->residual_[scvIdx][compIdx]);
                         }
                     }
 

dumux/multidomain/2cstokes2p2c/localoperator.hh

@@ -432,6 +432,11 @@ public:
         const Scalar normalMassFlux1 = boundaryVars1.normalVelocity()
                                        * cParams.elemVolVarsCur1[vertInElem1].density();
 
+        if (std::abs(bfNormal1[1]) < 1e-10)
+        {
+            DUNE_THROW(Dune::NotImplemented, "The coupling conditions are not implemented for vertical interfaces.");
+        }
+
         // MASS Balance
         // Neumann-like conditions
         if (cParams.boundaryTypes1.isCouplingNeumann(massBalanceIdx1))
@@ -517,8 +522,8 @@ public:
             // p*A as condition for free flow
             // pressure correction is done in stokeslocalresidual.hh
             couplingRes1.accumulate(lfsu1.child(momentumYIdx1), vertInElem1,
-                                    cParams.elemVolVarsCur2[vertInElem2].pressure(nPhaseIdx2) *
-                                    boundaryVars2.face().area);
+                                    cParams.elemVolVarsCur2[vertInElem2].pressure(nPhaseIdx2)
+                                    * bfNormal1[momentumYIdx1]);
         }
 
         // Dirichlet-like conditions

dumux/multidomain/2cstokes2p2c/stokesnccouplinglocalresidual.hh

@@ -157,10 +157,8 @@ public:
                     // if p.n comes from the pm
                     if (bcTypes.isCouplingNeumann(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx))
                     {
-                        DimVector pressureCorrection(intersection.centerUnitOuterNormal());
-                        pressureCorrection *= volVars.pressure();
-                        this->residual_[scvIdx][momentumYIdx] += pressureCorrection[momentumYIdx]
-                                                                 * boundaryVars.face().area;
+                        this->residual_[scvIdx][momentumYIdx] -= volVars.pressure()
+                                                                 * boundaryVars.face().normal[momentumYIdx];
                         Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]);
                     }
 
@@ -168,15 +166,13 @@ public:
                     for (int compIdx = 0; compIdx < numComponents; compIdx++)
                     {
                         int eqIdx = conti0EqIdx + compIdx;
-                        if (eqIdx != massBalanceIdx)
+                        if ((eqIdx != massBalanceIdx) && bcTypes.isCouplingDirichlet(eqIdx))
                         {
-                            if (bcTypes.isCouplingDirichlet(eqIdx))
-                            {
-                                if(useMoles)
-                                    this->residual_[scvIdx][eqIdx] = volVars.moleFraction(compIdx);
-                                else
-                                    this->residual_[scvIdx][eqIdx] = volVars.massFraction(compIdx);
-                            }
+                            if (useMoles)
+                                this->residual_[scvIdx][eqIdx] = volVars.moleFraction(compIdx);
+                            else
+                                this->residual_[scvIdx][eqIdx] = volVars.massFraction(compIdx);
+                            Valgrind::CheckDefined(this->residual_[scvIdx][compIdx]);
                         }
                     }
                 }

test/references/2cnistokes2p2cni-pm-reference.vtu

@@ -115,15 +115,15 @@
           -1.24484e-09 2.63388e-09 0 -1.40648e-09 2.47268e-09 0 -1.56187e-09 3.48588e-09 0
         </DataArray>
         <DataArray type="Float32" Name="velocityN" NumberOfComponents="3" format="ascii">
-          -9.11004e-16 9.11342e-10 0 -2.33277e-16 9.11378e-10 0 -1.30912e-14 2.85739e-09 0 -2.11362e-14 2.85745e-09 0
-          -7.93741e-16 9.11309e-10 0 -2.76293e-14 2.85753e-09 0 -1.73148e-15 9.11397e-10 0 -5.17278e-14 2.85765e-09 0
-          -1.84085e-15 9.11372e-10 0 -6.65425e-14 2.85791e-09 0 -7.77755e-16 9.1143e-10 0 -1.95408e-14 2.85812e-09 0
-          6.95172e-16 9.11414e-10 0 1.66251e-14 2.85806e-09 0 1.92997e-13 4.12376e-09 0 -4.26114e-14 4.1244e-09 0
-          -6.33254e-13 4.12507e-09 0 -1.23389e-12 4.12603e-09 0 -1.33455e-12 4.12705e-09 0 -6.84132e-13 4.1281e-09 0
-          -1.7865e-13 4.12905e-09 0 4.30742e-11 -1.28641e-09 0 3.31517e-11 -1.23795e-09 0 2.13677e-11 -1.2436e-09 0
-          2.02631e-11 -1.24388e-09 0 2.09812e-11 -1.23992e-09 0 3.13788e-11 -1.24544e-09 0 4.18153e-11 -1.19605e-09 0
-          3.73386e-09 -6.0169e-09 0 2.83077e-09 -5.92117e-09 0 1.91878e-09 -5.93359e-09 0 1.90508e-09 -5.93592e-09 0
-          1.9055e-09 -5.9295e-09 0 2.87766e-09 -5.94226e-09 0 3.8446e-09 -5.84549e-09 0
+          -4.27666e-16 9.11398e-10 0 5.35779e-17 9.114e-10 0 -1.65377e-14 2.85747e-09 0 -7.57308e-15 2.85753e-09 0
+          -5.31173e-16 9.11372e-10 0 -2.48814e-14 2.85751e-09 0 -1.10496e-15 9.11396e-10 0 -4.25772e-14 2.8577e-09 0
+          -1.02804e-15 9.11387e-10 0 -3.64365e-14 2.85781e-09 0 -5.53955e-16 9.1142e-10 0 -2.79711e-14 2.85795e-09 0
+          3.35419e-16 9.11376e-10 0 -1.70691e-14 2.85799e-09 0 -1.30459e-13 4.12427e-09 0 -2.49504e-13 4.12435e-09 0
+          -6.36421e-13 4.12502e-09 0 -1.14918e-12 4.1259e-09 0 -1.31946e-12 4.12705e-09 0 -9.10103e-13 4.12807e-09 0
+          -5.7535e-13 4.12836e-09 0 -1.08903e-12 -1.24603e-09 0 9.46607e-12 -1.25617e-09 0 2.07253e-11 -1.24637e-09 0
+          2.0015e-11 -1.24179e-09 0 1.881e-11 -1.24004e-09 0 9.07878e-12 -1.22727e-09 0 -8.61794e-13 -1.23507e-09 0
+          -1.39598e-10 -5.93706e-09 0 8.77525e-10 -5.95738e-09 0 1.89398e-09 -5.93912e-09 0 1.89665e-09 -5.93139e-09 0
+          1.89496e-09 -5.92995e-09 0 9.23398e-10 -5.90619e-09 0 -4.31519e-11 -5.92225e-09 0
         </DataArray>
       </PointData>
       <CellData Scalars="process rank">

test/references/2cnizeroeq2p2cni-ff-reference.vtu

@@ -92,24 +92,24 @@
           298.151 298.151 298.15 296.098 298.15 298.15 298.15 298.15 298.15 298.15 298.15 291.27
         </DataArray>
         <DataArray type="Float32" Name="v" NumberOfComponents="3" format="ascii">
-          0 0 0 0 0 0 3.5 0 0 3.71186 0.134918 0
-          -0.00881628 -0.00477543 0 3.15764 -0.102484 0 0.457328 0.00780399 0 3.54611 0.0340272 0
-          -0.0463155 -0.00621166 0 3.00943 -0.0107149 0 0.481567 0.00387028 0 3.46593 0.00676776 0
-          0 0 0 2.94966 -0.0348827 0 1.42325 -0.0157782 0 3.39995 -0.0513516 0
-          3.5 0 0 4.14566 0.118082 0 3.65518 0.0131766 0 4.05161 0.0621933 0
-          3.64544 0.0109064 0 4.03657 0.0275831 0 3.61726 -0.0207583 0 3.98511 -0.0304693 0
-          3.5 0 0 4.0414 0.0482092 0 3.61292 -0.0286125 0 4.10066 0.029847 0
-          3.68652 -0.0013876 0 4.12904 0.0188123 0 3.71514 -0.00667852 0 4.10585 -0.0112552 0
-          3.5 0 0 4.01777 -0.0052592 0 3.5756 -0.00309585 0 4.04829 -0.00490347 0
-          3.643 -0.0037801 0 4.08337 -3.19324e-05 0 3.67489 0.00479175 0 4.07276 0.00520353 0
-          3.5 0 0 4.04234 -0.0643566 0 3.62507 -0.0523174 0 4.11135 -0.0550917 0
-          3.70081 -0.0292999 0 4.14314 -0.0209893 0 3.72572 0.0145126 0 4.10934 0.0197014 0
-          3.5 0 0 4.07344 -0.133883 0 3.63522 -0.0369333 0 4.06249 -0.0694053 0
-          3.61398 -0.0197083 0 4.02474 -0.0248173 0 3.57513 0.0413183 0 3.93698 0.0524367 0
-          3.5 0 0 3.66584 -0.128329 0 3.11112 -0.0329058 0 3.48188 -0.0414777 0
-          2.91196 -0.00998838 0 3.39596 -0.00818624 0 2.86553 0.0615884 0 3.33484 0.0763954 0
+          0 0 0 0 0 0 3.5 0 0 3.71248 0.134867 0
+          -0.00814685 -0.00531453 0 3.15753 -0.102867 0 0.457021 0.00839452 0 3.54587 0.0345368 0
+          -0.0457672 -0.00642789 0 3.00958 -0.0108938 0 0.481019 0.00397879 0 3.4656 0.0069582 0
+          0 0 0 2.94988 -0.0348911 0 1.4228 -0.0158479 0 3.39987 -0.0514118 0
+          3.5 0 0 4.14603 0.118029 0 3.65505 0.0130453 0 4.05169 0.062583 0
+          3.64566 0.0107812 0 4.03656 0.0276851 0 3.61758 -0.0208556 0 3.98518 -0.0305992 0
+          3.5 0 0 4.04154 0.0481988 0 3.61285 -0.0287144 0 4.10077 0.0300385 0
+          3.6867 -0.0014391 0 4.12913 0.0189205 0 3.71539 -0.00673646 0 4.10593 -0.0113428 0
+          3.5 0 0 4.01787 -0.00525555 0 3.57559 -0.00311191 0 4.04834 -0.00482233 0
+          3.64307 -0.00381187 0 4.08341 -1.92334e-05 0 3.67501 0.00471022 0 4.07278 0.0051097 0
+          3.5 0 0 4.04239 -0.0643588 0 3.62507 -0.0523393 0 4.1114 -0.0550598 0
+          3.70086 -0.0293024 0 4.14318 -0.0209562 0 3.72576 0.0144946 0 4.10931 0.0196749 0
+          3.5 0 0 4.07347 -0.133885 0 3.63523 -0.0369328 0 4.06252 -0.0693878 0
+          3.61401 -0.0197134 0 4.02476 -0.0248193 0 3.57515 0.0412899 0 3.93694 0.0524073 0
+          3.5 0 0 3.66586 -0.128335 0 3.11111 -0.0329236 0 3.48191 -0.0414694 0
+          2.91199 -0.00997543 0 3.39598 -0.00815337 0 2.86555 0.061605 0 3.33479 0.0764045 0
           0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0 0 0 0 0 0 1.44882 0.00934465 0
+          0 0 0 0 0 0 0 0 0 1.44878 0.00934602 0
         </DataArray>
       </PointData>
       <CellData Scalars="mu_t">

test/references/2cnizeroeq2p2cni-pm-reference.vtu

@@ -104,24 +104,24 @@
           3 2 3 3 3 3
         </DataArray>
         <DataArray type="Float32" Name="velocityW" NumberOfComponents="3" format="ascii">
-          -2.00676e-08 -9.16679e-07 0 -3.55034e-08 -9.25937e-07 0 -7.2694e-09 -7.79487e-07 0 -1.85725e-08 -7.85804e-07 0
-          -5.72677e-08 -9.43895e-07 0 -3.2149e-08 -8.17695e-07 0 -5.81485e-08 -9.67281e-07 0 -2.82673e-08 -8.55201e-07 0
-          -3.6477e-08 -9.8846e-07 0 -1.5434e-08 -8.81287e-07 0 -2.0253e-08 -9.96464e-07 0 -8.75578e-09 -8.93124e-07 0
-          1.27988e-08 -3.91774e-07 0 -9.88897e-09 -3.67026e-07 0 -2.83273e-08 -4.37197e-07 0 -1.5594e-08 -5.0318e-07 0
-          -5.40254e-09 -5.31727e-07 0 -3.69487e-09 -5.49799e-07 0 8.71129e-08 -4.42703e-08 0 -6.62842e-09 1.27639e-07 0
-          -6.71148e-08 -3.88051e-08 0 -1.88233e-08 -1.38458e-07 0 -6.86887e-09 -1.55304e-07 0 -9.95107e-09 -1.84721e-07 0
-          1.92517e-07 5.27122e-08 0 7.64741e-08 3.43661e-07 0 -2.60563e-08 1.0529e-07 0 -6.81601e-09 -1.36785e-08 0
-          -7.77748e-09 -2.12678e-08 0 -1.44668e-08 -5.96281e-08 0
+          -2.06199e-08 -9.15899e-07 0 -3.64537e-08 -9.25417e-07 0 -7.208e-09 -7.78647e-07 0 -1.8506e-08 -7.85159e-07 0
+          -5.86978e-08 -9.43863e-07 0 -3.19732e-08 -8.17667e-07 0 -5.94499e-08 -9.67622e-07 0 -2.80051e-08 -8.55453e-07 0
+          -3.72174e-08 -9.88893e-07 0 -1.52651e-08 -8.81594e-07 0 -2.06431e-08 -9.96908e-07 0 -8.66232e-09 -8.93439e-07 0
+          1.29061e-08 -3.91165e-07 0 -9.79295e-09 -3.66427e-07 0 -2.8226e-08 -4.37129e-07 0 -1.55201e-08 -5.03201e-07 0
+          -5.3777e-09 -5.31872e-07 0 -3.67514e-09 -5.49948e-07 0 8.69343e-08 -4.42344e-08 0 -6.55548e-09 1.26961e-07 0
+          -6.705e-08 -3.83131e-08 0 -1.89759e-08 -1.38311e-07 0 -6.90394e-09 -1.55417e-07 0 -9.91087e-09 -1.84797e-07 0
+          1.91444e-07 5.24658e-08 0 7.62659e-08 3.41876e-07 0 -2.57961e-08 1.0616e-07 0 -6.92971e-09 -1.35043e-08 0
+          -7.78824e-09 -2.1386e-08 0 -1.43972e-08 -5.96684e-08 0
         </DataArray>
         <DataArray type="Float32" Name="velocityN" NumberOfComponents="3" format="ascii">
-          6.73309e-08 8.18601e-07 0 1.18843e-07 8.51098e-07 0 5.50954e-05 -8.92828e-05 0 9.02688e-05 -5.76987e-05 0
-          1.85973e-07 9.247e-07 0 0.000125914 -2.0172e-06 0 1.79421e-07 1.00211e-06 0 0.000113148 2.40232e-05 0
-          1.0809e-07 1.05695e-06 0 7.00806e-05 5.13838e-05 0 5.89289e-08 1.07687e-06 0 4.02512e-05 6.75515e-05 0
-          0.000152968 -0.000350157 0 0.000278122 -0.000245676 0 0.000265904 0.000129743 0 0.000197686 -1.48673e-05 0
-          0.000185526 0.000187731 0 0.000104212 0.000251117 0 1.68174e-05 -0.000562301 0 0.000808074 -0.000909317 0
-          0.000386486 0.00110686 0 0.000109459 -0.000708157 0 0.000527134 0.000631178 0 8.99139e-06 0.000394462 0
-          -0.00095173 -0.000603674 0 0.00158393 -0.00144353 0 0.000382565 0.00194928 0 -0.000225966 -0.00133953 0
-          0.00110891 0.000988605 0 -0.000684695 0.000420717 0
+          7.46646e-08 8.03738e-07 0 1.31524e-07 8.40301e-07 0 6.15559e-05 -9.98281e-05 0 0.000100778 -6.43549e-05 0
+          2.05061e-07 9.22636e-07 0 0.000139738 -1.04004e-06 0 1.96862e-07 1.00763e-06 0 0.000123602 2.83196e-05 0
+          1.18093e-07 1.06648e-06 0 7.52704e-05 5.58938e-05 0 6.42013e-08 1.08745e-06 0 4.28128e-05 7.19703e-05 0
+          0.000170439 -0.000390804 0 0.000311704 -0.000275673 0 0.000301279 0.000146191 0 0.000212835 5.43286e-07 0
+          0.000191522 0.000200466 0 0.000106961 0.000263158 0 1.67682e-05 -0.000625518 0 0.00089621 -0.0010165 0
+          0.000462243 0.00120823 0 0.000105831 -0.000699224 0 0.000534856 0.000654171 0 6.88585e-06 0.000408647 0
+          -0.00105676 -0.000669889 0 0.0017379 -0.0016112 0 0.000505535 0.00212108 0 -0.000278222 -0.0013439 0
+          0.00112951 0.00101813 0 -0.000706029 0.000433833 0
         </DataArray>
       </PointData>
       <CellData Scalars="process rank">

test/references/2cstokes2p2c-pm-reference.vtu

@@ -284,80 +284,80 @@
           3 3 3 3 3 3 3 3 3 3 3
         </DataArray>
         <DataArray type="Float32" Name="velocityW" NumberOfComponents="3" format="ascii">
-          -9.0203e-11 1.53631e-10 0 -1.75854e-10 1.45103e-10 0 -1.05722e-10 2.99748e-10 0 -2.06021e-10 2.81944e-10 0
-          -3.33987e-10 1.20423e-10 0 -3.91029e-10 2.30598e-10 0 -4.58662e-10 8.21357e-11 0 -5.36496e-10 1.51537e-10 0
-          -5.72792e-10 3.91221e-11 0 -6.44224e-10 5.319e-11 0 -5.67877e-10 1.06768e-11 0 -6.6833e-10 -5.77645e-11 0
-          -5.00886e-10 -6.86429e-11 0 -6.08142e-10 -1.7691e-10 0 -4.5174e-10 -1.16555e-10 0 -5.1904e-10 -2.583e-10 0
-          -3.30432e-10 -1.5171e-10 0 -3.83446e-10 -3.29822e-10 0 -1.73491e-10 -1.74642e-10 0 -2.01395e-10 -3.77277e-10 0
-          -8.88535e-11 -1.82557e-10 0 -1.0316e-10 -3.93617e-10 0 -1.48266e-10 6.14341e-10 0 -2.8851e-10 5.75879e-10 0
-          -5.4612e-10 4.65403e-10 0 -7.46307e-10 2.965e-10 0 -9.04521e-10 8.75193e-11 0 -9.38612e-10 -1.42219e-10 0
-          -8.27401e-10 -3.54203e-10 0 -6.99868e-10 -5.56911e-10 0 -5.26916e-10 -7.14371e-10 0 -2.76637e-10 -8.14266e-10 0
-          -1.41688e-10 -8.48619e-10 0 -2.14244e-10 9.99198e-10 0 -4.15901e-10 9.34476e-10 0 -7.8402e-10 7.49904e-10 0
-          -1.06541e-09 4.70789e-10 0 -1.25862e-09 1.34404e-10 0 -1.28962e-09 -2.07901e-10 0 -1.18826e-09 -5.62417e-10 0
-          -1.01894e-09 -8.98368e-10 0 -7.41166e-10 -1.15115e-09 0 -3.88921e-10 -1.3112e-09 0 -1.99211e-10 -1.36641e-09 0
-          -3.02511e-10 1.76028e-09 0 -5.84953e-10 1.66091e-09 0 -1.09568e-09 1.3813e-09 0 -1.47751e-09 9.67207e-10 0
-          -1.70716e-09 4.75901e-10 0 -1.75486e-09 -3.19085e-11 0 -1.6405e-09 -5.48727e-10 0 -1.39213e-09 -9.98713e-10 0
-          -1.0094e-09 -1.3538e-09 0 -5.30561e-10 -1.58382e-09 0 -2.72049e-10 -1.66366e-09 0 -4.13305e-10 2.91618e-09 0
-          -7.94082e-10 2.77035e-09 0 -1.47267e-09 2.36916e-09 0 -1.96234e-09 1.79396e-09 0 -2.23608e-09 1.12908e-09 0
-          -2.29163e-09 4.42374e-10 0 -2.14418e-09 -2.25128e-10 0 -1.8135e-09 -8.16932e-10 0 -1.31875e-09 -1.29099e-09 0
-          -6.95496e-10 -1.60228e-09 0 -3.57427e-10 -1.71142e-09 0 -5.48342e-10 4.24548e-09 0 -1.04268e-09 4.03608e-09 0
-          -1.9057e-09 3.48061e-09 0 -2.49905e-09 2.7192e-09 0 -2.81147e-09 1.86925e-09 0 -2.86236e-09 1.00697e-09 0
-          -2.67434e-09 1.83622e-10 0 -2.26632e-09 -5.55049e-10 0 -1.65664e-09 -1.15753e-09 0 -8.78535e-10 -1.56213e-09 0
-          -4.53242e-10 -1.70636e-09 0 -7.10581e-10 6.53185e-09 0 -1.32886e-09 6.23303e-09 0 -2.37951e-09 5.48524e-09 0
-          -3.05923e-09 4.51877e-09 0 -3.39418e-09 3.48272e-09 0 -3.43027e-09 2.45453e-09 0 -3.2014e-09 1.47993e-09 0
-          -2.7246e-09 5.98273e-10 0 -2.00751e-09 -1.38085e-10 0 -1.07352e-09 -6.48925e-10 0 -5.57201e-10 -8.35675e-10 0
-          -8.81613e-10 1.13828e-08 0 -1.60572e-09 1.09569e-08 0 -2.79651e-09 9.98018e-09 0 -3.51312e-09 8.80387e-09 0
-          -3.84175e-09 7.59393e-09 0 -3.85546e-09 6.41677e-09 0 -3.59544e-09 5.30651e-09 0 -3.07594e-09 4.29348e-09 0
-          -2.29031e-09 3.42385e-09 0 -1.23897e-09 2.795e-09 0 -6.488e-10 2.55698e-09 0 -1.01654e-09 1.92855e-08 0
-          -1.77818e-09 1.86843e-08 0 -2.98561e-09 1.74593e-08 0 -3.65553e-09 1.60908e-08 0 -3.94067e-09 1.47336e-08 0
-          -3.92905e-09 1.3434e-08 0 -3.66286e-09 1.2213e-08 0 -3.15274e-09 1.10887e-08 0 -2.37783e-09 1.00971e-08 0
-          -1.30595e-09 9.34497e-09 0 -6.92396e-10 9.0478e-09 0 -1.09977e-09 2.94864e-08 0 -1.81541e-09 2.86582e-08 0
-          -2.90824e-09 2.71954e-08 0 -3.46463e-09 2.56724e-08 0 -3.68364e-09 2.42034e-08 0 -3.65061e-09 2.28128e-08 0
-          -3.40318e-09 2.15095e-08 0 -2.94881e-09 2.03009e-08 0 -2.25787e-09 1.9209e-08 0 -1.26407e-09 1.83383e-08 0
-          -6.81455e-10 1.79773e-08 0 -1.24784e-09 4.11143e-08 0 -1.90961e-09 3.99989e-08 0 -2.88739e-09 3.8328e-08 0
-          -3.34613e-09 3.66908e-08 0 -3.51145e-09 3.5142e-08 0 -3.46041e-09 3.36871e-08 0 -3.2263e-09 3.23256e-08 0
-          -2.81492e-09 3.10568e-08 0 -2.19246e-09 2.98877e-08 0 -1.25639e-09 2.89069e-08 0 -6.91875e-10 2.84781e-08 0
-          -1.31967e-09 4.71612e-08 0 -1.89569e-09 4.5885e-08 0 -2.7324e-09 4.41182e-08 0 -3.10653e-09 4.24329e-08 0
-          -3.23116e-09 4.08507e-08 0 -3.17209e-09 3.93691e-08 0 -2.95795e-09 3.79835e-08 0 -2.59358e-09 3.66895e-08 0
-          -2.04706e-09 3.54869e-08 0 -1.19578e-09 3.44528e-08 0 -6.70507e-10 3.39888e-08 0
+          -9.02289e-11 1.535e-10 0 -1.76444e-10 1.45133e-10 0 -1.05676e-10 2.99626e-10 0 -2.06015e-10 2.81891e-10 0
+          -3.34647e-10 1.21161e-10 0 -3.91221e-10 2.30887e-10 0 -4.58701e-10 8.27743e-11 0 -5.38492e-10 1.50887e-10 0
+          -5.634e-10 3.86321e-11 0 -6.33856e-10 4.99678e-11 0 -5.63444e-10 8.54657e-12 0 -6.56601e-10 -4.43328e-11 0
+          -5.45136e-10 -6.23664e-11 0 -6.33702e-10 -1.56494e-10 0 -4.55145e-10 -9.25473e-11 0 -5.28766e-10 -2.37906e-10 0
+          -2.94785e-10 -1.48323e-10 0 -3.68133e-10 -3.2748e-10 0 -1.73855e-10 -1.75029e-10 0 -2.01475e-10 -3.77744e-10 0
+          -8.90109e-11 -1.82851e-10 0 -1.03198e-10 -3.93982e-10 0 -1.48255e-10 6.14267e-10 0 -2.88411e-10 5.75781e-10 0
+          -5.46293e-10 4.65292e-10 0 -7.57132e-10 2.94935e-10 0 -8.7572e-10 7.97767e-11 0 -9.00192e-10 -1.39138e-10 0
+          -8.58174e-10 -3.59407e-10 0 -7.28078e-10 -5.52489e-10 0 -5.25098e-10 -7.13821e-10 0 -2.76393e-10 -8.14608e-10 0
+          -1.41588e-10 -8.48915e-10 0 -2.14265e-10 9.9934e-10 0 -4.15828e-10 9.34631e-10 0 -7.8413e-10 7.50089e-10 0
+          -1.07988e-09 4.72729e-10 0 -1.23793e-09 1.40606e-10 0 -1.27138e-09 -2.27903e-10 0 -1.20986e-09 -5.86018e-10 0
+          -1.02339e-09 -9.01764e-10 0 -7.40569e-10 -1.15104e-09 0 -3.88702e-10 -1.31119e-09 0 -1.99109e-10 -1.36644e-09 0
+          -3.02463e-10 1.76056e-09 0 -5.84861e-10 1.66122e-09 0 -1.09549e-09 1.38169e-09 0 -1.47905e-09 9.70022e-10 0
+          -1.70071e-09 4.86671e-10 0 -1.75353e-09 -4.33412e-11 0 -1.64784e-09 -5.49549e-10 0 -1.39259e-09 -9.97512e-10 0
+          -1.00922e-09 -1.35331e-09 0 -5.30472e-10 -1.58361e-09 0 -2.72007e-10 -1.66352e-09 0 -4.13245e-10 2.91637e-09 0
+          -7.93956e-10 2.77055e-09 0 -1.47233e-09 2.3694e-09 0 -1.96153e-09 1.7942e-09 0 -2.23467e-09 1.12844e-09 0
+          -2.29207e-09 4.38789e-10 0 -2.14623e-09 -2.25155e-10 0 -1.81448e-09 -8.16156e-10 0 -1.319e-09 -1.29053e-09 0
+          -6.95552e-10 -1.60204e-09 0 -3.57446e-10 -1.71125e-09 0 -5.48285e-10 4.24559e-09 0 -1.04255e-09 4.03617e-09 0
+          -1.90539e-09 3.48067e-09 0 -2.49854e-09 2.71905e-09 0 -2.81106e-09 1.86852e-09 0 -2.86265e-09 1.00585e-09 0
+          -2.67522e-09 1.83092e-10 0 -2.26706e-09 -5.54904e-10 0 -1.65702e-09 -1.15728e-09 0 -8.78673e-10 -1.56197e-09 0
+          -4.53302e-10 -1.70624e-09 0 -7.10535e-10 6.53188e-09 0 -1.32875e-09 6.23304e-09 0 -2.37928e-09 5.4852e-09 0
+          -3.05893e-09 4.51855e-09 0 -3.39401e-09 3.48222e-09 0 -3.4305e-09 2.4539e-09 0 -3.20195e-09 1.47957e-09 0
+          -2.72515e-09 5.98222e-10 0 -2.00788e-09 -1.38008e-10 0 -1.07369e-09 -6.48853e-10 0 -5.57282e-10 -8.35621e-10 0
+          -8.81574e-10 1.13828e-08 0 -1.60563e-09 1.09568e-08 0 -2.79633e-09 9.98008e-09 0 -3.51292e-09 8.80366e-09 0
+          -3.84168e-09 7.59358e-09 0 -3.85566e-09 6.41638e-09 0 -3.59585e-09 5.30624e-09 0 -3.07638e-09 4.29338e-09 0
+          -2.29063e-09 3.42384e-09 0 -1.23914e-09 2.79501e-09 0 -6.48886e-10 2.55699e-09 0 -1.01651e-09 1.92855e-08 0
+          -1.77811e-09 1.86843e-08 0 -2.98548e-09 1.74593e-08 0 -3.65539e-09 1.60906e-08 0 -3.94064e-09 1.47333e-08 0
+          -3.9292e-09 1.34337e-08 0 -3.66317e-09 1.22128e-08 0 -3.15309e-09 1.10886e-08 0 -2.37811e-09 1.0097e-08 0
+          -1.30609e-09 9.34495e-09 0 -6.92472e-10 9.04777e-09 0 -1.09977e-09 2.94864e-08 0 -1.81537e-09 2.86582e-08 0
+          -2.90814e-09 2.71953e-08 0 -3.46454e-09 2.56723e-08 0 -3.68363e-09 2.42033e-08 0 -3.65073e-09 2.28126e-08 0
+          -3.40341e-09 2.15093e-08 0 -2.94909e-09 2.03008e-08 0 -2.25809e-09 1.9209e-08 0 -1.26419e-09 1.83383e-08 0
+          -6.81524e-10 1.79772e-08 0 -1.24783e-09 4.11143e-08 0 -1.90958e-09 3.99989e-08 0 -2.88731e-09 3.83279e-08 0
+          -3.34606e-09 3.66907e-08 0 -3.51145e-09 3.51419e-08 0 -3.46051e-09 3.3687e-08 0 -3.2265e-09 3.23256e-08 0
+          -2.81514e-09 3.10568e-08 0 -2.19264e-09 2.98877e-08 0 -1.25649e-09 2.89069e-08 0 -6.91938e-10 2.8478e-08 0
+          -1.31965e-09 4.71612e-08 0 -1.89565e-09 4.5885e-08 0 -2.73234e-09 4.41182e-08 0 -3.10646e-09 4.24329e-08 0
+          -3.23117e-09 4.08507e-08 0 -3.17218e-09 3.9369e-08 0 -2.95811e-09 3.79834e-08 0 -2.59376e-09 3.66895e-08 0
+          -2.04721e-09 3.54869e-08 0 -1.19586e-09 3.44528e-08 0 -6.70536e-10 3.39889e-08 0
         </DataArray>
         <DataArray type="Float32" Name="velocityN" NumberOfComponents="3" format="ascii">
-          8.09459e-17 3.42094e-11 0 9.22242e-17 3.41892e-11 0 6.39537e-18 8.76404e-11 0 -5.55732e-17 8.76337e-11 0
-          -3.27316e-16 3.416e-11 0 3.51537e-16 8.76329e-11 0 -2.34271e-16 3.43583e-11 0 -4.77571e-17 8.76595e-11 0
-          2.96387e-17 3.42827e-11 0 -4.80492e-16 8.76421e-11 0 -4.1242e-16 3.43325e-11 0 1.31783e-15 8.76724e-11 0
-          2.27562e-15 3.4604e-11 0 2.41035e-15 8.73849e-11 0 4.87143e-16 3.33305e-11 0 -9.04802e-16 8.70396e-11 0
-          -2.08885e-15 3.4256e-11 0 -1.92964e-15 8.76188e-11 0 -1.73708e-18 3.42585e-11 0 5.30731e-17 8.7616e-11 0
-          2.74304e-18 3.42579e-11 0 1.56311e-17 8.76149e-11 0 1.38463e-16 2.00051e-10 0 3.67234e-16 2.0005e-10 0
-          1.52674e-15 2.00043e-10 0 -3.09887e-15 1.99884e-10 0 -2.50972e-15 2.00251e-10 0 -1.60013e-14 2.00146e-10 0
-          -2.60095e-15 2.01444e-10 0 1.9909e-14 2.00163e-10 0 5.00245e-15 1.99918e-10 0 3.44731e-16 1.99898e-10 0
-          6.60268e-17 1.99895e-10 0 4.3126e-15 3.11059e-10 0 9.18451e-15 3.11046e-10 0 7.72043e-15 3.11001e-10 0
-          -2.73583e-15 3.10919e-10 0 2.9307e-14 3.11245e-10 0 1.62397e-13 3.10995e-10 0 -7.12411e-14 3.12082e-10 0
-          -1.52439e-13 3.11406e-10 0 5.68445e-14 3.10921e-10 0 9.40955e-15 3.10869e-10 0 2.0587e-15 3.10863e-10 0
-          2.16884e-13 1.17297e-10 0 3.55368e-13 1.17357e-10 0 4.38133e-13 1.17528e-10 0 5.44621e-13 1.17279e-10 0
-          1.42603e-12 1.16685e-10 0 3.37908e-12 1.16244e-10 0 -6.35881e-13 1.21656e-10 0 -2.86628e-12 1.1423e-10 0
-          4.23097e-13 1.17228e-10 0 4.43403e-13 1.1761e-10 0 1.92734e-13 1.1753e-10 0 7.51697e-13 -3.99639e-10 0
-          1.11589e-12 -3.99183e-10 0 1.67489e-12 -3.9866e-10 0 1.93377e-12 -3.99052e-10 0 1.72196e-12 -3.99915e-10 0
-          -1.38038e-13 -3.99322e-10 0 2.36498e-12 -3.92297e-10 0 3.32024e-12 -4.02513e-10 0 6.59382e-13 -3.99493e-10 0
-          9.58417e-13 -3.98072e-10 0 7.8646e-13 -3.97719e-10 0 2.76127e-12 -1.01393e-09 0 3.65395e-12 -1.01188e-09 0
-          4.53534e-12 -1.0106e-09 0 4.4785e-12 -1.01084e-09 0 4.07748e-12 -1.0112e-09 0 3.71857e-12 -1.00966e-09 0
-          4.54908e-12 -1.00832e-09 0 4.95792e-12 -1.01118e-09 0 4.3363e-12 -1.01155e-09 0 3.41942e-12 -1.00919e-09 0
-          2.69906e-12 -1.00698e-09 0 1.04902e-11 -2.52423e-09 0 1.21049e-11 -2.51582e-09 0 1.2977e-11 -2.51253e-09 0
-          1.20636e-11 -2.51305e-09 0 1.1765e-11 -2.51317e-09 0 1.16806e-11 -2.51197e-09 0 1.19772e-11 -2.51134e-09 0
-          1.22688e-11 -2.51312e-09 0 1.29316e-11 -2.51396e-09 0 1.19511e-11 -2.50995e-09 0 1.0346e-11 -2.50138e-09 0
-          9.02027e-11 -6.59633e-09 0 9.66638e-11 -6.54282e-09 0 9.60278e-11 -6.53157e-09 0 8.80769e-11 -6.53537e-09 0
-          8.71363e-11 -6.53584e-09 0 8.71253e-11 -6.53479e-09 0 8.74113e-11 -6.53433e-09 0 8.84258e-11 -6.53622e-09 0
-          9.60723e-11 -6.54025e-09 0 9.60604e-11 -6.5283e-09 0 8.92064e-11 -6.47537e-09 0 4.46088e-10 -1.38946e-08 0
-          4.39979e-10 -1.36212e-08 0 4.02099e-10 -1.35927e-08 0 3.6872e-10 -1.36115e-08 0 3.67102e-10 -1.3613e-08 0
-          3.67227e-10 -1.36121e-08 0 3.67539e-10 -1.36117e-08 0 3.69591e-10 -1.36145e-08 0 4.02625e-10 -1.36331e-08 0
-          4.37365e-10 -1.36031e-08 0 4.40948e-10 -1.3334e-08 0 1.54581e-09 -2.40176e-08 0 1.38018e-09 -2.29973e-08 0
-          1.13132e-09 -2.29842e-08 0 1.0491e-09 -2.30432e-08 0 1.04991e-09 -2.30449e-08 0 1.05001e-09 -2.30442e-08 0
-          1.05066e-09 -2.30438e-08 0 1.05121e-09 -2.30473e-08 0 1.13382e-09 -2.31051e-08 0 1.37136e-09 -2.30864e-08 0
-          1.52647e-09 -2.20834e-08 0 4.15941e-09 -3.65999e-08 0 3.32957e-09 -3.36353e-08 0 2.38314e-09 -3.38422e-08 0
-          2.27491e-09 -3.39647e-08 0 2.28228e-09 -3.39635e-08 0 2.28193e-09 -3.39633e-08 0 2.28349e-09 -3.3963e-08 0
-          2.27899e-09 -3.39647e-08 0 2.39151e-09 -3.40848e-08 0 3.30623e-09 -3.42736e-08 0 4.103e-09 -3.13608e-08 0
-          8.74859e-09 -4.34525e-08 0 6.44885e-09 -3.91106e-08 0 4.10099e-09 -3.95116e-08 0 4.06196e-09 -3.96695e-08 0
-          4.06986e-09 -3.96656e-08 0 4.06998e-09 -3.96659e-08 0 4.07172e-09 -3.96656e-08 0 4.06836e-09 -3.96655e-08 0
-          4.1209e-09 -3.98205e-08 0 6.40043e-09 -4.01939e-08 0 8.6236e-09 -3.59286e-08 0
+          -4.63284e-15 3.4408e-11 0 4.52443e-15 3.52512e-11 0 -1.47562e-14 8.77743e-11 0 3.74846e-15 8.88436e-11 0
+          5.21751e-15 3.24295e-11 0 1.07875e-14 8.66491e-11 0 -4.87018e-15 3.32149e-11 0 -4.60903e-15 8.69645e-11 0
+          -2.60043e-15 3.46028e-11 0 -3.46804e-15 8.79664e-11 0 5.64789e-17 3.43376e-11 0 5.74751e-16 8.77256e-11 0
+          1.25612e-14 3.45962e-11 0 -1.08728e-14 8.78854e-11 0 9.33912e-16 2.81314e-11 0 1.19118e-15 8.68778e-11 0
+          -1.22254e-14 3.42488e-11 0 1.19266e-14 8.76214e-11 0 -8.12516e-18 3.42581e-11 0 9.60099e-17 8.76158e-11 0
+          2.5571e-20 3.42581e-11 0 1.24467e-17 8.7615e-11 0 -1.70582e-14 2.00086e-10 0 -5.31465e-15 2.01225e-10 0
+          5.91334e-15 2.00296e-10 0 8.18902e-15 2.00042e-10 0 -1.60035e-17 1.99954e-10 0 -8.01105e-16 2.00305e-10 0
+          -1.89812e-14 1.99947e-10 0 -1.09009e-15 2.03454e-10 0 2.30758e-14 1.99935e-10 0 5.15086e-16 1.99898e-10 0
+          7.08141e-17 1.99895e-10 0 -1.05539e-13 3.11036e-10 0 -9.80052e-14 3.11687e-10 0 -1.89769e-14 3.11734e-10 0
+          1.32651e-13 3.11467e-10 0 9.06583e-14 3.10695e-10 0 -7.74487e-15 3.11159e-10 0 -1.21627e-13 3.10791e-10 0
+          -2.18598e-14 3.12444e-10 0 1.19336e-13 3.10948e-10 0 1.2679e-14 3.10869e-10 0 2.19919e-15 3.10862e-10 0
+          -1.23261e-12 1.18859e-10 0 -1.30894e-12 1.179e-10 0 -6.39038e-14 1.16706e-10 0 2.33608e-12 1.16758e-10 0
+          2.07595e-12 1.18773e-10 0 2.22139e-13 1.17846e-10 0 -1.28189e-12 1.18458e-10 0 -5.62927e-14 1.15487e-10 0
+          1.57929e-12 1.18152e-10 0 5.89105e-13 1.17769e-10 0 1.77885e-13 1.17512e-10 0 1.08924e-12 -3.96783e-10 0
+          1.54152e-12 -3.98463e-10 0 1.63131e-12 -4.00543e-10 0 6.55823e-13 -4.00032e-10 0 1.03109e-12 -3.96664e-10 0
+          1.7768e-12 -3.97422e-10 0 2.5413e-12 -3.97376e-10 0 1.65035e-12 -4.01481e-10 0 3.85789e-13 -3.97823e-10 0
+          7.41431e-13 -3.97666e-10 0 4.63634e-13 -3.98018e-10 0 1.1211e-12 -1.01028e-09 0 2.21282e-12 -1.01137e-09 0
+          3.54135e-12 -1.01216e-09 0 3.71788e-12 -1.01122e-09 0 3.90851e-12 -1.00949e-09 0 4.38503e-12 -1.00912e-09 0
+          4.65351e-12 -1.01008e-09 0 4.14384e-12 -1.01118e-09 0 3.08035e-12 -1.00982e-09 0 1.68816e-12 -1.0086e-09 0
+          8.06945e-13 -1.00884e-09 0 1.91336e-12 -2.51456e-09 0 5.29353e-12 -2.51644e-09 0 1.00146e-11 -2.51544e-09 0
+          1.14952e-11 -2.51359e-09 0 1.1732e-11 -2.51204e-09 0 1.19351e-11 -2.51163e-09 0 1.20298e-11 -2.51237e-09 0
+          1.16644e-11 -2.51281e-09 0 9.822e-12 -2.51083e-09 0 5.00763e-12 -2.50829e-09 0 1.68149e-12 -2.5095e-09 0
+          1.06862e-11 -6.54374e-09 0 3.80262e-11 -6.55438e-09 0 7.53784e-11 -6.5439e-09 0 8.63032e-11 -6.53695e-09 0
+          8.73109e-11 -6.53496e-09 0 8.75192e-11 -6.53452e-09 0 8.76032e-11 -6.53509e-09 0 8.64526e-11 -6.53475e-09 0
+          7.50445e-11 -6.5277e-09 0 3.72139e-11 -6.51584e-09 0 9.66985e-12 -6.52649e-09 0 2.8963e-11 -1.3638e-08 0
+          1.58451e-10 -1.36964e-08 0 3.25856e-10 -1.36412e-08 0 3.65859e-10 -1.36155e-08 0 3.68018e-10 -1.36123e-08 0
+          3.68206e-10 -1.36119e-08 0 3.68261e-10 -1.36123e-08 0 3.65918e-10 -1.36105e-08 0 3.25114e-10 -1.35844e-08 0
+          1.55416e-10 -1.35275e-08 0 2.4239e-11 -1.35891e-08 0 1.38539e-11 -2.30871e-08 0 4.45727e-10 -2.33237e-08 0
+          9.63752e-10 -2.31175e-08 0 1.05103e-09 -2.30488e-08 0 1.05229e-09 -2.30444e-08 0 1.05254e-09 -2.30439e-08 0
+          1.05239e-09 -2.30441e-08 0 1.05084e-09 -2.3041e-08 0 9.62573e-10 -2.29716e-08 0 4.36156e-10 -2.27611e-08 0
+          -3.26238e-12 -2.30115e-08 0 -1.85805e-10 -3.39613e-08 0 9.42816e-10 -3.46984e-08 0 2.18211e-09 -3.40968e-08 0
+          2.28974e-09 -3.39674e-08 0 2.28696e-09 -3.39634e-08 0 2.28743e-09 -3.39629e-08 0 2.28684e-09 -3.39626e-08 0
+          2.28888e-09 -3.39602e-08 0 2.18177e-09 -3.38301e-08 0 9.18587e-10 -3.32149e-08 0 -2.34637e-10 -3.39942e-08 0
+          -7.05201e-10 -3.96138e-08 0 1.63338e-09 -4.07131e-08 0 4.02863e-09 -3.98306e-08 0 4.08262e-09 -3.9669e-08 0
+          4.0795e-09 -3.96659e-08 0 4.07943e-09 -3.96654e-08 0 4.07916e-09 -3.96648e-08 0 4.08039e-09 -3.96633e-08 0
+          4.03166e-09 -3.95013e-08 0 1.58445e-09 -3.85982e-08 0 -8.12153e-10 -3.97601e-08 0
         </DataArray>
       </PointData>
       <CellData Scalars="process rank">

test/references/2czeroeq2p2c-ff-reference.vtu

@@ -130,39 +130,39 @@
           2.49337e-05 2.49399e-05 2.49341e-05 2.49395e-05 2.49347e-05 2.49389e-05 2.49354e-05 2.49383e-05 2.49361e-05 2.49377e-05 2.49368e-05
         </DataArray>
         <DataArray type="Float32" Name="v" NumberOfComponents="3" format="ascii">
-          0 0 0 0 0 0 3.5 0 0 3.42031 0.105684 0
-          0 0 0 2.99391 0.0104657 0 0 0 0 3.2173 0.050051 0
-          0 0 0 2.66765 0.0314799 0 0 0 0 2.9734 0.037838 0
-          -0.076672 -0.00203125 0 2.45282 -0.0378483 0 0.285243 0.00416202 0 2.82821 0.0195913 0
-          -0.10628 -0.00210128 0 2.33457 -0.00672857 0 0.313956 0.00202284 0 2.70291 0.0267442 0
-          -0.369267 -0.00125236 0 2.10005 0.0949269 0 3.5 0 0 3.953 0.116866 0
-          3.6054 0.0191981 0 3.91712 0.072603 0 3.55298 0.0410076 0 3.87468 0.0489078 0
-          3.49323 0.0146443 0 3.7544 0.0320538 0 3.40374 0.0197924 0 3.64359 0.049098 0
-          3.30567 0.127928 0 3.5 0 0 3.98608 0.0896666 0 3.65981 0.026136 0
-          4.01148 0.0712661 0 3.71644 0.0450585 0 4.0422 0.054771 0 3.74111 -0.00533968 0
-          4.03546 0.0395597 0 3.73293 0.0203734 0 4.0129 0.0611976 0 3.76148 0.073422 0
-          3.5 0 0 3.9465 0.0587035 0 3.61639 0.0228949 0 3.98955 0.0555594 0
-          3.69291 0.0340563 0 4.05209 0.0419241 0 3.74766 0.0227078 0 4.05047 0.0346005 0
-          3.77765 0.0293089 0 4.07153 0.0528636 0 3.82463 0.0594054 0 3.5 0 0
-          3.92744 0.0397325 0 3.59735 0.0204328 0 3.96913 0.0436784 0 3.67278 0.027576 0
-          4.0148 0.0354885 0 3.72553 -0.0016181 0 4.04589 0.0283468 0 3.76361 0.0170233 0
-          4.07856 0.0453794 0 3.82011 0.0527704 0 3.5 0 0 3.90187 0.0154156 0
-          3.56677 0.0141365 0 3.94613 0.0243739 0 3.63931 0.0188785 0 3.99716 0.02062 0
-          3.69538 0.0160513 0 4.01396 0.0197701 0 3.73658 0.0198135 0 4.0479 0.032187 0
-          3.79187 0.0358393 0 3.5 0 0 3.91037 0.00511844 0 3.57824 0.00767284 0
-          3.95003 0.0113931 0 3.64746 0.0105166 0 3.98934 0.0106819 0 3.70005 -0.00611094 0
-          4.02612 0.00981558 0 3.74363 0.00544337 0 4.064 0.022734 0 3.798 0.0269709 0
-          3.5 0 0 3.90672 -0.0254145 0 3.57208 0.0021926 0 3.94765 -0.0105211 0
-          3.64114 0.00423736 0 3.98995 -0.00503479 0 3.69632 0.00372801 0 4.01818 -0.000381241 0
-          3.74129 0.00509242 0 4.05577 0.00814756 0 3.79212 0.00886696 0 3.5 0 0
-          3.95969 -0.0429627 0 3.63299 -0.00450346 0 3.99108 -0.0264094 0 3.69718 -0.00555344 0
-          4.0272 -0.0184408 0 3.7477 -0.0177384 0 4.06737 -0.0134557 0 3.79265 -0.0111668 0
-          4.10707 -0.00197109 0 3.84013 -0.000461393 0 3.5 0 0 3.97012 -0.112655 0
-          3.64251 0.000689451 0 3.98523 -0.0503487 0 3.68891 -0.010443 0 4.0113 -0.0330436 0
-          3.72419 -0.0142373 0 4.03603 -0.0237837 0 3.75535 -0.0139811 0 4.06403 -0.0164696 0
-          3.7883 -0.0153316 0 3.5 0 0 3.63052 -0.117071 0 3.3024 0.00317627 0
-          3.58434 -0.048673 0 3.21489 -0.0144778 0 3.54019 -0.0329373 0 3.15204 -0.0276904 0
-          3.49769 -0.0263932 0 3.10308 -0.0225403 0 3.45986 -0.0200677 0 3.06467 0.00474018 0
+          0 0 0 0 0 0 3.5 0 0 3.42039 0.105685 0
+          0 0 0 2.99391 0.0104661 0 0 0 0 3.21739 0.0500521 0
+          0 0 0 2.66764 0.0314783 0 0 0 0 2.97359 0.0378259 0
+          -0.0763506 -0.00225326 0 2.45272 -0.0380239 0 0.285026 0.00443846 0 2.82785 0.0198553 0
+          -0.105987 -0.00227871 0 2.33453 -0.00687253 0 0.313534 0.00213547 0 2.70225 0.0268356 0
+          -0.368786 -0.00128892 0 2.09959 0.0950299 0 3.5 0 0 3.95307 0.116867 0
+          3.6054 0.0191937 0 3.9172 0.0726101 0 3.55299 0.0409807 0 3.87482 0.0488958 0
+          3.49317 0.0145875 0 3.75431 0.0322699 0 3.40378 0.0197108 0 3.6432 0.0492593 0
+          3.30538 0.128155 0 3.5 0 0 3.98614 0.0896676 0 3.65981 0.0261315 0
+          4.01155 0.0712719 0 3.71644 0.0450354 0 4.04225 0.054773 0 3.74103 -0.00541899 0
+          4.03545 0.039754 0 3.73295 0.0203309 0 4.01269 0.0613645 0 3.76142 0.0736453 0
+          3.5 0 0 3.94656 0.0587049 0 3.61639 0.0228875 0 3.98962 0.0555677 0
+          3.6929 0.0340277 0 4.05215 0.0419361 0 3.7476 0.0226901 0 4.05046 0.0347533 0
+          3.77771 0.0293125 0 4.07143 0.053051 0 3.82469 0.0596366 0 3.5 0 0
+          3.92749 0.0397335 0 3.59735 0.0204256 0 3.96918 0.0436871 0 3.67278 0.0275528 0
+          4.01482 0.0355096 0 3.72549 -0.0016625 0 4.0459 0.0284766 0 3.76365 0.0170372 0
+          4.07852 0.045538 0 3.82022 0.0529653 0 3.5 0 0 3.90192 0.0154167 0
+          3.56677 0.0141287 0 3.94618 0.0243829 0 3.63931 0.0188579 0 3.9972 0.0206442 0
+          3.69536 0.0160524 0 4.01396 0.0198738 0 3.73664 0.0198548 0 4.04789 0.0323454 0
+          3.792 0.0360274 0 3.5 0 0 3.91041 0.0051193 0 3.57824 0.00766547 0
+          3.95006 0.0114022 0 3.64746 0.0104989 0 3.98935 0.0107068 0 3.70005 -0.00613225 0
+          4.02613 0.00990023 0 3.74368 0.00547805 0 4.06402 0.0228593 0 3.79813 0.0271196 0
+          3.5 0 0 3.90675 -0.025414 0 3.57208 0.00218634 0 3.94768 -0.0105137 0
+          3.64114 0.00422612 0 3.98997 -0.00501142 0 3.69633 0.00373976 0 4.0182 -0.000310911 0
+          3.74134 0.00514388 0 4.0558 0.00826693 0 3.79225 0.00900438 0 3.5 0 0
+          3.95971 -0.0429623 0 3.63299 -0.00450921 0 3.9911 -0.0264027 0 3.69718 -0.00556516 0
+          4.02721 -0.0184224 0 3.74772 -0.0177478 0 4.0674 -0.0134027 0 3.79271 -0.0111331 0
+          4.10713 -0.00188518 0 3.84026 -0.000361623 0 3.5 0 0 3.97014 -0.112655 0
+          3.64251 0.000685823 0 3.98525 -0.0503451 0 3.68892 -0.0104478 0 4.01132 -0.0330299 0
+          3.72422 -0.0142259 0 4.03607 -0.0237427 0 3.75544 -0.0139416 0 4.06413 -0.0163948 0
+          3.78848 -0.015248 0 3.5 0 0 3.63054 -0.117071 0 3.3024 0.00317311 0
+          3.58435 -0.0486707 0 3.21489 -0.014486 0 3.5402 -0.0329333 0 3.15206 -0.0277051 0
+          3.49773 -0.0263782 0 3.10317 -0.0225365 0 3.45998 -0.0200415 0 3.0649 0.00477195 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

test/references/2czeroeq2p2c-pm-reference.vtu

@@ -13,13 +13,13 @@
           0.949676 1 1 1 0.967518 0.963144
         </DataArray>
         <DataArray type="Float32" Name="Sw" NumberOfComponents="1" format="ascii">
-          0.98762 0.987622 0.984158 0.984162 0.987626 0.98417 0.98763 0.984177 0.987632 0.984182 0.987633 0.984183
-          0.961085 0.961141 0.961252 0.961345 0.961405 0.961425 0.851987 0.852172 0.852542 0.852855 0.853064 0.853136
-          0.639811 0.640103 0.640716 0.64123 0.641585 0.64172 0.417826 0.418093 0.418699 0.419185 0.419562 0.41973
-          0.261103 0.26125 0.261715 0.262079 0.262484 0.262682 0.166881 0.166797 0.167095 0.167393 0.168011 0.168358
-          0.111476 0.110876 0.110968 0.111307 0.112536 0.113302 0.0775336 0.0757236 0.0755278 0.0760334 0.0786192 0.0803361
-          0.0553322 0.0506627 0.0501228 0.0508766 0.0562484 0.059865 0.0406149 0.0286217 0.0279776 0.0288008 0.0400897 0.0469303
-          0.050324 0 0 0 0.0324822 0.0368562
+          0.987614 0.987617 0.984148 0.984152 0.987621 0.984162 0.987626 0.984169 0.987628 0.984174 0.987629 0.984176
+          0.960951 0.961012 0.961135 0.961237 0.961303 0.961326 0.85155 0.851752 0.852158 0.852503 0.852733 0.852813
+          0.63913 0.63945 0.640122 0.640687 0.641075 0.641222 0.41721 0.417503 0.418168 0.4187 0.419107 0.419285
+          0.260671 0.260839 0.261351 0.261746 0.262172 0.262378 0.166599 0.16653 0.166866 0.16718 0.16781 0.168164
+          0.11129 0.110701 0.110824 0.111167 0.112399 0.113173 0.0774079 0.0756056 0.0754375 0.0759335 0.0785082 0.0802413
+          0.0552474 0.0505833 0.0500691 0.0508009 0.0561273 0.0597859 0.0405625 0.0285734 0.0279549 0.0287506 0.0398965 0.0468695
+          0.0502388 0 0 0 0.0318739 0.036944
         </DataArray>
         <DataArray type="Float32" Name="pn" NumberOfComponents="1" format="ascii">
           100004 100003 100003 100003 100003 100003 100002 100002 100002 100002 100002 100002
@@ -67,13 +67,13 @@
           1.15426 1.16059 1.15749 1.1547 1.1542 1.15424
         </DataArray>
         <DataArray type="Float32" Name="mobW" NumberOfComponents="1" format="ascii">
-          1123.09 1123.09 1120.77 1120.77 1123.1 1120.78 1123.1 1120.78 1123.1 1120.79 1123.1 1120.79
-          1064.13 1064.34 1064.75 1065.1 1065.32 1065.39 649.057 649.521 650.447 651.232 651.757 651.938
-          264.454 264.822 265.595 266.244 266.694 266.865 73.9109 74.0504 74.3685 74.6242 74.8227 74.9111
-          18.527 18.5577 18.655 18.7313 18.8164 18.8581 4.96724 4.95987 4.9861 5.01234 5.06711 5.09806
-          1.50112 1.47711 1.48076 1.49435 1.54414 1.57573 0.502756 0.467855 0.464177 0.473713 0.524481 0.560066
-          0.177732 0.134743 0.130262 0.136546 0.187089 0.22714 0.0664738 0.0206968 0.0191333 0.0211457 0.0637252 0.105731
-          0.13192 0 0 0 0.0318195 0.0484146
+          1123.09 1123.09 1120.76 1120.76 1123.09 1120.77 1123.1 1120.78 1123.1 1120.78 1123.1 1120.78
+          1063.63 1063.86 1064.32 1064.7 1064.94 1065.02 647.964 648.47 649.485 650.35 650.926 651.126
+          263.598 263.999 264.846 265.558 266.048 266.234 73.5889 73.7418 74.0896 74.3689 74.583 74.6765
+          18.4373 18.4722 18.5789 18.6615 18.7507 18.794 4.94253 4.93651 4.96594 4.99352 5.04923 5.08071
+          1.49365 1.47015 1.47503 1.48871 1.53852 1.57038 0.500281 0.465635 0.462485 0.471818 0.522232 0.558061
+          0.176882 0.134079 0.129822 0.135906 0.185835 0.22621 0.0661962 0.0205768 0.0190798 0.0210193 0.0627331 0.105296
+          0.131217 0 0 0 0.0298588 0.0487949
         </DataArray>
         <DataArray type="Float32" Name="mobN" NumberOfComponents="1" format="ascii">
           0.321915 0.321323 2.01535 2.01243 0.320115 2.00658 0.319084 2.00169 0.318428 1.99858 0.318203 1.99751
@@ -184,48 +184,48 @@
           3 2 2 2 3 3
         </DataArray>
         <DataArray type="Float32" Name="velocityW" NumberOfComponents="3" format="ascii">
-          -1.33926e-08 9.46302e-09 0 -2.82555e-08 1.0444e-08 0 -1.20347e-08 1.80772e-08 0 -2.75886e-08 2.12658e-08 0
-          -5.56285e-08 8.87458e-09 0 -5.8923e-08 1.97974e-08 0 -6.66631e-08 -4.55131e-10 0 -7.43491e-08 2.26234e-10 0
-          -4.60163e-08 -1.226e-08 0 -5.22395e-08 -2.51467e-08 0 -2.6845e-08 -1.73792e-08 0 -3.04844e-08 -3.58103e-08 0
-          -6.0865e-09 4.13938e-08 0 -2.25194e-08 5.26921e-08 0 -6.423e-08 5.63357e-08 0 -9.48689e-08 1.49852e-08 0
-          -7.04167e-08 -4.41364e-08 0 -4.06033e-08 -6.65735e-08 0 5.85259e-09 8.64961e-08 0 -5.76298e-09 1.12762e-07 0
-          -4.5284e-08 1.29865e-07 0 -9.01227e-08 7.11737e-08 0 -7.70688e-08 -3.26554e-08 0 -4.70815e-08 -7.01039e-08 0
-          1.47258e-08 1.57514e-07 0 8.26302e-09 2.02855e-07 0 -2.15366e-08 2.36426e-07 0 -6.33077e-08 1.69639e-07 0
-          -6.64819e-08 2.53262e-08 0 -5.12219e-08 -3.78634e-08 0 1.72582e-08 2.29816e-07 0 1.31746e-08 2.94321e-07 0
-          -7.33337e-09 3.41306e-07 0 -3.77508e-08 2.71267e-07 0 -4.41669e-08 1.01815e-07 0 -3.659e-08 9.25544e-09 0
-          1.9018e-08 2.7865e-07 0 1.47163e-08 3.61384e-07 0 -2.42294e-09 4.17195e-07 0 -2.78048e-08 3.46198e-07 0
-          -3.3547e-08 1.6012e-07 0 -2.6745e-08 4.73989e-08 0 2.26805e-08 3.02055e-07 0 1.66879e-08 4.03986e-07 0
-          -9.32722e-10 4.65351e-07 0 -2.57824e-08 3.95109e-07 0 -3.26179e-08 1.94767e-07 0 -2.62317e-08 6.72129e-08 0
-          2.88125e-08 3.07371e-07 0 1.97386e-08 4.31568e-07 0 -4.00844e-10 4.95857e-07 0 -2.69156e-08 4.28162e-07 0
-          -3.59431e-08 2.14178e-07 0 -2.95214e-08 7.1861e-08 0 3.83995e-08 3.00164e-07 0 2.40887e-08 4.51806e-07 0
-          -1.83282e-10 5.15852e-07 0 -2.95482e-08 4.53475e-07 0 -4.19855e-08 2.25426e-07 0 -3.50191e-08 6.65428e-08 0
-          5.63209e-08 2.81927e-07 0 3.16761e-08 4.70096e-07 0 -1.61565e-10 5.28744e-07 0 -3.43751e-08 4.766e-07 0
-          -5.19321e-08 2.32699e-07 0 -4.24685e-08 5.44669e-08 0 1.28271e-07 2.40243e-07 0 6.50952e-08 4.86886e-07 0
-          -7.34205e-11 5.33805e-07 0 -5.58924e-08 5.0434e-07 0 -7.64869e-08 2.34332e-07 0 -4.32555e-08 3.88595e-08 0
-          9.08253e-07 2.10528e-07 0 4.54126e-07 4.94424e-07 0 0 5.33961e-07 0 -9.40328e-08 5.20468e-07 0
-          -1.04016e-07 2.32963e-07 0 -1.99658e-08 3.07125e-08 0
+          -1.30732e-08 9.25e-09 0 -2.76715e-08 1.02675e-08 0 -1.16934e-08 1.76487e-08 0 -2.69584e-08 2.09073e-08 0
+          -5.468e-08 8.80799e-09 0 -5.78856e-08 1.96572e-08 0 -6.57344e-08 -3.75216e-10 0 -7.33385e-08 3.99453e-10 0
+          -4.54643e-08 -1.20885e-08 0 -5.16533e-08 -2.47966e-08 0 -2.65499e-08 -1.71823e-08 0 -3.01772e-08 -3.54168e-08 0
+          -5.70422e-09 4.05673e-08 0 -2.17915e-08 5.1989e-08 0 -6.29716e-08 5.60496e-08 0 -9.36149e-08 1.53931e-08 0
+          -6.97367e-08 -4.33569e-08 0 -4.03081e-08 -6.57824e-08 0 6.12228e-09 8.53975e-08 0 -5.2015e-09 1.11794e-07 0
+          -4.42252e-08 1.29445e-07 0 -8.88734e-08 7.1775e-08 0 -7.6402e-08 -3.12882e-08 0 -4.69825e-08 -6.90193e-08 0
+          1.48249e-08 1.56263e-07 0 8.55022e-09 2.01681e-07 0 -2.08957e-08 2.35854e-07 0 -6.23947e-08 1.70228e-07 0
+          -6.58023e-08 2.71827e-08 0 -5.0882e-08 -3.65409e-08 0 1.7249e-08 2.28441e-07 0 1.32712e-08 2.92935e-07 0
+          -7.01839e-09 3.40533e-07 0 -3.7279e-08 2.71716e-07 0 -4.36245e-08 1.03761e-07 0 -3.6021e-08 1.09273e-08 0
+          1.89648e-08 2.77176e-07 0 1.47417e-08 3.59802e-07 0 -2.2444e-09 4.1623e-07 0 -2.74893e-08 3.46493e-07 0
+          -3.32992e-08 1.62012e-07 0 -2.66271e-08 4.92883e-08 0 2.25984e-08 3.00537e-07 0 1.66836e-08 4.02254e-07 0
+          -7.98463e-10 4.6424e-07 0 -2.55019e-08 3.9525e-07 0 -3.24524e-08 1.96719e-07 0 -2.62668e-08 6.90945e-08 0
+          2.86971e-08 3.0586e-07 0 1.97102e-08 4.29721e-07 0 -2.8436e-10 4.94629e-07 0 -2.66125e-08 4.28124e-07 0
+          -3.57906e-08 2.16274e-07 0 -2.96482e-08 7.3702e-08 0 3.82424e-08 2.98702e-07 0 2.40288e-08 4.49868e-07 0
+          -8.45755e-11 5.14509e-07 0 -2.91761e-08 4.53166e-07 0 -4.18177e-08 2.27778e-07 0 -3.52677e-08 6.83556e-08 0
+          5.61098e-08 2.80548e-07 0 3.15696e-08 4.68082e-07 0 -9.76324e-11 5.27266e-07 0 -3.38361e-08 4.75781e-07 0
+          -5.17156e-08 2.35528e-07 0 -4.29838e-08 5.62569e-08 0 1.27756e-07 2.39026e-07 0 6.48367e-08 4.84221e-07 0
+          -5.76806e-11 5.32163e-07 0 -5.4557e-08 5.0224e-07 0 -7.59459e-08 2.38347e-07 0 -4.48105e-08 4.04185e-08 0
+          9.02879e-07 2.09426e-07 0 4.51439e-07 4.9114e-07 0 0 5.32231e-07 0 -8.75259e-08 5.17429e-07 0
+          -9.8865e-08 2.37849e-07 0 -2.26782e-08 3.20543e-08 0
         </DataArray>
         <DataArray type="Float32" Name="velocityN" NumberOfComponents="3" format="ascii">
-          4.51984e-10 -8.30227e-10 0 6.82785e-10 -6.8278e-10 0 2.91512e-09 -3.61333e-09 0 4.37192e-09 -2.53747e-09 0
-          8.39621e-10 -4.79636e-10 0 5.31818e-09 -1.08119e-09 0 6.23034e-10 -4.34058e-10 0 3.91285e-09 -7.9001e-10 0
-          3.22061e-10 -4.23635e-10 0 2.02356e-09 -7.15634e-10 0 1.63708e-10 -4.20078e-10 0 1.02908e-09 -6.89877e-10 0
-          7.80394e-08 -6.26795e-08 0 1.16185e-07 -3.75785e-08 0 1.39797e-07 -4.35256e-09 0 1.01822e-07 1.36007e-09 0
-          5.25102e-08 2.91738e-09 0 2.66389e-08 3.41499e-09 0 1.91031e-06 -1.15733e-06 0 2.84393e-06 -5.9512e-07 0
-          3.42151e-06 1.46742e-07 0 2.49579e-06 2.70225e-07 0 1.28962e-06 3.08417e-07 0 6.53135e-07 3.18149e-07 0
-          1.28935e-05 -8.40676e-06 0 1.95817e-05 -4.36193e-06 0 2.33501e-05 1.64359e-06 0 1.65717e-05 2.33447e-06 0
-          8.5075e-06 2.57393e-06 0 4.30187e-06 2.63869e-06 0 3.69309e-05 -3.06954e-05 0 5.72152e-05 -1.64458e-05 0
-          6.53694e-05 8.18748e-06 0 4.35587e-05 8.31039e-06 0 2.25871e-05 9.40634e-06 0 1.12962e-05 9.52287e-06 0