Merge from mainline.

git-svn-id: https://svn.dealii.org/branches/branch_port_the_testsuite@30918 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/examples/step-42/doc/intro.dox b/deal.II/examples/step-42/doc/intro.dox
index 7e1920895bd8b5d06aff4c00d358c66a1b22efa8..68f02e88b4313462bac00612ed08e690f7f4f239 100644 (file)
--- a/deal.II/examples/step-42/doc/intro.dox
+++ b/deal.II/examples/step-42/doc/intro.dox
@@ -392,5 +392,5 @@ Compared to step-41, the programs has two new classes:
   looks as follows (see also the picture at
   the top of this section on how the end result looks like):
 
-  <img src="http://www.dealii.org/images/steps/developer/step-42.character.png" alt="">
+  <img src="http://www.dealii.org/images/steps/developer/step-42.character.png" alt="" width="25%">
 </ul>

diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc
index 275c8c3623fc77171a91509712cf690723ac1943..2613a848b2d482a580a79fd4e9a48648d792527c 100644 (file)
--- a/deal.II/examples/step-42/step-42.cc
+++ b/deal.II/examples/step-42/step-42.cc
@@ -85,160 +85,63 @@ namespace Step42
 
 // @sect3{The <code>Input</code> class template}
 
-// This class has the the only purpose
-// to read in data from a picture file
-// stored in pbm ascii
-// format. This data will be bilinearly
-// interpolated and provides in this way
-// a function which describes an obstacle.
+// This class has the the only purpose to read in data from a picture file
+// stored in pbm ascii format. This data will be bilinearly interpolated and
+// provides in this way a function which describes an obstacle.
 //
-// The data which we read in by the
-// function read_obstacle () from the file
-// "obstacle_file.pbm" will be stored
-// in a double std::vector named
-// obstacle_data.
-// This vector composes the base
-// to calculate a piecewise bilinear
-// function as a polynomial interpolation.
-// This will be done by obstacle_function ().
+// The data which we read from the file will be stored in a double std::vector
+// named obstacle_data.  This vector composes the base to calculate a
+// piecewise bilinear function as a polynomial interpolation.  This will be
+// done by obstacle_function ().
 //
-// In the function run () of the class
-// PlasticityContactProblem we create
-// an object of the class Input which will
-// be used in the class Obstacle to
-// supply the obstacle function in
-// update_solution_and_constraints () of
-// the class PlasticityContactProblem.
-
+// In the function <code>run()</code> of the class
+// <code>PlasticityContactProblem</code> we create an object of the this class
+// which will be used in the class Obstacle to supply the obstacle function in
+// <code>update_solution_and_constraints()</code> of the class
+// <code>PlasticityContactProblem</code>.
+//
+// The <code>hx,hy</code> variables denote the spacing between pixels in $x$
+// and $y$ directions. <code>nx,ny</code> are the numbers of pixels in each of
+// these directions.  <code>get_value()</code> returns the value of the image
+// at a given location, interpolated from the adjacent pixel values.
   template <int dim>
   class Input
   {
   public:
-    Input (const std::string &name)
-      :
-      pcout(std::cout,
-            (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)),
-      obstacle_data(0),
-      hx(0),
-      hy(0),
-      nx(0),
-      ny(0)
-    {
-      read_obstacle(name);
-    }
-
-    double hv (const int i,
-               const int j);
+    Input (const std::string &name);
 
     double
-    obstacle_function (const double x,
-                       const double y);
-
-    void
-    read_obstacle (const std::string name);
+    get_value (const double x,
+               const double y);
 
   private:
-    ConditionalOStream pcout;
     std::vector<double> obstacle_data;
-    double hx, hy;
-    int nx, ny;
-  };
-
-// This function is used in obstacle_function ()
-// to provide the proper value of the obstacle.
-  template <int dim>
-  double
-  Input<dim>::hv (const int i,
-                  const int j)
-  {
-    assert(i >= 0 && i < nx);
-    assert(j >= 0 && j < ny);
-    return obstacle_data[nx * (ny - 1 - j) + i]; // i indiziert x-werte, j indiziert y-werte
-  }
-
-// obstacle_function () calculates the bilinear interpolated
-// value in the point (x,y).
-  template <int dim>
-  double
-  Input<dim>::obstacle_function (const double x,
-                                 const double y)
-  {
-    int ix = (int) (x / hx);
-    int iy = (int) (y / hy);
+    double              hx, hy;
+    int                 nx, ny;
 
-    if (ix < 0)
-      ix = 0;
-
-    if (iy < 0)
-      iy = 0;
-
-    if (ix >= nx - 1)
-      ix = nx - 2;
-
-    if (iy >= ny - 1)
-      iy = ny - 2;
-
-    double val = 0.0;
-    {
-      FullMatrix<double> H(4, 4);
-      Vector<double> X(4);
-      Vector<double> b(4);
-
-      double xx, yy;
-
-      xx = ix * hx;
-      yy = iy * hy;
-      H(0, 0) = xx;
-      H(0, 1) = yy;
-      H(0, 2) = xx * yy;
-      H(0, 3) = 1.0;
-      b(0) = hv(ix, iy);
-
-      xx = (ix + 1) * hx;
-      yy = iy * hy;
-      H(1, 0) = xx;
-      H(1, 1) = yy;
-      H(1, 2) = xx * yy;
-      H(1, 3) = 1.0;
-      b(1) = hv(ix + 1, iy);
-
-      xx = (ix + 1) * hx;
-      yy = (iy + 1) * hy;
-      H(2, 0) = xx;
-      H(2, 1) = yy;
-      H(2, 2) = xx * yy;
-      H(2, 3) = 1.0;
-      b(2) = hv(ix + 1, iy + 1);
-
-      xx = ix * hx;
-      yy = (iy + 1) * hy;
-      H(3, 0) = xx;
-      H(3, 1) = yy;
-      H(3, 2) = xx * yy;
-      H(3, 3) = 1.0;
-      b(3) = hv(ix, iy + 1);
-
-      H.gauss_jordan();
-      H.vmult(X, b);
-
-      val = X(0) * x + X(1) * y + X(2) * x * y + X(3);
-    }
+    double get_pixel_value (const int i,
+                            const int j);
+  };
 
-    return val;
-  }
 
-// As mentioned above this function reads in the
-// obstacle data and stores them in the std::vector
-// obstacle_data. It will be used only in run ().
-  template <int dim>
-  void
-  Input<dim>::read_obstacle (const std::string name)
+  // The constructor of this class reads in the data that describes
+  // the obstacle from the given file name.
+  template<int dim>
+  Input<dim>::Input(const std::string &name)
+    :
+    obstacle_data(0),
+    hx(0),
+    hy(0),
+    nx(0),
+    ny(0)
   {
     std::ifstream f(name.c_str());
 
     std::string temp;
     f >> temp >> nx >> ny;
-    assert(nx > 0 && ny > 0);
+
+    AssertThrow(nx > 0 && ny > 0,
+                ExcMessage ("Invalid file format."));
 
     for (int k = 0; k < nx * ny; k++)
       {
@@ -250,122 +153,184 @@ namespace Step42
     hx = 1.0 / (nx - 1);
     hy = 1.0 / (ny - 1);
 
-    pcout << "Resolution of the scanned obstacle picture: " << nx << " x "
-          << ny << std::endl;
+    if (Utilities::MPI::this_mpi_process(MPI_COMM_WORLD) == 0)
+      std::cout << "Resolution of the scanned obstacle picture: " << nx << " x " << ny
+                << std::endl;
+  }
+
+  template <int dim>
+  double
+  Input<dim>::get_pixel_value (const int i,
+                               const int j)
+  {
+    assert(i >= 0 && i < nx);
+    assert(j >= 0 && j < ny);
+    return obstacle_data[nx * (ny - 1 - j) + i];
+  }
+
+
+  template <int dim>
+  double
+  Input<dim>::get_value (const double x,
+                         const double y)
+  {
+    const int ix = std::min(std::max((int) (x / hx), 0), nx-2);
+    const int iy = std::min(std::max((int) (y / hy), 0), ny-2);
+
+    FullMatrix<double> H(4, 4);
+    Vector<double> X(4);
+    Vector<double> b(4);
+
+    double xx, yy;
+
+    xx = ix * hx;
+    yy = iy * hy;
+    H(0, 0) = xx;
+    H(0, 1) = yy;
+    H(0, 2) = xx * yy;
+    H(0, 3) = 1.0;
+    b(0) = get_value(ix, iy);
+
+    xx = (ix + 1) * hx;
+    yy = iy * hy;
+    H(1, 0) = xx;
+    H(1, 1) = yy;
+    H(1, 2) = xx * yy;
+    H(1, 3) = 1.0;
+    b(1) = get_pixel_value(ix + 1, iy);
+
+    xx = (ix + 1) * hx;
+    yy = (iy + 1) * hy;
+    H(2, 0) = xx;
+    H(2, 1) = yy;
+    H(2, 2) = xx * yy;
+    H(2, 3) = 1.0;
+    b(2) = get_pixel_value(ix + 1, iy + 1);
+
+    xx = ix * hx;
+    yy = (iy + 1) * hy;
+    H(3, 0) = xx;
+    H(3, 1) = yy;
+    H(3, 2) = xx * yy;
+    H(3, 3) = 1.0;
+    b(3) = get_pixel_value(ix, iy + 1);
+
+    H.gauss_jordan();
+    H.vmult(X, b);
+
+    return (X(0) * x + X(1) * y + X(2) * x * y + X(3));
   }
 
-// @sect3{The <code>ConstitutiveLaw</code> class template}
 
-// This class provides an interface
-// for a constitutive law. In this
-// example we are using an elastoplastic
-// material behavior with linear,
-// isotropic hardening.
-// For $\gamma = 0$ we obtain perfect elastoplastic
-// behavior.
+  // @sect3{The <code>ConstitutiveLaw</code> class template}
+
+  // This class provides an interface for a constitutive law, i.e., for the
+  // relationship between strain $\varepsilon(\mathbf u)$ and stress
+  // $\sigma$. In this example we are using an elastoplastic material behavior
+  // with linear, isotropic hardening. Such materials are characterized by
+  // Young's modulus $E$, Poisson's ratio $\nu$, the initial yield stress
+  // $\sigma_0$ and the isotropic hardening parameter $\gamma$.  For $\gamma =
+  // 0$ we obtain perfect elastoplastic behavior.
+  //
+  // As explained in the paper that describes this program, the first Newton
+  // steps are solved with a completely elastic material model to avoid having
+  // to deal with both nonlinearities (plasticity and contact) at once. To this
+  // end, this class has a function <code>set_sigma_0()</code> that we use later
+  // on to simply set $\sigma_0$ to a very large value -- essentially
+  // guaranteeing that the actual stress will not exceed it, and thereby
+  // producing an elastic material. When we are ready to use a plastic model, we
+  // set $\sigma_0$ back to its proper value, using the same function.  As a
+  // result of this approach, we need to leave <code>sigma_0</code> as the only
+  // non-const member variable of this class.
   template <int dim>
   class ConstitutiveLaw
   {
   public:
-    ConstitutiveLaw (const double _E,
-                     const double _nu,
-                     const double _sigma_0,
-                     const double _gamma);
-
-    bool
-    plast_linear_hardening (SymmetricTensor<4, dim> &stress_strain_tensor,
-                            const SymmetricTensor<2, dim> &strain_tensor) const;
+    ConstitutiveLaw (const double E,
+                     const double nu,
+                     const double sigma_0,
+                     const double gamma);
 
     void
-    linearized_plast_linear_hardening (SymmetricTensor<4, dim> &stress_strain_tensor_linearized,
-      SymmetricTensor<4, dim> &stress_strain_tensor,
-      const SymmetricTensor<2, dim> &strain_tensor) const;
+    set_sigma_0 (double sigma_zero);
 
-    SymmetricTensor<2, dim>
-    get_strain (const FEValues<dim> &fe_values,
-      const unsigned int shape_func,
-      const unsigned int q_point) const;
+    bool
+    get_stress_strain_tensor (const SymmetricTensor<2, dim> &strain_tensor,
+                              SymmetricTensor<4, dim> &stress_strain_tensor) const;
 
     void
-    set_sigma_0 (double sigma_zero)
-    {
-      sigma_0 = sigma_zero;
-    }
+    get_linearized_stress_strain_tensors (const SymmetricTensor<2, dim> &strain_tensor,
+                                          SymmetricTensor<4, dim> &stress_strain_tensor_linearized,
+                                          SymmetricTensor<4, dim> &stress_strain_tensor) const;
 
   private:
-    const double E;
-    const double nu;
+    const double kappa;
+    const double mu;
     double       sigma_0;
     const double gamma;
-    const double mu;
-    const double kappa;
 
     const SymmetricTensor<4, dim> stress_strain_tensor_kappa;
     const SymmetricTensor<4, dim> stress_strain_tensor_mu;
   };
 
-// The constructor of the ConstitutiveLaw class sets the
-// required material parameter for our deformable body:
-// E -> elastic modulus
-// nu -> Passion's number
-// sigma_0 -> yield stress
-// gamma -> hardening parameter.
-// Also it supplies the stress strain tensor of forth order
-// of the volumetric and deviator part. For further details
-// see the documentation above.
+  // The constructor of the ConstitutiveLaw class sets the required material
+  // parameter for our deformable body. Material parameters for elastic
+  // isotropic media can be defined in a variety of ways, such as the pair $E,
+  // \nu$ (elastic modulus and Poisson's number), using the Lame parameters
+  // $\lambda,mu$ or several other commonly used conventions. Here, the
+  // constructor takes a description of material parameters in the form of
+  // $E,\nu$, but since this turns out to these are not the coefficients that
+  // appear in the equations of the plastic projector, we immediately convert
+  // them into the more suitable set $\kappa,\mu$ of bulk and shear moduli.  In
+  // addition, the constructor takes $\sigma_0$ (the yield stress absent any
+  // plastic strain) and $\gamma$ (the hardening parameter) as arguments. In
+  // this constructor, we also compute the two principal components of the
+  // stress-strain relation and its linearization.
   template <int dim>
-  ConstitutiveLaw<dim>::ConstitutiveLaw (double _E,
-                                         double _nu,
-                                         double _sigma_0,
-                                         double _gamma)
+  ConstitutiveLaw<dim>::ConstitutiveLaw (double E,
+                                         double nu,
+                                         double sigma_0,
+                                         double gamma)
     :
-    E(_E),
-    nu(_nu),
-    sigma_0(_sigma_0),
-    gamma(_gamma),
-    mu (E / (2 * (1 + nu))),
     kappa (E / (3 * (1 - 2 * nu))),
+    mu (E / (2 * (1 + nu))),
+    sigma_0(sigma_0),
+    gamma(gamma),
     stress_strain_tensor_kappa (kappa
-                                 * outer_product(unit_symmetric_tensor<dim>(),
-                                                 unit_symmetric_tensor<dim>())),
+                                * outer_product(unit_symmetric_tensor<dim>(),
+                                                unit_symmetric_tensor<dim>())),
     stress_strain_tensor_mu (2 * mu
-                              * (identity_tensor<dim>()
-                                 - outer_product(unit_symmetric_tensor<dim>(),
-                                                 unit_symmetric_tensor<dim>()) / 3.0))
+                             * (identity_tensor<dim>()
+                                - outer_product(unit_symmetric_tensor<dim>(),
+                                                unit_symmetric_tensor<dim>()) / 3.0))
   {}
 
