From: Denis Davydov <davydden@gmail.com>
Date: Wed, 2 Sep 2015 20:51:22 +0000 (+0200)
Subject: test make_hanging_node_constraints() with hp-ref and neither_element_dominates
X-Git-Tag: v8.4.0-rc2~484^2~1
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=599b12b1549383f2cbe995b6d3dd37ca0f3974c8;p=dealii.git

test make_hanging_node_constraints() with hp-ref and neither_element_dominates
---

diff --git a/tests/hp/hp_constraints_neither_dominate_02.cc b/tests/hp/hp_constraints_neither_dominate_02.cc
new file mode 100644
index 0000000000..a9b6626cc4
--- /dev/null
+++ b/tests/hp/hp_constraints_neither_dominate_02.cc
@@ -0,0 +1,456 @@
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2005 - 2015 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE at
+// the top level of the deal.II distribution.
+//
+// ---------------------------------------------------------------------
+
+
+
+// test hp_constraints when neither_element_dominates.
+// 2D and 3D hp-refinement.
+// Check continuity across the face by evaluating the
+// field from both sides.
+
+#include <deal.II/base/logstream.h>
+#include <deal.II/base/utilities.h>
+#include <deal.II/base/index_set.h>
+#include <deal.II/base/conditional_ostream.h>
+#include <deal.II/base/function.h>
+#include <deal.II/base/function_parser.h>
+
+#include <deal.II/numerics/data_postprocessor.h>
+#include <deal.II/numerics/data_out.h>
+
+#include <deal.II/hp/dof_handler.h>
+#include <deal.II/hp/fe_values.h>
+#include <deal.II/hp/q_collection.h>
+#include <deal.II/hp/fe_collection.h>
+
+
+#include <deal.II/dofs/dof_renumbering.h>
+#include <deal.II/dofs/dof_tools.h>
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_nothing.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/fe_tools.h>
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/lac/sparsity_tools.h>
+#include <deal.II/lac/compressed_simple_sparsity_pattern.h>
+#include <deal.II/hp/fe_values.h>
+
+#include <deal.II/lac/petsc_parallel_vector.h>
+#include <deal.II/lac/petsc_parallel_sparse_matrix.h>
+#include <deal.II/lac/petsc_solver.h>
+#include <deal.II/lac/petsc_precondition.h>
+
+#include <deal.II/lac/slepc_solver.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <deal.II/grid/grid_refinement.h>
+#include <deal.II/numerics/error_estimator.h>
+
+#include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_in.h>
+#include <deal.II/grid/grid_tools.h>
+#include <deal.II/lac/vector.h>
+
+#include "../tests.h"
+
+#include <fstream>
+#include <iostream>
+
+//#define DEBUG_OUTPUT_VTK
+unsigned int counter=0;
+
+const double eps = 1e-10;
+
+using namespace dealii;
+
+template <int dim>
+struct less_than_key
+{
+  inline bool operator() (const std::pair<Point<dim>,Vector<double>> &pair1, const std::pair<Point<dim>,Vector<double>> &pair2)
+  {
+    const double precision = 1e-3;
+    const Point<dim> &p1 = pair1.first;
+    const Point<dim> &p2 = pair2.first;
+
+    for (unsigned int d = 0; d < dim; d++)
+      {
+        const bool is_equal = ( std::abs(p1[d] - p2[d]) < precision );
+        if (!is_equal)
+          return p1[d] < p2[d];
+      }
+
+    return false; //< They are equal
+  }
+};
+
+//  ------------
+//  |0 |0 |    |
+//  |-----| 2  |
+//  |1 |1 |    |
+//  ------------
+template<int dim>
+void test2cells(const FiniteElement<dim> &fe_0,
+                const FiniteElement<dim> &fe_1,
+                const FiniteElement<dim> &fe_2,
+                const FiniteElement<dim> &fe_common)
+{
+  const unsigned int n_comp = fe_0.n_components ();
+  Triangulation<dim>   triangulation;
+  {
+    Point<dim> p1,p2;
+    for (unsigned int d=0; d<dim; d++)
+      p2[d] = 1.0;
+    p1[0] = -1.0;
+    std::vector<unsigned int> repetitoins(dim,1);
+    repetitoins[0]=2;
+    GridGenerator::subdivided_hyper_rectangle(triangulation,
+                                              repetitoins,
+                                              p1,
+                                              p2);
+  }
+
+  hp::DoFHandler<dim> dof_handler(triangulation);
+
+  hp::FECollection<dim> fe_collection;
+  fe_collection.push_back(fe_0);
+  fe_collection.push_back(fe_1);
+  fe_collection.push_back(fe_2);
+  fe_collection.push_back(fe_common);
+
+  const QIterated<dim-1> quad_formula(QTrapez<1>(),2);
+
+  hp::QCollection< dim-1 > q_face_collection;
+  q_face_collection.push_back(quad_formula);
+  q_face_collection.push_back(quad_formula);
+  q_face_collection.push_back(quad_formula);
+  q_face_collection.push_back(quad_formula);
+
+  deallog << "hp-ref: "<<fe_0.get_name()<< " vs "<<fe_1.get_name()<<" vs " << fe_2.get_name() <<" with "<< fe_common.get_name()<<std::endl;
+
+  // refine once:
+  triangulation.begin_active()->set_refine_flag();
+  triangulation.execute_coarsening_and_refinement();
+
+  for (typename hp::DoFHandler<dim>::active_cell_iterator
+       cell = dof_handler.begin_active();
+       cell != dof_handler.end(); ++cell)
+    {
+      const Point<dim> &center = cell->center();
+      if (center[0] > 0.0 ) // right cell
+        cell->set_active_fe_index(2);
+      else if (center[1] > 0.5)
+        cell->set_active_fe_index(0);
+      else
+        cell->set_active_fe_index(1);
+    }
+  dof_handler.distribute_dofs(fe_collection);
+
+  ConstraintMatrix constraints;
+  constraints.clear();
+  dealii::DoFTools::make_hanging_node_constraints  (dof_handler, constraints);
+  constraints.close ();
+
+  constraints.print(deallog.get_file_stream());
+
+#ifdef DEBUG_OUTPUT_VTK
+  // output to check if all is good:
+  counter++;
+  std::vector<Vector<double>> shape_functions;
+  std::vector<std::string> names;
+  for (unsigned int s=0; s < dof_handler.n_dofs(); s++)
+    {
+      Vector<double> shape_function;
+      shape_function.reinit(dof_handler.n_dofs());
+      shape_function[s] = 1.0;
+
+      // if the dof is constrained, first output unconstrained vector
+      if (constraints.is_constrained(s))
+        {
+          names.push_back(std::string("UN_") +
+                          dealii::Utilities::int_to_string(s,2));
+          shape_functions.push_back(shape_function);
+        }
+
+      names.push_back(std::string("N_") +
+                      dealii::Utilities::int_to_string(s,2));
+
+      // make continuous/constrain:
+      constraints.distribute(shape_function);
+      shape_functions.push_back(shape_function);
+    }
+
+  DataOut<dim,hp::DoFHandler<dim>> data_out;
+  data_out.attach_dof_handler (dof_handler);
+
+  // get material ids:
+  Vector<float> fe_index(triangulation.n_active_cells());
+  typename hp::DoFHandler<dim>::active_cell_iterator
+  cell = dof_handler.begin_active (),
+  endc = dof_handler.end ();
+  for (unsigned int index=0; cell!=endc; ++cell,++index)
+    {
+      fe_index[index] = cell->active_fe_index();
+    }
+  data_out.add_data_vector(fe_index, "fe_index");
+
+  for (unsigned int i = 0; i < shape_functions.size(); i++)
+    data_out.add_data_vector (shape_functions[i], names[i]);
+
+  data_out.build_patches(0);
+
+  std::string filename = "shape_functions_"+dealii::Utilities::int_to_string(counter,1)
+                         +"_"+dealii::Utilities::int_to_string(dim)+"D.vtu";
+  std::ofstream output (filename.c_str ());
+  data_out.write_vtu (output);
+#endif
+
+  // fill some vector
+  Vector<double> solution(dof_handler.n_dofs());
+  for (unsigned int dof=0; dof < dof_handler.n_dofs(); dof++)
+    solution[dof] = 1.5*(dof%7) + 0.5*dim + 2.0*(dof%3);
+
+  constraints.distribute(solution);
+
+
+  // evaluate field at the interface:
+  hp::FEFaceValues< dim > fe_face_values_hp(fe_collection,
+                                            q_face_collection,
+                                            update_values | update_quadrature_points);
+
+  std::vector<unsigned int> local_face_dof_indices;
+  std::vector<Vector<double> > values;
+
+  std::vector<std::pair<Point<dim>,Vector<double>> > pairs_point_value;
+
+  for (typename hp::DoFHandler<dim>::active_cell_iterator
+       cell = dof_handler.begin_active();
+       cell != dof_handler.end(); ++cell)
+    {
+      const unsigned int fe_index = cell->active_fe_index();
+      local_face_dof_indices.resize(fe_collection[fe_index].dofs_per_face);
+      for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
+        if (std::abs(cell->face(f)->center()[0]) < 0.1)
+          {
+            // deallog << "cell="<<cell<<" face="<<f<<std::endl;
+            fe_face_values_hp.reinit(cell,f);
+            const FEFaceValues<dim> &fe_face_values = fe_face_values_hp.get_present_fe_values();
+
+            const unsigned int n_q_points = fe_face_values.n_quadrature_points;
+            values.resize (n_q_points,
+                           Vector<double>(n_comp));
+
+            fe_face_values.get_function_values (solution,values);
+
+            const std::vector<dealii::Point<dim> > &q_points = fe_face_values.get_quadrature_points();
+
+            for (unsigned int q = 0; q < n_q_points; q++)
+              {
+                // since our face is [0,1]^{dim-1}, the quadrature rule
+                // will coincide with quadrature points on mother face.
+                // Use that to limit output of sub-faces at the same quadrature
+                // points only.
+                Point<dim-1> qpt;
+                for (unsigned int d=0; d < dim-1; d++)
+                  qpt[d] = q_points[q][d+1];
+
+                unsigned int q_found = 0;
+                for (; q_found < quad_formula.size(); q_found++)
+                  if (quad_formula.point(q_found).distance(qpt) < 1e-5)
+                    break;
+
+                if (q_found < quad_formula.size())
+                  {
+                    pairs_point_value.push_back(
+                      std::make_pair(q_points[q],
+                                     values[q]));
+//                      deallog << "@"<<q_points[q]<<" u = {"<<values[q][0];
+//                      for (unsigned int c = 1; c < n_comp; c++)
+//                        deallog << ","<<values[q][c];
+//                      deallog << "}"<<std::endl;
+                  }
+              }
+          }
+    }
+
+  // sort
+  std::sort(pairs_point_value.begin(),
+            pairs_point_value.end(),
+            less_than_key<dim>());
+  for (unsigned int p=0; p < pairs_point_value.size(); p++)
+    {
+      const Point<dim> &pt      = pairs_point_value[p].first;
+      const Vector<double> &val = pairs_point_value[p].second;
+
+      Assert(val.size() == n_comp,
+             ExcInternalError());
+      deallog << "@"<<pt<<" u = {"<<val[0];
+      for (unsigned int c = 1; c < n_comp; c++)
+        deallog << ","<<val[c];
+      deallog << "}"<<std::endl;
+    }
+
+  dof_handler.clear();
+}
+int main (int argc,char **argv)
+{
+  std::ofstream logfile ("output");
+  deallog << std::setprecision(4);
+  deallog << std::fixed;
+  deallog.attach(logfile);
+  deallog.depth_console(0);
+  deallog.threshold_double(1.e-10);
+
+
+  try
+    {
+      {
+        const unsigned int dim = 2;
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+      }
+
+      {
+        const unsigned int dim = 3;
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+        test2cells<dim>(FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(2),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(2),1),
+                        FESystem<dim>(FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1,
+                                      FE_Q<dim>(1),1)
+                       );
+      }
+    }
+  catch (std::exception &exc)
+    {
+      std::cerr << std::endl << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      std::cerr << "Exception on processing: " << std::endl
+                << exc.what() << std::endl
+                << "Aborting!" << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+
+      return 1;
+    }
+  catch (...)
+    {
+      std::cerr << std::endl << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      std::cerr << "Unknown exception!" << std::endl
+                << "Aborting!" << std::endl
+                << "----------------------------------------------------"
+                << std::endl;
+      return 1;
+    };
+}
diff --git a/tests/hp/hp_constraints_neither_dominate_02.output b/tests/hp/hp_constraints_neither_dominate_02.output
new file mode 100644
index 0000000000..66d61ec0e2
--- /dev/null
+++ b/tests/hp/hp_constraints_neither_dominate_02.output
@@ -0,0 +1,869 @@
+
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+    6 24:  0.500000
+    6 51:  0.500000
+    20 30:  0.500000
+    20 33:  0.500000
+    25 23:  0.750000
+    25 50:  0.250000
+    27 23:  0.250000
+    27 33:  0.500000
+    27 50:  0.250000
+    43 29:  0.500000
+    43 32:  0.500000
+    46 22:  0.500000
+    46 49:  0.500000
+    47 23:  0.500000
+    47 50:  0.500000
+    48 24:  0.500000
+    48 51:  0.500000
+    52 22:  0.250000
+    52 49:  0.750000
+    53 22:  0.250000
+    53 32:  0.500000
+    53 49:  0.250000
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.5000 u = {3.7500,7.2500,4.7500}
+DEAL::@0.0000 0.5000 u = {3.7500,7.2500,4.7500}
+DEAL::@0.0000 0.5000 u = {3.7500,7.2500,4.7500}
+DEAL::@0.0000 1.0000 u = {3.0000,6.5000,4.0000}
+DEAL::@0.0000 1.0000 u = {3.0000,6.5000,4.0000}
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(2)-FE_Q<2>(2)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+    6 32:  0.500000
+    6 59:  0.500000
+    7 33:  0.500000
+    7 60:  0.500000
+    8 34:  0.500000
+    8 61:  0.500000
+    30 40:  0.500000
+    30 43:  0.500000
+    35 33:  0.750000
+    35 60:  0.250000
+    37 33:  0.250000
+    37 43:  0.500000
+    37 60:  0.250000
+    53 39:  0.500000
+    53 42:  0.500000
+    56 32:  0.500000
+    56 59:  0.500000
+    57 33:  0.500000
+    57 60:  0.500000
+    58 34:  0.500000
+    58 61:  0.500000
+    62 32:  0.250000
+    62 59:  0.750000
+    63 32:  0.250000
+    63 42:  0.500000
+    63 59:  0.250000
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.5000 u = {10.2500,7.7500,11.2500}
+DEAL::@0.0000 0.5000 u = {10.2500,7.7500,11.2500}
+DEAL::@0.0000 0.5000 u = {10.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 u = {9.5000,7.0000,10.5000}
+DEAL::@0.0000 1.0000 u = {9.5000,7.0000,10.5000}
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(2)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(2)-FE_Q<2>(2)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+    6 32:  0.500000
+    6 54:  0.500000
+    7 33:  0.500000
+    7 55:  0.500000
+    8 34:  0.500000
+    8 56:  0.500000
+    30 40:  0.500000
+    30 43:  0.500000
+    35 33:  0.750000
+    35 55:  0.250000
+    37 33:  0.250000
+    37 43:  0.500000
+    37 55:  0.250000
+    51 32:  0.500000
+    51 54:  0.500000
+    52 33:  0.500000
+    52 55:  0.500000
+    53 34:  0.500000
+    53 56:  0.500000
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.0000 u = {11.0000,8.5000,12.0000}
+DEAL::@0.0000 0.5000 u = {9.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 u = {9.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 u = {9.7500,10.2500,8.5000}
+DEAL::@0.0000 1.0000 u = {8.5000,12.0000,5.0000}
+DEAL::@0.0000 1.0000 u = {8.5000,12.0000,5.0000}
+DEAL::hp-ref: FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(2)-FE_Q<2>(1)-FE_Q<2>(1)] vs FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(2)] with FESystem<2>[FE_Q<2>(1)-FE_Q<2>(1)-FE_Q<2>(1)]
+    6 24:  0.500000
+    6 50:  0.500000
+    25 22:  0.750000
+    25 48:  0.250000
+    45 22:  0.500000
+    45 48:  0.500000
+    46 23:  0.500000
+    46 49:  0.500000
+    47 24:  0.500000
+    47 50:  0.500000
+    51 22:  0.250000
+    51 48:  0.750000
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.0000 u = {4.5000,8.0000,5.5000}
+DEAL::@0.0000 0.5000 u = {7.2500,5.5000,6.0000}
+DEAL::@0.0000 0.5000 u = {7.2500,5.5000,6.0000}
+DEAL::@0.0000 0.5000 u = {7.2500,5.5000,6.0000}
+DEAL::@0.0000 1.0000 u = {10.0000,3.0000,6.5000}
+DEAL::@0.0000 1.0000 u = {10.0000,3.0000,6.5000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+    12 64:  0.500000
+    12 130:  0.500000
+    16 173:  0.500000
+    16 215:  0.500000
+    20 64:  0.500000
+    20 173:  0.500000
+    22 130:  0.500000
+    22 215:  0.500000
+    24 64:  0.250000
+    24 130:  0.250000
+    24 173:  0.250000
+    24 215:  0.250000
+    46 83:  0.500000
+    46 86:  0.500000
+    50 95:  0.500000
+    50 98:  0.500000
+    53 83:  0.500000
+    53 95:  0.500000
+    54 86:  0.500000
+    54 98:  0.500000
+    58 83:  0.250000
+    58 86:  0.250000
+    58 95:  0.250000
+    58 98:  0.250000
+    65 62:  0.500000
+    65 171:  0.500000
+    66 63:  0.500000
+    66 172:  0.500000
+    67 64:  0.500000
+    67 173:  0.500000
+    68 63:  0.750000
+    68 129:  0.250000
+    70 63:  0.250000
+    70 86:  0.500000
+    70 129:  0.250000
+    71 63:  0.375000
+    71 129:  0.125000
+    71 172:  0.375000
+    71 214:  0.125000
+    73 63:  0.125000
+    73 98:  0.500000
+    73 129:  0.125000
+    73 172:  0.125000
+    73 214:  0.125000
+    74 63:  0.750000
+    74 172:  0.250000
+    75 63:  0.375000
+    75 129:  0.375000
+    75 172:  0.125000
+    75 214:  0.125000
+    76 63:  0.562500
+    76 129:  0.187500
+    76 172:  0.187500
+    76 214:  0.0625000
+    78 63:  0.187500
+    78 86:  0.250000
+    78 98:  0.250000
+    78 129:  0.187500
+    78 172:  0.0625000
+    78 214:  0.0625000
+    108 82:  0.500000
+    108 85:  0.500000
+    112 94:  0.500000
+    112 97:  0.500000
+    114 82:  0.500000
+    114 94:  0.500000
+    115 85:  0.500000
+    115 97:  0.500000
+    120 82:  0.250000
+    120 85:  0.250000
+    120 94:  0.250000
+    120 97:  0.250000
+    125 62:  0.500000
+    125 128:  0.500000
+    126 63:  0.500000
+    126 129:  0.500000
+    127 64:  0.500000
+    127 130:  0.500000
+    131 62:  0.250000
+    131 128:  0.250000
+    131 171:  0.250000
+    131 213:  0.250000
+    132 63:  0.250000
+    132 129:  0.250000
+    132 172:  0.250000
+    132 214:  0.250000
+    133 64:  0.250000
+    133 130:  0.250000
+    133 173:  0.250000
+    133 215:  0.250000
+    134 128:  0.500000
+    134 213:  0.500000
+    135 129:  0.500000
+    135 214:  0.500000
+    136 130:  0.500000
+    136 215:  0.500000
+    137 62:  0.250000
+    137 128:  0.750000
+    138 62:  0.250000
+    138 85:  0.500000
+    138 128:  0.250000
+    140 62:  0.125000
+    140 128:  0.375000
+    140 171:  0.125000
+    140 213:  0.375000
+    141 62:  0.125000
+    141 97:  0.500000
+    141 128:  0.125000
+    141 171:  0.125000
+    141 213:  0.125000
+    143 62:  0.375000
+    143 128:  0.375000
+    143 171:  0.125000
+    143 213:  0.125000
+    144 128:  0.750000
+    144 213:  0.250000
+    145 62:  0.187500
+    145 128:  0.562500
+    145 171:  0.0625000
+    145 213:  0.187500
+    146 62:  0.187500
+    146 85:  0.250000
+    146 97:  0.250000
+    146 128:  0.187500
+    146 171:  0.0625000
+    146 213:  0.0625000
+    160 185:  0.500000
+    160 188:  0.500000
+    163 95:  0.500000
+    163 185:  0.500000
+    164 98:  0.500000
+    164 188:  0.500000
+    168 95:  0.250000
+    168 98:  0.250000
+    168 185:  0.250000
+    168 188:  0.250000
+    174 172:  0.750000
+    174 214:  0.250000
+    176 172:  0.250000
+    176 188:  0.500000
+    176 214:  0.250000
+    177 63:  0.250000
+    177 172:  0.750000
+    178 63:  0.125000
+    178 129:  0.125000
+    178 172:  0.375000
+    178 214:  0.375000
+    179 63:  0.187500
+    179 129:  0.0625000
+    179 172:  0.562500
+    179 214:  0.187500
+    181 63:  0.0625000
+    181 98:  0.250000
+    181 129:  0.0625000
+    181 172:  0.187500
+    181 188:  0.250000
+    181 214:  0.187500
+    198 184:  0.500000
+    198 187:  0.500000
+    200 94:  0.500000
+    200 184:  0.500000
+    201 97:  0.500000
+    201 187:  0.500000
+    206 94:  0.250000
+    206 97:  0.250000
+    206 184:  0.250000
+    206 187:  0.250000
+    210 171:  0.500000
+    210 213:  0.500000
+    211 172:  0.500000
+    211 214:  0.500000
+    212 173:  0.500000
+    212 215:  0.500000
+    216 171:  0.250000
+    216 213:  0.750000
+    217 171:  0.250000
+    217 187:  0.500000
+    217 213:  0.250000
+    219 62:  0.125000
+    219 128:  0.125000
+    219 171:  0.375000
+    219 213:  0.375000
+    220 128:  0.250000
+    220 213:  0.750000
+    221 62:  0.0625000
+    221 128:  0.187500
+    221 171:  0.187500
+    221 213:  0.562500
+    222 62:  0.0625000
+    222 97:  0.250000
+    222 128:  0.0625000
+    222 171:  0.187500
+    222 187:  0.250000
+    222 213:  0.187500
+DEAL::@0.0000 0.0000 0.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.0000 0.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.2500,5.5000,9.0000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.2500,5.5000,9.0000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.2500,5.5000,9.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {11.5000,3.7500,7.2500}
+DEAL::@0.0000 0.5000 0.0000 u = {11.5000,3.7500,7.2500}
+DEAL::@0.0000 0.5000 0.0000 u = {11.5000,3.7500,7.2500}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 0.5000 u = {8.7500,6.6250,10.1250}
+DEAL::@0.0000 0.5000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {8.5000,6.0000,9.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {8.5000,6.0000,9.5000}
+DEAL::@0.0000 1.0000 0.5000 u = {7.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {7.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {7.2500,7.7500,11.2500}
+DEAL::@0.0000 1.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::@0.0000 1.0000 1.0000 u = {6.0000,9.5000,13.0000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(2)-FE_Q<3>(2)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+    12 100:  0.500000
+    12 166:  0.500000
+    13 101:  0.500000
+    13 167:  0.500000
+    14 102:  0.500000
+    14 168:  0.500000
+    24 209:  0.500000
+    24 251:  0.500000
+    25 210:  0.500000
+    25 252:  0.500000
+    26 211:  0.500000
+    26 253:  0.500000
+    36 100:  0.500000
+    36 209:  0.500000
+    37 101:  0.500000
+    37 210:  0.500000
+    38 102:  0.500000
+    38 211:  0.500000
+    42 166:  0.500000
+    42 251:  0.500000
+    43 167:  0.500000
+    43 252:  0.500000
+    44 168:  0.500000
+    44 253:  0.500000
+    48 100:  0.250000
+    48 166:  0.250000
+    48 209:  0.250000
+    48 251:  0.250000
+    49 101:  0.250000
+    49 167:  0.250000
+    49 210:  0.250000
+    49 252:  0.250000
+    50 102:  0.250000
+    50 168:  0.250000
+    50 211:  0.250000
+    50 253:  0.250000
+    84 121:  0.500000
+    84 124:  0.500000
+    88 133:  0.500000
+    88 136:  0.500000
+    91 121:  0.500000
+    91 133:  0.500000
+    92 124:  0.500000
+    92 136:  0.500000
+    96 121:  0.250000
+    96 124:  0.250000
+    96 133:  0.250000
+    96 136:  0.250000
+    103 100:  0.500000
+    103 209:  0.500000
+    104 101:  0.500000
+    104 210:  0.500000
+    105 102:  0.500000
+    105 211:  0.500000
+    106 101:  0.750000
+    106 167:  0.250000
+    108 101:  0.250000
+    108 124:  0.500000
+    108 167:  0.250000
+    109 101:  0.375000
+    109 167:  0.125000
+    109 210:  0.375000
+    109 252:  0.125000
+    111 101:  0.125000
+    111 136:  0.500000
+    111 167:  0.125000
+    111 210:  0.125000
+    111 252:  0.125000
+    112 101:  0.750000
+    112 210:  0.250000
+    113 101:  0.375000
+    113 167:  0.375000
+    113 210:  0.125000
+    113 252:  0.125000
+    114 101:  0.562500
+    114 167:  0.187500
+    114 210:  0.187500
+    114 252:  0.0625000
+    116 101:  0.187500
+    116 124:  0.250000
+    116 136:  0.250000
+    116 167:  0.187500
+    116 210:  0.0625000
+    116 252:  0.0625000
+    146 120:  0.500000
+    146 123:  0.500000
+    150 132:  0.500000
+    150 135:  0.500000
+    152 120:  0.500000
+    152 132:  0.500000
+    153 123:  0.500000
+    153 135:  0.500000
+    158 120:  0.250000
+    158 123:  0.250000
+    158 132:  0.250000
+    158 135:  0.250000
+    163 100:  0.500000
+    163 166:  0.500000
+    164 101:  0.500000
+    164 167:  0.500000
+    165 102:  0.500000
+    165 168:  0.500000
+    169 100:  0.250000
+    169 166:  0.250000
+    169 209:  0.250000
+    169 251:  0.250000
+    170 101:  0.250000
+    170 167:  0.250000
+    170 210:  0.250000
+    170 252:  0.250000
+    171 102:  0.250000
+    171 168:  0.250000
+    171 211:  0.250000
+    171 253:  0.250000
+    172 166:  0.500000
+    172 251:  0.500000
+    173 167:  0.500000
+    173 252:  0.500000
+    174 168:  0.500000
+    174 253:  0.500000
+    175 100:  0.250000
+    175 166:  0.750000
+    176 100:  0.250000
+    176 123:  0.500000
+    176 166:  0.250000
+    178 100:  0.125000
+    178 166:  0.375000
+    178 209:  0.125000
+    178 251:  0.375000
+    179 100:  0.125000
+    179 135:  0.500000
+    179 166:  0.125000
+    179 209:  0.125000
+    179 251:  0.125000
+    181 100:  0.375000
+    181 166:  0.375000
+    181 209:  0.125000
+    181 251:  0.125000
+    182 166:  0.750000
+    182 251:  0.250000
+    183 100:  0.187500
+    183 166:  0.562500
+    183 209:  0.0625000
+    183 251:  0.187500
+    184 100:  0.187500
+    184 123:  0.250000
+    184 135:  0.250000
+    184 166:  0.187500
+    184 209:  0.0625000
+    184 251:  0.0625000
+    198 223:  0.500000
+    198 226:  0.500000
+    201 133:  0.500000
+    201 223:  0.500000
+    202 136:  0.500000
+    202 226:  0.500000
+    206 133:  0.250000
+    206 136:  0.250000
+    206 223:  0.250000
+    206 226:  0.250000
+    212 210:  0.750000
+    212 252:  0.250000
+    214 210:  0.250000
+    214 226:  0.500000
+    214 252:  0.250000
+    215 101:  0.250000
+    215 210:  0.750000
+    216 101:  0.125000
+    216 167:  0.125000
+    216 210:  0.375000
+    216 252:  0.375000
+    217 101:  0.187500
+    217 167:  0.0625000
+    217 210:  0.562500
+    217 252:  0.187500
+    219 101:  0.0625000
+    219 136:  0.250000
+    219 167:  0.0625000
+    219 210:  0.187500
+    219 226:  0.250000
+    219 252:  0.187500
+    236 222:  0.500000
+    236 225:  0.500000
+    238 132:  0.500000
+    238 222:  0.500000
+    239 135:  0.500000
+    239 225:  0.500000
+    244 132:  0.250000
+    244 135:  0.250000
+    244 222:  0.250000
+    244 225:  0.250000
+    248 209:  0.500000
+    248 251:  0.500000
+    249 210:  0.500000
+    249 252:  0.500000
+    250 211:  0.500000
+    250 253:  0.500000
+    254 209:  0.250000
+    254 251:  0.750000
+    255 209:  0.250000
+    255 225:  0.500000
+    255 251:  0.250000
+    257 100:  0.125000
+    257 166:  0.125000
+    257 209:  0.375000
+    257 251:  0.375000
+    258 166:  0.250000
+    258 251:  0.750000
+    259 100:  0.0625000
+    259 166:  0.187500
+    259 209:  0.187500
+    259 251:  0.562500
+    260 100:  0.0625000
+    260 135:  0.250000
+    260 166:  0.0625000
+    260 209:  0.187500
+    260 225:  0.250000
+    260 251:  0.187500
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {10.5000,5.7500,6.2500}
+DEAL::@0.0000 0.0000 0.5000 u = {10.5000,5.7500,6.2500}
+DEAL::@0.0000 0.0000 0.5000 u = {10.5000,5.7500,6.2500}
+DEAL::@0.0000 0.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {8.7500,12.2500,4.5000}
+DEAL::@0.0000 0.5000 0.0000 u = {8.7500,12.2500,4.5000}
+DEAL::@0.0000 0.5000 0.0000 u = {8.7500,12.2500,4.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {11.6250,6.8750,4.7500}
+DEAL::@0.0000 0.5000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {11.0000,14.5000,1.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {11.0000,14.5000,1.5000}
+DEAL::@0.0000 1.0000 0.5000 u = {12.7500,8.0000,3.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {12.7500,8.0000,3.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {12.7500,8.0000,3.2500}
+DEAL::@0.0000 1.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::@0.0000 1.0000 1.0000 u = {14.5000,1.5000,5.0000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(2)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(2)-FE_Q<3>(2)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+    12 100:  0.500000
+    12 147:  0.500000
+    13 101:  0.500000
+    13 148:  0.500000
+    14 102:  0.500000
+    14 149:  0.500000
+    24 176:  0.500000
+    24 204:  0.500000
+    25 177:  0.500000
+    25 205:  0.500000
+    26 178:  0.500000
+    26 206:  0.500000
+    36 100:  0.500000
+    36 176:  0.500000
+    37 101:  0.500000
+    37 177:  0.500000
+    38 102:  0.500000
+    38 178:  0.500000
+    42 147:  0.500000
+    42 204:  0.500000
+    43 148:  0.500000
+    43 205:  0.500000
+    44 149:  0.500000
+    44 206:  0.500000
+    48 100:  0.250000
+    48 147:  0.250000
+    48 176:  0.250000
+    48 204:  0.250000
+    49 101:  0.250000
+    49 148:  0.250000
+    49 177:  0.250000
+    49 205:  0.250000
+    50 102:  0.250000
+    50 149:  0.250000
+    50 178:  0.250000
+    50 206:  0.250000
+    84 121:  0.500000
+    84 124:  0.500000
+    88 133:  0.500000
+    88 136:  0.500000
+    91 121:  0.500000
+    91 133:  0.500000
+    92 124:  0.500000
+    92 136:  0.500000
+    96 121:  0.250000
+    96 124:  0.250000
+    96 133:  0.250000
+    96 136:  0.250000
+    103 100:  0.500000
+    103 176:  0.500000
+    104 101:  0.500000
+    104 177:  0.500000
+    105 102:  0.500000
+    105 178:  0.500000
+    106 101:  0.750000
+    106 148:  0.250000
+    108 101:  0.250000
+    108 124:  0.500000
+    108 148:  0.250000
+    109 101:  0.375000
+    109 148:  0.125000
+    109 177:  0.375000
+    109 205:  0.125000
+    111 101:  0.125000
+    111 136:  0.500000
+    111 148:  0.125000
+    111 177:  0.125000
+    111 205:  0.125000
+    112 101:  0.750000
+    112 177:  0.250000
+    113 101:  0.375000
+    113 148:  0.375000
+    113 177:  0.125000
+    113 205:  0.125000
+    114 101:  0.562500
+    114 148:  0.187500
+    114 177:  0.187500
+    114 205:  0.0625000
+    116 101:  0.187500
+    116 124:  0.250000
+    116 136:  0.250000
+    116 148:  0.187500
+    116 177:  0.0625000
+    116 205:  0.0625000
+    144 100:  0.500000
+    144 147:  0.500000
+    145 101:  0.500000
+    145 148:  0.500000
+    146 102:  0.500000
+    146 149:  0.500000
+    150 100:  0.250000
+    150 147:  0.250000
+    150 176:  0.250000
+    150 204:  0.250000
+    151 101:  0.250000
+    151 148:  0.250000
+    151 177:  0.250000
+    151 205:  0.250000
+    152 102:  0.250000
+    152 149:  0.250000
+    152 178:  0.250000
+    152 206:  0.250000
+    153 147:  0.500000
+    153 204:  0.500000
+    154 148:  0.500000
+    154 205:  0.500000
+    155 149:  0.500000
+    155 206:  0.500000
+    165 190:  0.500000
+    165 193:  0.500000
+    168 133:  0.500000
+    168 190:  0.500000
+    169 136:  0.500000
+    169 193:  0.500000
+    173 133:  0.250000
+    173 136:  0.250000
+    173 190:  0.250000
+    173 193:  0.250000
+    179 177:  0.750000
+    179 205:  0.250000
+    181 177:  0.250000
+    181 193:  0.500000
+    181 205:  0.250000
+    182 101:  0.250000
+    182 177:  0.750000
+    183 101:  0.125000
+    183 148:  0.125000
+    183 177:  0.375000
+    183 205:  0.375000
+    184 101:  0.187500
+    184 148:  0.0625000
+    184 177:  0.562500
+    184 205:  0.187500
+    186 101:  0.0625000
+    186 136:  0.250000
+    186 148:  0.0625000
+    186 177:  0.187500
+    186 193:  0.250000
+    186 205:  0.187500
+    201 176:  0.500000
+    201 204:  0.500000
+    202 177:  0.500000
+    202 205:  0.500000
+    203 178:  0.500000
+    203 206:  0.500000
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {6.7500,7.2500,7.7500}
+DEAL::@0.0000 0.0000 0.5000 u = {6.7500,7.2500,7.7500}
+DEAL::@0.0000 0.0000 0.5000 u = {6.7500,7.2500,7.7500}
+DEAL::@0.0000 0.0000 1.0000 u = {7.0000,4.5000,8.0000}
+DEAL::@0.0000 0.0000 1.0000 u = {7.0000,4.5000,8.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.0000,7.5000,8.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.0000,7.5000,8.0000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.0000,7.5000,8.0000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 0.5000 u = {4.5000,6.5000,8.5000}
+DEAL::@0.0000 0.5000 1.0000 u = {5.0000,5.5000,9.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {5.0000,5.5000,9.0000}
+DEAL::@0.0000 0.5000 1.0000 u = {5.0000,5.5000,9.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {1.5000,5.0000,8.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {1.5000,5.0000,8.5000}
+DEAL::@0.0000 1.0000 0.5000 u = {2.2500,5.7500,9.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {2.2500,5.7500,9.2500}
+DEAL::@0.0000 1.0000 0.5000 u = {2.2500,5.7500,9.2500}
+DEAL::@0.0000 1.0000 1.0000 u = {3.0000,6.5000,10.0000}
+DEAL::@0.0000 1.0000 1.0000 u = {3.0000,6.5000,10.0000}
+DEAL::hp-ref: FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(2)-FE_Q<3>(1)-FE_Q<3>(1)] vs FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(2)] with FESystem<3>[FE_Q<3>(1)-FE_Q<3>(1)-FE_Q<3>(1)]
+    12 60:  0.500000
+    12 122:  0.500000
+    16 161:  0.500000
+    16 200:  0.500000
+    20 60:  0.500000
+    20 161:  0.500000
+    22 122:  0.500000
+    22 200:  0.500000
+    24 60:  0.250000
+    24 122:  0.250000
+    24 161:  0.250000
+    24 200:  0.250000
+    61 58:  0.500000
+    61 159:  0.500000
+    62 59:  0.500000
+    62 160:  0.500000
+    63 60:  0.500000
+    63 161:  0.500000
+    64 58:  0.750000
+    64 120:  0.250000
+    66 58:  0.375000
+    66 120:  0.125000
+    66 159:  0.375000
+    66 198:  0.125000
+    68 58:  0.750000
+    68 159:  0.250000
+    69 58:  0.562500
+    69 120:  0.187500
+    69 159:  0.187500
+    69 198:  0.0625000
+    117 58:  0.500000
+    117 120:  0.500000
+    118 59:  0.500000
+    118 121:  0.500000
+    119 60:  0.500000
+    119 122:  0.500000
+    123 58:  0.250000
+    123 120:  0.250000
+    123 159:  0.250000
+    123 198:  0.250000
+    124 59:  0.250000
+    124 121:  0.250000
+    124 160:  0.250000
+    124 199:  0.250000
+    125 60:  0.250000
+    125 122:  0.250000
+    125 161:  0.250000
+    125 200:  0.250000
+    126 120:  0.500000
+    126 198:  0.500000
+    127 121:  0.500000
+    127 199:  0.500000
+    128 122:  0.500000
+    128 200:  0.500000
+    129 58:  0.250000
+    129 120:  0.750000
+    132 58:  0.125000
+    132 120:  0.375000
+    132 159:  0.125000
+    132 198:  0.375000
+    135 58:  0.375000
+    135 120:  0.375000
+    135 159:  0.125000
+    135 198:  0.125000
+    136 120:  0.750000
+    136 198:  0.250000
+    137 58:  0.187500
+    137 120:  0.562500
+    137 159:  0.0625000
+    137 198:  0.187500
+    162 159:  0.750000
+    162 198:  0.250000
+    164 58:  0.250000
+    164 159:  0.750000
+    165 58:  0.187500
+    165 120:  0.0625000
+    165 159:  0.562500
+    165 198:  0.187500
+    195 159:  0.500000
+    195 198:  0.500000
+    196 160:  0.500000
+    196 199:  0.500000
+    197 161:  0.500000
+    197 200:  0.500000
+    201 159:  0.250000
+    201 198:  0.750000
+    204 58:  0.125000
+    204 120:  0.125000
+    204 159:  0.375000
+    204 198:  0.375000
+    205 120:  0.250000
+    205 198:  0.750000
+    206 58:  0.0625000
+    206 120:  0.187500
+    206 159:  0.187500
+    206 198:  0.562500
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.0000 u = {6.5000,10.0000,7.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {7.7500,11.2500,6.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {7.7500,11.2500,6.5000}
+DEAL::@0.0000 0.0000 0.5000 u = {7.7500,11.2500,6.5000}
+DEAL::@0.0000 0.0000 1.0000 u = {9.0000,12.5000,5.5000}
+DEAL::@0.0000 0.0000 1.0000 u = {9.0000,12.5000,5.5000}
+DEAL::@0.0000 0.5000 0.0000 u = {4.7500,8.2500,8.7500}
+DEAL::@0.0000 0.5000 0.0000 u = {4.7500,8.2500,8.7500}
+DEAL::@0.0000 0.5000 0.0000 u = {4.7500,8.2500,8.7500}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 0.5000 u = {5.7500,9.2500,8.6250}
+DEAL::@0.0000 0.5000 1.0000 u = {6.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.7500,10.2500,8.5000}
+DEAL::@0.0000 0.5000 1.0000 u = {6.7500,10.2500,8.5000}
+DEAL::@0.0000 1.0000 0.0000 u = {3.0000,6.5000,10.0000}
+DEAL::@0.0000 1.0000 0.0000 u = {3.0000,6.5000,10.0000}
+DEAL::@0.0000 1.0000 0.5000 u = {3.7500,7.2500,10.7500}
+DEAL::@0.0000 1.0000 0.5000 u = {3.7500,7.2500,10.7500}
+DEAL::@0.0000 1.0000 0.5000 u = {3.7500,7.2500,10.7500}
+DEAL::@0.0000 1.0000 1.0000 u = {4.5000,8.0000,11.5000}
+DEAL::@0.0000 1.0000 1.0000 u = {4.5000,8.0000,11.5000}