]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Add new test case
authorMartin Kronbichler <martin.kronbichler@rub.de>
Wed, 20 Mar 2024 15:23:56 +0000 (16:23 +0100)
committerMartin Kronbichler <martin.kronbichler@rub.de>
Thu, 21 Mar 2024 13:50:36 +0000 (14:50 +0100)
tests/matrix_free/point_evaluation_30.cc [new file with mode: 0644]
tests/matrix_free/point_evaluation_30.output [new file with mode: 0644]
tests/matrix_free/point_evaluation_31.cc [new file with mode: 0644]
tests/matrix_free/point_evaluation_31.output [new file with mode: 0644]
tests/matrix_free/point_evaluation_32.cc [new file with mode: 0644]
tests/matrix_free/point_evaluation_32.output [new file with mode: 0644]

diff --git a/tests/matrix_free/point_evaluation_30.cc b/tests/matrix_free/point_evaluation_30.cc
new file mode 100644 (file)
index 0000000..36420fc
--- /dev/null
@@ -0,0 +1,152 @@
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2024 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// Part of the source code is dual licensed under Apache-2.0 WITH
+// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
+// governing the source code and code contributions can be found in
+// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
+//
+// ------------------------------------------------------------------------
+
+
+// check FEPointEvaluation for dim != spacedim, otherwise same test as
+// point_evaluation_01
+
+#include <deal.II/base/function_lib.h>
+
+#include <deal.II/dofs/dof_handler.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/mapping_q.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/lac/vector.h>
+
+#include <deal.II/matrix_free/fe_point_evaluation.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <iostream>
+
+#include "../tests.h"
+
+
+
+template <int dim, int spacedim, typename Number = double>
+void
+test(const unsigned int degree)
+{
+  Triangulation<dim, spacedim> tria;
+
+  if constexpr (dim + 1 == spacedim)
+    GridGenerator::hyper_sphere(tria, Point<spacedim>(), 1.);
+  else if constexpr (dim == 1)
+    GridGenerator::hyper_rectangle(tria, Point<dim>(), Point<dim>(1.2));
+  else
+    static_assert("Case not implemented");
+
+  MappingQ<dim, spacedim> mapping(degree);
+  deallog << "Mapping of degree " << degree << std::endl;
+
+  std::vector<Point<dim>> unit_points;
+  for (unsigned int i = 0; i < 13; ++i)
+    {
+      Point<dim> p;
+      for (unsigned int d = 0; d < dim; ++d)
+        p[d] = static_cast<Number>(i) / 17. + 0.015625 * d;
+      unit_points.push_back(p);
+    }
+
+  FE_Q<dim, spacedim>     fe(degree);
+  FEValues<dim, spacedim> fe_values(mapping,
+                                    fe,
+                                    Quadrature<dim>(unit_points),
+                                    update_values | update_gradients);
+
+  DoFHandler<dim, spacedim> dof_handler(tria);
+  dof_handler.distribute_dofs(fe);
+  Vector<Number> vector(dof_handler.n_dofs());
+
+  FEPointEvaluation<1, dim, spacedim, Number> evaluator(
+    mapping, fe, update_values | update_gradients);
+
+  Tensor<1, spacedim, Number> exponents;
+  exponents[0] = 1.;
+  VectorTools::interpolate(mapping,
+                           dof_handler,
+                           Functions::Monomial<spacedim, Number>(exponents),
+                           vector);
+
+  std::vector<Number>                      solution_values(fe.dofs_per_cell);
+  std::vector<Number>                      function_values(unit_points.size());
+  std::vector<Tensor<1, spacedim, Number>> function_gradients(
+    unit_points.size());
+
+  // For float numbers that are sensitive to roundoff in the numdiff
+  // tolerances (absolute 1e-8), we multiply by 1e-3 to ensure that the test
+  // remains robust
+  const double factor_float = std::is_same_v<Number, float> ? 0.001 : 1.;
+
+  for (const auto &cell : dof_handler.active_cell_iterators())
+    {
+      fe_values.reinit(cell);
+      fe_values.get_function_values(vector, function_values);
+      fe_values.get_function_gradients(vector, function_gradients);
+
+      cell->get_dof_values(vector,
+                           solution_values.begin(),
+                           solution_values.end());
+
+      evaluator.reinit(cell, unit_points);
+      evaluator.evaluate(solution_values,
+                         EvaluationFlags::values | EvaluationFlags::gradients);
+
+      deallog << "Cell with center " << cell->center(true) << std::endl;
+      for (unsigned int i = 0; i < function_values.size(); ++i)
+        deallog << mapping.transform_unit_to_real_cell(cell, unit_points[i])
+                << ": " << factor_float * evaluator.get_value(i)
+                << " error value "
+                << factor_float * (function_values[i] - evaluator.get_value(i))
+                << " error grad "
+                << factor_float *
+                     (evaluator.get_gradient(i) - function_gradients[i]).norm()
+                << std::endl;
+      deallog << std::endl;
+
+      for (unsigned int i = 0; i < unit_points.size(); ++i)
+        {
+          evaluator.submit_value(evaluator.get_value(i), i);
+          evaluator.submit_gradient(evaluator.get_gradient(i), i);
+        }
+
+      evaluator.test_and_sum(solution_values,
+                             EvaluationFlags::values |
+                               EvaluationFlags::gradients);
+
+      for (const auto i : solution_values)
+        deallog << factor_float * i << ' ';
+      deallog << std::endl;
+    }
+}
+
+
+
+int
+main()
+{
+  initlog();
+  deallog << std::setprecision(10);
+
+  test<1, 2>(3);
+  test<2, 3>(2);
+  test<1, 3>(3);
+
+  test<2, 3, float>(5);
+}
diff --git a/tests/matrix_free/point_evaluation_30.output b/tests/matrix_free/point_evaluation_30.output
new file mode 100644 (file)
index 0000000..4246ff2
--- /dev/null
@@ -0,0 +1,277 @@
+
+DEAL::Mapping of degree 3
+DEAL::Cell with center 1.110223025e-16 -1.000000000
+DEAL::0.7071067812 -0.7071067812: 0.7071067812 error value 0.000000000 error grad 0.000000000
+DEAL::0.6391202566 -0.7712838782: 0.6391202566 error value 0.000000000 error grad 0.000000000
+DEAL::0.5654327913 -0.8274388380: 0.5654327913 error value 0.000000000 error grad 0.000000000
+DEAL::0.4868045108 -0.8755716608: 0.4868045108 error value 0.000000000 error grad 0.000000000
+DEAL::0.4039955404 -0.9156823464: 0.4039955404 error value 0.000000000 error grad 0.000000000
+DEAL::0.3177660056 -0.9477708949: 0.3177660056 error value 0.000000000 error grad 0.000000000
+DEAL::0.2288760319 -0.9718373063: 0.2288760319 error value 0.000000000 error grad 0.000000000
+DEAL::0.1380857445 -0.9878815805: 0.1380857445 error value 0.000000000 error grad 0.000000000
+DEAL::0.04615526908 -0.9959037176: 0.04615526908 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04615526908 -0.9959037176: -0.04615526908 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1380857445 -0.9878815805: -0.1380857445 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2288760319 -0.9718373063: -0.2288760319 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3177660056 -0.9477708949: -0.3177660056 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::11.45242453 -0.1762428574 2.726194254 -11.19991605 
+DEAL::Cell with center -1.000000000 -1.110223025e-16
+DEAL::-0.7071067812 -0.7071067812: -0.7071067812 error value 0.000000000 error grad 0.000000000
+DEAL::-0.7712838782 -0.6391202566: -0.7712838782 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8274388380 -0.5654327913: -0.8274388380 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8755716608 -0.4868045108: -0.8755716608 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9156823464 -0.4039955404: -0.9156823464 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9477708949 -0.3177660056: -0.9477708949 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9718373063 -0.2288760319: -0.9718373063 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9878815805 -0.1380857445: -0.9878815805 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9959037176 -0.04615526908: -0.9959037176 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9959037176 0.04615526908: -0.9959037176 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9878815805 0.1380857445: -0.9878815805 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9718373063 0.2288760319: -0.9718373063 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9477708949 0.3177660056: -0.9477708949 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::6.989633632 0.7389983457 -16.26413810 -3.368364381 
+DEAL::Cell with center 1.000000000 1.110223025e-16
+DEAL::0.7071067812 0.7071067812: 0.7071067812 error value 0.000000000 error grad 0.000000000
+DEAL::0.7712838782 0.6391202566: 0.7712838782 error value 0.000000000 error grad 0.000000000
+DEAL::0.8274388380 0.5654327913: 0.8274388380 error value 0.000000000 error grad 0.000000000
+DEAL::0.8755716608 0.4868045108: 0.8755716608 error value 0.000000000 error grad 0.000000000
+DEAL::0.9156823464 0.4039955404: 0.9156823464 error value 0.000000000 error grad 0.000000000
+DEAL::0.9477708949 0.3177660056: 0.9477708949 error value 0.000000000 error grad 0.000000000
+DEAL::0.9718373063 0.2288760319: 0.9718373063 error value 0.000000000 error grad 0.000000000
+DEAL::0.9878815805 0.1380857445: 0.9878815805 error value 0.000000000 error grad 0.000000000
+DEAL::0.9959037176 0.04615526908: 0.9959037176 error value 0.000000000 error grad 0.000000000
+DEAL::0.9959037176 -0.04615526908: 0.9959037176 error value 0.000000000 error grad 0.000000000
+DEAL::0.9878815805 -0.1380857445: 0.9878815805 error value 0.000000000 error grad 0.000000000
+DEAL::0.9718373063 -0.2288760319: 0.9718373063 error value 0.000000000 error grad 0.000000000
+DEAL::0.9477708949 -0.3177660056: 0.9477708949 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-6.989633632 -0.7389983457 16.26413810 3.368364381 
+DEAL::Cell with center -1.110223025e-16 1.000000000
+DEAL::-0.7071067812 0.7071067812: -0.7071067812 error value 0.000000000 error grad 0.000000000
+DEAL::-0.6391202566 0.7712838782: -0.6391202566 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5654327913 0.8274388380: -0.5654327913 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4868045108 0.8755716608: -0.4868045108 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4039955404 0.9156823464: -0.4039955404 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3177660056 0.9477708949: -0.3177660056 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2288760319 0.9718373063: -0.2288760319 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1380857445 0.9878815805: -0.1380857445 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04615526908 0.9959037176: -0.04615526908 error value 0.000000000 error grad 0.000000000
+DEAL::0.04615526908 0.9959037176: 0.04615526908 error value 0.000000000 error grad 0.000000000
+DEAL::0.1380857445 0.9878815805: 0.1380857445 error value 0.000000000 error grad 0.000000000
+DEAL::0.2288760319 0.9718373063: 0.2288760319 error value 0.000000000 error grad 0.000000000
+DEAL::0.3177660056 0.9477708949: 0.3177660056 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-11.45242453 0.1762428574 -2.726194254 11.19991605 
+DEAL::Mapping of degree 2
+DEAL::Cell with center 0.000000000 0.000000000 -1.000000000
+DEAL::-0.5593080733 -0.5853333358 -0.5853333358: -0.5593080733 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5158409604 -0.5409831617 -0.6518067606: -0.5158409604 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4629781872 -0.4873495734 -0.7224794539: -0.4629781872 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4019874746 -0.4257002918 -0.7922007171: -0.4019874746 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3341365434 -0.3573030376 -0.8565698965: -0.3341365434 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2606931144 -0.2834255319 -0.9119363833: -0.2606931144 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1829249085 -0.2053354953 -0.9553996133: -0.1829249085 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1020996466 -0.1243006489 -0.9848090674: -0.1020996466 error value 0.000000000 error grad 0.000000000
+DEAL::-0.01948504954 -0.04158871337 -0.9987642713: -0.01948504954 error value 0.000000000 error grad 0.000000000
+DEAL::0.06365116184 0.04153259032 -0.9966147953: 0.06365116184 error value 0.000000000 error grad 0.000000000
+DEAL::0.1460412667 0.1237955413 -0.9784602549: 0.1460412667 error value 0.000000000 error grad 0.000000000
+DEAL::0.2264175441 0.2039324188 -0.9451503101: 0.2264175441 error value 0.000000000 error grad 0.000000000
+DEAL::0.3035122733 0.2806755020 -0.8982846661: 0.3035122733 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-7.395100711 0.8334063071 -1.342717800 0.6437129783 7.114003954 -1.223724496 -8.167443670 5.414419780 2.023611946 
+DEAL::Cell with center 1.000000000 0.000000000 0.000000000
+DEAL::0.5853333358 -0.5853333358 -0.5593080733: 0.5853333358 error value 0.000000000 error grad 0.000000000
+DEAL::0.6518067606 -0.5409831617 -0.5158409604: 0.6518067606 error value 0.000000000 error grad 0.000000000
+DEAL::0.7224794539 -0.4873495734 -0.4629781872: 0.7224794539 error value 0.000000000 error grad 0.000000000
+DEAL::0.7922007171 -0.4257002918 -0.4019874746: 0.7922007171 error value 0.000000000 error grad 0.000000000
+DEAL::0.8565698965 -0.3573030376 -0.3341365434: 0.8565698965 error value 0.000000000 error grad 0.000000000
+DEAL::0.9119363833 -0.2834255319 -0.2606931144: 0.9119363833 error value 0.000000000 error grad 0.000000000
+DEAL::0.9553996133 -0.2053354953 -0.1829249085: 0.9553996133 error value 0.000000000 error grad 0.000000000
+DEAL::0.9848090674 -0.1243006489 -0.1020996466: 0.9848090674 error value 0.000000000 error grad 0.000000000
+DEAL::0.9987642713 -0.04158871337 -0.01948504954: 0.9987642713 error value 0.000000000 error grad 0.000000000
+DEAL::0.9966147953 0.04153259032 0.06365116184: 0.9966147953 error value 0.000000000 error grad 0.000000000
+DEAL::0.9784602549 0.1237955413 0.1460412667: 0.9784602549 error value 0.000000000 error grad 0.000000000
+DEAL::0.9451503101 0.2039324188 0.2264175441: 0.9451503101 error value 0.000000000 error grad 0.000000000
+DEAL::0.8982846661 0.2806755020 0.3035122733: 0.8982846661 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-5.832306001 -0.7080816281 -0.6117074595 -0.1384243340 3.525869726 -0.9668013027 3.683918369 -0.9730667899 13.29840895 
+DEAL::Cell with center 0.000000000 0.000000000 1.000000000
+DEAL::-0.5853333358 -0.5593080733 0.5853333358: -0.5853333358 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5409831617 -0.5158409604 0.6518067606: -0.5409831617 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4873495734 -0.4629781872 0.7224794539: -0.4873495734 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4257002918 -0.4019874746 0.7922007171: -0.4257002918 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3573030376 -0.3341365434 0.8565698965: -0.3573030376 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2834255319 -0.2606931144 0.9119363833: -0.2834255319 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2053354953 -0.1829249085 0.9553996133: -0.2053354953 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1243006489 -0.1020996466 0.9848090674: -0.1243006489 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04158871337 -0.01948504954 0.9987642713: -0.04158871337 error value 0.000000000 error grad 0.000000000
+DEAL::0.04153259032 0.06365116184 0.9966147953: 0.04153259032 error value 0.000000000 error grad 0.000000000
+DEAL::0.1237955413 0.1460412667 0.9784602549: 0.1237955413 error value 0.000000000 error grad 0.000000000
+DEAL::0.2039324188 0.2264175441 0.9451503101: 0.2039324188 error value 0.000000000 error grad 0.000000000
+DEAL::0.2806755020 0.3035122733 0.8982846661: 0.2806755020 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-6.890811517 -1.452971756 0.8592313911 0.7857593783 -9.106658205 4.919866875 6.809982664 -1.343083854 3.017301286 
+DEAL::Cell with center -1.000000000 0.000000000 0.000000000
+DEAL::-0.5853333358 -0.5593080733 -0.5853333358: -0.5853333358 error value 0.000000000 error grad 0.000000000
+DEAL::-0.6518067606 -0.5158409604 -0.5409831617: -0.6518067606 error value 0.000000000 error grad 0.000000000
+DEAL::-0.7224794539 -0.4629781872 -0.4873495734: -0.7224794539 error value 0.000000000 error grad 0.000000000
+DEAL::-0.7922007171 -0.4019874746 -0.4257002918: -0.7922007171 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8565698965 -0.3341365434 -0.3573030376: -0.8565698965 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9119363833 -0.2606931144 -0.2834255319: -0.9119363833 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9553996133 -0.1829249085 -0.2053354953: -0.9553996133 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9848090674 -0.1020996466 -0.1243006489: -0.9848090674 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9987642713 -0.01948504954 -0.04158871337: -0.9987642713 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9966147953 0.06365116184 0.04153259032: -0.9966147953 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9784602549 0.1460412667 0.1237955413: -0.9784602549 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9451503101 0.2264175441 0.2039324188: -0.9451503101 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8982846661 0.3035122733 0.2806755020: -0.8982846661 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::5.832306001 0.7080816281 0.6117074595 0.1384243340 -3.525869726 0.9668013027 -3.683918369 0.9730667899 -13.29840895 
+DEAL::Cell with center 0.000000000 -1.000000000 0.000000000
+DEAL::-0.5853333358 -0.5853333358 -0.5593080733: -0.5853333358 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5409831617 -0.6518067606 -0.5158409604: -0.5409831617 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4873495734 -0.7224794539 -0.4629781872: -0.4873495734 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4257002918 -0.7922007171 -0.4019874746: -0.4257002918 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3573030376 -0.8565698965 -0.3341365434: -0.3573030376 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2834255319 -0.9119363833 -0.2606931144: -0.2834255319 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2053354953 -0.9553996133 -0.1829249085: -0.2053354953 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1243006489 -0.9848090674 -0.1020996466: -0.1243006489 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04158871337 -0.9987642713 -0.01948504954: -0.04158871337 error value 0.000000000 error grad 0.000000000
+DEAL::0.04153259032 -0.9966147953 0.06365116184: 0.04153259032 error value 0.000000000 error grad 0.000000000
+DEAL::0.1237955413 -0.9784602549 0.1460412667: 0.1237955413 error value 0.000000000 error grad 0.000000000
+DEAL::0.2039324188 -0.9451503101 0.2264175441: 0.2039324188 error value 0.000000000 error grad 0.000000000
+DEAL::0.2806755020 -0.8982846661 0.3035122733: 0.2806755020 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-6.890811517 -1.452971756 0.8592313911 0.7857593783 -9.106658205 4.919866875 6.809982664 -1.343083854 3.017301286 
+DEAL::Cell with center 0.000000000 1.000000000 0.000000000
+DEAL::-0.5593080733 0.5853333358 -0.5853333358: -0.5593080733 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5158409604 0.6518067606 -0.5409831617: -0.5158409604 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4629781872 0.7224794539 -0.4873495734: -0.4629781872 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4019874746 0.7922007171 -0.4257002918: -0.4019874746 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3341365434 0.8565698965 -0.3573030376: -0.3341365434 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2606931144 0.9119363833 -0.2834255319: -0.2606931144 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1829249085 0.9553996133 -0.2053354953: -0.1829249085 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1020996466 0.9848090674 -0.1243006489: -0.1020996466 error value 0.000000000 error grad 0.000000000
+DEAL::-0.01948504954 0.9987642713 -0.04158871337: -0.01948504954 error value 0.000000000 error grad 0.000000000
+DEAL::0.06365116184 0.9966147953 0.04153259032: 0.06365116184 error value 0.000000000 error grad 0.000000000
+DEAL::0.1460412667 0.9784602549 0.1237955413: 0.1460412667 error value 0.000000000 error grad 0.000000000
+DEAL::0.2264175441 0.9451503101 0.2039324188: 0.2264175441 error value 0.000000000 error grad 0.000000000
+DEAL::0.3035122733 0.8982846661 0.2806755020: 0.3035122733 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-7.395100711 0.8334063071 -1.342717800 0.6437129783 7.114003954 -1.223724496 -8.167443670 5.414419780 2.023611946 
+DEAL::Mapping of degree 3
+DEAL::Cell with center 0.6000000000 0.000000000 0.000000000
+DEAL::0.000000000 0.000000000 0.000000000: 0.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.07058823529 0.000000000 0.000000000: 0.07058823529 error value 0.000000000 error grad 0.000000000
+DEAL::0.1411764706 0.000000000 0.000000000: 0.1411764706 error value 0.000000000 error grad 0.000000000
+DEAL::0.2117647059 0.000000000 0.000000000: 0.2117647059 error value 0.000000000 error grad 0.000000000
+DEAL::0.2823529412 0.000000000 0.000000000: 0.2823529412 error value 0.000000000 error grad 0.000000000
+DEAL::0.3529411765 0.000000000 0.000000000: 0.3529411765 error value 0.000000000 error grad 0.000000000
+DEAL::0.4235294118 0.000000000 0.000000000: 0.4235294118 error value 0.000000000 error grad 0.000000000
+DEAL::0.4941176471 0.000000000 0.000000000: 0.4941176471 error value 0.000000000 error grad 0.000000000
+DEAL::0.5647058824 0.000000000 0.000000000: 0.5647058824 error value 0.000000000 error grad 0.000000000
+DEAL::0.6352941176 0.000000000 0.000000000: 0.6352941176 error value 0.000000000 error grad 0.000000000
+DEAL::0.7058823529 0.000000000 0.000000000: 0.7058823529 error value 0.000000000 error grad 0.000000000
+DEAL::0.7764705882 0.000000000 0.000000000: 0.7764705882 error value 0.000000000 error grad 0.000000000
+DEAL::0.8470588235 0.000000000 0.000000000: 0.8470588235 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-16.83261934 0.3424029086 6.096939517 15.89915927 
+DEAL::Mapping of degree 5
+DEAL::Cell with center 0.000000000 0.000000000 -1.000000000
+DEAL::-0.5615400474 -0.5850941618 -0.5850941618: -0.0005615400672 error value 0.000000000 error grad 7.633119026e-10
+DEAL::-0.5266213958 -0.5502723541 -0.6479689802: -0.0005266214013 error value 0.000000000 error grad 1.908279756e-10
+DEAL::-0.4823717789 -0.5062307415 -0.7148766286: -0.0004823718071 error value 0.000000000 error grad 2.980232239e-11
+DEAL::-0.4279792629 -0.4521105898 -0.7825765235: -0.0004279792309 error value 5.960464478e-11 error grad 2.980232239e-11
+DEAL::-0.3632308363 -0.3876570904 -0.8472126133: -0.0003632308245 error value 2.980232239e-11 error grad 1.264405398e-10
+DEAL::-0.2886473745 -0.3133543858 -0.9046955109: -0.0002886474133 error value 0.000000000 error grad 1.812804555e-10
+DEAL::-0.2055209229 -0.2304635804 -0.9511304479: -0.0002055209130 error value 0.000000000 error grad 3.003425064e-11
+DEAL::-0.1158540190 -0.1409630743 -0.9832290326: -0.0001158540398 error value 0.000000000 error grad 3.839941343e-11
+DEAL::-0.02221219975 -0.04740207253 -0.9986555384: -2.221219428e-05 error value 0.000000000 error grad 4.190951586e-11
+DEAL::0.07249034961 0.04731300103 -0.9962711890: 7.249034941e-05 error value 7.450580597e-12 error grad 1.362606383e-11
+DEAL::0.1652458760 0.1401731972 -0.9762526478: 0.0001652459055 error value -2.980232239e-11 error grad 1.202668898e-10
+DEAL::0.2532191833 0.2283332267 -0.9400736644: 0.0002532191873 error value 0.000000000 error grad 6.064333036e-11
+DEAL::0.3339941925 0.3093585177 -0.8903515675: 0.0003339942098 error value 0.000000000 error grad 1.862645149e-10
+DEAL::
+DEAL::-0.01070370674 -0.0001433241516 0.0003489109278 2.614885755e-05 0.01308790112 -0.004225476265 0.002049915552 -0.001000132084 0.0004281401634 -0.0006316589713 0.0005374653339 -0.0001766962856 -0.004902071476 0.002495401859 -0.001557525635 0.0003954216838 -0.0005633962154 0.0009477905035 -0.001705689430 -0.0001163646355 -0.004820720196 -0.009116593361 0.004755000114 -0.001203763604 0.01133630943 -0.002579056978 -0.009948725700 0.001889498115 -0.004190996647 0.01017806911 0.002595032215 -0.001814017773 0.001710586548 -0.003080873489 0.006801867008 0.0007283001542 
+DEAL::Cell with center 1.000000000 0.000000000 0.000000000
+DEAL::0.5850941618 -0.5850941618 -0.5615400474: 0.0005850941539 error value 0.000000000 error grad 3.397987598e-10
+DEAL::0.6479689802 -0.5502723541 -0.5266213958: 0.0006479689479 error value 0.000000000 error grad 4.470348358e-10
+DEAL::0.7148766286 -0.5062307415 -0.4823717789: 0.0007148766518 error value 0.000000000 error grad 2.421152203e-10
+DEAL::0.7825765235 -0.4521105898 -0.4279792629: 0.0007825765610 error value 0.000000000 error grad 5.144678425e-10
+DEAL::0.8472126133 -0.3876570904 -0.3632308363: 0.0008472126722 error value 1.192092896e-10 error grad 1.365714013e-10
+DEAL::0.9046955109 -0.3133543858 -0.2886473745: 0.0009046955109 error value 0.000000000 error grad 3.919883511e-10
+DEAL::0.9511304479 -0.2304635804 -0.2055209229: 0.0009511304498 error value 0.000000000 error grad 1.373819174e-10
+DEAL::0.9832290326 -0.1409630743 -0.1158540190: 0.0009832290411 error value 0.000000000 error grad 3.857789466e-10
+DEAL::0.9986555384 -0.04740207253 -0.02221219975: 0.0009986555576 error value 0.000000000 error grad 2.343133616e-10
+DEAL::0.9962711890 0.04731300103 0.07249034961: 0.0009962711930 error value 0.000000000 error grad 2.977178326e-10
+DEAL::0.9762526478 0.1401731972 0.1652458760: 0.0009762526751 error value 0.000000000 error grad 3.417738128e-10
+DEAL::0.9400736644 0.2283332267 0.2532191833: 0.0009400737882 error value 0.000000000 error grad 5.960464478e-11
+DEAL::0.8903515675 0.3093585177 0.3339941925: 0.0008903515935 error value 0.000000000 error grad 1.264405398e-10
+DEAL::
+DEAL::-0.01736172295 0.0004503556788 0.0001789980233 -5.144322291e-06 0.003745887756 -0.001515596509 0.001099653721 -0.0004737645090 -4.256209731e-05 -1.399764698e-05 0.0004909711778 6.559645385e-05 0.007672265053 -0.003499441385 0.002271631718 -0.001207180738 -7.368709147e-05 1.825724170e-05 0.0005129659772 5.779013783e-05 0.006666445732 0.0003004362285 -0.0009794752598 9.978847951e-05 0.0009141848087 0.008405218124 0.001603288889 6.987068802e-05 -0.001206738114 0.002039612055 0.004898586273 -0.001774510145 0.0001739098728 -8.566421457e-07 -0.001979302645 -0.0003833430707 
+DEAL::Cell with center 0.000000000 0.000000000 1.000000000
+DEAL::-0.5850941618 -0.5615400474 0.5850941618: -0.0005850941539 error value 0.000000000 error grad 2.421152203e-10
+DEAL::-0.5502723541 -0.5266213958 0.6479689802: -0.0005502722859 error value 0.000000000 error grad 3.686343746e-10
+DEAL::-0.5062307415 -0.4823717789 0.7148766286: -0.0005062307119 error value 0.000000000 error grad 1.794337408e-10
+DEAL::-0.4521105898 -0.4279792629 0.7825765235: -0.0004521105289 error value 0.000000000 error grad 1.032382713e-10
+DEAL::-0.3876570904 -0.3632308363 0.8472126133: -0.0003876571059 error value 0.000000000 error grad 8.940696716e-11
+DEAL::-0.3133543858 -0.2886473745 0.9046955109: -0.0003133543432 error value 0.000000000 error grad 7.028863536e-11
+DEAL::-0.2304635804 -0.2055209229 0.9511304479: -0.0002304636240 error value 0.000000000 error grad 1.535976502e-11
+DEAL::-0.1409630743 -0.1158540190 0.9832290326: -0.0001409631073 error value 1.490116119e-11 error grad 1.237783351e-10
+DEAL::-0.04740207253 -0.02221219975 0.9986555384: -4.740206152e-05 error value -1.490116119e-11 error grad 3.072747035e-11
+DEAL::0.04731300103 0.07249034961 0.9962711890: 4.731298983e-05 error value 7.450580597e-12 error grad 1.192674830e-10
+DEAL::0.1401731972 0.1652458760 0.9762526478: 0.0001401732117 error value 0.000000000 error grad 1.862645149e-12
+DEAL::0.2283332267 0.2532191833 0.9400736644: 0.0002283332497 error value -1.490116119e-11 error grad 6.705522537e-11
+DEAL::0.3093585177 0.3339941925 0.8903515675: 0.0003093585074 error value -2.980232239e-11 error grad 6.828570065e-11
+DEAL::
+DEAL::-0.009020513535 0.0005723311901 -0.0002013786137 9.273748100e-05 -0.01016140079 0.003486890078 -0.001876659513 0.0005260861516 -0.0004842967987 0.0008831225038 -0.001359152317 -0.0003440904319 0.01230634785 -0.005418019772 0.002946031094 -0.001612691998 0.0004967076778 -0.0006768767834 0.0006906953454 -0.0003812220991 0.002883476019 0.008129294395 -0.003512819529 0.001526075363 -0.01236182308 0.001248167872 0.008252410889 -0.002913965225 0.005315990925 -0.01169214058 0.005116312504 0.005297445774 -0.001320580840 0.001920366287 -0.002288476944 0.001447246909 
+DEAL::Cell with center -1.000000000 0.000000000 0.000000000
+DEAL::-0.5850941618 -0.5615400474 -0.5850941618: -0.0005850941539 error value 0.000000000 error grad 3.397987598e-10
+DEAL::-0.6479689802 -0.5266213958 -0.5502723541: -0.0006479689479 error value 0.000000000 error grad 4.470348358e-10
+DEAL::-0.7148766286 -0.4823717789 -0.5062307415: -0.0007148766518 error value 0.000000000 error grad 2.421152203e-10
+DEAL::-0.7825765235 -0.4279792629 -0.4521105898: -0.0007825765610 error value 0.000000000 error grad 5.144678425e-10
+DEAL::-0.8472126133 -0.3632308363 -0.3876570904: -0.0008472126722 error value -1.192092896e-10 error grad 1.365714013e-10
+DEAL::-0.9046955109 -0.2886473745 -0.3133543858: -0.0009046955109 error value 0.000000000 error grad 3.919883511e-10
+DEAL::-0.9511304479 -0.2055209229 -0.2304635804: -0.0009511304498 error value 0.000000000 error grad 1.373819174e-10
+DEAL::-0.9832290326 -0.1158540190 -0.1409630743: -0.0009832290411 error value 0.000000000 error grad 3.857789466e-10
+DEAL::-0.9986555384 -0.02221219975 -0.04740207253: -0.0009986555576 error value 0.000000000 error grad 2.343133616e-10
+DEAL::-0.9962711890 0.07249034961 0.04731300103: -0.0009962711930 error value 0.000000000 error grad 2.977178326e-10
+DEAL::-0.9762526478 0.1652458760 0.1401731972: -0.0009762526751 error value 0.000000000 error grad 3.417738128e-10
+DEAL::-0.9400736644 0.2532191833 0.2283332267: -0.0009400737882 error value 0.000000000 error grad 5.960464478e-11
+DEAL::-0.8903515675 0.3339941925 0.3093585177: -0.0008903515935 error value 0.000000000 error grad 1.264405398e-10
+DEAL::
+DEAL::0.01736172295 -0.0004503556788 -0.0001789980233 5.144322291e-06 -0.003745887756 0.001515596509 -0.001099653721 0.0004737645090 4.256209731e-05 1.399764698e-05 -0.0004909711778 -6.559645385e-05 -0.007672265053 0.003499441385 -0.002271631718 0.001207180738 7.368709147e-05 -1.825724170e-05 -0.0005129659772 -5.779013783e-05 -0.006666445732 -0.0003004362285 0.0009794752598 -9.978847951e-05 -0.0009141848087 -0.008405218124 -0.001603288889 -6.987068802e-05 0.001206738114 -0.002039612055 -0.004898586273 0.001774510145 -0.0001739098728 8.566421457e-07 0.001979302645 0.0003833430707 
+DEAL::Cell with center 0.000000000 -1.000000000 0.000000000
+DEAL::-0.5850941618 -0.5850941618 -0.5615400474: -0.0005850941539 error value 0.000000000 error grad 2.421152203e-10
+DEAL::-0.5502723541 -0.6479689802 -0.5266213958: -0.0005502722859 error value 0.000000000 error grad 3.686343746e-10
+DEAL::-0.5062307415 -0.7148766286 -0.4823717789: -0.0005062307119 error value 0.000000000 error grad 1.794337408e-10
+DEAL::-0.4521105898 -0.7825765235 -0.4279792629: -0.0004521105289 error value 0.000000000 error grad 1.032382713e-10
+DEAL::-0.3876570904 -0.8472126133 -0.3632308363: -0.0003876571059 error value 0.000000000 error grad 8.940696716e-11
+DEAL::-0.3133543858 -0.9046955109 -0.2886473745: -0.0003133543432 error value 0.000000000 error grad 7.028863536e-11
+DEAL::-0.2304635804 -0.9511304479 -0.2055209229: -0.0002304636240 error value 0.000000000 error grad 1.535976502e-11
+DEAL::-0.1409630743 -0.9832290326 -0.1158540190: -0.0001409631073 error value 1.490116119e-11 error grad 1.237783351e-10
+DEAL::-0.04740207253 -0.9986555384 -0.02221219975: -4.740206152e-05 error value -1.490116119e-11 error grad 3.072747035e-11
+DEAL::0.04731300103 -0.9962711890 0.07249034961: 4.731298983e-05 error value 7.450580597e-12 error grad 1.192674830e-10
+DEAL::0.1401731972 -0.9762526478 0.1652458760: 0.0001401732117 error value 0.000000000 error grad 1.862645149e-12
+DEAL::0.2283332267 -0.9400736644 0.2532191833: 0.0002283332497 error value -1.490116119e-11 error grad 6.705522537e-11
+DEAL::0.3093585177 -0.8903515675 0.3339941925: 0.0003093585074 error value -2.980232239e-11 error grad 6.828570065e-11
+DEAL::
+DEAL::-0.009020513535 0.0005723311901 -0.0002013786137 9.273748100e-05 -0.01016140079 0.003486890078 -0.001876659513 0.0005260861516 -0.0004842967987 0.0008831225038 -0.001359152317 -0.0003440904319 0.01230634785 -0.005418019772 0.002946031094 -0.001612691998 0.0004967076778 -0.0006768767834 0.0006906953454 -0.0003812220991 0.002883476019 0.008129294395 -0.003512819529 0.001526075363 -0.01236182308 0.001248167872 0.008252410889 -0.002913965225 0.005315990925 -0.01169214058 0.005116312504 0.005297445774 -0.001320580840 0.001920366287 -0.002288476944 0.001447246909 
+DEAL::Cell with center 0.000000000 1.000000000 0.000000000
+DEAL::-0.5615400474 0.5850941618 -0.5850941618: -0.0005615400672 error value 0.000000000 error grad 7.633119026e-10
+DEAL::-0.5266213958 0.6479689802 -0.5502723541: -0.0005266214013 error value 0.000000000 error grad 1.908279756e-10
+DEAL::-0.4823717789 0.7148766286 -0.5062307415: -0.0004823718071 error value 0.000000000 error grad 2.980232239e-11
+DEAL::-0.4279792629 0.7825765235 -0.4521105898: -0.0004279792309 error value 5.960464478e-11 error grad 2.980232239e-11
+DEAL::-0.3632308363 0.8472126133 -0.3876570904: -0.0003632308245 error value 2.980232239e-11 error grad 1.264405398e-10
+DEAL::-0.2886473745 0.9046955109 -0.3133543858: -0.0002886474133 error value 0.000000000 error grad 1.812804555e-10
+DEAL::-0.2055209229 0.9511304479 -0.2304635804: -0.0002055209130 error value 0.000000000 error grad 3.003425064e-11
+DEAL::-0.1158540190 0.9832290326 -0.1409630743: -0.0001158540398 error value 0.000000000 error grad 3.839941343e-11
+DEAL::-0.02221219975 0.9986555384 -0.04740207253: -2.221219428e-05 error value 0.000000000 error grad 4.190951586e-11
+DEAL::0.07249034961 0.9962711890 0.04731300103: 7.249034941e-05 error value 7.450580597e-12 error grad 1.362606383e-11
+DEAL::0.1652458760 0.9762526478 0.1401731972: 0.0001652459055 error value -2.980232239e-11 error grad 1.202668898e-10
+DEAL::0.2532191833 0.9400736644 0.2283332267: 0.0002532191873 error value 0.000000000 error grad 6.064333036e-11
+DEAL::0.3339941925 0.8903515675 0.3093585177: 0.0003339942098 error value 0.000000000 error grad 1.862645149e-10
+DEAL::
+DEAL::-0.01070370674 -0.0001433241516 0.0003489109278 2.614885755e-05 0.01308790112 -0.004225476265 0.002049915552 -0.001000132084 0.0004281401634 -0.0006316589713 0.0005374653339 -0.0001766962856 -0.004902071476 0.002495401859 -0.001557525635 0.0003954216838 -0.0005633962154 0.0009477905035 -0.001705689430 -0.0001163646355 -0.004820720196 -0.009116593361 0.004755000114 -0.001203763604 0.01133630943 -0.002579056978 -0.009948725700 0.001889498115 -0.004190996647 0.01017806911 0.002595032215 -0.001814017773 0.001710586548 -0.003080873489 0.006801867008 0.0007283001542 
diff --git a/tests/matrix_free/point_evaluation_31.cc b/tests/matrix_free/point_evaluation_31.cc
new file mode 100644 (file)
index 0000000..0a97d92
--- /dev/null
@@ -0,0 +1,160 @@
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2024 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// Part of the source code is dual licensed under Apache-2.0 WITH
+// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
+// governing the source code and code contributions can be found in
+// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
+//
+// ------------------------------------------------------------------------
+
+
+// check FEPointEvaluation for dim != spacedim, otherwise same test as
+// point_evaluation_03
+
+#include <deal.II/base/function_lib.h>
+
+#include <deal.II/dofs/dof_handler.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/fe_values.h>
+#include <deal.II/fe/mapping_q.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/lac/vector.h>
+
+#include <deal.II/matrix_free/fe_point_evaluation.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <iostream>
+
+#include "../tests.h"
+
+
+
+template <int dim>
+class MyFunction : public Function<dim>
+{
+public:
+  MyFunction()
+    : Function<dim>(dim)
+  {}
+
+  double
+  value(const Point<dim> &p, const unsigned int component) const override
+  {
+    return component + p[component];
+  }
+};
+
+
+
+template <int dim, int spacedim>
+void
+test(const unsigned int degree)
+{
+  Triangulation<dim, spacedim> tria;
+
+  if constexpr (dim + 1 == spacedim)
+    GridGenerator::hyper_sphere(tria, Point<spacedim>(), 1.);
+  else if constexpr (dim == 1)
+    GridGenerator::hyper_rectangle(tria, Point<dim>(), Point<dim>(1.2));
+  else
+    static_assert("Case not implemented");
+
+  MappingQ<dim, spacedim> mapping(degree);
+  deallog << "Mapping of degree " << degree << std::endl;
+
+  std::vector<Point<dim>> unit_points;
+  for (unsigned int i = 0; i < 13; ++i)
+    {
+      Point<dim> p;
+      for (unsigned int d = 0; d < dim; ++d)
+        p[d] = static_cast<double>(i) / 17. + 0.015625 * d;
+      unit_points.push_back(p);
+    }
+
+  FESystem<dim, spacedim> fe(FE_Q<dim, spacedim>(degree), spacedim);
+  FEValues<dim, spacedim> fe_values(mapping,
+                                    fe,
+                                    Quadrature<dim>(unit_points),
+                                    update_values | update_gradients);
+
+  DoFHandler<dim, spacedim> dof_handler(tria);
+  dof_handler.distribute_dofs(fe);
+  Vector<double> vector(dof_handler.n_dofs());
+
+  FEPointEvaluation<spacedim, dim, spacedim> evaluator(
+    mapping, fe, update_values | update_gradients);
+
+  VectorTools::interpolate(mapping,
+                           dof_handler,
+                           MyFunction<spacedim>(),
+                           vector);
+
+  std::vector<double>              solution_values(fe.dofs_per_cell);
+  std::vector<Tensor<1, spacedim>> function_values(unit_points.size());
+  std::vector<Tensor<2, spacedim>> function_gradients(unit_points.size());
+
+  FEValuesExtractors::Vector extractor(0);
+
+  for (const auto &cell : dof_handler.active_cell_iterators())
+    {
+      fe_values.reinit(cell);
+      fe_values[extractor].get_function_values(vector, function_values);
+      fe_values[extractor].get_function_gradients(vector, function_gradients);
+
+      cell->get_dof_values(vector,
+                           solution_values.begin(),
+                           solution_values.end());
+
+      evaluator.reinit(cell, unit_points);
+      evaluator.evaluate(solution_values,
+                         EvaluationFlags::values | EvaluationFlags::gradients);
+
+      deallog << "Cell with center " << cell->center(true) << std::endl;
+      for (unsigned int i = 0; i < function_values.size(); ++i)
+        deallog << mapping.transform_unit_to_real_cell(cell, unit_points[i])
+                << ": " << evaluator.get_value(i) << " error value "
+                << (function_values[i] - evaluator.get_value(i)).norm()
+                << " error grad "
+                << (evaluator.get_gradient(i) - function_gradients[i]).norm()
+                << std::endl;
+      deallog << std::endl;
+
+      for (unsigned int i = 0; i < unit_points.size(); ++i)
+        {
+          evaluator.submit_value(evaluator.get_value(i), i);
+          evaluator.submit_gradient(evaluator.get_gradient(i), i);
+        }
+
+      evaluator.test_and_sum(solution_values,
+                             EvaluationFlags::values |
+                               EvaluationFlags::gradients);
+
+      for (const auto i : solution_values)
+        deallog << i << ' ';
+      deallog << std::endl;
+    }
+}
+
+
+
+int
+main()
+{
+  initlog();
+  deallog << std::setprecision(10);
+
+  test<1, 2>(2);
+  test<1, 3>(6);
+  test<2, 3>(5);
+}
diff --git a/tests/matrix_free/point_evaluation_31.output b/tests/matrix_free/point_evaluation_31.output
new file mode 100644 (file)
index 0000000..28850f7
--- /dev/null
@@ -0,0 +1,180 @@
+
+DEAL::Mapping of degree 2
+DEAL::Cell with center 1.110223025e-16 -1.000000000
+DEAL::0.7071067812 -0.7071067812: 0.7071067812 0.2928932188 error value 0.000000000 error grad 0.000000000
+DEAL::0.6239177481 -0.7719689473: 0.6239177481 0.2280310527 error value 0.000000000 error grad 0.000000000
+DEAL::0.5407287150 -0.8287233426: 0.5407287150 0.1712766574 error value 0.000000000 error grad 0.000000000
+DEAL::0.4575396819 -0.8773699672: 0.4575396819 0.1226300328 error value 0.000000000 error grad 0.000000000
+DEAL::0.3743506489 -0.9179088210: 0.3743506489 0.08209117898 error value 0.000000000 error grad 0.000000000
+DEAL::0.2911616158 -0.9503399041: 0.2911616158 0.04966009592 error value 0.000000000 error grad 0.000000000
+DEAL::0.2079725827 -0.9746632164: 0.2079725827 0.02533678363 error value 0.000000000 error grad 0.000000000
+DEAL::0.1247835496 -0.9908787579: 0.1247835496 0.009121242108 error value 0.000000000 error grad 0.000000000
+DEAL::0.04159451654 -0.9989865287: 0.04159451654 0.001013471345 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04159451654 -0.9989865287: -0.04159451654 0.001013471345 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1247835496 -0.9908787579: -0.1247835496 0.009121242108 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2079725827 -0.9746632164: -0.2079725827 0.02533678363 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2911616158 -0.9503399041: -0.2911616158 0.04966009592 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::13.87776080 5.327068422 -4.547355687 1.462288083 -6.626761534 -5.722171179 
+DEAL::Cell with center -1.000000000 -1.110223025e-16
+DEAL::-0.7071067812 -0.7071067812: -0.7071067812 0.2928932188 error value 0.000000000 error grad 0.000000000
+DEAL::-0.7719689473 -0.6239177481: -0.7719689473 0.3760822519 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8287233426 -0.5407287150: -0.8287233426 0.4592712850 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8773699672 -0.4575396819: -0.8773699672 0.5424603181 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9179088210 -0.3743506489: -0.9179088210 0.6256493511 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9503399041 -0.2911616158: -0.9503399041 0.7088383842 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9746632164 -0.2079725827: -0.9746632164 0.7920274173 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9908787579 -0.1247835496: -0.9908787579 0.8752164504 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9989865287 -0.04159451654: -0.9989865287 0.9584054835 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9989865287 0.04159451654: -0.9989865287 1.041594517 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9908787579 0.1247835496: -0.9908787579 1.124783550 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9746632164 0.2079725827: -0.9746632164 1.207972583 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9503399041 0.2911616158: -0.9503399041 1.291161616 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::1.593504408 -10.14419678 1.552253481 4.457390289 -15.07857256 15.98316292 
+DEAL::Cell with center 1.000000000 1.110223025e-16
+DEAL::0.7071067812 0.7071067812: 0.7071067812 1.707106781 error value 0.000000000 error grad 0.000000000
+DEAL::0.7719689473 0.6239177481: 0.7719689473 1.623917748 error value 0.000000000 error grad 0.000000000
+DEAL::0.8287233426 0.5407287150: 0.8287233426 1.540728715 error value 0.000000000 error grad 0.000000000
+DEAL::0.8773699672 0.4575396819: 0.8773699672 1.457539682 error value 0.000000000 error grad 0.000000000
+DEAL::0.9179088210 0.3743506489: 0.9179088210 1.374350649 error value 0.000000000 error grad 0.000000000
+DEAL::0.9503399041 0.2911616158: 0.9503399041 1.291161616 error value 0.000000000 error grad 0.000000000
+DEAL::0.9746632164 0.2079725827: 0.9746632164 1.207972583 error value 0.000000000 error grad 0.000000000
+DEAL::0.9908787579 0.1247835496: 0.9908787579 1.124783550 error value 0.000000000 error grad 0.000000000
+DEAL::0.9989865287 0.04159451654: 0.9989865287 1.041594517 error value 0.000000000 error grad 0.000000000
+DEAL::0.9989865287 -0.04159451654: 0.9989865287 0.9584054835 error value 0.000000000 error grad 0.000000000
+DEAL::0.9908787579 -0.1247835496: 0.9908787579 0.8752164504 error value 0.000000000 error grad 0.000000000
+DEAL::0.9746632164 -0.2079725827: 0.9746632164 0.7920274173 error value 0.000000000 error grad 0.000000000
+DEAL::0.9503399041 -0.2911616158: 0.9503399041 0.7088383842 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-1.593504408 17.61132481 -1.552253481 -4.637321085 15.07857256 2.729639850 
+DEAL::Cell with center -1.110223025e-16 1.000000000
+DEAL::-0.7071067812 0.7071067812: -0.7071067812 1.707106781 error value 0.000000000 error grad 0.000000000
+DEAL::-0.6239177481 0.7719689473: -0.6239177481 1.771968947 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5407287150 0.8287233426: -0.5407287150 1.828723343 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4575396819 0.8773699672: -0.4575396819 1.877369967 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3743506489 0.9179088210: -0.3743506489 1.917908821 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2911616158 0.9503399041: -0.2911616158 1.950339904 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2079725827 0.9746632164: -0.2079725827 1.974663216 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1247835496 0.9908787579: -0.1247835496 1.990878758 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04159451654 0.9989865287: -0.04159451654 1.998986529 error value 0.000000000 error grad 0.000000000
+DEAL::0.04159451654 0.9989865287: 0.04159451654 1.998986529 error value 0.000000000 error grad 0.000000000
+DEAL::0.1247835496 0.9908787579: 0.1247835496 1.990878758 error value 0.000000000 error grad 0.000000000
+DEAL::0.2079725827 0.9746632164: 0.2079725827 1.974663216 error value 0.000000000 error grad 0.000000000
+DEAL::0.2911616158 0.9503399041: 0.2911616158 1.950339904 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-13.87776080 2.140059606 4.547355687 -1.642218879 6.626761534 24.43497395 
+DEAL::Mapping of degree 6
+DEAL::Cell with center 0.6000000000 0.000000000 0.000000000
+DEAL::0.000000000 0.000000000 0.000000000: 0.000000000 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.07058823529 0.000000000 0.000000000: 0.07058823529 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.1411764706 0.000000000 0.000000000: 0.1411764706 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.2117647059 0.000000000 0.000000000: 0.2117647059 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.2823529412 0.000000000 0.000000000: 0.2823529412 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.3529411765 0.000000000 0.000000000: 0.3529411765 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.4235294118 0.000000000 0.000000000: 0.4235294118 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.4941176471 0.000000000 0.000000000: 0.4941176471 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.5647058824 0.000000000 0.000000000: 0.5647058824 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.6352941176 0.000000000 0.000000000: 0.6352941176 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.7058823529 0.000000000 0.000000000: 0.7058823529 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.7764705882 0.000000000 0.000000000: 0.7764705882 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::0.8470588235 0.000000000 0.000000000: 0.8470588235 1.000000000 2.000000000 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-24.08827623 0.9701664240 1.940332848 -0.4858920148 0.1065407208 0.2130814416 14.27400475 -5.252146187 6.306755451 13.56627423 1.185162351 2.314958554 3.555859998 4.439742677 1.963991039 -0.3512594117 4.629917107 7.111719996 8.879485353 3.927982077 -0.7025188235 
+DEAL::Mapping of degree 5
+DEAL::Cell with center 0.000000000 0.000000000 -1.000000000
+DEAL::-0.5615400474 -0.5850941618 -0.5850941618: -0.5615400474 0.4149058382 1.414905838 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5266213958 -0.5502723541 -0.6479689802: -0.5266213958 0.4497276459 1.352031020 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4823717789 -0.5062307415 -0.7148766286: -0.4823717789 0.4937692585 1.285123371 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4279792629 -0.4521105898 -0.7825765235: -0.4279792629 0.5478894102 1.217423477 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3632308363 -0.3876570904 -0.8472126133: -0.3632308363 0.6123429096 1.152787387 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2886473745 -0.3133543858 -0.9046955109: -0.2886473745 0.6866456142 1.095304489 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2055209229 -0.2304635804 -0.9511304479: -0.2055209229 0.7695364196 1.048869552 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1158540190 -0.1409630743 -0.9832290326: -0.1158540190 0.8590369257 1.016770967 error value 0.000000000 error grad 0.000000000
+DEAL::-0.02221219975 -0.04740207253 -0.9986555384: -0.02221219975 0.9525979275 1.001344462 error value 0.000000000 error grad 0.000000000
+DEAL::0.07249034961 0.04731300103 -0.9962711890: 0.07249034961 1.047313001 1.003728811 error value 0.000000000 error grad 0.000000000
+DEAL::0.1652458760 0.1401731972 -0.9762526478: 0.1652458760 1.140173197 1.023747352 error value 0.000000000 error grad 0.000000000
+DEAL::0.2532191833 0.2283332267 -0.9400736644: 0.2532191833 1.228333227 1.059926336 error value 0.000000000 error grad 0.000000000
+DEAL::0.3339941925 0.3093585177 -0.8903515675: 0.3339941925 1.309358518 1.109648432 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-10.70371237 -8.122658156 19.15740149 -0.1433240303 0.5972630792 -0.4004905652 0.3489112012 -0.1611056471 -0.09845159279 0.02614886614 0.1115009155 0.04267123379 13.08791002 -4.225478682 2.049916726 -1.000132943 -9.781366612 3.231256214 -1.643273662 0.4081279801 -2.985803869 1.004325273 -0.6328797264 0.2378468053 0.4281401656 -0.6316589999 0.5374651994 -0.1766963013 -0.5455935299 0.9228626320 -1.271200321 -0.4079450041 -0.08003103416 0.09347832520 -0.3150671469 -0.1933057582 -4.902069661 2.495401110 -1.557525053 0.3954213599 12.20631835 -5.482816834 3.076073867 -1.687742861 -7.872304524 3.369838395 -2.011541150 1.057075704 -0.5633963802 0.9477903655 -1.705689444 -0.1163646624 0.4109480283 -0.5975662584 0.7249374710 -0.4438234038 -0.09783196767 0.1403642153 -0.4444810152 -0.1829929382 -4.820722504 -9.116591790 4.754999932 -1.203763635 11.33631086 -2.579056086 -9.948725546 1.889498014 -4.190996521 10.17806803 2.595031753 -1.814017208 1.710586628 -3.080873557 6.801866671 0.7283002287 5.718738768 8.295059088 -3.861814768 1.715693331 -11.69241558 5.443572338 8.623575194 -3.066962809 4.745580737 -10.81803110 8.856762526 5.071415346 -1.059414309 1.643660833 -2.312463637 1.664477689 -0.9959065666 0.03108901506 0.2814856958 0.2794464910 0.4246215590 -0.01441023211 -0.8609560551 -0.3758672815 0.06591907327 -0.2913930944 2.582314109 1.322449834 0.3484222708 -0.5525554637 1.931326786 0.8178051974 
+DEAL::Cell with center 1.000000000 0.000000000 0.000000000
+DEAL::0.5850941618 -0.5850941618 -0.5615400474: 0.5850941618 0.4149058382 1.438459953 error value 0.000000000 error grad 0.000000000
+DEAL::0.6479689802 -0.5502723541 -0.5266213958: 0.6479689802 0.4497276459 1.473378604 error value 0.000000000 error grad 0.000000000
+DEAL::0.7148766286 -0.5062307415 -0.4823717789: 0.7148766286 0.4937692585 1.517628221 error value 0.000000000 error grad 0.000000000
+DEAL::0.7825765235 -0.4521105898 -0.4279792629: 0.7825765235 0.5478894102 1.572020737 error value 0.000000000 error grad 0.000000000
+DEAL::0.8472126133 -0.3876570904 -0.3632308363: 0.8472126133 0.6123429096 1.636769164 error value 0.000000000 error grad 0.000000000
+DEAL::0.9046955109 -0.3133543858 -0.2886473745: 0.9046955109 0.6866456142 1.711352626 error value 0.000000000 error grad 0.000000000
+DEAL::0.9511304479 -0.2304635804 -0.2055209229: 0.9511304479 0.7695364196 1.794479077 error value 0.000000000 error grad 0.000000000
+DEAL::0.9832290326 -0.1409630743 -0.1158540190: 0.9832290326 0.8590369257 1.884145981 error value 0.000000000 error grad 0.000000000
+DEAL::0.9986555384 -0.04740207253 -0.02221219975: 0.9986555384 0.9525979275 1.977787800 error value 0.000000000 error grad 0.000000000
+DEAL::0.9962711890 0.04731300103 0.07249034961: 0.9962711890 1.047313001 2.072490350 error value 0.000000000 error grad 0.000000000
+DEAL::0.9762526478 0.1401731972 0.1652458760: 0.9762526478 1.140173197 2.165245876 error value 0.000000000 error grad 0.000000000
+DEAL::0.9400736644 0.2283332267 0.2532191833: 0.9400736644 1.228333227 2.253219183 error value 0.000000000 error grad 0.000000000
+DEAL::0.8903515675 0.3093585177 0.3339941925: 0.8903515675 1.309358518 2.333994192 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-17.36172704 -8.122658156 -8.908037911 0.4503556671 0.5972630792 -0.09345892828 0.1789981479 -0.1611056471 0.4294577564 -0.005144343710 0.1115009155 0.06367575623 3.745891545 -1.515598339 1.099654314 -0.4737648367 -9.781366612 3.231256214 -1.643273662 0.4081279801 13.84799770 -4.736751748 2.516691314 -1.236050975 -0.04256204942 -0.01399795100 0.4909711954 0.06559652817 -0.5455935299 0.9228626320 -1.271200321 -0.4079450041 0.3055470820 -0.5521786257 0.7133692479 -0.3044055313 7.672264922 -3.499441755 2.271632032 -1.207180772 12.20631835 -5.482816834 3.076073867 -1.687742861 -5.102109263 2.365797750 -1.297434172 0.2453162918 -0.07368722284 0.01825683652 0.5129655849 0.05779023369 0.4109480283 -0.5975662584 0.7249374710 -0.4438234038 -0.7349155707 1.106411417 -1.637204875 -0.2415673670 6.666442640 0.3004379645 -0.9794747511 0.09978833582 0.9141887537 8.405218410 1.603286274 0.06987162062 -1.206738694 2.039608771 4.898586820 -1.774510572 0.1739101565 -0.0008554078016 -1.979301178 -0.3833433380 5.718738768 8.295059088 -3.861814768 1.715693331 -11.69241558 5.443572338 8.623575194 -3.066962809 4.745580737 -10.81803110 8.856762526 5.071415346 -1.059414309 1.643660833 -2.312463637 1.664477689 0.8498135689 -8.785064811 4.057010877 -0.8245288086 12.67512117 5.811752092 -9.206395327 1.583502353 -5.331816142 11.92628371 10.07593268 -2.266077945 2.232919055 -3.634284429 6.753892279 1.162762088 
+DEAL::Cell with center 0.000000000 0.000000000 1.000000000
+DEAL::-0.5850941618 -0.5615400474 0.5850941618: -0.5850941618 0.4384599526 2.585094162 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5502723541 -0.5266213958 0.6479689802: -0.5502723541 0.4733786042 2.647968980 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5062307415 -0.4823717789 0.7148766286: -0.5062307415 0.5176282211 2.714876629 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4521105898 -0.4279792629 0.7825765235: -0.4521105898 0.5720207371 2.782576523 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3876570904 -0.3632308363 0.8472126133: -0.3876570904 0.6367691637 2.847212613 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3133543858 -0.2886473745 0.9046955109: -0.3133543858 0.7113526255 2.904695511 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2304635804 -0.2055209229 0.9511304479: -0.2304635804 0.7944790771 2.951130448 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1409630743 -0.1158540190 0.9832290326: -0.1409630743 0.8841459810 2.983229033 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04740207253 -0.02221219975 0.9986555384: -0.04740207253 0.9777878002 2.998655538 error value 0.000000000 error grad 0.000000000
+DEAL::0.04731300103 0.07249034961 0.9962711890: 0.04731300103 1.072490350 2.996271189 error value 0.000000000 error grad 0.000000000
+DEAL::0.1401731972 0.1652458760 0.9762526478: 0.1401731972 1.165245876 2.976252648 error value 0.000000000 error grad 0.000000000
+DEAL::0.2283332267 0.2532191833 0.9400736644: 0.2283332267 1.253219183 2.940073664 error value 0.000000000 error grad 0.000000000
+DEAL::0.3093585177 0.3339941925 0.8903515675: 0.3093585177 1.333994192 2.890351568 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-9.020495384 -9.805875139 -15.56605258 0.5723305282 -0.1183914793 0.5002207691 -0.2013789246 0.3891844788 0.2595447031 0.09273747043 0.04491231118 0.03238254637 -10.16141045 3.486892747 -1.876660956 0.5260869957 13.46795386 -4.481115215 2.283304020 -1.118091959 4.505979221 -2.026871405 1.566428902 -0.7096828680 -0.4842969881 0.8831224449 -1.359152345 -0.3440903891 0.3668436238 -0.5919188128 0.6254172237 -0.2405509163 -0.1651551330 0.06548242319 0.6668752440 -0.06211270183 12.30633815 -5.418015154 2.946028426 -1.612690327 -5.002089462 2.430599430 -1.427479613 0.3203688258 7.472225319 -3.629045115 2.531722913 -1.357285840 0.4967076236 -0.6768767844 0.6906951862 -0.3812220515 -0.6491559754 1.027100891 -1.671447159 -0.1789660147 -0.2452064134 0.1768778884 0.5814501545 -0.06741247083 2.883470732 8.129295598 -3.512820241 1.526075918 -12.36182074 1.248168249 8.252410084 -2.913964979 5.315990547 -11.69213894 5.116312061 5.297445715 -1.320580523 1.920366269 -2.288476441 1.447246760 -1.985454468 -8.950828301 4.406005405 -1.014146222 12.00571602 1.616348003 -9.577560437 1.736500184 -4.761406332 11.05217587 6.335482217 -2.040047576 1.971752842 -3.357578993 6.777879475 0.9455311583 12.33697871 0.6319649440 -1.677463806 0.4790231627 2.252999066 16.79602659 2.345616494 -0.2361240403 -2.347558315 3.787824448 12.37948775 -2.226571309 0.6962425838 -0.5542662793 -2.027275571 0.05111852134 
+DEAL::Cell with center -1.000000000 0.000000000 0.000000000
+DEAL::-0.5850941618 -0.5615400474 -0.5850941618: -0.5850941618 0.4384599526 1.414905838 error value 0.000000000 error grad 0.000000000
+DEAL::-0.6479689802 -0.5266213958 -0.5502723541: -0.6479689802 0.4733786042 1.449727646 error value 0.000000000 error grad 0.000000000
+DEAL::-0.7148766286 -0.4823717789 -0.5062307415: -0.7148766286 0.5176282211 1.493769258 error value 0.000000000 error grad 0.000000000
+DEAL::-0.7825765235 -0.4279792629 -0.4521105898: -0.7825765235 0.5720207371 1.547889410 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8472126133 -0.3632308363 -0.3876570904: -0.8472126133 0.6367691637 1.612342910 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9046955109 -0.2886473745 -0.3133543858: -0.9046955109 0.7113526255 1.686645614 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9511304479 -0.2055209229 -0.2304635804: -0.9511304479 0.7944790771 1.769536420 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9832290326 -0.1158540190 -0.1409630743: -0.9832290326 0.8841459810 1.859036926 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9986555384 -0.02221219975 -0.04740207253: -0.9986555384 0.9777878002 1.952597927 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9962711890 0.07249034961 0.04731300103: -0.9962711890 1.072490350 2.047313001 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9762526478 0.1652458760 0.1401731972: -0.9762526478 1.165245876 2.140173197 error value 0.000000000 error grad 0.000000000
+DEAL::-0.9400736644 0.2532191833 0.2283332267: -0.9400736644 1.253219183 2.228333227 error value 0.000000000 error grad 0.000000000
+DEAL::-0.8903515675 0.3339941925 0.3093585177: -0.8903515675 1.333994192 2.309358518 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::17.36172704 -9.805875139 -7.224820928 -0.4503556671 -0.1183914793 0.6221956302 -0.1789981479 0.3891844788 -0.1208323695 0.005144343710 0.04491231118 0.1302643605 -3.745891545 1.515598339 -1.099654314 0.4737648367 13.46795386 -4.481115215 2.283304020 -1.118091959 -9.401322774 2.975619681 -1.409886368 0.2901689644 0.04256204942 0.01399795100 -0.4909711954 -0.06559652817 0.3668436238 -0.5919188128 0.6254172237 -0.2405509163 -0.6068900717 0.9626028191 -1.183248296 -0.4717996191 -7.672264922 3.499441755 -2.271632032 1.207180772 -5.002089462 2.430599430 -1.427479613 0.3203688258 12.10629855 -5.547618513 3.206119308 -1.762795395 0.07368722284 -0.01825683652 -0.5129655849 -0.05779023369 -0.6491559754 1.027100891 -1.671447159 -0.1789660147 0.3251884331 -0.5182557325 0.7591797558 -0.5064247560 -6.666442640 -0.3004379645 0.9794747511 -0.09978833582 -0.9141887537 -8.405218410 -1.603286274 -0.06987162062 1.206738694 -2.039608771 -4.898586820 1.774510572 -0.1739101565 0.0008554078016 1.979301178 0.3833433380 -1.985454468 -8.950828301 4.406005405 -1.014146222 12.00571602 1.616348003 -9.577560437 1.736500184 -4.761406332 11.05217587 6.335482217 -2.040047576 1.971752842 -3.357578993 6.777879475 0.9455311583 8.554006805 8.460822578 -4.210809296 1.905310745 -11.02301043 9.638976426 8.994740303 -3.219960639 4.175170926 -9.943923266 12.59721299 4.845384977 -0.7982480956 1.366955397 -2.336450833 1.881708619 
+DEAL::Cell with center 0.000000000 -1.000000000 0.000000000
+DEAL::-0.5850941618 -0.5850941618 -0.5615400474: -0.5850941618 0.4149058382 1.438459953 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5502723541 -0.6479689802 -0.5266213958: -0.5502723541 0.3520310198 1.473378604 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5062307415 -0.7148766286 -0.4823717789: -0.5062307415 0.2851233714 1.517628221 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4521105898 -0.7825765235 -0.4279792629: -0.4521105898 0.2174234765 1.572020737 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3876570904 -0.8472126133 -0.3632308363: -0.3876570904 0.1527873867 1.636769164 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3133543858 -0.9046955109 -0.2886473745: -0.3133543858 0.09530448905 1.711352626 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2304635804 -0.9511304479 -0.2055209229: -0.2304635804 0.04886955206 1.794479077 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1409630743 -0.9832290326 -0.1158540190: -0.1409630743 0.01677096741 1.884145981 error value 0.000000000 error grad 0.000000000
+DEAL::-0.04740207253 -0.9986555384 -0.02221219975: -0.04740207253 0.001344461580 1.977787800 error value 0.000000000 error grad 0.000000000
+DEAL::0.04731300103 -0.9962711890 0.07249034961: 0.04731300103 0.003728811032 2.072490350 error value 0.000000000 error grad 0.000000000
+DEAL::0.1401731972 -0.9762526478 0.1652458760: 0.1401731972 0.02374735218 2.165245876 error value 0.000000000 error grad 0.000000000
+DEAL::0.2283332267 -0.9400736644 0.2532191833: 0.2283332267 0.05992633563 2.253219183 error value 0.000000000 error grad 0.000000000
+DEAL::0.3093585177 -0.8903515675 0.3339941925: 0.3093585177 0.1096484325 2.333994192 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-9.020495384 18.25956426 -8.908037911 0.5723305282 -0.4254231162 -0.09345892828 -0.2013789246 -0.1387248704 0.4294577564 0.09273747043 0.02390778875 0.06367575623 -10.16141045 3.486892747 -1.876660956 0.5260869957 -3.365847707 1.259961806 -0.8662670204 0.3558058210 13.84799770 -4.736751748 2.516691314 -1.236050975 -0.4842969881 0.8831224449 -1.359152345 -0.3440903891 -0.01873449237 0.05373813810 -0.4030191711 -0.1294511432 0.3055470820 -0.5521786257 0.7133692479 -0.3044055313 12.30633815 -5.418015154 2.946028426 -1.612690327 -7.772284723 3.434640075 -2.141586591 1.132128238 -5.102109263 2.365797750 -1.297434172 0.2453162918 0.4967076236 -0.6768767844 0.6906951862 -0.3812220515 -0.01207237241 0.06105368941 -0.4787233000 -0.1203915859 -0.7349155707 1.106411417 -1.637204875 -0.2415673670 2.883470732 8.129295598 -3.512820241 1.526075918 -12.36182074 1.248168249 8.252410084 -2.913964979 5.315990547 -11.69213894 5.116312061 5.297445715 -1.320580523 1.920366269 -2.288476441 1.447246760 -3.831174603 -0.1346744747 0.6304802234 0.08982907761 -0.2447835973 -4.209814321 -1.232121165 -0.2228694511 0.6363288837 -1.165500933 -1.158136355 1.548480203 0.08725605713 -0.2758500280 1.955313982 0.6005742677 0.8498135689 -8.785064811 4.057010877 -0.8245288086 12.67512117 5.811752092 -9.206395327 1.583502353 -5.331816142 11.92628371 10.07593268 -2.266077945 2.232919055 -3.634284429 6.753892279 1.162762088 
+DEAL::Cell with center 0.000000000 1.000000000 0.000000000
+DEAL::-0.5615400474 0.5850941618 -0.5850941618: -0.5615400474 1.585094162 1.414905838 error value 0.000000000 error grad 0.000000000
+DEAL::-0.5266213958 0.6479689802 -0.5502723541: -0.5266213958 1.647968980 1.449727646 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4823717789 0.7148766286 -0.5062307415: -0.4823717789 1.714876629 1.493769258 error value 0.000000000 error grad 0.000000000
+DEAL::-0.4279792629 0.7825765235 -0.4521105898: -0.4279792629 1.782576523 1.547889410 error value 0.000000000 error grad 0.000000000
+DEAL::-0.3632308363 0.8472126133 -0.3876570904: -0.3632308363 1.847212613 1.612342910 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2886473745 0.9046955109 -0.3133543858: -0.2886473745 1.904695511 1.686645614 error value 0.000000000 error grad 0.000000000
+DEAL::-0.2055209229 0.9511304479 -0.2304635804: -0.2055209229 1.951130448 1.769536420 error value 0.000000000 error grad 0.000000000
+DEAL::-0.1158540190 0.9832290326 -0.1409630743: -0.1158540190 1.983229033 1.859036926 error value 0.000000000 error grad 0.000000000
+DEAL::-0.02221219975 0.9986555384 -0.04740207253: -0.02221219975 1.998655538 1.952597927 error value 0.000000000 error grad 0.000000000
+DEAL::0.07249034961 0.9962711890 0.04731300103: 0.07249034961 1.996271189 2.047313001 error value 0.000000000 error grad 0.000000000
+DEAL::0.1652458760 0.9762526478 0.1401731972: 0.1652458760 1.976252648 2.140173197 error value 0.000000000 error grad 0.000000000
+DEAL::0.2532191833 0.9400736644 0.2283332267: 0.2532191833 1.940073664 2.228333227 error value 0.000000000 error grad 0.000000000
+DEAL::0.3339941925 0.8903515675 0.3093585177: 0.3339941925 1.890351568 2.309358518 error value 0.000000000 error grad 0.000000000
+DEAL::
+DEAL::-10.70371237 -16.46388981 -7.224820928 -0.1433240303 0.4752882181 0.6221956302 0.3489112012 0.2192714255 -0.1208323695 0.02614886614 0.01361910133 0.1302643605 13.08791002 -4.225478682 2.049916726 -1.000132943 4.125935383 -1.771234872 1.333041608 -0.5917238523 -9.401322774 2.975619681 -1.409886368 0.2901689644 0.4281401656 -0.6316589999 0.5374651994 -0.1766963013 -0.1038585912 0.02574223609 0.5789232197 0.001741913173 -0.6068900717 0.9626028191 -1.183248296 -0.4717996191 -4.902069661 2.495401110 -1.557525053 0.3954213599 7.572245120 -3.564243435 2.401677473 -1.282233306 12.10629855 -5.547618513 3.206119308 -1.762795395 -0.5633963802 0.9477903655 -1.705689444 -0.1163646624 -0.1594468181 0.09756736246 0.5472078697 -0.004811118571 0.3251884331 -0.5182557325 0.7591797558 -0.5064247560 -4.820722504 -9.116591790 4.754999932 -1.203763635 11.33631086 -2.579056086 -9.948725546 1.889498014 -4.190996521 10.17806803 2.595031753 -1.814017208 1.710586628 -3.080873557 6.801866671 0.7283002287 9.501710676 0.4662014542 -1.328469279 0.2894057493 1.583593910 12.60062250 1.974451384 -0.08312620984 -1.777148505 2.913716609 8.639037284 -2.000540940 0.4350763702 -0.2775608436 -2.003288375 -0.1661124083 8.554006805 8.460822578 -4.210809296 1.905310745 -11.02301043 9.638976426 8.994740303 -3.219960639 4.175170926 -9.943923266 12.59721299 4.845384977 -0.7982480956 1.366955397 -2.336450833 1.881708619 
diff --git a/tests/matrix_free/point_evaluation_32.cc b/tests/matrix_free/point_evaluation_32.cc
new file mode 100644 (file)
index 0000000..b69106e
--- /dev/null
@@ -0,0 +1,167 @@
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2024 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// Part of the source code is dual licensed under Apache-2.0 WITH
+// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
+// governing the source code and code contributions can be found in
+// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
+//
+// ------------------------------------------------------------------------
+
+
+// Check face path of FEPointEvaluation with dim < spacedim, otherwise similar
+// to point_evaluation_23
+
+
+#include <deal.II/base/function_lib.h>
+
+#include <deal.II/dofs/dof_handler.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/fe_system.h>
+#include <deal.II/fe/mapping_q.h>
+
+#include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/tria.h>
+
+#include <deal.II/lac/vector.h>
+
+#include <deal.II/matrix_free/fe_point_evaluation.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <iostream>
+
+#include "../tests.h"
+
+
+
+template <int dim>
+class MyFunction : public Function<dim>
+{
+public:
+  MyFunction()
+    : Function<dim>()
+  {}
+
+  double
+  value(const Point<dim> &p, const unsigned int component) const override
+  {
+    return component + p[component];
+  }
+};
+
+
+
+template <int dim, int spacedim>
+void
+test(const unsigned int degree)
+{
+  Triangulation<dim, spacedim> tria;
+
+
+  if constexpr (dim + 1 == spacedim)
+    GridGenerator::hyper_sphere(tria, Point<spacedim>(), 1.);
+  else if constexpr (dim == 1)
+    GridGenerator::hyper_rectangle(tria, Point<dim>(), Point<dim>(1.2));
+  else
+    static_assert("Case not implemented");
+
+  MappingQ<dim, spacedim> mapping(degree);
+  deallog << "Mapping of degree " << degree << std::endl;
+
+  std::vector<Point<dim - 1>> unit_points;
+  for (unsigned int i = 0; i < (dim > 1 ? 13 : 1); ++i)
+    {
+      Point<dim - 1> p;
+      for (unsigned int d = 0; d < dim - 1; ++d)
+        p[d] = static_cast<double>(i) / 17. + 0.015625 * d;
+      unit_points.push_back(p);
+    }
+
+  FE_Q<dim, spacedim>       fe(degree);
+  DoFHandler<dim, spacedim> dof_handler(tria);
+  dof_handler.distribute_dofs(fe);
+  Vector<double> vector(dof_handler.n_dofs());
+
+  NonMatching::MappingInfo<dim, spacedim> mapping_info(mapping,
+                                                       update_values |
+                                                         update_gradients);
+
+  std::vector<std::vector<Quadrature<dim - 1>>> quad_vec_faces(
+    tria.n_active_cells());
+
+  for (const auto &cell : tria.active_cell_iterators())
+    {
+      for (const auto f : cell->face_indices())
+        {
+          quad_vec_faces[cell->active_cell_index()].push_back(
+            Quadrature<dim - 1>(unit_points));
+        }
+    }
+
+  mapping_info.reinit_faces(tria.active_cell_iterators(), quad_vec_faces);
+
+  FEFacePointEvaluation<1, dim, spacedim> evaluator(mapping_info, fe);
+
+  VectorTools::interpolate(mapping,
+                           dof_handler,
+                           MyFunction<spacedim>(),
+                           vector);
+
+  std::vector<double> solution_values_in(fe.dofs_per_cell);
+  std::vector<double> solution_values_out(fe.dofs_per_cell);
+
+  for (const auto &cell : dof_handler.active_cell_iterators())
+    {
+      cell->get_dof_values(vector,
+                           solution_values_in.begin(),
+                           solution_values_in.end());
+
+      deallog << "Cell with center " << cell->center(true) << std::endl;
+      for (const auto f : cell->face_indices())
+        {
+          evaluator.reinit(cell->active_cell_index(), f);
+
+          evaluator.evaluate(solution_values_in,
+                             EvaluationFlags::values |
+                               EvaluationFlags::gradients);
+          deallog << "face number " << f << std::endl;
+          for (unsigned int i = 0; i < unit_points.size(); ++i)
+            deallog << "value " << evaluator.get_value(i) << " gradient "
+                    << evaluator.get_gradient(i) << std::endl;
+          deallog << std::endl;
+
+          for (unsigned int i = 0; i < unit_points.size(); ++i)
+            {
+              evaluator.submit_value(evaluator.get_value(i), i);
+              evaluator.submit_gradient(evaluator.get_gradient(i), i);
+            }
+
+          evaluator.test_and_sum(solution_values_out,
+                                 EvaluationFlags::values |
+                                   EvaluationFlags::gradients);
+
+          for (double i : solution_values_out)
+            deallog << i << ' ';
+          deallog << std::endl;
+        }
+    }
+}
+
+
+
+int
+main()
+{
+  initlog();
+  deallog << std::setprecision(10);
+
+  test<1, 2>(2);
+  test<1, 3>(2);
+  test<2, 3>(3);
+}
diff --git a/tests/matrix_free/point_evaluation_32.output b/tests/matrix_free/point_evaluation_32.output
new file mode 100644 (file)
index 0000000..ce5302d
--- /dev/null
@@ -0,0 +1,439 @@
+
+DEAL::Mapping of degree 2
+DEAL::Cell with center 1.110223025e-16 -1.000000000
+DEAL::face number 0
+DEAL::value 0.7071067812 gradient 0.5930172805 0.4912716006
+DEAL::
+DEAL::1.965086402 0.4193265404 -1.677306162 
+DEAL::face number 1
+DEAL::value -0.7071067812 gradient 0.5930172805 -0.4912716006
+DEAL::
+DEAL::-0.4193265404 -1.965086402 1.677306162 
+DEAL::Cell with center -1.000000000 -1.110223025e-16
+DEAL::face number 0
+DEAL::value -0.7071067812 gradient 0.4069827195 -0.4912716006
+DEAL::
+DEAL::0.3350376593 0.3473814802 -1.389525921 
+DEAL::face number 1
+DEAL::value -0.7071067812 gradient 0.4069827195 0.4912716006
+DEAL::
+DEAL::0.3473814802 0.3350376593 -1.389525921 
+DEAL::Cell with center 1.000000000 1.110223025e-16
+DEAL::face number 0
+DEAL::value 0.7071067812 gradient 0.4069827195 -0.4912716006
+DEAL::
+DEAL::-0.3350376593 -0.3473814802 1.389525921 
+DEAL::face number 1
+DEAL::value 0.7071067812 gradient 0.4069827195 0.4912716006
+DEAL::
+DEAL::-0.3473814802 -0.3350376593 1.389525921 
+DEAL::Cell with center -1.110223025e-16 1.000000000
+DEAL::face number 0
+DEAL::value -0.7071067812 gradient 0.5930172805 0.4912716006
+DEAL::
+DEAL::-1.965086402 -0.4193265404 1.677306162 
+DEAL::face number 1
+DEAL::value 0.7071067812 gradient 0.5930172805 -0.4912716006
+DEAL::
+DEAL::0.4193265404 1.965086402 -1.677306162 
+DEAL::Mapping of degree 2
+DEAL::Cell with center 0.6000000000 0.000000000 0.000000000
+DEAL::face number 0
+DEAL::value 0.000000000 gradient 1.000000000 0.000000000 0.000000000
+DEAL::
+DEAL::-2.500000000 -0.8333333333 3.333333333 
+DEAL::face number 1
+DEAL::value 1.200000000 gradient 1.000000000 0.000000000 0.000000000
+DEAL::
+DEAL::0.8333333333 3.700000000 -3.333333333 
+DEAL::Mapping of degree 3
+DEAL::Cell with center 0.000000000 0.000000000 -1.000000000
+DEAL::face number 0
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3457731317 -0.3265807932
+DEAL::value -0.5168487451 gradient 0.7285722781 -0.3281813615 -0.3000861670
+DEAL::value -0.4535639572 gradient 0.7951660259 -0.2995682560 -0.2704364567
+DEAL::value -0.3878670072 gradient 0.8530112273 -0.2631442345 -0.2369349811
+DEAL::value -0.3201289971 gradient 0.9015397880 -0.2211799013 -0.1996127500
+DEAL::value -0.2507210286 gradient 0.9404568577 -0.1753423628 -0.1589113348
+DEAL::value -0.1800142036 gradient 0.9696432327 -0.1268797053 -0.1154849532
+DEAL::value -0.1083796240 gradient 0.9890795498 -0.07674256641 -0.07008403854
+DEAL::value -0.03618839165 gradient 0.9987871482 -0.02567941921 -0.02349357912
+DEAL::value 0.03618839165 gradient 0.9987871482 0.02567941921 0.02349357912
+DEAL::value 0.1083796240 gradient 0.9890795498 0.07674256641 0.07008403854
+DEAL::value 0.1800142036 gradient 0.9696432327 0.1268797053 0.1154849532
+DEAL::value 0.2507210286 gradient 0.9404568577 0.1753423628 0.1589113348
+DEAL::
+DEAL::-14.93776624 -0.03458205383 0.4157039597 -0.003056132278 -0.7651666446 13.45578877 -0.05494126836 0.02185965012 -0.2797746925 0.1068644233 -0.02472462950 0.009443968110 -0.4444841979 0.1697778561 0.1768482844 -0.06755003380 
+DEAL::face number 1
+DEAL::value -0.5773502692 gradient 0.6542268683 0.3457731317 -0.3265807932
+DEAL::value -0.5168487451 gradient 0.7285722781 0.3281813615 -0.3000861670
+DEAL::value -0.4535639572 gradient 0.7951660259 0.2995682560 -0.2704364567
+DEAL::value -0.3878670072 gradient 0.8530112273 0.2631442345 -0.2369349811
+DEAL::value -0.3201289971 gradient 0.9015397880 0.2211799013 -0.1996127500
+DEAL::value -0.2507210286 gradient 0.9404568577 0.1753423628 -0.1589113348
+DEAL::value -0.1800142036 gradient 0.9696432327 0.1268797053 -0.1154849532
+DEAL::value -0.1083796240 gradient 0.9890795498 0.07674256641 -0.07008403854
+DEAL::value -0.03618839165 gradient 0.9987871482 0.02567941921 -0.02349357912
+DEAL::value 0.03618839165 gradient 0.9987871482 -0.02567941921 0.02349357912
+DEAL::value 0.1083796240 gradient 0.9890795498 -0.07674256641 0.07008403854
+DEAL::value 0.1800142036 gradient 0.9696432327 -0.1268797053 0.1154849532
+DEAL::value 0.2507210286 gradient 0.9404568577 -0.1753423628 0.1589113348
+DEAL::
+DEAL::-0.03458205383 -14.93776624 -0.003056132278 0.4157039597 -0.05494126836 0.02185965012 -0.7651666446 13.45578877 0.1068644233 -0.2797746925 0.009443968110 -0.02472462950 0.1697778561 -0.4444841979 -0.06755003380 0.1768482844 
+DEAL::face number 2
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3457731317 -0.3265807932
+DEAL::value -0.6059008348 gradient 0.6031982095 -0.3281813615 -0.3628320872
+DEAL::value -0.6308825797 gradient 0.5618835185 -0.2995682560 -0.3955114285
+DEAL::value -0.6522955039 gradient 0.5289103591 -0.2631442345 -0.4241689557
+DEAL::value -0.6701396074 gradient 0.5031439836 -0.2211799013 -0.4484078129
+DEAL::value -0.6844148903 gradient 0.4836526427 -0.1753423628 -0.4679613443
+DEAL::value -0.6951213524 gradient 0.4696912413 -0.1268797053 -0.4826830425
+DEAL::value -0.7022589938 gradient 0.4606979203 -0.07674256641 -0.4925098223
+DEAL::value -0.7058278145 gradient 0.4562958566 -0.02567941921 -0.4974238789
+DEAL::value -0.7058278145 gradient 0.4562958566 0.02567941921 -0.4974238789
+DEAL::value -0.7022589938 gradient 0.4606979203 0.07674256641 -0.4925098223
+DEAL::value -0.6951213524 gradient 0.4696912413 0.1268797053 -0.4826830425
+DEAL::value -0.6844148903 gradient 0.4836526427 0.1753423628 -0.4679613443
+DEAL::
+DEAL::-7.069933216 1.627626224 1.135550187 -0.1943807830 9.186793990 -3.509043057 -1.572573568 0.6006696533 -27.71278762 -13.77501731 3.633545959 1.795000808 29.39600430 14.52186159 -11.22827451 -5.546857546 
+DEAL::face number 3
+DEAL::value 0.5773502692 gradient 0.6542268683 0.3457731317 0.3265807932
+DEAL::value 0.6059008348 gradient 0.6031982095 0.3281813615 0.3628320872
+DEAL::value 0.6308825797 gradient 0.5618835185 0.2995682560 0.3955114285
+DEAL::value 0.6522955039 gradient 0.5289103591 0.2631442345 0.4241689557
+DEAL::value 0.6701396074 gradient 0.5031439836 0.2211799013 0.4484078129
+DEAL::value 0.6844148903 gradient 0.4836526427 0.1753423628 0.4679613443
+DEAL::value 0.6951213524 gradient 0.4696912413 0.1268797053 0.4826830425
+DEAL::value 0.7022589938 gradient 0.4606979203 0.07674256641 0.4925098223
+DEAL::value 0.7058278145 gradient 0.4562958566 0.02567941921 0.4974238789
+DEAL::value 0.7058278145 gradient 0.4562958566 -0.02567941921 0.4974238789
+DEAL::value 0.7022589938 gradient 0.4606979203 -0.07674256641 0.4925098223
+DEAL::value 0.6951213524 gradient 0.4696912413 -0.1268797053 0.4826830425
+DEAL::value 0.6844148903 gradient 0.4836526427 -0.1753423628 0.4679613443
+DEAL::
+DEAL::-1.135550187 0.1943807830 7.069933216 -1.627626224 3.509043057 -9.186793990 -0.6006696533 1.572573568 -3.633545959 -1.795000808 27.71278762 13.77501731 11.22827451 5.546857546 -29.39600430 -14.52186159 
+DEAL::Cell with center 1.000000000 0.000000000 0.000000000
+DEAL::face number 0
+DEAL::value 0.5773502692 gradient 0.6915462634 0.3265807932 0.3265807932
+DEAL::value 0.6059008348 gradient 0.6682295124 0.3628320872 0.3000861670
+DEAL::value 0.6308825797 gradient 0.6429504556 0.3955114285 0.2704364567
+DEAL::value 0.6522955039 gradient 0.6180784136 0.4241689557 0.2369349811
+DEAL::value 0.6701396074 gradient 0.5953162285 0.4484078129 0.1996127500
+DEAL::value 0.6844148903 gradient 0.5758904996 0.4679613443 0.1589113348
+DEAL::value 0.6951213524 gradient 0.5606655261 0.4826830425 0.1154849532
+DEAL::value 0.7022589938 gradient 0.5502225299 0.4925098223 0.07008403854
+DEAL::value 0.7058278145 gradient 0.5449169952 0.4974238789 0.02349357912
+DEAL::value 0.7058278145 gradient 0.5449169952 0.4974238789 -0.02349357912
+DEAL::value 0.7022589938 gradient 0.5502225299 0.4925098223 -0.07008403854
+DEAL::value 0.6951213524 gradient 0.5606655261 0.4826830425 -0.1154849532
+DEAL::value 0.6844148903 gradient 0.5758904996 0.4679613443 -0.1589113348
+DEAL::
+DEAL::-17.74341922 1.193482859 0.6931112378 -0.2135999875 -5.616698021 -10.11729690 3.967904584 1.968232178 9.655479154 -3.688064859 -1.728060199 0.6600602614 32.10102240 -12.26149948 15.92333281 -6.082171919 
+DEAL::face number 1
+DEAL::value 0.5773502692 gradient 0.6915462634 -0.3265807932 0.3265807932
+DEAL::value 0.6059008348 gradient 0.6682295124 -0.3628320872 0.3000861670
+DEAL::value 0.6308825797 gradient 0.6429504556 -0.3955114285 0.2704364567
+DEAL::value 0.6522955039 gradient 0.6180784136 -0.4241689557 0.2369349811
+DEAL::value 0.6701396074 gradient 0.5953162285 -0.4484078129 0.1996127500
+DEAL::value 0.6844148903 gradient 0.5758904996 -0.4679613443 0.1589113348
+DEAL::value 0.6951213524 gradient 0.5606655261 -0.4826830425 0.1154849532
+DEAL::value 0.7022589938 gradient 0.5502225299 -0.4925098223 0.07008403854
+DEAL::value 0.7058278145 gradient 0.5449169952 -0.4974238789 0.02349357912
+DEAL::value 0.7058278145 gradient 0.5449169952 -0.4974238789 -0.02349357912
+DEAL::value 0.7022589938 gradient 0.5502225299 -0.4925098223 -0.07008403854
+DEAL::value 0.6951213524 gradient 0.5606655261 -0.4826830425 -0.1154849532
+DEAL::value 0.6844148903 gradient 0.5758904996 -0.4679613443 -0.1589113348
+DEAL::
+DEAL::1.193482859 -17.74341922 -0.2135999875 0.6931112378 3.967904584 1.968232178 -5.616698021 -10.11729690 -3.688064859 9.655479154 0.6600602614 -1.728060199 -12.26149948 32.10102240 -6.082171919 15.92333281 
+DEAL::face number 2
+DEAL::value 0.5773502692 gradient 0.6915462634 0.3265807932 0.3265807932
+DEAL::value 0.6059008348 gradient 0.6682295124 0.3000861670 0.3628320872
+DEAL::value 0.6308825797 gradient 0.6429504556 0.2704364567 0.3955114285
+DEAL::value 0.6522955039 gradient 0.6180784136 0.2369349811 0.4241689557
+DEAL::value 0.6701396074 gradient 0.5953162285 0.1996127500 0.4484078129
+DEAL::value 0.6844148903 gradient 0.5758904996 0.1589113348 0.4679613443
+DEAL::value 0.6951213524 gradient 0.5606655261 0.1154849532 0.4826830425
+DEAL::value 0.7022589938 gradient 0.5502225299 0.07008403854 0.4925098223
+DEAL::value 0.7058278145 gradient 0.5449169952 0.02349357912 0.4974238789
+DEAL::value 0.7058278145 gradient 0.5449169952 -0.02349357912 0.4974238789
+DEAL::value 0.7022589938 gradient 0.5502225299 -0.07008403854 0.4925098223
+DEAL::value 0.6951213524 gradient 0.5606655261 -0.1154849532 0.4826830425
+DEAL::value 0.6844148903 gradient 0.5758904996 -0.1589113348 0.4679613443
+DEAL::
+DEAL::-17.74341922 0.6931112378 1.193482859 -0.2135999875 9.655479154 -3.688064859 -1.728060199 0.6600602614 -5.616698021 -10.11729690 3.967904584 1.968232178 32.10102240 15.92333281 -12.26149948 -6.082171919 
+DEAL::face number 3
+DEAL::value 0.5773502692 gradient 0.6915462634 0.3265807932 -0.3265807932
+DEAL::value 0.6059008348 gradient 0.6682295124 0.3000861670 -0.3628320872
+DEAL::value 0.6308825797 gradient 0.6429504556 0.2704364567 -0.3955114285
+DEAL::value 0.6522955039 gradient 0.6180784136 0.2369349811 -0.4241689557
+DEAL::value 0.6701396074 gradient 0.5953162285 0.1996127500 -0.4484078129
+DEAL::value 0.6844148903 gradient 0.5758904996 0.1589113348 -0.4679613443
+DEAL::value 0.6951213524 gradient 0.5606655261 0.1154849532 -0.4826830425
+DEAL::value 0.7022589938 gradient 0.5502225299 0.07008403854 -0.4925098223
+DEAL::value 0.7058278145 gradient 0.5449169952 0.02349357912 -0.4974238789
+DEAL::value 0.7058278145 gradient 0.5449169952 -0.02349357912 -0.4974238789
+DEAL::value 0.7022589938 gradient 0.5502225299 -0.07008403854 -0.4925098223
+DEAL::value 0.6951213524 gradient 0.5606655261 -0.1154849532 -0.4826830425
+DEAL::value 0.6844148903 gradient 0.5758904996 -0.1589113348 -0.4679613443
+DEAL::
+DEAL::1.193482859 -0.2135999875 -17.74341922 0.6931112378 -3.688064859 9.655479154 0.6600602614 -1.728060199 3.967904584 1.968232178 -5.616698021 -10.11729690 -12.26149948 -6.082171919 32.10102240 15.92333281 
+DEAL::Cell with center 0.000000000 0.000000000 1.000000000
+DEAL::face number 0
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3457731317 0.3265807932
+DEAL::value -0.6059008348 gradient 0.6031982095 -0.3281813615 0.3628320872
+DEAL::value -0.6308825797 gradient 0.5618835185 -0.2995682560 0.3955114285
+DEAL::value -0.6522955039 gradient 0.5289103591 -0.2631442345 0.4241689557
+DEAL::value -0.6701396074 gradient 0.5031439836 -0.2211799013 0.4484078129
+DEAL::value -0.6844148903 gradient 0.4836526427 -0.1753423628 0.4679613443
+DEAL::value -0.6951213524 gradient 0.4696912413 -0.1268797053 0.4826830425
+DEAL::value -0.7022589938 gradient 0.4606979203 -0.07674256641 0.4925098223
+DEAL::value -0.7058278145 gradient 0.4562958566 -0.02567941921 0.4974238789
+DEAL::value -0.7058278145 gradient 0.4562958566 0.02567941921 0.4974238789
+DEAL::value -0.7022589938 gradient 0.4606979203 0.07674256641 0.4925098223
+DEAL::value -0.6951213524 gradient 0.4696912413 0.1268797053 0.4826830425
+DEAL::value -0.6844148903 gradient 0.4836526427 0.1753423628 0.4679613443
+DEAL::
+DEAL::-7.069933216 1.135550187 1.627626224 -0.1943807830 -27.71278762 -13.77501731 3.633545959 1.795000808 9.186793990 -3.509043057 -1.572573568 0.6006696533 29.39600430 -11.22827451 14.52186159 -5.546857546 
+DEAL::face number 1
+DEAL::value 0.5773502692 gradient 0.6542268683 0.3457731317 -0.3265807932
+DEAL::value 0.6059008348 gradient 0.6031982095 0.3281813615 -0.3628320872
+DEAL::value 0.6308825797 gradient 0.5618835185 0.2995682560 -0.3955114285
+DEAL::value 0.6522955039 gradient 0.5289103591 0.2631442345 -0.4241689557
+DEAL::value 0.6701396074 gradient 0.5031439836 0.2211799013 -0.4484078129
+DEAL::value 0.6844148903 gradient 0.4836526427 0.1753423628 -0.4679613443
+DEAL::value 0.6951213524 gradient 0.4696912413 0.1268797053 -0.4826830425
+DEAL::value 0.7022589938 gradient 0.4606979203 0.07674256641 -0.4925098223
+DEAL::value 0.7058278145 gradient 0.4562958566 0.02567941921 -0.4974238789
+DEAL::value 0.7058278145 gradient 0.4562958566 -0.02567941921 -0.4974238789
+DEAL::value 0.7022589938 gradient 0.4606979203 -0.07674256641 -0.4925098223
+DEAL::value 0.6951213524 gradient 0.4696912413 -0.1268797053 -0.4826830425
+DEAL::value 0.6844148903 gradient 0.4836526427 -0.1753423628 -0.4679613443
+DEAL::
+DEAL::-1.135550187 7.069933216 0.1943807830 -1.627626224 -3.633545959 -1.795000808 27.71278762 13.77501731 3.509043057 -9.186793990 -0.6006696533 1.572573568 11.22827451 -29.39600430 5.546857546 -14.52186159 
+DEAL::face number 2
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3457731317 0.3265807932
+DEAL::value -0.5168487451 gradient 0.7285722781 -0.3281813615 0.3000861670
+DEAL::value -0.4535639572 gradient 0.7951660259 -0.2995682560 0.2704364567
+DEAL::value -0.3878670072 gradient 0.8530112273 -0.2631442345 0.2369349811
+DEAL::value -0.3201289971 gradient 0.9015397880 -0.2211799013 0.1996127500
+DEAL::value -0.2507210286 gradient 0.9404568577 -0.1753423628 0.1589113348
+DEAL::value -0.1800142036 gradient 0.9696432327 -0.1268797053 0.1154849532
+DEAL::value -0.1083796240 gradient 0.9890795498 -0.07674256641 0.07008403854
+DEAL::value -0.03618839165 gradient 0.9987871482 -0.02567941921 0.02349357912
+DEAL::value 0.03618839165 gradient 0.9987871482 0.02567941921 -0.02349357912
+DEAL::value 0.1083796240 gradient 0.9890795498 0.07674256641 -0.07008403854
+DEAL::value 0.1800142036 gradient 0.9696432327 0.1268797053 -0.1154849532
+DEAL::value 0.2507210286 gradient 0.9404568577 0.1753423628 -0.1589113348
+DEAL::
+DEAL::-14.93776624 0.4157039597 -0.03458205383 -0.003056132278 -0.2797746925 0.1068644233 -0.02472462950 0.009443968110 -0.7651666446 13.45578877 -0.05494126836 0.02185965012 -0.4444841979 0.1768482844 0.1697778561 -0.06755003380 
+DEAL::face number 3
+DEAL::value -0.5773502692 gradient 0.6542268683 0.3457731317 0.3265807932
+DEAL::value -0.5168487451 gradient 0.7285722781 0.3281813615 0.3000861670
+DEAL::value -0.4535639572 gradient 0.7951660259 0.2995682560 0.2704364567
+DEAL::value -0.3878670072 gradient 0.8530112273 0.2631442345 0.2369349811
+DEAL::value -0.3201289971 gradient 0.9015397880 0.2211799013 0.1996127500
+DEAL::value -0.2507210286 gradient 0.9404568577 0.1753423628 0.1589113348
+DEAL::value -0.1800142036 gradient 0.9696432327 0.1268797053 0.1154849532
+DEAL::value -0.1083796240 gradient 0.9890795498 0.07674256641 0.07008403854
+DEAL::value -0.03618839165 gradient 0.9987871482 0.02567941921 0.02349357912
+DEAL::value 0.03618839165 gradient 0.9987871482 -0.02567941921 -0.02349357912
+DEAL::value 0.1083796240 gradient 0.9890795498 -0.07674256641 -0.07008403854
+DEAL::value 0.1800142036 gradient 0.9696432327 -0.1268797053 -0.1154849532
+DEAL::value 0.2507210286 gradient 0.9404568577 -0.1753423628 -0.1589113348
+DEAL::
+DEAL::-0.03458205383 -0.003056132278 -14.93776624 0.4157039597 0.1068644233 -0.2797746925 0.009443968110 -0.02472462950 -0.05494126836 0.02185965012 -0.7651666446 13.45578877 0.1697778561 -0.06755003380 -0.4444841979 0.1768482844 
+DEAL::Cell with center -1.000000000 0.000000000 0.000000000
+DEAL::face number 0
+DEAL::value -0.5773502692 gradient 0.6915462634 -0.3265807932 -0.3265807932
+DEAL::value -0.6059008348 gradient 0.6682295124 -0.3000861670 -0.3628320872
+DEAL::value -0.6308825797 gradient 0.6429504556 -0.2704364567 -0.3955114285
+DEAL::value -0.6522955039 gradient 0.6180784136 -0.2369349811 -0.4241689557
+DEAL::value -0.6701396074 gradient 0.5953162285 -0.1996127500 -0.4484078129
+DEAL::value -0.6844148903 gradient 0.5758904996 -0.1589113348 -0.4679613443
+DEAL::value -0.6951213524 gradient 0.5606655261 -0.1154849532 -0.4826830425
+DEAL::value -0.7022589938 gradient 0.5502225299 -0.07008403854 -0.4925098223
+DEAL::value -0.7058278145 gradient 0.5449169952 -0.02349357912 -0.4974238789
+DEAL::value -0.7058278145 gradient 0.5449169952 0.02349357912 -0.4974238789
+DEAL::value -0.7022589938 gradient 0.5502225299 0.07008403854 -0.4925098223
+DEAL::value -0.6951213524 gradient 0.5606655261 0.1154849532 -0.4826830425
+DEAL::value -0.6844148903 gradient 0.5758904996 0.1589113348 -0.4679613443
+DEAL::
+DEAL::17.74341922 -1.193482859 -0.6931112378 0.2135999875 5.616698021 10.11729690 -3.967904584 -1.968232178 -9.655479154 3.688064859 1.728060199 -0.6600602614 -32.10102240 12.26149948 -15.92333281 6.082171919 
+DEAL::face number 1
+DEAL::value -0.5773502692 gradient 0.6915462634 -0.3265807932 0.3265807932
+DEAL::value -0.6059008348 gradient 0.6682295124 -0.3000861670 0.3628320872
+DEAL::value -0.6308825797 gradient 0.6429504556 -0.2704364567 0.3955114285
+DEAL::value -0.6522955039 gradient 0.6180784136 -0.2369349811 0.4241689557
+DEAL::value -0.6701396074 gradient 0.5953162285 -0.1996127500 0.4484078129
+DEAL::value -0.6844148903 gradient 0.5758904996 -0.1589113348 0.4679613443
+DEAL::value -0.6951213524 gradient 0.5606655261 -0.1154849532 0.4826830425
+DEAL::value -0.7022589938 gradient 0.5502225299 -0.07008403854 0.4925098223
+DEAL::value -0.7058278145 gradient 0.5449169952 -0.02349357912 0.4974238789
+DEAL::value -0.7058278145 gradient 0.5449169952 0.02349357912 0.4974238789
+DEAL::value -0.7022589938 gradient 0.5502225299 0.07008403854 0.4925098223
+DEAL::value -0.6951213524 gradient 0.5606655261 0.1154849532 0.4826830425
+DEAL::value -0.6844148903 gradient 0.5758904996 0.1589113348 0.4679613443
+DEAL::
+DEAL::-1.193482859 17.74341922 0.2135999875 -0.6931112378 -3.967904584 -1.968232178 5.616698021 10.11729690 3.688064859 -9.655479154 -0.6600602614 1.728060199 12.26149948 -32.10102240 6.082171919 -15.92333281 
+DEAL::face number 2
+DEAL::value -0.5773502692 gradient 0.6915462634 -0.3265807932 -0.3265807932
+DEAL::value -0.6059008348 gradient 0.6682295124 -0.3628320872 -0.3000861670
+DEAL::value -0.6308825797 gradient 0.6429504556 -0.3955114285 -0.2704364567
+DEAL::value -0.6522955039 gradient 0.6180784136 -0.4241689557 -0.2369349811
+DEAL::value -0.6701396074 gradient 0.5953162285 -0.4484078129 -0.1996127500
+DEAL::value -0.6844148903 gradient 0.5758904996 -0.4679613443 -0.1589113348
+DEAL::value -0.6951213524 gradient 0.5606655261 -0.4826830425 -0.1154849532
+DEAL::value -0.7022589938 gradient 0.5502225299 -0.4925098223 -0.07008403854
+DEAL::value -0.7058278145 gradient 0.5449169952 -0.4974238789 -0.02349357912
+DEAL::value -0.7058278145 gradient 0.5449169952 -0.4974238789 0.02349357912
+DEAL::value -0.7022589938 gradient 0.5502225299 -0.4925098223 0.07008403854
+DEAL::value -0.6951213524 gradient 0.5606655261 -0.4826830425 0.1154849532
+DEAL::value -0.6844148903 gradient 0.5758904996 -0.4679613443 0.1589113348
+DEAL::
+DEAL::17.74341922 -0.6931112378 -1.193482859 0.2135999875 -9.655479154 3.688064859 1.728060199 -0.6600602614 5.616698021 10.11729690 -3.967904584 -1.968232178 -32.10102240 -15.92333281 12.26149948 6.082171919 
+DEAL::face number 3
+DEAL::value -0.5773502692 gradient 0.6915462634 0.3265807932 -0.3265807932
+DEAL::value -0.6059008348 gradient 0.6682295124 0.3628320872 -0.3000861670
+DEAL::value -0.6308825797 gradient 0.6429504556 0.3955114285 -0.2704364567
+DEAL::value -0.6522955039 gradient 0.6180784136 0.4241689557 -0.2369349811
+DEAL::value -0.6701396074 gradient 0.5953162285 0.4484078129 -0.1996127500
+DEAL::value -0.6844148903 gradient 0.5758904996 0.4679613443 -0.1589113348
+DEAL::value -0.6951213524 gradient 0.5606655261 0.4826830425 -0.1154849532
+DEAL::value -0.7022589938 gradient 0.5502225299 0.4925098223 -0.07008403854
+DEAL::value -0.7058278145 gradient 0.5449169952 0.4974238789 -0.02349357912
+DEAL::value -0.7058278145 gradient 0.5449169952 0.4974238789 0.02349357912
+DEAL::value -0.7022589938 gradient 0.5502225299 0.4925098223 0.07008403854
+DEAL::value -0.6951213524 gradient 0.5606655261 0.4826830425 0.1154849532
+DEAL::value -0.6844148903 gradient 0.5758904996 0.4679613443 0.1589113348
+DEAL::
+DEAL::-1.193482859 0.2135999875 17.74341922 -0.6931112378 3.688064859 -9.655479154 -0.6600602614 1.728060199 -3.967904584 -1.968232178 5.616698021 10.11729690 12.26149948 6.082171919 -32.10102240 -15.92333281 
+DEAL::Cell with center 0.000000000 -1.000000000 0.000000000
+DEAL::face number 0
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3265807932 -0.3457731317
+DEAL::value -0.6059008348 gradient 0.6031982095 -0.3628320872 -0.3281813615
+DEAL::value -0.6308825797 gradient 0.5618835185 -0.3955114285 -0.2995682560
+DEAL::value -0.6522955039 gradient 0.5289103591 -0.4241689557 -0.2631442345
+DEAL::value -0.6701396074 gradient 0.5031439836 -0.4484078129 -0.2211799013
+DEAL::value -0.6844148903 gradient 0.4836526427 -0.4679613443 -0.1753423628
+DEAL::value -0.6951213524 gradient 0.4696912413 -0.4826830425 -0.1268797053
+DEAL::value -0.7022589938 gradient 0.4606979203 -0.4925098223 -0.07674256641
+DEAL::value -0.7058278145 gradient 0.4562958566 -0.4974238789 -0.02567941921
+DEAL::value -0.7058278145 gradient 0.4562958566 -0.4974238789 0.02567941921
+DEAL::value -0.7022589938 gradient 0.4606979203 -0.4925098223 0.07674256641
+DEAL::value -0.6951213524 gradient 0.4696912413 -0.4826830425 0.1268797053
+DEAL::value -0.6844148903 gradient 0.4836526427 -0.4679613443 0.1753423628
+DEAL::
+DEAL::-7.069933216 1.135550187 1.627626224 -0.1943807830 -27.71278762 -13.77501731 3.633545959 1.795000808 9.186793990 -3.509043057 -1.572573568 0.6006696533 29.39600430 -11.22827451 14.52186159 -5.546857546 
+DEAL::face number 1
+DEAL::value 0.5773502692 gradient 0.6542268683 0.3265807932 0.3457731317
+DEAL::value 0.6059008348 gradient 0.6031982095 0.3628320872 0.3281813615
+DEAL::value 0.6308825797 gradient 0.5618835185 0.3955114285 0.2995682560
+DEAL::value 0.6522955039 gradient 0.5289103591 0.4241689557 0.2631442345
+DEAL::value 0.6701396074 gradient 0.5031439836 0.4484078129 0.2211799013
+DEAL::value 0.6844148903 gradient 0.4836526427 0.4679613443 0.1753423628
+DEAL::value 0.6951213524 gradient 0.4696912413 0.4826830425 0.1268797053
+DEAL::value 0.7022589938 gradient 0.4606979203 0.4925098223 0.07674256641
+DEAL::value 0.7058278145 gradient 0.4562958566 0.4974238789 0.02567941921
+DEAL::value 0.7058278145 gradient 0.4562958566 0.4974238789 -0.02567941921
+DEAL::value 0.7022589938 gradient 0.4606979203 0.4925098223 -0.07674256641
+DEAL::value 0.6951213524 gradient 0.4696912413 0.4826830425 -0.1268797053
+DEAL::value 0.6844148903 gradient 0.4836526427 0.4679613443 -0.1753423628
+DEAL::
+DEAL::-1.135550187 7.069933216 0.1943807830 -1.627626224 -3.633545959 -1.795000808 27.71278762 13.77501731 3.509043057 -9.186793990 -0.6006696533 1.572573568 11.22827451 -29.39600430 5.546857546 -14.52186159 
+DEAL::face number 2
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3265807932 -0.3457731317
+DEAL::value -0.5168487451 gradient 0.7285722781 -0.3000861670 -0.3281813615
+DEAL::value -0.4535639572 gradient 0.7951660259 -0.2704364567 -0.2995682560
+DEAL::value -0.3878670072 gradient 0.8530112273 -0.2369349811 -0.2631442345
+DEAL::value -0.3201289971 gradient 0.9015397880 -0.1996127500 -0.2211799013
+DEAL::value -0.2507210286 gradient 0.9404568577 -0.1589113348 -0.1753423628
+DEAL::value -0.1800142036 gradient 0.9696432327 -0.1154849532 -0.1268797053
+DEAL::value -0.1083796240 gradient 0.9890795498 -0.07008403854 -0.07674256641
+DEAL::value -0.03618839165 gradient 0.9987871482 -0.02349357912 -0.02567941921
+DEAL::value 0.03618839165 gradient 0.9987871482 0.02349357912 0.02567941921
+DEAL::value 0.1083796240 gradient 0.9890795498 0.07008403854 0.07674256641
+DEAL::value 0.1800142036 gradient 0.9696432327 0.1154849532 0.1268797053
+DEAL::value 0.2507210286 gradient 0.9404568577 0.1589113348 0.1753423628
+DEAL::
+DEAL::-14.93776624 0.4157039597 -0.03458205383 -0.003056132278 -0.2797746925 0.1068644233 -0.02472462950 0.009443968110 -0.7651666446 13.45578877 -0.05494126836 0.02185965012 -0.4444841979 0.1768482844 0.1697778561 -0.06755003380 
+DEAL::face number 3
+DEAL::value -0.5773502692 gradient 0.6542268683 -0.3265807932 0.3457731317
+DEAL::value -0.5168487451 gradient 0.7285722781 -0.3000861670 0.3281813615
+DEAL::value -0.4535639572 gradient 0.7951660259 -0.2704364567 0.2995682560
+DEAL::value -0.3878670072 gradient 0.8530112273 -0.2369349811 0.2631442345
+DEAL::value -0.3201289971 gradient 0.9015397880 -0.1996127500 0.2211799013
+DEAL::value -0.2507210286 gradient 0.9404568577 -0.1589113348 0.1753423628
+DEAL::value -0.1800142036 gradient 0.9696432327 -0.1154849532 0.1268797053
+DEAL::value -0.1083796240 gradient 0.9890795498 -0.07008403854 0.07674256641
+DEAL::value -0.03618839165 gradient 0.9987871482 -0.02349357912 0.02567941921
+DEAL::value 0.03618839165 gradient 0.9987871482 0.02349357912 -0.02567941921
+DEAL::value 0.1083796240 gradient 0.9890795498 0.07008403854 -0.07674256641
+DEAL::value 0.1800142036 gradient 0.9696432327 0.1154849532 -0.1268797053
+DEAL::value 0.2507210286 gradient 0.9404568577 0.1589113348 -0.1753423628
+DEAL::
+DEAL::-0.03458205383 -0.003056132278 -14.93776624 0.4157039597 0.1068644233 -0.2797746925 0.009443968110 -0.02472462950 -0.05494126836 0.02185965012 -0.7651666446 13.45578877 0.1697778561 -0.06755003380 -0.4444841979 0.1768482844 
+DEAL::Cell with center 0.000000000 1.000000000 0.000000000
+DEAL::face number 0
+DEAL::value -0.5773502692 gradient 0.6542268683 0.3265807932 -0.3457731317
+DEAL::value -0.5168487451 gradient 0.7285722781 0.3000861670 -0.3281813615
+DEAL::value -0.4535639572 gradient 0.7951660259 0.2704364567 -0.2995682560
+DEAL::value -0.3878670072 gradient 0.8530112273 0.2369349811 -0.2631442345
+DEAL::value -0.3201289971 gradient 0.9015397880 0.1996127500 -0.2211799013
+DEAL::value -0.2507210286 gradient 0.9404568577 0.1589113348 -0.1753423628
+DEAL::value -0.1800142036 gradient 0.9696432327 0.1154849532 -0.1268797053
+DEAL::value -0.1083796240 gradient 0.9890795498 0.07008403854 -0.07674256641
+DEAL::value -0.03618839165 gradient 0.9987871482 0.02349357912 -0.02567941921
+DEAL::value 0.03618839165 gradient 0.9987871482 -0.02349357912 0.02567941921
+DEAL::value 0.1083796240 gradient 0.9890795498 -0.07008403854 0.07674256641
+DEAL::value 0.1800142036 gradient 0.9696432327 -0.1154849532 0.1268797053
+DEAL::value 0.2507210286 gradient 0.9404568577 -0.1589113348 0.1753423628
+DEAL::
+DEAL::-14.93776624 -0.03458205383 0.4157039597 -0.003056132278 -0.7651666446 13.45578877 -0.05494126836 0.02185965012 -0.2797746925 0.1068644233 -0.02472462950 0.009443968110 -0.4444841979 0.1697778561 0.1768482844 -0.06755003380 
+DEAL::face number 1
+DEAL::value -0.5773502692 gradient 0.6542268683 0.3265807932 0.3457731317
+DEAL::value -0.5168487451 gradient 0.7285722781 0.3000861670 0.3281813615
+DEAL::value -0.4535639572 gradient 0.7951660259 0.2704364567 0.2995682560
+DEAL::value -0.3878670072 gradient 0.8530112273 0.2369349811 0.2631442345
+DEAL::value -0.3201289971 gradient 0.9015397880 0.1996127500 0.2211799013
+DEAL::value -0.2507210286 gradient 0.9404568577 0.1589113348 0.1753423628
+DEAL::value -0.1800142036 gradient 0.9696432327 0.1154849532 0.1268797053
+DEAL::value -0.1083796240 gradient 0.9890795498 0.07008403854 0.07674256641
+DEAL::value -0.03618839165 gradient 0.9987871482 0.02349357912 0.02567941921
+DEAL::value 0.03618839165 gradient 0.9987871482 -0.02349357912 -0.02567941921
+DEAL::value 0.1083796240 gradient 0.9890795498 -0.07008403854 -0.07674256641
+DEAL::value 0.1800142036 gradient 0.9696432327 -0.1154849532 -0.1268797053
+DEAL::value 0.2507210286 gradient 0.9404568577 -0.1589113348 -0.1753423628
+DEAL::
+DEAL::-0.03458205383 -14.93776624 -0.003056132278 0.4157039597 -0.05494126836 0.02185965012 -0.7651666446 13.45578877 0.1068644233 -0.2797746925 0.009443968110 -0.02472462950 0.1697778561 -0.4444841979 -0.06755003380 0.1768482844 
+DEAL::face number 2
+DEAL::value -0.5773502692 gradient 0.6542268683 0.3265807932 -0.3457731317
+DEAL::value -0.6059008348 gradient 0.6031982095 0.3628320872 -0.3281813615
+DEAL::value -0.6308825797 gradient 0.5618835185 0.3955114285 -0.2995682560
+DEAL::value -0.6522955039 gradient 0.5289103591 0.4241689557 -0.2631442345
+DEAL::value -0.6701396074 gradient 0.5031439836 0.4484078129 -0.2211799013
+DEAL::value -0.6844148903 gradient 0.4836526427 0.4679613443 -0.1753423628
+DEAL::value -0.6951213524 gradient 0.4696912413 0.4826830425 -0.1268797053
+DEAL::value -0.7022589938 gradient 0.4606979203 0.4925098223 -0.07674256641
+DEAL::value -0.7058278145 gradient 0.4562958566 0.4974238789 -0.02567941921
+DEAL::value -0.7058278145 gradient 0.4562958566 0.4974238789 0.02567941921
+DEAL::value -0.7022589938 gradient 0.4606979203 0.4925098223 0.07674256641
+DEAL::value -0.6951213524 gradient 0.4696912413 0.4826830425 0.1268797053
+DEAL::value -0.6844148903 gradient 0.4836526427 0.4679613443 0.1753423628
+DEAL::
+DEAL::-7.069933216 1.627626224 1.135550187 -0.1943807830 9.186793990 -3.509043057 -1.572573568 0.6006696533 -27.71278762 -13.77501731 3.633545959 1.795000808 29.39600430 14.52186159 -11.22827451 -5.546857546 
+DEAL::face number 3
+DEAL::value 0.5773502692 gradient 0.6542268683 -0.3265807932 0.3457731317
+DEAL::value 0.6059008348 gradient 0.6031982095 -0.3628320872 0.3281813615
+DEAL::value 0.6308825797 gradient 0.5618835185 -0.3955114285 0.2995682560
+DEAL::value 0.6522955039 gradient 0.5289103591 -0.4241689557 0.2631442345
+DEAL::value 0.6701396074 gradient 0.5031439836 -0.4484078129 0.2211799013
+DEAL::value 0.6844148903 gradient 0.4836526427 -0.4679613443 0.1753423628
+DEAL::value 0.6951213524 gradient 0.4696912413 -0.4826830425 0.1268797053
+DEAL::value 0.7022589938 gradient 0.4606979203 -0.4925098223 0.07674256641
+DEAL::value 0.7058278145 gradient 0.4562958566 -0.4974238789 0.02567941921
+DEAL::value 0.7058278145 gradient 0.4562958566 -0.4974238789 -0.02567941921
+DEAL::value 0.7022589938 gradient 0.4606979203 -0.4925098223 -0.07674256641
+DEAL::value 0.6951213524 gradient 0.4696912413 -0.4826830425 -0.1268797053
+DEAL::value 0.6844148903 gradient 0.4836526427 -0.4679613443 -0.1753423628
+DEAL::
+DEAL::-1.135550187 0.1943807830 7.069933216 -1.627626224 3.509043057 -9.186793990 -0.6006696533 1.572573568 -3.633545959 -1.795000808 27.71278762 13.77501731 11.22827451 5.546857546 -29.39600430 -14.52186159 

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.