-// @sect4{ConstitutiveLaw::ConstitutiveLaw}
 
-// Calculates the strain $\varepsilon(\varphi)=\dfrac{1}{2}\left(\nabla\varphi + \nabla\varphi^T$
-// for the shape functions $\varphi$.
   template <int dim>
-  SymmetricTensor<2, dim>
-  ConstitutiveLaw<dim>::get_strain (const FEValues<dim> &fe_values,
-                                    const unsigned int shape_func,
-                                    const unsigned int q_point) const
+  void
+  ConstitutiveLaw<dim>::set_sigma_0 (double sigma_zero)
   {
-    const FEValuesExtractors::Vector displacement(0);
-    return fe_values[displacement].symmetric_gradient(shape_func, q_point);
+    sigma_0 = sigma_zero;
   }
 
-// @sect4{ConstitutiveLaw::plast_linear_hardening}
 
-// This is the implemented constitutive law. It projects the
-// deviatoric part of the stresses in a quadrature point back to
-// the yield stress (i.e., the original yield stress $\sigma_0$ plus
-// the term that describes linear isotropic hardening).
-// We need this function to calculate the nonlinear
-// residual in
-// PlasticityContactProblem::residual_nl_system(TrilinosWrappers::MPI::Vector &u).
-//
-// The function returns whether the quadrature point is plastic to allow for
-// some statistics downstream on how many of the quadrature points are
-// plastic and how many are elastic.
+  // @sect4{ConstitutiveLaw::get_stress_strain_tensor}
+
+  // This is the principal component of the constitutive law. It projects the
+  // deviatoric part of the stresses in a quadrature point back to the yield
+  // stress (i.e., the original yield stress $\sigma_0$ plus the term that
+  // describes linear isotropic hardening).  We need this function to calculate
+  // the nonlinear residual in PlasticityContactProblem::residual_nl_system. The
+  // computations follow the formulas laid out in the introduction.
+  //
+  // The function returns whether the quadrature point is plastic to allow for
+  // some statistics downstream on how many of the quadrature points are
+  // plastic and how many are elastic.
   template <int dim>
   bool
   ConstitutiveLaw<dim>::
-  plast_linear_hardening (SymmetricTensor<4, dim> &stress_strain_tensor,
-                                 const SymmetricTensor<2, dim> &strain_tensor) const
+  get_stress_strain_tensor (const SymmetricTensor<2, dim> &strain_tensor,
+                            SymmetricTensor<4, dim> &stress_strain_tensor) const
   {
     Assert (dim == 3, ExcNotImplemented());
 
@@ -388,24 +353,24 @@ namespace Step42
     return (deviator_stress_tensor_norm > sigma_0);
   }
 
-// @sect4{ConstitutiveLaw::linearized_plast_linear_hardening}
-
-// This function returns the linearized stress strain tensor, linearized
-// around the solution $u^{i-1}$ of the previous Newton step $i-1$.
-// The parameter strain_tensor $\varepsilon(u^{i-1})$ is calculated
-// by $u^{i-1}$. It contains the derivative of the nonlinear
-// constitutive law. As the result this function returns
-// the stress_strain_tensor of the nonlinear problem as well as
-// the stress_strain_tensor_linearized of the linearized problem.
-// See
-// PlasticityContactProblem::assemble_nl_system(TrilinosWrappers::MPI::Vector &u)
-// where this function is used.
+
+  // @sect4{ConstitutiveLaw::get_linearized_stress_strain_tensors}
+
+  // This function returns the linearized stress strain tensor, linearized
+  // around the solution $u^{i-1}$ of the previous Newton step $i-1$.  The
+  // parameter <code>strain_tensor</code> (commonly denoted
+  // $\varepsilon(u^{i-1})$) must be passed as an argument, and serves as the
+  // linearization point. The function returns the derivative of the nonlinear
+  // constitutive law in the variable stress_strain_tensor, as well as the
+  // stress-strain tensor of the linearized problem in
+  // stress_strain_tensor_linearized.  See
+  // PlasticityContactProblem::assemble_nl_system where this function is used.
   template <int dim>
   void
   ConstitutiveLaw<dim>::
-  linearized_plast_linear_hardening (SymmetricTensor<4, dim> &stress_strain_tensor_linearized,
-                                     SymmetricTensor<4, dim> &stress_strain_tensor,
-                                     const SymmetricTensor<2, dim> &strain_tensor) const
+  get_linearized_stress_strain_tensors (const SymmetricTensor<2, dim> &strain_tensor,
+                                        SymmetricTensor<4, dim> &stress_strain_tensor_linearized,
+                                        SymmetricTensor<4, dim> &stress_strain_tensor) const
   {
     Assert (dim == 3, ExcNotImplemented());
 
@@ -434,142 +399,133 @@ namespace Step42
     stress_strain_tensor_linearized += stress_strain_tensor_kappa;
   }
 
+  // <h3>Equation data: right hand side, boundary values, obstacles</h3>
+  //
+  // The following should be relatively standard. We need classes for
+  // the right hand side forcing term (which we here choose to be zero)
+  // and boundary values on those part of the boundary that are not part
+  // of the contact surface (also chosen to be zero here).
   namespace EquationData
   {
-// It possible to apply an additional body force
-// but in here it is set to zero.
     template <int dim>
     class RightHandSide : public Function<dim>
     {
     public:
-      RightHandSide ()
-        :
-        Function<dim>(dim)
-      {
-      }
+      RightHandSide ();
 
-      virtual double
-      value (
-        const Point<dim> &p, const unsigned int component = 0) const;
+      virtual
+      double value (const Point<dim> &p,
+                    const unsigned int component = 0) const;
 
-      virtual void
-      vector_value (
-        const Point<dim> &p, Vector<double> &values) const;
+      virtual
+      void vector_value (const Point<dim> &p,
+                         Vector<double> &values) const;
     };
 
     template <int dim>
-    double
-    RightHandSide<dim>::value (
-      const Point<dim> &p, const unsigned int component) const
-    {
-      double return_value = 0.0;
+    RightHandSide<dim>::RightHandSide ()
+      :
+      Function<dim>(dim)
+    {}
 
-      if (component == 0)
-        return_value = 0.0;
-      if (component == 1)
-        return_value = 0.0;
-      if (component == 2)
-        return_value = 0.0;
 
-      return return_value;
+    template <int dim>
+    double
+    RightHandSide<dim>::value (const Point<dim> &,
+                               const unsigned int) const
+    {
+      return 0.;
     }
 
     template <int dim>
     void
-    RightHandSide<dim>::vector_value (
-      const Point<dim> &p, Vector<double> &values) const
+    RightHandSide<dim>::vector_value (const Point<dim> &p,
+                                      Vector<double> &values) const
     {
       for (unsigned int c = 0; c < this->n_components; ++c)
         values(c) = RightHandSide<dim>::value(p, c);
     }
 
-// This function class is used to describe the prescribed displacements
-// at the boundary. But again we set this to zero.
+
+
     template <int dim>
     class BoundaryValues : public Function<dim>
     {
     public:
-      BoundaryValues ()
-        :
-        Function<dim>(dim)
-      {
-      }
-      ;
+      BoundaryValues ();
 
-      virtual double
-      value (
-        const Point<dim> &p, const unsigned int component = 0) const;
+      virtual double value (const Point<dim> &p,
+                            const unsigned int component = 0) const;
 
-      virtual void
-      vector_value (
-        const Point<dim> &p, Vector<double> &values) const;
+      virtual
+      void vector_value (const Point<dim> &p,
+                         Vector<double> &values) const;
     };
 
+
     template <int dim>
-    double
-    BoundaryValues<dim>::value (
-      const Point<dim> &p, const unsigned int component) const
-    {
-      double return_value = 0;
+    BoundaryValues<dim>::BoundaryValues ()
+      :
+      Function<dim>(dim)
+    {}
 
-      if (component == 0)
-        return_value = 0.0;
-      if (component == 1)
-        return_value = 0.0;
-      if (component == 2)
-        return_value = 0.0;
 
-      return return_value;
+    template <int dim>
+    double
+    BoundaryValues<dim>::value (const Point<dim> &,
+                                const unsigned int) const
+    {
+      return 0.;
     }
 
     template <int dim>
     void
-    BoundaryValues<dim>::vector_value (
-      const Point<dim> &p, Vector<double> &values) const
+    BoundaryValues<dim>::vector_value (const Point<dim> &p,
+                                       Vector<double> &values) const
     {
       for (unsigned int c = 0; c < this->n_components; ++c)
         values(c) = BoundaryValues<dim>::value(p, c);
     }
 
-// This function is obviously implemented to
-// define the obstacle that penetrates our deformable
-// body. You can choose between two ways to define
-// your obstacle: to read it from a file or to use
-// a function (here a ball).
-// z_max_domain is the z value of the surface of the work piece
+    // This function is obviously implemented to
+    // define the obstacle that penetrates our deformable
+    // body. You can choose between two ways to define
+    // your obstacle: to read it from a file or to use
+    // a function (here a ball).
+    // z_max_domain is the z value of the surface of the work piece
     template <int dim>
     class Obstacle : public Function<dim>
     {
     public:
-      Obstacle (
-        std_cxx1x::shared_ptr<Input<dim> > const &_input,
-        bool _use_read_obstacle, double z_max_domain)
+      Obstacle (const std_cxx1x::shared_ptr<Input<dim> > &input,
+                const bool use_read_obstacle,
+                const double z_max_domain)
         :
         Function<dim>(dim),
-        input_obstacle_copy(_input),
-        use_read_obstacle(_use_read_obstacle),
+        input_obstacle(input),
+        use_read_obstacle(use_read_obstacle),
         z_max_domain(z_max_domain)
-      {
-      }
+      {}
 
-      virtual double
-      value (
-        const Point<dim> &p, const unsigned int component = 0) const;
+      virtual
+      double value (const Point<dim> &p,
+                    const unsigned int component = 0) const;
 
-      virtual void
-      vector_value (
-        const Point<dim> &p, Vector<double> &values) const;
+      virtual
+      void vector_value (const Point<dim> &p,
+                         Vector<double> &values) const;
 
     private:
-      std_cxx1x::shared_ptr<Input<dim> > const &input_obstacle_copy;
+      const std_cxx1x::shared_ptr<Input<dim> > &input_obstacle;
       bool use_read_obstacle;
       double z_max_domain;
     };
 
+
     template <int dim>
     double
-    Obstacle<dim>::value (
-      const Point<dim> &p, const unsigned int component) const
+    Obstacle<dim>::value (const Point<dim> &p,
+                          const unsigned int component) const
     {
       if (component == 0)
         return p(0);
@@ -581,92 +537,68 @@ namespace Step42
         {
           if (p(0) >= 0.0 && p(0) <= 1.0 && p(1) >= 0.0 && p(1) <= 1.0)
             return z_max_domain + 0.999
-                   - input_obstacle_copy->obstacle_function(p(0), p(1));
+                   - input_obstacle->get_value(p(0), p(1));
           else
             return 10000.0;
         }
       else
         {
           //sphere:
-          return -std::sqrt(
-                   0.36 - (p(0) - 0.5) * (p(0) - 0.5)
-                   - (p(1) - 0.5) * (p(1) - 0.5)) + z_max_domain + 0.59;
+          return -std::sqrt(0.36 - (p(0) - 0.5) * (p(0) - 0.5)
+                            - (p(1) - 0.5) * (p(1) - 0.5)) + z_max_domain + 0.59;
         }
     }
 
     template <int dim>
     void
-    Obstacle<dim>::vector_value (
-      const Point<dim> &p, Vector<double> &values) const
+    Obstacle<dim>::vector_value (const Point<dim> &p,
+                                 Vector<double> &values) const
     {
       for (unsigned int c = 0; c < this->n_components; ++c)
         values(c) = Obstacle<dim>::value(p, c);
     }
   }
 
-// @sect3{The <code>PlasticityContactProblem</code> class template}
-
-// This class supplies all function
-// and variables needed to describe
-// the nonlinear contact problem. It is
-// close to step-41 but with some additional
-// features like: handling hanging nodes,
-// a newton method, using Trilinos and p4est
-// for parallel distributed computing.
-// To deal with hanging nodes makes
-// life a bit more complicated since
-// we need an other ConstraintMatrix now.
-// We create a newton method for the
-// active set method for the contact
-// situation and to handle the nonlinear
-// operator for the constitutive law.
-
+  // @sect3{The <code>PlasticityContactProblem</code> class template}
+
+  // This class supplies all function
+  // and variables needed to describe
+  // the nonlinear contact problem. It is
+  // close to step-41 but with some additional
+  // features like: handling hanging nodes,
+  // a newton method, using Trilinos and p4est
+  // for parallel distributed computing.
+  // To deal with hanging nodes makes
+  // life a bit more complicated since
+  // we need an other ConstraintMatrix now.
+  // We create a newton method for the
+  // active set method for the contact
+  // situation and to handle the nonlinear
+  // operator for the constitutive law.
   template <int dim>
   class PlasticityContactProblem
   {
   public:
-    PlasticityContactProblem (
-      const ParameterHandler &prm);
-    void
-    run ();
+    PlasticityContactProblem (const ParameterHandler &prm);
 
-    static void
-    declare (
-      ParameterHandler &prm);
+    void run ();
+
+    static void declare (ParameterHandler &prm);
 
   private:
-    void
-    make_grid ();
-    void
-    setup_system ();
-    void
-    assemble_nl_system (
-      TrilinosWrappers::MPI::Vector &u);
-    void
-    residual_nl_system (
-      TrilinosWrappers::MPI::Vector &u);
-    void
-    assemble_mass_matrix_diagonal (
-      TrilinosWrappers::SparseMatrix &mass_matrix);
-    void
-    update_solution_and_constraints ();
-    void
-    dirichlet_constraints ();
-    void
-    solve ();
-    void
-    solve_newton ();
-    void
-    refine_grid ();
-    void
-    move_mesh (
-      const TrilinosWrappers::MPI::Vector &_complete_displacement) const;
-    void
-    output_results (
-      const std::string &title);
-    void
-    output_contact_force (
-      const unsigned int cycle);
+    void make_grid ();
+    void setup_system ();
+    void assemble_nl_system (TrilinosWrappers::MPI::Vector &u);
+    void residual_nl_system (TrilinosWrappers::MPI::Vector &u);
+    void assemble_mass_matrix_diagonal (TrilinosWrappers::SparseMatrix &mass_matrix);
+    void update_solution_and_constraints ();
+    void dirichlet_constraints ();
+    void solve ();
+    void solve_newton ();
+    void refine_grid ();
+    void move_mesh (const TrilinosWrappers::MPI::Vector &_complete_displacement) const;
+    void output_results (const std::string &title);
+    void output_contact_force (const unsigned int cycle);
 
     double to_refine_factor;
     double to_coarsen_factor;
@@ -716,7 +648,7 @@ namespace Step42
 
     double sigma_0; // Yield stress
     double gamma; // Parameter for the linear isotropic hardening
-    double e_modul; // E-Modul
+    double e_modulus; // E-Modul
     double nu; // Poisson ratio
 
     TimerOutput computing_timer;
@@ -759,13 +691,13 @@ namespace Step42
           (Utilities::MPI::this_mpi_process(mpi_communicator) == 0)),
     sigma_0(400.0),
     gamma(0.01),
-    e_modul(2.0e+5),
+    e_modulus(2.0e+5),
     nu(0.3),
     computing_timer(MPI_COMM_WORLD, pcout, TimerOutput::never,
                     TimerOutput::wall_times)
   {
     // double _E, double _nu, double _sigma_0, double _gamma
-    plast_lin_hard.reset(new ConstitutiveLaw<dim>(e_modul, nu, sigma_0, gamma));
+    plast_lin_hard.reset(new ConstitutiveLaw<dim>(e_modulus, nu, sigma_0, gamma));
 
     degree = prm.get_integer("polynomial degree");
     n_initial_refinements = prm.get_integer("number of initial refinements");
@@ -1039,19 +971,19 @@ namespace Step42
               SymmetricTensor<4, dim> stress_strain_tensor;
               SymmetricTensor<2, dim> stress_tensor;
 
-              plast_lin_hard->linearized_plast_linear_hardening(
-                stress_strain_tensor_linearized, stress_strain_tensor,
-                strain_tensor[q_point]);
+              plast_lin_hard->get_linearized_stress_strain_tensors(strain_tensor[q_point],
+                                                                   stress_strain_tensor_linearized,
+                                                                   stress_strain_tensor);
 
               for (unsigned int i = 0; i < dofs_per_cell; ++i)
                 {
                   stress_tensor = stress_strain_tensor_linearized
-                                  * plast_lin_hard->get_strain(fe_values, i, q_point);
+                                  * fe_values[displacement].symmetric_gradient(i, q_point);
 
                   for (unsigned int j = 0; j < dofs_per_cell; ++j)
                     {
                       cell_matrix(i, j) += (stress_tensor
-                                            * plast_lin_hard->get_strain(fe_values, j, q_point)
+                                            * fe_values[displacement].symmetric_gradient(i, q_point)
                                             * fe_values.JxW(q_point));
                     }
 
@@ -1062,7 +994,7 @@ namespace Step42
                   // the residual part a(v^i;v) of the rhs
                   cell_rhs(i) -= (strain_tensor[q_point]
                                   * stress_strain_tensor
-                                  * plast_lin_hard->get_strain(fe_values, i, q_point)
+                                  * fe_values[displacement].symmetric_gradient(i, q_point)
                                   * fe_values.JxW(q_point));
 
                   // the residual part F(v) of the rhs
@@ -1169,22 +1101,21 @@ namespace Step42
               SymmetricTensor<2, dim> stress_tensor;
 
               const bool q_point_is_plastic
-              = plast_lin_hard->plast_linear_hardening(stress_strain_tensor,
-                                                     strain_tensor[q_point]);
+                = plast_lin_hard->get_stress_strain_tensor(strain_tensor[q_point],
+                                                           stress_strain_tensor);
               if (q_point_is_plastic)
-              {
-                 ++plast_points;
+                {
+                  ++plast_points;
                   ++cell_constitution(cell_number);
-              }
+                }
               else
-                 ++elast_points;
+                ++elast_points;
 
               for (unsigned int i = 0; i < dofs_per_cell; ++i)
                 {
                   cell_rhs(i) -= (strain_tensor[q_point]
                                   * stress_strain_tensor
-                                  * //(stress_tensor) *
-                                  plast_lin_hard->get_strain(fe_values, i, q_point)
+                                  * fe_values[displacement].symmetric_gradient(i, q_point)
                                   * fe_values.JxW(q_point));
 
                   Tensor<1, dim> rhs_values;
@@ -1203,7 +1134,7 @@ namespace Step42
                   fe_values_face.reinit(cell, face);
 
                   right_hand_side.vector_value_list(fe_values_face.get_quadrature_points(),
-                                 right_hand_side_values_face);
+                                                    right_hand_side_values_face);
 
                   for (unsigned int q_point = 0; q_point < n_face_q_points;
                        ++q_point)
@@ -1241,9 +1172,9 @@ namespace Step42
 //    constraints_hanging_nodes.condense(system_rhs_lambda);
 
     const unsigned int sum_elast_points = Utilities::MPI::sum(elast_points,
-                                                        mpi_communicator);
+                                                              mpi_communicator);
     const unsigned int sum_plast_points = Utilities::MPI::sum(plast_points,
-                                                        mpi_communicator);
+                                                              mpi_communicator);
     pcout << "      Number of elastic quadrature points: " << sum_elast_points
           << " and plastic quadrature points: " << sum_plast_points
           << std::endl;
@@ -1328,7 +1259,7 @@ namespace Step42
     active_set.clear();
     IndexSet active_set_locally_owned;
     active_set_locally_owned.set_size(locally_owned_dofs.size());
-    const double c = 100.0 * e_modul;
+    const double c = 100.0 * e_modulus;
 
     Quadrature<dim - 1> face_quadrature(fe.get_unit_face_support_points());
     FEFaceValues<dim> fe_values_face(fe, face_quadrature,

diff --git a/deal.II/include/deal.II/integrators/elasticity.h b/deal.II/include/deal.II/integrators/elasticity.h
index d772cf695cc84853ab2bef8aba7b6cf6136e8ba1..6623ba93903fcd5bf30613060cde023ffde37b2d 100644 (file)
--- a/deal.II/include/deal.II/integrators/elasticity.h
+++ b/deal.II/include/deal.II/integrators/elasticity.h
@@ -214,6 +214,60 @@ namespace LocalIntegrators
         }
     }
 
+    /**
+     * Homogeneous weak boundary condition for the elasticity operator by Nitsche,
+     * namely on the face <i>F</i> the vector
+     * @f[
+     * \int_F \Bigl(\gamma u \cdot v - n^T \epsilon(u) v - u \epsilon(v) n^T\Bigr)\;ds.
+     * @f]
+     *
+     * Here, <i>u</i> is the finite element function whose values and
+     * gradient are given in the arguments <tt>input</tt> and
+     * <tt>Dinput</tt>, respectively. $n$ is the outer
+     * normal vector and $\gamma$ is the usual penalty parameter.
+     *
+     * @author Guido Kanschat
+     * @date 2013
+     */
+    template <int dim, typename number>
+    void nitsche_residual_homogeneous (
+      Vector<number> &result,
+      const FEValuesBase<dim> &fe,
+      const VectorSlice<const std::vector<std::vector<double> > > &input,
+      const VectorSlice<const std::vector<std::vector<Tensor<1,dim> > > > &Dinput,
+      double penalty,
+      double factor = 1.)
+    {
+      const unsigned int n_dofs = fe.dofs_per_cell;
+      AssertVectorVectorDimension(input, dim, fe.n_quadrature_points);
+      AssertVectorVectorDimension(Dinput, dim, fe.n_quadrature_points);
+
+      for (unsigned int k=0; k<fe.n_quadrature_points; ++k)
+        {
+          const double dx = factor * fe.JxW(k);
+          const Point<dim> &n = fe.normal_vector(k);
+          for (unsigned int i=0; i<n_dofs; ++i)
+            for (unsigned int d1=0; d1<dim; ++d1)
+              {
+                const double u= input[d1][k];
+                const double v= fe.shape_value_component(i,k,d1);
+                result(i) += dx + 2.*penalty * u * v;
+
+                for (unsigned int d2=0; d2<dim; ++d2)
+                  {
+                    // v . nabla u n
+                    result(i) -= .5*dx* v * Dinput[d1][k][d2] * n(d2);
+                    // v . (nabla u)^T n
+                    result(i) -= .5*dx* v * Dinput[d2][k][d1] * n(d2);
+                    // u  nabla v n
+                    result(i) -= .5*dx * u * fe.shape_grad_component(i,k,d1)[d2] * n(d2);
+                    // u  (nabla v)^T n
+                    result(i) -= .5*dx * u * fe.shape_grad_component(i,k,d2)[d1] * n(d2);
+                  }
+              }
+        }
+    }
+
     /**
      * The interior penalty flux
      * for symmetric gradients.

diff --git a/deal.II/source/distributed/tria.cc b/deal.II/source/distributed/tria.cc
index 60f62492495516886c0011beec0f5d2fa8019ca2..5b0caa979f73fc30374920e6ad0f3ac390345da8 100644 (file)
--- a/deal.II/source/distributed/tria.cc
+++ b/deal.II/source/distributed/tria.cc
@@ -3511,6 +3511,18 @@ namespace parallel
                        /* user_data_size = */ 0,
                        /* user_data_constructor = */ NULL,
                        /* user_pointer */ this);
+
+
+      try
+      {
+        copy_local_forest_to_triangulation ();
+      }
+      catch (const typename Triangulation<dim>::DistortedCellList &)
+      {
+        // the underlying triangulation should not be checking for distorted
+        // cells
+        AssertThrow (false, ExcInternalError());
+      }
       
       #else
       Assert(false, ExcMessage ("Need p4est version >= 0.3.4.1!"));

diff --git a/deal.II/source/lac/petsc_matrix_base.cc b/deal.II/source/lac/petsc_matrix_base.cc
index fb96b311eb4ce3e7b7a6d3dec09dd3e3a4ff7541..81780da99d92e45529d987efede62f95c9e2e87c 100644 (file)
--- a/deal.II/source/lac/petsc_matrix_base.cc
+++ b/deal.II/source/lac/petsc_matrix_base.cc
@@ -358,13 +358,16 @@ namespace PETScWrappers
     ierr = MatGetRow(*this, row, &ncols, &colnums, &values);
     AssertThrow (ierr == 0, MatrixBase::ExcPETScError(ierr));
 
-    // then restore the matrix and return the
-    // number of columns in this row as
-    // queried previously
+    // then restore the matrix and return the number of columns in this row as
+    // queried previously. Starting with PETSc 3.4, MatRestoreRow actually
+    // resets the last three arguments to zero/NULL, to avoid abuse of pointers
+    // now dangling. as a consequence, we need to save the size of the array
+    // and return the saved value.
+    const PetscInt ncols_saved = ncols;
     ierr = MatRestoreRow(*this, row, &ncols, &colnums, &values);
     AssertThrow (ierr == 0, MatrixBase::ExcPETScError(ierr));
 
-    return ncols;
+    return ncols_saved;
   }
 
 

diff --git a/deal.II/source/meshworker/mesh_worker_vector_selector.inst.in b/deal.II/source/meshworker/mesh_worker_vector_selector.inst.in
index 797b59f8ea60001c1f280c4e8bd29e7b23e9d006..4de6186fc85d229041ef1fe3855661451e11a825 100644 (file)
--- a/deal.II/source/meshworker/mesh_worker_vector_selector.inst.in
+++ b/deal.II/source/meshworker/mesh_worker_vector_selector.inst.in
@@ -29,6 +29,8 @@ for (VECTOR : SERIAL_VECTORS ; deal_II_dimension : DIMENSIONS)
   namespace MeshWorker
   \{
     template class VectorData<VECTOR,deal_II_dimension>;
+    template class VectorData<const VECTOR,deal_II_dimension>;
     template class MGVectorData<VECTOR,deal_II_dimension>;
+    template class MGVectorData<const VECTOR,deal_II_dimension>;
   \}
 }

diff --git a/tests/fe/nedelec_2.cc b/tests/fe/nedelec_2.cc
index 2f9fdfff5e3f1da49b9e1ebcaf888b9855fb914c..bc305f06a588ccb45affac9f3971a1ce1ff4e32f 100644 (file)
--- a/tests/fe/nedelec_2.cc
+++ b/tests/fe/nedelec_2.cc
@@ -47,7 +47,7 @@ plot_diff(const Vector<double> &v1, const Vector<double> &v2, const Tensor<1,2>
   double p1 = v1(0)*n[1] - v1(1)*n[0];
   double p2 = v2(0)*n[1] - v2(1)*n[0];
 
-  deallog << " tangential  diff " << p1-p2 << " (" << p1 << '-' << p2 << ')' << std::endl;
+  deallog << " tangential  diff " << p1-p2 << " (" << p1 << " - " << p2 << ')' << std::endl;
 }
 
 
@@ -67,7 +67,7 @@ plot_diff(const Vector<double> &v1, const Vector<double> &v2, const Tensor<1,3>
       p2[d] = v2(d1)*n[d2] - v2(d2)*n[d1];
     }
 
-  deallog << " tangential  diff " << p1-p2 << " (" << p1 << '-' << p2 << ')' << std::endl;
+  deallog << " tangential  diff " << p1-p2 << " (" << p1 << " - " << p2 << ')' << std::endl;
 }
 
 

diff --git a/tests/fe/nedelec_2.output b/tests/fe/nedelec_2.output
index b89319def1aa3fdae7e4c6fb5ae8572e208ae435..bfdce69d4629774e9b4f735ee85e6a8f86e81b8d 100644 (file)
--- a/tests/fe/nedelec_2.output
+++ b/tests/fe/nedelec_2.output
@@ -5,735 +5,735 @@ DEAL::cell 1.1 face 0 neighbor 1.0 refined
 DEAL::cell 1.1 face 1 boundary
 DEAL::cell 1.1 face 2 boundary
 DEAL::cell 1.1 face 3 neighbor 1.3
-DEAL::0.50 0.50 tangential  diff 0 (6.00-6.00)
-DEAL::1.00 0.50 tangential  diff 0 (6.00-6.00)
+DEAL::0.50 0.50 tangential  diff 0 (6.00 - 6.00)
+DEAL::1.00 0.50 tangential  diff 0 (6.00 - 6.00)
 DEAL::cell 1.2 face 0 boundary
 DEAL::cell 1.2 face 1 neighbor 1.3
-DEAL::0.50 0.50 tangential  diff 0 (-10.00--10.00)
-DEAL::0.50 1.00 tangential  diff 0 (-10.00--10.00)
+DEAL::0.50 0.50 tangential  diff 0 (-10.00 - -10.00)
+DEAL::0.50 1.00 tangential  diff 0 (-10.00 - -10.00)
 DEAL::cell 1.2 face 2 neighbor 1.0 refined
 DEAL::cell 1.2 face 3 boundary
 DEAL::cell 1.3 face 0 neighbor 1.2
-DEAL::0.50 0.50 tangential  diff 0 (10.00-10.00)
-DEAL::0.50 1.00 tangential  diff 0 (10.00-10.00)
+DEAL::0.50 0.50 tangential  diff 0 (10.00 - 10.00)
+DEAL::0.50 1.00 tangential  diff 0 (10.00 - 10.00)
 DEAL::cell 1.3 face 1 boundary
 DEAL::cell 1.3 face 2 neighbor 1.1
-DEAL::0.50 0.50 tangential  diff 0 (-6.00--6.00)
-DEAL::1.00 0.50 tangential  diff 0 (-6.00--6.00)
+DEAL::0.50 0.50 tangential  diff 0 (-6.00 - -6.00)
+DEAL::1.00 0.50 tangential  diff 0 (-6.00 - -6.00)
 DEAL::cell 1.3 face 3 boundary
 DEAL::cell 2.0 face 0 boundary
 DEAL::cell 2.0 face 1 neighbor 2.1
-DEAL::0.25 0.00 tangential  diff 0 (-44.00--44.00)
-DEAL::0.25 0.25 tangential  diff 0 (-44.00--44.00)
+DEAL::0.25 0.00 tangential  diff 0 (-44.00 - -44.00)
+DEAL::0.25 0.25 tangential  diff 0 (-44.00 - -44.00)
 DEAL::cell 2.0 face 2 boundary
 DEAL::cell 2.0 face 3 neighbor 2.2
-DEAL::0.00 0.25 tangential  diff 0 (52.00-52.00)
-DEAL::0.25 0.25 tangential  diff 0 (52.00-52.00)
+DEAL::0.00 0.25 tangential  diff 0 (52.00 - 52.00)
+DEAL::0.25 0.25 tangential  diff 0 (52.00 - 52.00)
 DEAL::cell 2.1 face 0 neighbor 2.0
-DEAL::0.25 0.00 tangential  diff 0 (44.00-44.00)
-DEAL::0.25 0.25 tangential  diff 0 (44.00-44.00)
+DEAL::0.25 0.00 tangential  diff 0 (44.00 - 44.00)
+DEAL::0.25 0.25 tangential  diff 0 (44.00 - 44.00)
 DEAL::cell 2.1 face 1 neighbor 1.1
-DEAL::0.50 0.00 tangential  diff 0 (0-0)
-DEAL::0.50 0.25 tangential  diff 0 (0-0)
+DEAL::0.50 0.00 tangential  diff 0 (0 - 0)
+DEAL::0.50 0.25 tangential  diff 0 (0 - 0)
 DEAL::cell 2.1 face 2 boundary
 DEAL::cell 2.1 face 3 neighbor 2.3
-DEAL::0.25 0.25 tangential  diff 0 (64.00-64.00)
-DEAL::0.50 0.25 tangential  diff 0 (64.00-64.00)
+DEAL::0.25 0.25 tangential  diff 0 (64.00 - 64.00)
+DEAL::0.50 0.25 tangential  diff 0 (64.00 - 64.00)
 DEAL::cell 2.2 face 0 boundary
 DEAL::cell 2.2 face 1 neighbor 2.3
-DEAL::0.25 0.25 tangential  diff 0 (-72.00--72.00)
-DEAL::0.25 0.50 tangential  diff 0 (-72.00--72.00)
+DEAL::0.25 0.25 tangential  diff 0 (-72.00 - -72.00)
+DEAL::0.25 0.50 tangential  diff 0 (-72.00 - -72.00)
 DEAL::cell 2.2 face 2 neighbor 2.0
-DEAL::0.00 0.25 tangential  diff 0 (-52.00--52.00)
-DEAL::0.25 0.25 tangential  diff 0 (-52.00--52.00)
+DEAL::0.00 0.25 tangential  diff 0 (-52.00 - -52.00)
+DEAL::0.25 0.25 tangential  diff 0 (-52.00 - -52.00)
 DEAL::cell 2.2 face 3 neighbor 1.2
-DEAL::0.00 0.50 tangential  diff 0 (12.00-12.00)
-DEAL::0.25 0.50 tangential  diff 0 (12.00-12.00)
+DEAL::0.00 0.50 tangential  diff 0 (12.00 - 12.00)
+DEAL::0.25 0.50 tangential  diff 0 (12.00 - 12.00)
 DEAL::cell 2.3 face 0 neighbor 2.2
-DEAL::0.25 0.25 tangential  diff 0 (72.00-72.00)
-DEAL::0.25 0.50 tangential  diff 0 (72.00-72.00)
+DEAL::0.25 0.25 tangential  diff 0 (72.00 - 72.00)
+DEAL::0.25 0.50 tangential  diff 0 (72.00 - 72.00)
 DEAL::cell 2.3 face 1 neighbor 1.1
-DEAL::0.50 0.25 tangential  diff 0 (0-0)
-DEAL::0.50 0.50 tangential  diff 0 (0-0)
+DEAL::0.50 0.25 tangential  diff 0 (0 - 0)
+DEAL::0.50 0.50 tangential  diff 0 (0 - 0)
 DEAL::cell 2.3 face 2 neighbor 2.1
-DEAL::0.25 0.25 tangential  diff 0 (-64.00--64.00)
-DEAL::0.50 0.25 tangential  diff 0 (-64.00--64.00)
+DEAL::0.25 0.25 tangential  diff 0 (-64.00 - -64.00)
+DEAL::0.50 0.25 tangential  diff 0 (-64.00 - -64.00)
 DEAL::cell 2.3 face 3 neighbor 1.2
-DEAL::0.25 0.50 tangential  diff 0 (12.00-12.00)
-DEAL::0.50 0.50 tangential  diff 0 (12.00-12.00)
+DEAL::0.25 0.50 tangential  diff 0 (12.00 - 12.00)
+DEAL::0.50 0.50 tangential  diff 0 (12.00 - 12.00)
 DEAL::3d, 15 CELLS
 DEAL::cell 1.1 face 0 neighbor 1.0 refined
 DEAL::cell 1.1 face 1 boundary
 DEAL::cell 1.1 face 2 boundary
 DEAL::cell 1.1 face 3 neighbor 1.3
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 6.00--20.00 0.00 6.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 14.00--20.00 0.00 14.00)
-DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-22.00 0.00 6.00--22.00 0.00 6.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-22.00 0.00 14.00--22.00 0.00 14.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 6.00 - -20.00 0.00 6.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 14.00 - -20.00 0.00 14.00)
+DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-22.00 0.00 6.00 - -22.00 0.00 6.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-22.00 0.00 14.00 - -22.00 0.00 14.00)
 DEAL::cell 1.1 face 4 boundary
 DEAL::cell 1.1 face 5 neighbor 1.5
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -12.00 0.00-8.00 -12.00 0.00)
-DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (10.00 -12.00 0.00-10.00 -12.00 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -14.00 0.00-8.00 -14.00 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (10.00 -14.00 0.00-10.00 -14.00 0.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -12.00 0.00 - 8.00 -12.00 0.00)
+DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (10.00 -12.00 0.00 - 10.00 -12.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -14.00 0.00 - 8.00 -14.00 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (10.00 -14.00 0.00 - 10.00 -14.00 0.00)
 DEAL::cell 1.2 face 0 boundary
 DEAL::cell 1.2 face 1 neighbor 1.3
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -26.000.00 20.00 -26.00)
-DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 44.00 -26.000.00 44.00 -26.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -34.000.00 20.00 -34.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 44.00 -34.000.00 44.00 -34.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -26.00 - 0.00 20.00 -26.00)
+DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 44.00 -26.00 - 0.00 44.00 -26.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -34.00 - 0.00 20.00 -34.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 44.00 -34.00 - 0.00 44.00 -34.00)
 DEAL::cell 1.2 face 2 neighbor 1.0 refined
 DEAL::cell 1.2 face 3 boundary
 DEAL::cell 1.2 face 4 boundary
 DEAL::cell 1.2 face 5 neighbor 1.6
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -36.00 0.00-32.00 -36.00 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -36.00 0.00-34.00 -36.00 0.00)
-DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -38.00 0.00-32.00 -38.00 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -38.00 0.00-34.00 -38.00 0.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -36.00 0.00 - 32.00 -36.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -36.00 0.00 - 34.00 -36.00 0.00)
+DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -38.00 0.00 - 32.00 -38.00 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -38.00 0.00 - 34.00 -38.00 0.00)
 DEAL::cell 1.3 face 0 neighbor 1.2
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -20.00 26.000.00 -20.00 26.00)
-DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -44.00 26.000.00 -44.00 26.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -20.00 34.000.00 -20.00 34.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -44.00 34.000.00 -44.00 34.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -20.00 26.00 - 0.00 -20.00 26.00)
+DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -44.00 26.00 - 0.00 -44.00 26.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -20.00 34.00 - 0.00 -20.00 34.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -44.00 34.00 - 0.00 -44.00 34.00)
 DEAL::cell 1.3 face 1 boundary
 DEAL::cell 1.3 face 2 neighbor 1.1
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (20.00 0.00 -6.00-20.00 0.00 -6.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (20.00 0.00 -14.00-20.00 0.00 -14.00)
-DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (22.00 0.00 -6.00-22.00 0.00 -6.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (22.00 0.00 -14.00-22.00 0.00 -14.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (20.00 0.00 -6.00 - 20.00 0.00 -6.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (20.00 0.00 -14.00 - 20.00 0.00 -14.00)
+DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (22.00 0.00 -6.00 - 22.00 0.00 -6.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (22.00 0.00 -14.00 - 22.00 0.00 -14.00)
 DEAL::cell 1.3 face 3 boundary
 DEAL::cell 1.3 face 4 boundary
 DEAL::cell 1.3 face 5 neighbor 1.7
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -14.00 0.00-34.00 -14.00 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (50.00 -14.00 0.00-50.00 -14.00 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -52.00 0.00-34.00 -52.00 0.00)
-DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (50.00 -52.00 0.00-50.00 -52.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -14.00 0.00 - 34.00 -14.00 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (50.00 -14.00 0.00 - 50.00 -14.00 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (34.00 -52.00 0.00 - 34.00 -52.00 0.00)
+DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (50.00 -52.00 0.00 - 50.00 -52.00 0.00)
 DEAL::cell 1.4 face 0 boundary
 DEAL::cell 1.4 face 1 neighbor 1.5
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 70.00 -8.000.00 70.00 -8.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -8.000.00 74.00 -8.00)
-DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 70.00 -62.000.00 70.00 -62.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -62.000.00 74.00 -62.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 70.00 -8.00 - 0.00 70.00 -8.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -8.00 - 0.00 74.00 -8.00)
+DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 70.00 -62.00 - 0.00 70.00 -62.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -62.00 - 0.00 74.00 -62.00)
 DEAL::cell 1.4 face 2 boundary
 DEAL::cell 1.4 face 3 neighbor 1.6
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-72.00 0.00 36.00--72.00 0.00 36.00)
-DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-72.00 0.00 66.00--72.00 0.00 66.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 36.00--74.00 0.00 36.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 66.00--74.00 0.00 66.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-72.00 0.00 36.00 - -72.00 0.00 36.00)
+DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-72.00 0.00 66.00 - -72.00 0.00 66.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 36.00 - -74.00 0.00 36.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 66.00 - -74.00 0.00 66.00)
 DEAL::cell 1.4 face 4 neighbor 1.0 refined
 DEAL::cell 1.4 face 5 boundary
 DEAL::cell 1.5 face 0 neighbor 1.4
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -70.00 8.000.00 -70.00 8.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 8.000.00 -74.00 8.00)
-DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -70.00 62.000.00 -70.00 62.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 62.000.00 -74.00 62.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -70.00 8.00 - 0.00 -70.00 8.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 8.00 - 0.00 -74.00 8.00)
+DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -70.00 62.00 - 0.00 -70.00 62.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 62.00 - 0.00 -74.00 62.00)
 DEAL::cell 1.5 face 1 boundary
 DEAL::cell 1.5 face 2 boundary
 DEAL::cell 1.5 face 3 neighbor 1.7
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 14.00--74.00 0.00 14.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 80.00--74.00 0.00 80.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-84.00 0.00 14.00--84.00 0.00 14.00)
-DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-84.00 0.00 80.00--84.00 0.00 80.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 14.00 - -74.00 0.00 14.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-74.00 0.00 80.00 - -74.00 0.00 80.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-84.00 0.00 14.00 - -84.00 0.00 14.00)
+DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-84.00 0.00 80.00 - -84.00 0.00 80.00)
 DEAL::cell 1.5 face 4 neighbor 1.1
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-8.00 12.00 0.00--8.00 12.00 0.00)
-DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-10.00 12.00 0.00--10.00 12.00 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-8.00 14.00 0.00--8.00 14.00 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-10.00 14.00 0.00--10.00 14.00 0.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-8.00 12.00 0.00 - -8.00 12.00 0.00)
+DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-10.00 12.00 0.00 - -10.00 12.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-8.00 14.00 0.00 - -8.00 14.00 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-10.00 14.00 0.00 - -10.00 14.00 0.00)
 DEAL::cell 1.5 face 5 boundary
 DEAL::cell 1.6 face 0 boundary
 DEAL::cell 1.6 face 1 neighbor 1.7
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -34.000.00 74.00 -34.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 94.00 -34.000.00 94.00 -34.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -88.000.00 74.00 -88.00)
-DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 94.00 -88.000.00 94.00 -88.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -34.00 - 0.00 74.00 -34.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 94.00 -34.00 - 0.00 94.00 -34.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 74.00 -88.00 - 0.00 74.00 -88.00)
+DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 94.00 -88.00 - 0.00 94.00 -88.00)
 DEAL::cell 1.6 face 2 neighbor 1.4
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (72.00 0.00 -36.00-72.00 0.00 -36.00)
-DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (72.00 0.00 -66.00-72.00 0.00 -66.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -36.00-74.00 0.00 -36.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -66.00-74.00 0.00 -66.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (72.00 0.00 -36.00 - 72.00 0.00 -36.00)
+DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (72.00 0.00 -66.00 - 72.00 0.00 -66.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -36.00 - 74.00 0.00 -36.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -66.00 - 74.00 0.00 -66.00)
 DEAL::cell 1.6 face 3 boundary
 DEAL::cell 1.6 face 4 neighbor 1.2
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-32.00 36.00 0.00--32.00 36.00 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 36.00 0.00--34.00 36.00 0.00)
-DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-32.00 38.00 0.00--32.00 38.00 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 38.00 0.00--34.00 38.00 0.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-32.00 36.00 0.00 - -32.00 36.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 36.00 0.00 - -34.00 36.00 0.00)
+DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-32.00 38.00 0.00 - -32.00 38.00 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 38.00 0.00 - -34.00 38.00 0.00)
 DEAL::cell 1.6 face 5 boundary
 DEAL::cell 1.7 face 0 neighbor 1.6
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 34.000.00 -74.00 34.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -94.00 34.000.00 -94.00 34.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 88.000.00 -74.00 88.00)
-DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -94.00 88.000.00 -94.00 88.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 34.00 - 0.00 -74.00 34.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -94.00 34.00 - 0.00 -94.00 34.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -74.00 88.00 - 0.00 -74.00 88.00)
+DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -94.00 88.00 - 0.00 -94.00 88.00)
 DEAL::cell 1.7 face 1 boundary
 DEAL::cell 1.7 face 2 neighbor 1.5
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -14.00-74.00 0.00 -14.00)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -80.00-74.00 0.00 -80.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (84.00 0.00 -14.00-84.00 0.00 -14.00)
-DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (84.00 0.00 -80.00-84.00 0.00 -80.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -14.00 - 74.00 0.00 -14.00)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (74.00 0.00 -80.00 - 74.00 0.00 -80.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (84.00 0.00 -14.00 - 84.00 0.00 -14.00)
+DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (84.00 0.00 -80.00 - 84.00 0.00 -80.00)
 DEAL::cell 1.7 face 3 boundary
 DEAL::cell 1.7 face 4 neighbor 1.3
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 14.00 0.00--34.00 14.00 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-50.00 14.00 0.00--50.00 14.00 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 52.00 0.00--34.00 52.00 0.00)
-DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-50.00 52.00 0.00--50.00 52.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 14.00 0.00 - -34.00 14.00 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-50.00 14.00 0.00 - -50.00 14.00 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-34.00 52.00 0.00 - -34.00 52.00 0.00)
+DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-50.00 52.00 0.00 - -50.00 52.00 0.00)
 DEAL::cell 1.7 face 5 boundary
 DEAL::cell 2.0 face 0 boundary
 DEAL::cell 2.0 face 1 neighbor 2.1
-DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 240.00 -208.000.00 240.00 -208.00)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -208.000.00 248.00 -208.00)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 240.00 -224.000.00 240.00 -224.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -224.000.00 248.00 -224.00)
+DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 240.00 -208.00 - 0.00 240.00 -208.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -208.00 - 0.00 248.00 -208.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 240.00 -224.00 - 0.00 240.00 -224.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -224.00 - 0.00 248.00 -224.00)
 DEAL::cell 2.0 face 2 boundary
 DEAL::cell 2.0 face 3 neighbor 2.2
-DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-244.00 0.00 216.00--244.00 0.00 216.00)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-244.00 0.00 232.00--244.00 0.00 232.00)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 216.00--248.00 0.00 216.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 232.00--248.00 0.00 232.00)
+DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-244.00 0.00 216.00 - -244.00 0.00 216.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-244.00 0.00 232.00 - -244.00 0.00 232.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 216.00 - -248.00 0.00 216.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 232.00 - -248.00 0.00 232.00)
 DEAL::cell 2.0 face 4 boundary
 DEAL::cell 2.0 face 5 neighbor 2.4
-DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (220.00 -228.00 0.00-220.00 -228.00 0.00)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -228.00 0.00-224.00 -228.00 0.00)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (220.00 -232.00 0.00-220.00 -232.00 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -232.00 0.00-224.00 -232.00 0.00)
+DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (220.00 -228.00 0.00 - 220.00 -228.00 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -228.00 0.00 - 224.00 -228.00 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (220.00 -232.00 0.00 - 220.00 -232.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -232.00 0.00 - 224.00 -232.00 0.00)
 DEAL::cell 2.1 face 0 neighbor 2.0
-DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -240.00 208.000.00 -240.00 208.00)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 208.000.00 -248.00 208.00)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -240.00 224.000.00 -240.00 224.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 224.000.00 -248.00 224.00)
+DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -240.00 208.00 - 0.00 -240.00 208.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 208.00 - 0.00 -248.00 208.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -240.00 224.00 - 0.00 -240.00 224.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 224.00 - 0.00 -248.00 224.00)
 DEAL::cell 2.1 face 1 neighbor 1.1
-DEAL::0.50 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 16.00 0.000.00 16.00 0.00)
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 18.00 0.000.00 18.00 0.00)
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 16.00 -4.000.00 16.00 -4.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.000.00 18.00 -4.00)
+DEAL::0.50 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 16.00 0.00 - 0.00 16.00 0.00)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 18.00 0.00 - 0.00 18.00 0.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 16.00 -4.00 - 0.00 16.00 -4.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.00 - 0.00 18.00 -4.00)
 DEAL::cell 2.1 face 2 boundary
 DEAL::cell 2.1 face 3 neighbor 2.3
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 260.00--248.00 0.00 260.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 272.00--248.00 0.00 272.00)
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 260.00--18.00 0.00 260.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 272.00--18.00 0.00 272.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 260.00 - -248.00 0.00 260.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-248.00 0.00 272.00 - -248.00 0.00 272.00)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 260.00 - -18.00 0.00 260.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 272.00 - -18.00 0.00 272.00)
 DEAL::cell 2.1 face 4 boundary
 DEAL::cell 2.1 face 5 neighbor 2.5
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -268.00 0.00-224.00 -268.00 0.00)
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -268.00 0.00-4.00 -268.00 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -272.00 0.00-224.00 -272.00 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -272.00 0.00-4.00 -272.00 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -268.00 0.00 - 224.00 -268.00 0.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -268.00 0.00 - 4.00 -268.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (224.00 -272.00 0.00 - 224.00 -272.00 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -272.00 0.00 - 4.00 -272.00 0.00)
 DEAL::cell 2.2 face 0 boundary
 DEAL::cell 2.2 face 1 neighbor 2.3
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -288.000.00 248.00 -288.00)
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -288.000.00 30.00 -288.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -300.000.00 248.00 -300.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -300.000.00 30.00 -300.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -288.00 - 0.00 248.00 -288.00)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -288.00 - 0.00 30.00 -288.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 248.00 -300.00 - 0.00 248.00 -300.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -300.00 - 0.00 30.00 -300.00)
 DEAL::cell 2.2 face 2 neighbor 2.0
-DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (244.00 0.00 -216.00-244.00 0.00 -216.00)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (244.00 0.00 -232.00-244.00 0.00 -232.00)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -216.00-248.00 0.00 -216.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -232.00-248.00 0.00 -232.00)
+DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (244.00 0.00 -216.00 - 244.00 0.00 -216.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (244.00 0.00 -232.00 - 244.00 0.00 -232.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -216.00 - 248.00 0.00 -216.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -232.00 - 248.00 0.00 -232.00)
 DEAL::cell 2.2 face 3 neighbor 1.2
-DEAL::0.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 28.00--40.00 0.00 28.00)
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 32.00--40.00 0.00 32.00)
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 28.00--30.00 0.00 28.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00--30.00 0.00 32.00)
+DEAL::0.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 28.00 - -40.00 0.00 28.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 32.00 - -40.00 0.00 32.00)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 28.00 - -30.00 0.00 28.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00 - -30.00 0.00 32.00)
 DEAL::cell 2.2 face 4 boundary
 DEAL::cell 2.2 face 5 neighbor 2.6
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (296.00 -232.00 0.00-296.00 -232.00 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -232.00 0.00-300.00 -232.00 0.00)
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (296.00 -32.00 0.00-296.00 -32.00 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -32.00 0.00-300.00 -32.00 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (296.00 -232.00 0.00 - 296.00 -232.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -232.00 0.00 - 300.00 -232.00 0.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (296.00 -32.00 0.00 - 296.00 -32.00 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -32.00 0.00 - 300.00 -32.00 0.00)
 DEAL::cell 2.3 face 0 neighbor 2.2
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 288.000.00 -248.00 288.00)
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 288.000.00 -30.00 288.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 300.000.00 -248.00 300.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 300.000.00 -30.00 300.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 288.00 - 0.00 -248.00 288.00)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 288.00 - 0.00 -30.00 288.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -248.00 300.00 - 0.00 -248.00 300.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 300.00 - 0.00 -30.00 300.00)
 DEAL::cell 2.3 face 1 neighbor 1.1
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 18.00 0.000.00 18.00 0.00)
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 20.00 0.000.00 20.00 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.000.00 18.00 -4.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -4.000.00 20.00 -4.00)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 18.00 0.00 - 0.00 18.00 0.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 20.00 0.00 - 0.00 20.00 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.00 - 0.00 18.00 -4.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -4.00 - 0.00 20.00 -4.00)
 DEAL::cell 2.3 face 2 neighbor 2.1
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -260.00-248.00 0.00 -260.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -272.00-248.00 0.00 -272.00)
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -260.00-18.00 0.00 -260.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -272.00-18.00 0.00 -272.00)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -260.00 - 248.00 0.00 -260.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (248.00 0.00 -272.00 - 248.00 0.00 -272.00)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -260.00 - 18.00 0.00 -260.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -272.00 - 18.00 0.00 -272.00)
 DEAL::cell 2.3 face 3 neighbor 1.2
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 28.00--30.00 0.00 28.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00--30.00 0.00 32.00)
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 28.00--20.00 0.00 28.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 32.00--20.00 0.00 32.00)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 28.00 - -30.00 0.00 28.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00 - -30.00 0.00 32.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 28.00 - -20.00 0.00 28.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 32.00 - -20.00 0.00 32.00)
 DEAL::cell 2.3 face 4 boundary
 DEAL::cell 2.3 face 5 neighbor 2.7
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -272.00 0.00-300.00 -272.00 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -272.00 0.00-4.00 -272.00 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -32.00 0.00-300.00 -32.00 0.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -32.00 0.00-4.00 -32.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -272.00 0.00 - 300.00 -272.00 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -272.00 0.00 - 4.00 -272.00 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (300.00 -32.00 0.00 - 300.00 -32.00 0.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4.00 -32.00 0.00 - 4.00 -32.00 0.00)
 DEAL::cell 2.4 face 0 boundary
 DEAL::cell 2.4 face 1 neighbor 2.5
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 356.00 -224.000.00 356.00 -224.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -224.000.00 364.00 -224.00)
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 356.00 -32.000.00 356.00 -32.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -32.000.00 364.00 -32.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 356.00 -224.00 - 0.00 356.00 -224.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -224.00 - 0.00 364.00 -224.00)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 356.00 -32.00 - 0.00 356.00 -32.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -32.00 - 0.00 364.00 -32.00)
 DEAL::cell 2.4 face 2 boundary
 DEAL::cell 2.4 face 3 neighbor 2.6
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-360.00 0.00 232.00--360.00 0.00 232.00)
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-360.00 0.00 47.00--360.00 0.00 47.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 232.00--364.00 0.00 232.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 47.00--364.00 0.00 47.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-360.00 0.00 232.00 - -360.00 0.00 232.00)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-360.00 0.00 47.00 - -360.00 0.00 47.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 232.00 - -364.00 0.00 232.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 47.00 - -364.00 0.00 47.00)
 DEAL::cell 2.4 face 4 neighbor 2.0
-DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-220.00 228.00 0.00--220.00 228.00 0.00)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 228.00 0.00--224.00 228.00 0.00)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-220.00 232.00 0.00--220.00 232.00 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 232.00 0.00--224.00 232.00 0.00)
+DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-220.00 228.00 0.00 - -220.00 228.00 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 228.00 0.00 - -224.00 228.00 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-220.00 232.00 0.00 - -220.00 232.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 232.00 0.00 - -224.00 232.00 0.00)
 DEAL::cell 2.4 face 5 neighbor 1.4
-DEAL::0.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -58.00 0.00-56.00 -58.00 0.00)
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -58.00 0.00-32.00 -58.00 0.00)
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -47.00 0.00-56.00 -47.00 0.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00-32.00 -47.00 0.00)
+DEAL::0.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -58.00 0.00 - 56.00 -58.00 0.00)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -58.00 0.00 - 32.00 -58.00 0.00)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -47.00 0.00 - 56.00 -47.00 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00 - 32.00 -47.00 0.00)
 DEAL::cell 2.5 face 0 neighbor 2.4
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -356.00 224.000.00 -356.00 224.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 224.000.00 -364.00 224.00)
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -356.00 32.000.00 -356.00 32.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 32.000.00 -364.00 32.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -356.00 224.00 - 0.00 -356.00 224.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 224.00 - 0.00 -364.00 224.00)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -356.00 32.00 - 0.00 -356.00 32.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 32.00 - 0.00 -364.00 32.00)
 DEAL::cell 2.5 face 1 neighbor 1.1
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 16.00 -4.000.00 16.00 -4.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.000.00 18.00 -4.00)
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 16.00 -8.000.00 16.00 -8.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -8.000.00 18.00 -8.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 16.00 -4.00 - 0.00 16.00 -4.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.00 - 0.00 18.00 -4.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 16.00 -8.00 - 0.00 16.00 -8.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -8.00 - 0.00 18.00 -8.00)
 DEAL::cell 2.5 face 2 boundary
 DEAL::cell 2.5 face 3 neighbor 2.7
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 272.00--364.00 0.00 272.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 47.00--364.00 0.00 47.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 272.00--18.00 0.00 272.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 47.00--18.00 0.00 47.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 272.00 - -364.00 0.00 272.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-364.00 0.00 47.00 - -364.00 0.00 47.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 272.00 - -18.00 0.00 272.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-18.00 0.00 47.00 - -18.00 0.00 47.00)
 DEAL::cell 2.5 face 4 neighbor 2.1
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 268.00 0.00--224.00 268.00 0.00)
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 268.00 0.00--4.00 268.00 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 272.00 0.00--224.00 272.00 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 272.00 0.00--4.00 272.00 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 268.00 0.00 - -224.00 268.00 0.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 268.00 0.00 - -4.00 268.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-224.00 272.00 0.00 - -224.00 272.00 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 272.00 0.00 - -4.00 272.00 0.00)
 DEAL::cell 2.5 face 5 neighbor 1.4
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -58.00 0.00-32.00 -58.00 0.00)
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -58.00 0.00-8.00 -58.00 0.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00-32.00 -47.00 0.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -47.00 0.00-8.00 -47.00 0.00)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -58.00 0.00 - 32.00 -58.00 0.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -58.00 0.00 - 8.00 -58.00 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00 - 32.00 -47.00 0.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -47.00 0.00 - 8.00 -47.00 0.00)
 DEAL::cell 2.6 face 0 boundary
 DEAL::cell 2.6 face 1 neighbor 2.7
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -300.000.00 364.00 -300.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -300.000.00 30.00 -300.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -32.000.00 364.00 -32.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -32.000.00 30.00 -32.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -300.00 - 0.00 364.00 -300.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -300.00 - 0.00 30.00 -300.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 364.00 -32.00 - 0.00 364.00 -32.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 30.00 -32.00 - 0.00 30.00 -32.00)
 DEAL::cell 2.6 face 2 neighbor 2.4
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (360.00 0.00 -232.00-360.00 0.00 -232.00)
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (360.00 0.00 -47.00-360.00 0.00 -47.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -232.00-364.00 0.00 -232.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -47.00-364.00 0.00 -47.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (360.00 0.00 -232.00 - 360.00 0.00 -232.00)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (360.00 0.00 -47.00 - 360.00 0.00 -47.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -232.00 - 364.00 0.00 -232.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -47.00 - 364.00 0.00 -47.00)
 DEAL::cell 2.6 face 3 neighbor 1.2
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 32.00--40.00 0.00 32.00)
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 36.00--40.00 0.00 36.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00--30.00 0.00 32.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 36.00--30.00 0.00 36.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 32.00 - -40.00 0.00 32.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 36.00 - -40.00 0.00 36.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00 - -30.00 0.00 32.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 36.00 - -30.00 0.00 36.00)
 DEAL::cell 2.6 face 4 neighbor 2.2
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-296.00 232.00 0.00--296.00 232.00 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 232.00 0.00--300.00 232.00 0.00)
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-296.00 32.00 0.00--296.00 32.00 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 32.00 0.00--300.00 32.00 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-296.00 232.00 0.00 - -296.00 232.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 232.00 0.00 - -300.00 232.00 0.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-296.00 32.00 0.00 - -296.00 32.00 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 32.00 0.00 - -300.00 32.00 0.00)
 DEAL::cell 2.6 face 5 neighbor 1.4
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -47.00 0.00-56.00 -47.00 0.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00-32.00 -47.00 0.00)
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -36.00 0.00-56.00 -36.00 0.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -36.00 0.00-32.00 -36.00 0.00)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -47.00 0.00 - 56.00 -47.00 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00 - 32.00 -47.00 0.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (56.00 -36.00 0.00 - 56.00 -36.00 0.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -36.00 0.00 - 32.00 -36.00 0.00)
 DEAL::cell 2.7 face 0 neighbor 2.6
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 300.000.00 -364.00 300.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 300.000.00 -30.00 300.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 32.000.00 -364.00 32.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 32.000.00 -30.00 32.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 300.00 - 0.00 -364.00 300.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 300.00 - 0.00 -30.00 300.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -364.00 32.00 - 0.00 -364.00 32.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -30.00 32.00 - 0.00 -30.00 32.00)
 DEAL::cell 2.7 face 1 neighbor 1.1
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.000.00 18.00 -4.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -4.000.00 20.00 -4.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -8.000.00 18.00 -8.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -8.000.00 20.00 -8.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -4.00 - 0.00 18.00 -4.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -4.00 - 0.00 20.00 -4.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 18.00 -8.00 - 0.00 18.00 -8.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 20.00 -8.00 - 0.00 20.00 -8.00)
 DEAL::cell 2.7 face 2 neighbor 2.5
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -272.00-364.00 0.00 -272.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -47.00-364.00 0.00 -47.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -272.00-18.00 0.00 -272.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -47.00-18.00 0.00 -47.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -272.00 - 364.00 0.00 -272.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (364.00 0.00 -47.00 - 364.00 0.00 -47.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -272.00 - 18.00 0.00 -272.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (18.00 0.00 -47.00 - 18.00 0.00 -47.00)
 DEAL::cell 2.7 face 3 neighbor 1.2
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00--30.00 0.00 32.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 36.00--30.00 0.00 36.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 32.00--20.00 0.00 32.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 36.00--20.00 0.00 36.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 32.00 - -30.00 0.00 32.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-30.00 0.00 36.00 - -30.00 0.00 36.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 32.00 - -20.00 0.00 32.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-20.00 0.00 36.00 - -20.00 0.00 36.00)
 DEAL::cell 2.7 face 4 neighbor 2.3
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 272.00 0.00--300.00 272.00 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 272.00 0.00--4.00 272.00 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 32.00 0.00--300.00 32.00 0.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 32.00 0.00--4.00 32.00 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 272.00 0.00 - -300.00 272.00 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 272.00 0.00 - -4.00 272.00 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-300.00 32.00 0.00 - -300.00 32.00 0.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4.00 32.00 0.00 - -4.00 32.00 0.00)
 DEAL::cell 2.7 face 5 neighbor 1.4
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00-32.00 -47.00 0.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -47.00 0.00-8.00 -47.00 0.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -36.00 0.00-32.00 -36.00 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -36.00 0.00-8.00 -36.00 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -47.00 0.00 - 32.00 -47.00 0.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -47.00 0.00 - 8.00 -47.00 0.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (32.00 -36.00 0.00 - 32.00 -36.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (8.00 -36.00 0.00 - 8.00 -36.00 0.00)
 DEAL::Degree 1:
 DEAL::2d, 7 CELLS
 DEAL::cell 1.1 face 0 neighbor 1.0 refined
 DEAL::cell 1.1 face 1 boundary
 DEAL::cell 1.1 face 2 boundary
 DEAL::cell 1.1 face 3 neighbor 1.3
-DEAL::0.50 0.50 tangential  diff 0 (-12.25--12.25)
-DEAL::1.00 0.50 tangential  diff 0 (36.25-36.25)
+DEAL::0.50 0.50 tangential  diff 0 (-12.25 - -12.25)
+DEAL::1.00 0.50 tangential  diff 0 (36.25 - 36.25)
 DEAL::cell 1.2 face 0 boundary
 DEAL::cell 1.2 face 1 neighbor 1.3
-DEAL::0.50 0.50 tangential  diff 0 (23.96-23.96)
-DEAL::0.50 1.00 tangential  diff 0 (-79.96--79.96)
+DEAL::0.50 0.50 tangential  diff 0 (23.96 - 23.96)
+DEAL::0.50 1.00 tangential  diff 0 (-79.96 - -79.96)
 DEAL::cell 1.2 face 2 neighbor 1.0 refined
 DEAL::cell 1.2 face 3 boundary
 DEAL::cell 1.3 face 0 neighbor 1.2
-DEAL::0.50 0.50 tangential  diff 0 (-23.96--23.96)
-DEAL::0.50 1.00 tangential  diff 0 (79.96-79.96)
+DEAL::0.50 0.50 tangential  diff 0 (-23.96 - -23.96)
+DEAL::0.50 1.00 tangential  diff 0 (79.96 - 79.96)
 DEAL::cell 1.3 face 1 boundary
 DEAL::cell 1.3 face 2 neighbor 1.1
-DEAL::0.50 0.50 tangential  diff 0 (12.25-12.25)
-DEAL::1.00 0.50 tangential  diff 0 (-36.25--36.25)
+DEAL::0.50 0.50 tangential  diff 0 (12.25 - 12.25)
+DEAL::1.00 0.50 tangential  diff 0 (-36.25 - -36.25)
 DEAL::cell 1.3 face 3 boundary
 DEAL::cell 2.0 face 0 boundary
 DEAL::cell 2.0 face 1 neighbor 2.1
-DEAL::0.25 0.00 tangential  diff 0 (106.49-106.49)
-DEAL::0.25 0.25 tangential  diff 0 (-378.49--378.49)
+DEAL::0.25 0.00 tangential  diff 0 (106.49 - 106.49)
+DEAL::0.25 0.25 tangential  diff 0 (-378.49 - -378.49)
 DEAL::cell 2.0 face 2 boundary
 DEAL::cell 2.0 face 3 neighbor 2.2
-DEAL::0.00 0.25 tangential  diff 0 (-118.20--118.20)
-DEAL::0.25 0.25 tangential  diff 0 (422.20-422.20)
+DEAL::0.00 0.25 tangential  diff 0 (-118.20 - -118.20)
+DEAL::0.25 0.25 tangential  diff 0 (422.20 - 422.20)
 DEAL::cell 2.1 face 0 neighbor 2.0
-DEAL::0.25 0.00 tangential  diff 0 (-106.49--106.49)
-DEAL::0.25 0.25 tangential  diff 0 (378.49-378.49)
+DEAL::0.25 0.00 tangential  diff 0 (-106.49 - -106.49)
+DEAL::0.25 0.25 tangential  diff 0 (378.49 - 378.49)
 DEAL::cell 2.1 face 1 neighbor 1.1
-DEAL::0.50 0.00 tangential  diff 0 (3.46-3.46)
-DEAL::0.50 0.25 tangential  diff 0 (0-0)
+DEAL::0.50 0.00 tangential  diff 0 (3.46 - 3.46)
+DEAL::0.50 0.25 tangential  diff 0 (0 - 0)
 DEAL::cell 2.1 face 2 boundary
 DEAL::cell 2.1 face 3 neighbor 2.3
-DEAL::0.25 0.25 tangential  diff 0 (-147.48--147.48)
-DEAL::0.50 0.25 tangential  diff 0 (531.48-531.48)
+DEAL::0.25 0.25 tangential  diff 0 (-147.48 - -147.48)
+DEAL::0.50 0.25 tangential  diff 0 (531.48 - 531.48)
 DEAL::cell 2.2 face 0 boundary
 DEAL::cell 2.2 face 1 neighbor 2.3
-DEAL::0.25 0.25 tangential  diff 0 (170.91-170.91)
-DEAL::0.25 0.50 tangential  diff 0 (-618.91--618.91)
+DEAL::0.25 0.25 tangential  diff 0 (170.91 - 170.91)
+DEAL::0.25 0.50 tangential  diff 0 (-618.91 - -618.91)
 DEAL::cell 2.2 face 2 neighbor 2.0
-DEAL::0.00 0.25 tangential  diff 0 (118.20-118.20)
-DEAL::0.25 0.25 tangential  diff 0 (-422.20--422.20)
+DEAL::0.00 0.25 tangential  diff 0 (118.20 - 118.20)
+DEAL::0.25 0.25 tangential  diff 0 (-422.20 - -422.20)
 DEAL::cell 2.2 face 3 neighbor 1.2
-DEAL::0.00 0.50 tangential  diff 0 (-26.89--26.89)
-DEAL::0.25 0.50 tangential  diff 0 (32.00-32.00)
+DEAL::0.00 0.50 tangential  diff 0 (-26.89 - -26.89)
+DEAL::0.25 0.50 tangential  diff 0 (32.00 - 32.00)
 DEAL::cell 2.3 face 0 neighbor 2.2
-DEAL::0.25 0.25 tangential  diff 0 (-170.91--170.91)
-DEAL::0.25 0.50 tangential  diff 0 (618.91-618.91)
+DEAL::0.25 0.25 tangential  diff 0 (-170.91 - -170.91)
+DEAL::0.25 0.50 tangential  diff 0 (618.91 - 618.91)
 DEAL::cell 2.3 face 1 neighbor 1.1
-DEAL::0.50 0.25 tangential  diff 0 (0-0)
-DEAL::0.50 0.50 tangential  diff 0 (-3.46--3.46)
+DEAL::0.50 0.25 tangential  diff 0 (0 - 0)
+DEAL::0.50 0.50 tangential  diff 0 (-3.46 - -3.46)
 DEAL::cell 2.3 face 2 neighbor 2.1
-DEAL::0.25 0.25 tangential  diff 0 (147.48-147.48)
-DEAL::0.50 0.25 tangential  diff 0 (-531.48--531.48)
+DEAL::0.25 0.25 tangential  diff 0 (147.48 - 147.48)
+DEAL::0.50 0.25 tangential  diff 0 (-531.48 - -531.48)
 DEAL::cell 2.3 face 3 neighbor 1.2
-DEAL::0.25 0.50 tangential  diff 0 (32.00-32.00)
-DEAL::0.50 0.50 tangential  diff 0 (90.89-90.89)
+DEAL::0.25 0.50 tangential  diff 0 (32.00 - 32.00)
+DEAL::0.50 0.50 tangential  diff 0 (90.89 - 90.89)
 DEAL::3d, 15 CELLS
 DEAL::cell 1.1 face 0 neighbor 1.0 refined
 DEAL::cell 1.1 face 1 boundary
 DEAL::cell 1.1 face 2 boundary
 DEAL::cell 1.1 face 3 neighbor 1.3
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (32.75 0.00 -12.25-32.75 0.00 -12.25)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-112.75 0.00 -23.96--112.75 0.00 -23.96)
-DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (35.67 0.00 36.25-35.67 0.00 36.25)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-123.67 0.00 79.96--123.67 0.00 79.96)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (32.75 0.00 -12.25 - 32.75 0.00 -12.25)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-112.75 0.00 -23.96 - -112.75 0.00 -23.96)
+DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (35.67 0.00 36.25 - 35.67 0.00 36.25)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-123.67 0.00 79.96 - -123.67 0.00 79.96)
 DEAL::cell 1.1 face 4 boundary
 DEAL::cell 1.1 face 5 neighbor 1.5
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-15.18 21.03 0.00--15.18 21.03 0.00)
-DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-18.11 -69.03 0.00--18.11 -69.03 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (47.18 23.96 0.00-47.18 23.96 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (58.11 -79.96 0.00-58.11 -79.96 0.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-15.18 21.03 0.00 - -15.18 21.03 0.00)
+DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-18.11 -69.03 0.00 - -18.11 -69.03 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (47.18 23.96 0.00 - 47.18 23.96 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (58.11 -79.96 0.00 - 58.11 -79.96 0.00)
 DEAL::cell 1.2 face 0 boundary
 DEAL::cell 1.2 face 1 neighbor 1.3
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -32.75 85.450.00 -32.75 85.45)
-DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -111.81 -309.450.00 -111.81 -309.45)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 112.75 97.170.00 112.75 97.17)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 407.81 -353.170.00 407.81 -353.17)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -32.75 85.45 - 0.00 -32.75 85.45)
+DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -111.81 -309.45 - 0.00 -111.81 -309.45)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 112.75 97.17 - 0.00 112.75 97.17)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 407.81 -353.17 - 0.00 407.81 -353.17)
 DEAL::cell 1.2 face 2 neighbor 1.0 refined
 DEAL::cell 1.2 face 3 boundary
 DEAL::cell 1.2 face 4 boundary
 DEAL::cell 1.2 face 5 neighbor 1.6
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-94.24 100.09 0.00--94.24 100.09 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-97.17 -364.09 0.00--97.17 -364.09 0.00)
-DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (342.24 103.02 0.00-342.24 103.02 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (353.17 -375.02 0.00-353.17 -375.02 0.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-94.24 100.09 0.00 - -94.24 100.09 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-97.17 -364.09 0.00 - -97.17 -364.09 0.00)
+DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (342.24 103.02 0.00 - 342.24 103.02 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (353.17 -375.02 0.00 - 353.17 -375.02 0.00)
 DEAL::cell 1.3 face 0 neighbor 1.2
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 32.75 -85.450.00 32.75 -85.45)
-DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 111.81 309.450.00 111.81 309.45)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -112.75 -97.170.00 -112.75 -97.17)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -407.81 353.170.00 -407.81 353.17)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 32.75 -85.45 - 0.00 32.75 -85.45)
+DEAL::0.50 1.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 111.81 309.45 - 0.00 111.81 309.45)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -112.75 -97.17 - 0.00 -112.75 -97.17)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -407.81 353.17 - 0.00 -407.81 353.17)
 DEAL::cell 1.3 face 1 boundary
 DEAL::cell 1.3 face 2 neighbor 1.1
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-32.75 0.00 12.25--32.75 0.00 12.25)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (112.75 0.00 23.96-112.75 0.00 23.96)
-DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-35.67 0.00 -36.25--35.67 0.00 -36.25)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (123.67 0.00 -79.96-123.67 0.00 -79.96)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-32.75 0.00 12.25 - -32.75 0.00 12.25)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (112.75 0.00 23.96 - 112.75 0.00 23.96)
+DEAL::1.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (-35.67 0.00 -36.25 - -35.67 0.00 -36.25)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (123.67 0.00 -79.96 - 123.67 0.00 -79.96)
 DEAL::cell 1.3 face 3 boundary
 DEAL::cell 1.3 face 4 boundary
 DEAL::cell 1.3 face 5 neighbor 1.7
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-97.17 23.96 0.00--97.17 23.96 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-164.52 -79.96 0.00--164.52 -79.96 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (353.17 167.44 0.00-353.17 167.44 0.00)
-DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (604.52 -615.44 0.00-604.52 -615.44 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-97.17 23.96 0.00 - -97.17 23.96 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-164.52 -79.96 0.00 - -164.52 -79.96 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (353.17 167.44 0.00 - 353.17 167.44 0.00)
+DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (604.52 -615.44 0.00 - 604.52 -615.44 0.00)
 DEAL::cell 1.4 face 0 boundary
 DEAL::cell 1.4 face 1 neighbor 1.5
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -226.01 15.180.00 -226.01 15.18)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -231.86 -47.180.00 -231.86 -47.18)
-DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 834.01 214.290.00 834.01 214.29)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 855.86 -790.290.00 855.86 -790.29)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -226.01 15.18 - 0.00 -226.01 15.18)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -231.86 -47.18 - 0.00 -231.86 -47.18)
+DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 834.01 214.29 - 0.00 834.01 214.29)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 855.86 -790.29 - 0.00 855.86 -790.29)
 DEAL::cell 1.4 face 2 boundary
 DEAL::cell 1.4 face 3 neighbor 1.6
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (228.94 0.00 -100.09-228.94 0.00 -100.09)
-DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-844.94 0.00 -220.15--844.94 0.00 -220.15)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (231.86 0.00 364.09-231.86 0.00 364.09)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-855.86 0.00 812.15--855.86 0.00 812.15)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (228.94 0.00 -100.09 - 228.94 0.00 -100.09)
+DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-844.94 0.00 -220.15 - -844.94 0.00 -220.15)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (231.86 0.00 364.09 - 231.86 0.00 364.09)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-855.86 0.00 812.15 - -855.86 0.00 812.15)
 DEAL::cell 1.4 face 4 neighbor 1.0 refined
 DEAL::cell 1.4 face 5 boundary
 DEAL::cell 1.5 face 0 neighbor 1.4
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 226.01 -15.180.00 226.01 -15.18)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 231.86 47.180.00 231.86 47.18)
-DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -834.01 -214.290.00 -834.01 -214.29)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -855.86 790.290.00 -855.86 790.29)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 226.01 -15.18 - 0.00 226.01 -15.18)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 231.86 47.18 - 0.00 231.86 47.18)
+DEAL::0.50 0.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -834.01 -214.29 - 0.00 -834.01 -214.29)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -855.86 790.29 - 0.00 -855.86 790.29)
 DEAL::cell 1.5 face 1 boundary
 DEAL::cell 1.5 face 2 boundary
 DEAL::cell 1.5 face 3 neighbor 1.7
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (231.86 0.00 -23.96-231.86 0.00 -23.96)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-855.86 0.00 -284.57--855.86 0.00 -284.57)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (290.43 0.00 79.96-290.43 0.00 79.96)
-DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-1074.43 0.00 1052.57--1074.43 0.00 1052.57)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (231.86 0.00 -23.96 - 231.86 0.00 -23.96)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-855.86 0.00 -284.57 - -855.86 0.00 -284.57)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (290.43 0.00 79.96 - 290.43 0.00 79.96)
+DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (-1074.43 0.00 1052.57 - -1074.43 0.00 1052.57)
 DEAL::cell 1.5 face 4 neighbor 1.1
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (15.18 -21.03 0.00-15.18 -21.03 0.00)
-DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (18.11 69.03 0.00-18.11 69.03 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-47.18 -23.96 0.00--47.18 -23.96 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-58.11 79.96 0.00--58.11 79.96 0.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (15.18 -21.03 0.00 - 15.18 -21.03 0.00)
+DEAL::1.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (18.11 69.03 0.00 - 18.11 69.03 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-47.18 -23.96 0.00 - -47.18 -23.96 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-58.11 79.96 0.00 - -58.11 79.96 0.00)
 DEAL::cell 1.5 face 5 boundary
 DEAL::cell 1.6 face 0 boundary
 DEAL::cell 1.6 face 1 neighbor 1.7
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -231.86 97.170.00 -231.86 97.17)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -337.28 -353.170.00 -337.28 -353.17)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 855.86 328.490.00 855.86 328.49)
-DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 1249.28 -1216.490.00 1249.28 -1216.49)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -231.86 97.17 - 0.00 -231.86 97.17)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -337.28 -353.17 - 0.00 -337.28 -353.17)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 855.86 328.49 - 0.00 855.86 328.49)
+DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 1249.28 -1216.49 - 0.00 1249.28 -1216.49)
 DEAL::cell 1.6 face 2 neighbor 1.4
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-228.94 0.00 100.09--228.94 0.00 100.09)
-DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (844.94 0.00 220.15-844.94 0.00 220.15)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-231.86 0.00 -364.09--231.86 0.00 -364.09)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (855.86 0.00 -812.15-855.86 0.00 -812.15)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-228.94 0.00 100.09 - -228.94 0.00 100.09)
+DEAL::0.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (844.94 0.00 220.15 - 844.94 0.00 220.15)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-231.86 0.00 -364.09 - -231.86 0.00 -364.09)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (855.86 0.00 -812.15 - 855.86 0.00 -812.15)
 DEAL::cell 1.6 face 3 boundary
 DEAL::cell 1.6 face 4 neighbor 1.2
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (94.24 -100.09 0.00-94.24 -100.09 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (97.17 364.09 0.00-97.17 364.09 0.00)
-DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-342.24 -103.02 0.00--342.24 -103.02 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-353.17 375.02 0.00--353.17 375.02 0.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (94.24 -100.09 0.00 - 94.24 -100.09 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (97.17 364.09 0.00 - 97.17 364.09 0.00)
+DEAL::0.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-342.24 -103.02 0.00 - -342.24 -103.02 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-353.17 375.02 0.00 - -353.17 375.02 0.00)
 DEAL::cell 1.6 face 5 boundary
 DEAL::cell 1.7 face 0 neighbor 1.6
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 231.86 -97.170.00 231.86 -97.17)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 337.28 353.170.00 337.28 353.17)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -855.86 -328.490.00 -855.86 -328.49)
-DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -1249.28 1216.490.00 -1249.28 1216.49)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 231.86 -97.17 - 0.00 231.86 -97.17)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 337.28 353.17 - 0.00 337.28 353.17)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -855.86 -328.49 - 0.00 -855.86 -328.49)
+DEAL::0.50 1.00 1.00 tangential  diff 0.00 0.00 0.00 (0.00 -1249.28 1216.49 - 0.00 -1249.28 1216.49)
 DEAL::cell 1.7 face 1 boundary
 DEAL::cell 1.7 face 2 neighbor 1.5
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-231.86 0.00 23.96--231.86 0.00 23.96)
-DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (855.86 0.00 284.57-855.86 0.00 284.57)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-290.43 0.00 -79.96--290.43 0.00 -79.96)
-DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (1074.43 0.00 -1052.57-1074.43 0.00 -1052.57)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-231.86 0.00 23.96 - -231.86 0.00 23.96)
+DEAL::0.50 0.50 1.00 tangential  diff 0.00 0.00 0.00 (855.86 0.00 284.57 - 855.86 0.00 284.57)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-290.43 0.00 -79.96 - -290.43 0.00 -79.96)
+DEAL::1.00 0.50 1.00 tangential  diff 0.00 0.00 0.00 (1074.43 0.00 -1052.57 - 1074.43 0.00 -1052.57)
 DEAL::cell 1.7 face 3 boundary
 DEAL::cell 1.7 face 4 neighbor 1.3
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (97.17 -23.96 0.00-97.17 -23.96 0.00)
-DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (164.52 79.96 0.00-164.52 79.96 0.00)
-DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-353.17 -167.44 0.00--353.17 -167.44 0.00)
-DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-604.52 615.44 0.00--604.52 615.44 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (97.17 -23.96 0.00 - 97.17 -23.96 0.00)
+DEAL::1.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (164.52 79.96 0.00 - 164.52 79.96 0.00)
+DEAL::0.50 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-353.17 -167.44 0.00 - -353.17 -167.44 0.00)
+DEAL::1.00 1.00 0.50 tangential  diff 0.00 0.00 0.00 (-604.52 615.44 0.00 - -604.52 615.44 0.00)
 DEAL::cell 1.7 face 5 boundary
 DEAL::cell 2.0 face 0 boundary
 DEAL::cell 2.0 face 1 neighbor 2.1
-DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -867.82 820.970.00 -867.82 820.97)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -879.53 -3044.970.00 -879.53 -3044.97)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 3219.82 844.390.00 3219.82 844.39)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 3263.53 -3132.390.00 3263.53 -3132.39)
+DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -867.82 820.97 - 0.00 -867.82 820.97)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -879.53 -3044.97 - 0.00 -879.53 -3044.97)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 3219.82 844.39 - 0.00 3219.82 844.39)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 3263.53 -3132.39 - 0.00 3263.53 -3132.39)
 DEAL::cell 2.0 face 2 boundary
 DEAL::cell 2.0 face 3 neighbor 2.2
-DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (873.68 0.00 -832.68-873.68 0.00 -832.68)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3241.68 0.00 -856.11--3241.68 0.00 -856.11)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (879.53 0.00 3088.68-879.53 0.00 3088.68)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3263.53 0.00 3176.11--3263.53 0.00 3176.11)
+DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (873.68 0.00 -832.68 - 873.68 0.00 -832.68)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3241.68 0.00 -856.11 - -3241.68 0.00 -856.11)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (879.53 0.00 3088.68 - 879.53 0.00 3088.68)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3263.53 0.00 3176.11 - -3263.53 0.00 3176.11)
 DEAL::cell 2.0 face 4 boundary
 DEAL::cell 2.0 face 5 neighbor 2.4
-DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-838.54 850.25 0.00--838.54 850.25 0.00)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-844.39 -3154.25 0.00--844.39 -3154.25 0.00)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3110.54 856.11 0.00-3110.54 856.11 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3132.39 -3176.11 0.00-3132.39 -3176.11 0.00)
+DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-838.54 850.25 0.00 - -838.54 850.25 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-844.39 -3154.25 0.00 - -844.39 -3154.25 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3110.54 856.11 0.00 - 3110.54 856.11 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3132.39 -3176.11 0.00 - 3132.39 -3176.11 0.00)
 DEAL::cell 2.1 face 0 neighbor 2.0
-DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 867.82 -820.970.00 867.82 -820.97)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 879.53 3044.970.00 879.53 3044.97)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -3219.82 -844.390.00 -3219.82 -844.39)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -3263.53 3132.390.00 -3263.53 3132.39)
+DEAL::0.25 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 867.82 -820.97 - 0.00 867.82 -820.97)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 879.53 3044.97 - 0.00 879.53 3044.97)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -3219.82 -844.39 - 0.00 -3219.82 -844.39)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -3263.53 3132.39 - 0.00 -3263.53 3132.39)
 DEAL::cell 2.1 face 1 neighbor 1.1
-DEAL::0.50 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -26.89 3.460.00 -26.89 3.46)
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -11.83 0.000.00 -11.83 0.00)
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 32.00 -7.390.00 32.00 -7.39)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.780.00 13.48 12.78)
+DEAL::0.50 0.00 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -26.89 3.46 - 0.00 -26.89 3.46)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -11.83 0.00 - 0.00 -11.83 0.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 32.00 -7.39 - 0.00 32.00 -7.39)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.78 - 0.00 13.48 12.78)
 DEAL::cell 2.1 face 2 boundary
 DEAL::cell 2.1 face 3 neighbor 2.3
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (879.53 0.00 -984.95-879.53 0.00 -984.95)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3263.53 0.00 -1002.52--3263.53 0.00 -1002.52)
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (11.83 0.00 3656.95-11.83 0.00 3656.95)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-13.48 0.00 3722.52--13.48 0.00 3722.52)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (879.53 0.00 -984.95 - 879.53 0.00 -984.95)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3263.53 0.00 -1002.52 - -3263.53 0.00 -1002.52)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (11.83 0.00 3656.95 - 11.83 0.00 3656.95)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-13.48 0.00 3722.52 - -13.48 0.00 3722.52)
 DEAL::cell 2.1 face 4 boundary
 DEAL::cell 2.1 face 5 neighbor 2.5
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-844.39 996.66 0.00--844.39 996.66 0.00)
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (7.39 -3700.66 0.00-7.39 -3700.66 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3132.39 1002.52 0.00-3132.39 1002.52 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-12.78 -3722.52 0.00--12.78 -3722.52 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-844.39 996.66 0.00 - -844.39 996.66 0.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (7.39 -3700.66 0.00 - 7.39 -3700.66 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3132.39 1002.52 0.00 - 3132.39 1002.52 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-12.78 -3722.52 0.00 - -12.78 -3722.52 0.00)
 DEAL::cell 2.2 face 0 boundary
 DEAL::cell 2.2 face 1 neighbor 2.3
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -879.53 1102.080.00 -879.53 1102.08)
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -14.59 -4094.080.00 -14.59 -4094.08)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 3263.53 1119.650.00 3263.53 1119.65)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 17.25 -4159.650.00 17.25 -4159.65)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -879.53 1102.08 - 0.00 -879.53 1102.08)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -14.59 -4094.08 - 0.00 -14.59 -4094.08)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 3263.53 1119.65 - 0.00 3263.53 1119.65)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 17.25 -4159.65 - 0.00 17.25 -4159.65)
 DEAL::cell 2.2 face 2 neighbor 2.0
-DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-873.68 0.00 832.68--873.68 0.00 832.68)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3241.68 0.00 856.11-3241.68 0.00 856.11)
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-879.53 0.00 -3088.68--879.53 0.00 -3088.68)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3263.53 0.00 -3176.11-3263.53 0.00 -3176.11)
+DEAL::0.00 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-873.68 0.00 832.68 - -873.68 0.00 832.68)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3241.68 0.00 856.11 - 3241.68 0.00 856.11)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-879.53 0.00 -3088.68 - -879.53 0.00 -3088.68)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3263.53 0.00 -3176.11 - 3263.53 0.00 -3176.11)
 DEAL::cell 2.2 face 3 neighbor 1.2
-DEAL::0.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (105.95 0.00 -88.38-105.95 0.00 -88.38)
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-140.00 0.00 -38.22--140.00 0.00 -38.22)
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (14.59 0.00 116.00-14.59 0.00 116.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52--17.25 0.00 49.52)
+DEAL::0.00 0.50 0.00 tangential  diff 0.00 0.00 0.00 (105.95 0.00 -88.38 - 105.95 0.00 -88.38)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-140.00 0.00 -38.22 - -140.00 0.00 -38.22)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (14.59 0.00 116.00 - 14.59 0.00 116.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52 - -17.25 0.00 49.52)
 DEAL::cell 2.2 face 4 boundary
 DEAL::cell 2.2 face 5 neighbor 2.6
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1113.79 856.11 0.00--1113.79 856.11 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1119.65 -3176.11 0.00--1119.65 -3176.11 0.00)
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4137.79 38.22 0.00-4137.79 38.22 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4159.65 -49.52 0.00-4159.65 -49.52 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1113.79 856.11 0.00 - -1113.79 856.11 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1119.65 -3176.11 0.00 - -1119.65 -3176.11 0.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4137.79 38.22 0.00 - 4137.79 38.22 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4159.65 -49.52 0.00 - 4159.65 -49.52 0.00)
 DEAL::cell 2.3 face 0 neighbor 2.2
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 879.53 -1102.080.00 879.53 -1102.08)
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 14.59 4094.080.00 14.59 4094.08)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -3263.53 -1119.650.00 -3263.53 -1119.65)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -17.25 4159.650.00 -17.25 4159.65)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 879.53 -1102.08 - 0.00 879.53 -1102.08)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 14.59 4094.08 - 0.00 14.59 4094.08)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -3263.53 -1119.65 - 0.00 -3263.53 -1119.65)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -17.25 4159.65 - 0.00 -17.25 4159.65)
 DEAL::cell 2.3 face 1 neighbor 1.1
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -11.83 0.000.00 -11.83 0.00)
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -32.75 -3.460.00 -32.75 -3.46)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.780.00 13.48 12.78)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 40.00 32.960.00 40.00 32.96)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -11.83 0.00 - 0.00 -11.83 0.00)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (0.00 -32.75 -3.46 - 0.00 -32.75 -3.46)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.78 - 0.00 13.48 12.78)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 40.00 32.96 - 0.00 40.00 32.96)
 DEAL::cell 2.3 face 2 neighbor 2.1
-DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-879.53 0.00 984.95--879.53 0.00 984.95)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3263.53 0.00 1002.52-3263.53 0.00 1002.52)
-DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-11.83 0.00 -3656.95--11.83 0.00 -3656.95)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (13.48 0.00 -3722.52-13.48 0.00 -3722.52)
+DEAL::0.25 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-879.53 0.00 984.95 - -879.53 0.00 984.95)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (3263.53 0.00 1002.52 - 3263.53 0.00 1002.52)
+DEAL::0.50 0.25 0.00 tangential  diff 0.00 0.00 0.00 (-11.83 0.00 -3656.95 - -11.83 0.00 -3656.95)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (13.48 0.00 -3722.52 - 13.48 0.00 -3722.52)
 DEAL::cell 2.3 face 3 neighbor 1.2
-DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (14.59 0.00 116.00-14.59 0.00 116.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52--17.25 0.00 49.52)
-DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (32.75 0.00 320.38-32.75 0.00 320.38)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 137.26--40.00 0.00 137.26)
+DEAL::0.25 0.50 0.00 tangential  diff 0.00 0.00 0.00 (14.59 0.00 116.00 - 14.59 0.00 116.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52 - -17.25 0.00 49.52)
+DEAL::0.50 0.50 0.00 tangential  diff 0.00 0.00 0.00 (32.75 0.00 320.38 - 32.75 0.00 320.38)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 137.26 - -40.00 0.00 137.26)
 DEAL::cell 2.3 face 4 boundary
 DEAL::cell 2.3 face 5 neighbor 2.7
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1119.65 1002.52 0.00--1119.65 1002.52 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-12.78 -3722.52 0.00--12.78 -3722.52 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4159.65 -49.52 0.00-4159.65 -49.52 0.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-32.96 -137.26 0.00--32.96 -137.26 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1119.65 1002.52 0.00 - -1119.65 1002.52 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-12.78 -3722.52 0.00 - -12.78 -3722.52 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (4159.65 -49.52 0.00 - 4159.65 -49.52 0.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-32.96 -137.26 0.00 - -32.96 -137.26 0.00)
 DEAL::cell 2.4 face 0 boundary
 DEAL::cell 2.4 face 1 neighbor 2.5
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -1342.19 844.390.00 -1342.19 844.39)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -1353.90 -3132.390.00 -1353.90 -3132.39)
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 4990.19 -2.740.00 4990.19 -2.74)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 5033.90 6.420.00 5033.90 6.42)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -1342.19 844.39 - 0.00 -1342.19 844.39)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -1353.90 -3132.39 - 0.00 -1353.90 -3132.39)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 4990.19 -2.74 - 0.00 4990.19 -2.74)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 5033.90 6.42 - 0.00 5033.90 6.42)
 DEAL::cell 2.4 face 2 boundary
 DEAL::cell 2.4 face 3 neighbor 2.6
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1348.05 0.00 -856.11-1348.05 0.00 -856.11)
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-5012.05 0.00 -42.45--5012.05 0.00 -42.45)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1353.90 0.00 3176.11-1353.90 0.00 3176.11)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-5033.90 0.00 55.31--5033.90 0.00 55.31)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1348.05 0.00 -856.11 - 1348.05 0.00 -856.11)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-5012.05 0.00 -42.45 - -5012.05 0.00 -42.45)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1353.90 0.00 3176.11 - 1353.90 0.00 3176.11)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-5033.90 0.00 55.31 - -5033.90 0.00 55.31)
 DEAL::cell 2.4 face 4 neighbor 2.0
-DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (838.54 -850.25 0.00-838.54 -850.25 0.00)
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (844.39 3154.25 0.00-844.39 3154.25 0.00)
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3110.54 -856.11 0.00--3110.54 -856.11 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3132.39 3176.11 0.00--3132.39 3176.11 0.00)
+DEAL::0.00 0.00 0.25 tangential  diff 0.00 0.00 0.00 (838.54 -850.25 0.00 - 838.54 -850.25 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (844.39 3154.25 0.00 - 844.39 3154.25 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3110.54 -856.11 0.00 - -3110.54 -856.11 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3132.39 3176.11 0.00 - -3132.39 3176.11 0.00)
 DEAL::cell 2.4 face 5 neighbor 1.4
-DEAL::0.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-205.51 208.44 0.00--205.51 208.44 0.00)
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (2.74 -280.00 0.00-2.74 -280.00 0.00)
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (276.00 42.45 0.00-276.00 42.45 0.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00--6.42 -55.31 0.00)
+DEAL::0.00 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-205.51 208.44 0.00 - -205.51 208.44 0.00)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (2.74 -280.00 0.00 - 2.74 -280.00 0.00)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (276.00 42.45 0.00 - 276.00 42.45 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00 - -6.42 -55.31 0.00)
 DEAL::cell 2.5 face 0 neighbor 2.4
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 1342.19 -844.390.00 1342.19 -844.39)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 1353.90 3132.390.00 1353.90 3132.39)
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -4990.19 2.740.00 -4990.19 2.74)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -5033.90 -6.420.00 -5033.90 -6.42)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 1342.19 -844.39 - 0.00 1342.19 -844.39)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 1353.90 3132.39 - 0.00 1353.90 3132.39)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -4990.19 2.74 - 0.00 -4990.19 2.74)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -5033.90 -6.42 - 0.00 -5033.90 -6.42)
 DEAL::cell 2.5 face 1 neighbor 1.1
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 32.00 -7.390.00 32.00 -7.39)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.780.00 13.48 12.78)
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 90.89 15.180.00 90.89 15.18)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 38.80 -16.000.00 38.80 -16.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (0.00 32.00 -7.39 - 0.00 32.00 -7.39)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.78 - 0.00 13.48 12.78)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (0.00 90.89 15.18 - 0.00 90.89 15.18)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 38.80 -16.00 - 0.00 38.80 -16.00)
 DEAL::cell 2.5 face 2 boundary
 DEAL::cell 2.5 face 3 neighbor 2.7
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1353.90 0.00 -1002.52-1353.90 0.00 -1002.52)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-5033.90 0.00 55.31--5033.90 0.00 55.31)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-13.48 0.00 3722.52--13.48 0.00 3722.52)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-38.80 0.00 153.08--38.80 0.00 153.08)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1353.90 0.00 -1002.52 - 1353.90 0.00 -1002.52)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-5033.90 0.00 55.31 - -5033.90 0.00 55.31)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-13.48 0.00 3722.52 - -13.48 0.00 3722.52)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-38.80 0.00 153.08 - -38.80 0.00 153.08)
 DEAL::cell 2.5 face 4 neighbor 2.1
-DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (844.39 -996.66 0.00-844.39 -996.66 0.00)
-DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-7.39 3700.66 0.00--7.39 3700.66 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3132.39 -1002.52 0.00--3132.39 -1002.52 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (12.78 3722.52 0.00-12.78 3722.52 0.00)
+DEAL::0.25 0.00 0.25 tangential  diff 0.00 0.00 0.00 (844.39 -996.66 0.00 - 844.39 -996.66 0.00)
+DEAL::0.50 0.00 0.25 tangential  diff 0.00 0.00 0.00 (-7.39 3700.66 0.00 - -7.39 3700.66 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-3132.39 -1002.52 0.00 - -3132.39 -1002.52 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (12.78 3722.52 0.00 - 12.78 3722.52 0.00)
 DEAL::cell 2.5 face 5 neighbor 1.4
-DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (2.74 -280.00 0.00-2.74 -280.00 0.00)
-DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-15.18 -768.44 0.00--15.18 -768.44 0.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00--6.42 -55.31 0.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (16.00 -153.08 0.00-16.00 -153.08 0.00)
+DEAL::0.25 0.00 0.50 tangential  diff 0.00 0.00 0.00 (2.74 -280.00 0.00 - 2.74 -280.00 0.00)
+DEAL::0.50 0.00 0.50 tangential  diff 0.00 0.00 0.00 (-15.18 -768.44 0.00 - -15.18 -768.44 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00 - -6.42 -55.31 0.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (16.00 -153.08 0.00 - 16.00 -153.08 0.00)
 DEAL::cell 2.6 face 0 boundary
 DEAL::cell 2.6 face 1 neighbor 2.7
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -1353.90 1119.650.00 -1353.90 1119.65)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 17.25 -4159.650.00 17.25 -4159.65)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 5033.90 6.420.00 5033.90 6.42)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 49.10 15.580.00 49.10 15.58)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -1353.90 1119.65 - 0.00 -1353.90 1119.65)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 17.25 -4159.65 - 0.00 17.25 -4159.65)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 5033.90 6.42 - 0.00 5033.90 6.42)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 49.10 15.58 - 0.00 49.10 15.58)
 DEAL::cell 2.6 face 2 neighbor 2.4
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1348.05 0.00 856.11--1348.05 0.00 856.11)
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (5012.05 0.00 42.45-5012.05 0.00 42.45)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1353.90 0.00 -3176.11--1353.90 0.00 -3176.11)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (5033.90 0.00 -55.31-5033.90 0.00 -55.31)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1348.05 0.00 856.11 - -1348.05 0.00 856.11)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (5012.05 0.00 42.45 - 5012.05 0.00 42.45)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1353.90 0.00 -3176.11 - -1353.90 0.00 -3176.11)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (5033.90 0.00 -55.31 - 5033.90 0.00 -55.31)
 DEAL::cell 2.6 face 3 neighbor 1.2
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-140.00 0.00 -38.22--140.00 0.00 -38.22)
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-385.95 0.00 -100.09--385.95 0.00 -100.09)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52--17.25 0.00 49.52)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-49.10 0.00 132.00--49.10 0.00 132.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-140.00 0.00 -38.22 - -140.00 0.00 -38.22)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-385.95 0.00 -100.09 - -385.95 0.00 -100.09)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52 - -17.25 0.00 49.52)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-49.10 0.00 132.00 - -49.10 0.00 132.00)
 DEAL::cell 2.6 face 4 neighbor 2.2
-DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1113.79 -856.11 0.00-1113.79 -856.11 0.00)
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1119.65 3176.11 0.00-1119.65 3176.11 0.00)
-DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4137.79 -38.22 0.00--4137.79 -38.22 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4159.65 49.52 0.00--4159.65 49.52 0.00)
+DEAL::0.00 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1113.79 -856.11 0.00 - 1113.79 -856.11 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1119.65 3176.11 0.00 - 1119.65 3176.11 0.00)
+DEAL::0.00 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4137.79 -38.22 0.00 - -4137.79 -38.22 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4159.65 49.52 0.00 - -4159.65 49.52 0.00)
 DEAL::cell 2.6 face 5 neighbor 1.4
-DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (276.00 42.45 0.00-276.00 42.45 0.00)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00--6.42 -55.31 0.00)
-DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (757.51 100.09 0.00-757.51 100.09 0.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-15.58 -132.00 0.00--15.58 -132.00 0.00)
+DEAL::0.00 0.25 0.50 tangential  diff 0.00 0.00 0.00 (276.00 42.45 0.00 - 276.00 42.45 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00 - -6.42 -55.31 0.00)
+DEAL::0.00 0.50 0.50 tangential  diff 0.00 0.00 0.00 (757.51 100.09 0.00 - 757.51 100.09 0.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-15.58 -132.00 0.00 - -15.58 -132.00 0.00)
 DEAL::cell 2.7 face 0 neighbor 2.6
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 1353.90 -1119.650.00 1353.90 -1119.65)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -17.25 4159.650.00 -17.25 4159.65)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -5033.90 -6.420.00 -5033.90 -6.42)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -49.10 -15.580.00 -49.10 -15.58)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 1353.90 -1119.65 - 0.00 1353.90 -1119.65)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 -17.25 4159.65 - 0.00 -17.25 4159.65)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -5033.90 -6.42 - 0.00 -5033.90 -6.42)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 -49.10 -15.58 - 0.00 -49.10 -15.58)
 DEAL::cell 2.7 face 1 neighbor 1.1
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.780.00 13.48 12.78)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 40.00 32.960.00 40.00 32.96)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 38.80 -16.000.00 38.80 -16.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 112.75 -47.180.00 112.75 -47.18)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (0.00 13.48 12.78 - 0.00 13.48 12.78)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (0.00 40.00 32.96 - 0.00 40.00 32.96)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (0.00 38.80 -16.00 - 0.00 38.80 -16.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (0.00 112.75 -47.18 - 0.00 112.75 -47.18)
 DEAL::cell 2.7 face 2 neighbor 2.5
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1353.90 0.00 1002.52--1353.90 0.00 1002.52)
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (5033.90 0.00 -55.31-5033.90 0.00 -55.31)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (13.48 0.00 -3722.52-13.48 0.00 -3722.52)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (38.80 0.00 -153.08-38.80 0.00 -153.08)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (-1353.90 0.00 1002.52 - -1353.90 0.00 1002.52)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (5033.90 0.00 -55.31 - 5033.90 0.00 -55.31)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (13.48 0.00 -3722.52 - 13.48 0.00 -3722.52)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (38.80 0.00 -153.08 - 38.80 0.00 -153.08)
 DEAL::cell 2.7 face 3 neighbor 1.2
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52--17.25 0.00 49.52)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-49.10 0.00 132.00--49.10 0.00 132.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 137.26--40.00 0.00 137.26)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-112.75 0.00 364.09--112.75 0.00 364.09)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-17.25 0.00 49.52 - -17.25 0.00 49.52)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-49.10 0.00 132.00 - -49.10 0.00 132.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-40.00 0.00 137.26 - -40.00 0.00 137.26)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-112.75 0.00 364.09 - -112.75 0.00 364.09)
 DEAL::cell 2.7 face 4 neighbor 2.3
-DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1119.65 -1002.52 0.00-1119.65 -1002.52 0.00)
-DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (12.78 3722.52 0.00-12.78 3722.52 0.00)
-DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4159.65 49.52 0.00--4159.65 49.52 0.00)
-DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (32.96 137.26 0.00-32.96 137.26 0.00)
+DEAL::0.25 0.25 0.25 tangential  diff 0.00 0.00 0.00 (1119.65 -1002.52 0.00 - 1119.65 -1002.52 0.00)
+DEAL::0.50 0.25 0.25 tangential  diff 0.00 0.00 0.00 (12.78 3722.52 0.00 - 12.78 3722.52 0.00)
+DEAL::0.25 0.50 0.25 tangential  diff 0.00 0.00 0.00 (-4159.65 49.52 0.00 - -4159.65 49.52 0.00)
+DEAL::0.50 0.50 0.25 tangential  diff 0.00 0.00 0.00 (32.96 137.26 0.00 - 32.96 137.26 0.00)
 DEAL::cell 2.7 face 5 neighbor 1.4
-DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00--6.42 -55.31 0.00)
-DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (16.00 -153.08 0.00-16.00 -153.08 0.00)
-DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-15.58 -132.00 0.00--15.58 -132.00 0.00)
-DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (47.18 -364.09 0.00-47.18 -364.09 0.00)
+DEAL::0.25 0.25 0.50 tangential  diff 0.00 0.00 0.00 (-6.42 -55.31 0.00 - -6.42 -55.31 0.00)
+DEAL::0.50 0.25 0.50 tangential  diff 0.00 0.00 0.00 (16.00 -153.08 0.00 - 16.00 -153.08 0.00)
+DEAL::0.25 0.50 0.50 tangential  diff 0.00 0.00 0.00 (-15.58 -132.00 0.00 - -15.58 -132.00 0.00)
+DEAL::0.50 0.50 0.50 tangential  diff 0.00 0.00 0.00 (47.18 -364.09 0.00 - 47.18 -364.09 0.00)