Add test

diff --git a/tests/matrix_free/point_evaluation_35.cc b/tests/matrix_free/point_evaluation_35.cc
new file mode 100644 (file)
index 0000000..b5cccad
--- /dev/null
+++ b/tests/matrix_free/point_evaluation_35.cc
@@ -0,0 +1,137 @@
+// ------------------------------------------------------------------------
+//
+// SPDX-License-Identifier: LGPL-2.1-or-later
+// Copyright (C) 2020 - 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 vector-valued FE_Q and MappingQ by
+// comparing to the output of FEValues with the same settings for curl
+
+#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>
+void
+test(const unsigned int degree)
+{
+  Triangulation<dim> tria;
+
+  if (dim > 1)
+    GridGenerator::hyper_shell(tria, Point<dim>(), 0.5, 1, 6);
+  else
+    GridGenerator::subdivided_hyper_cube(tria, 2, 0, 1);
+
+  MappingQ<dim> 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> fe(FE_Q<dim>(degree), dim);
+  FEValues<dim> fe_values(mapping,
+                          fe,
+                          Quadrature<dim>(unit_points),
+                          update_values | update_gradients);
+
+  DoFHandler<dim> dof_handler(tria);
+  dof_handler.distribute_dofs(fe);
+  Vector<double> vector(dof_handler.n_dofs());
+
+  FEPointEvaluation<dim, dim> evaluator(mapping,
+                                        fe,
+                                        update_values | update_gradients);
+
+  VectorTools::interpolate(mapping, dof_handler, MyFunction<dim>(), vector);
+  const FEValuesExtractors::Vector extractor(0);
+
+  std::vector<double> solution_values(fe.dofs_per_cell);
+  std::vector<typename FEValuesViews::Vector<dim>::curl_type> function_curls(
+    unit_points.size());
+
+  for (const auto &cell : dof_handler.active_cell_iterators())
+    {
+      fe_values.reinit(cell);
+      fe_values[extractor].get_function_curls(vector, function_curls);
+
+      cell->get_dof_values(vector,
+                           solution_values.begin(),
+                           solution_values.end());
+
+      evaluator.reinit(cell, unit_points);
+      evaluator.evaluate(solution_values, EvaluationFlags::gradients);
+
+      deallog << "Cell with center " << cell->center(true) << std::endl;
+      for (unsigned int i = 0; i < function_curls.size(); ++i)
+        deallog << mapping.transform_unit_to_real_cell(cell, unit_points[i])
+                << ": "
+                << " error curl " << (evaluator.get_curl(i) - function_curls[i])
+                << std::endl;
+      deallog << std::endl;
+    }
+}
+
+
+
+int
+main()
+{
+  initlog();
+  deallog << std::setprecision(10);
+
+  test<2>(2);
+  test<2>(6);
+  test<3>(5);
+}

diff --git a/tests/matrix_free/point_evaluation_35.output b/tests/matrix_free/point_evaluation_35.output
new file mode 100644 (file)
index 0000000..2d3c796
--- /dev/null
+++ b/tests/matrix_free/point_evaluation_35.output
@@ -0,0 +1,274 @@
+
+DEAL::Mapping of degree 2
+DEAL::Cell with center 0.6495190528 0.3750000000
+DEAL::0.9921875000 0.000000000:  error curl 1.970645869e-15
+DEAL::0.9591965610 0.06332871922:  error curl 2.914335440e-16
+DEAL::0.9234265460 0.1210574346:  error curl 1.026956298e-15
+DEAL::0.8851608469 0.1733497625:  error curl -1.110223025e-16
+DEAL::0.8446828558 0.2203693193:  error curl 4.440892099e-16
+DEAL::0.8022759646 0.2622797217:  error curl 4.996003611e-16
+DEAL::0.7582235654 0.2992445859:  error curl 0.000000000
+DEAL::0.7128090500 0.3314275284:  error curl -2.775557562e-16
+DEAL::0.6663158106 0.3589921657:  error curl 3.330669074e-16
+DEAL::0.6190272391 0.3821021142:  error curl 9.436895709e-16
+DEAL::0.5712267275 0.4009209903:  error curl -7.771561172e-16
+DEAL::0.5231976677 0.4156124105:  error curl -5.551115123e-16
+DEAL::0.4752234518 0.4263399913:  error curl -9.992007222e-16
+DEAL::
+DEAL::Cell with center 4.592425497e-17 0.7500000000
+DEAL::0.4960937500 0.8592595803:  error curl 3.330669074e-16
+DEAL::0.4247540009 0.8623529487:  error curl -4.996003611e-16
+DEAL::0.3568744593 0.8602395646:  error curl 8.326672685e-16
+DEAL::0.2924551254 0.8532466611:  error curl -3.330669074e-16
+DEAL::0.2314959991 0.8417014709:  error curl -6.106226635e-16
+DEAL::0.1739970804 0.8259312270:  error curl -9.436895709e-16
+DEAL::0.1199583694 0.8062631623:  error curl -3.885780586e-16
+DEAL::0.06937986592 0.7830245096:  error curl -3.191891196e-16
+DEAL::0.02226157007 0.7565425018:  error curl 1.734723476e-17
+DEAL::-0.02139651817 0.7271443718:  error curl 3.816391647e-17
+DEAL::-0.06159439879 0.6951573525:  error curl 2.914335440e-16
+DEAL::-0.09833207180 0.6609086767:  error curl 1.387778781e-16
+DEAL::-0.1316095372 0.6247255774:  error curl -4.163336342e-16
+DEAL::
+DEAL::Cell with center -0.6495190528 0.3750000000
+DEAL::-0.4960937500 0.8592595803:  error curl 4.440892099e-16
+DEAL::-0.5344425601 0.7990242294:  error curl 9.992007222e-16
+DEAL::-0.5665520866 0.7391821301:  error curl -1.054711873e-15
+DEAL::-0.5927057215 0.6798968986:  error curl -1.387778781e-15
+DEAL::-0.6131868567 0.6213321516:  error curl 9.436895709e-16
+DEAL::-0.6282788842 0.5636515054:  error curl 8.326672685e-16
+DEAL::-0.6382651960 0.5070185764:  error curl 4.440892099e-16
+DEAL::-0.6434291841 0.4515969812:  error curl -3.330669074e-16
+DEAL::-0.6440542406 0.3975503361:  error curl -7.632783294e-17
+DEAL::-0.6404237573 0.3450422576:  error curl 5.551115123e-16
+DEAL::-0.6328211262 0.2942363621:  error curl 2.220446049e-16
+DEAL::-0.6215297395 0.2452962662:  error curl 1.498801083e-15
+DEAL::-0.6068329890 0.1983855861:  error curl 1.776356839e-15
+DEAL::
+DEAL::Cell with center -0.6495190528 -0.3750000000
+DEAL::-0.9921875000 1.215079246e-16:  error curl 1.610251688e-15
+DEAL::-0.9591965610 -0.06332871922:  error curl 4.163336342e-16
+DEAL::-0.9234265460 -0.1210574346:  error curl 2.498001805e-16
+DEAL::-0.8851608469 -0.1733497625:  error curl 8.326672685e-17
+DEAL::-0.8446828558 -0.2203693193:  error curl -9.992007222e-16
+DEAL::-0.8022759646 -0.2622797217:  error curl -8.326672685e-16
+DEAL::-0.7582235654 -0.2992445859:  error curl -1.110223025e-16
+DEAL::-0.7128090500 -0.3314275284:  error curl -6.661338148e-16
+DEAL::-0.6663158106 -0.3589921657:  error curl -2.775557562e-17
+DEAL::-0.6190272391 -0.3821021142:  error curl -1.665334537e-16
+DEAL::-0.5712267275 -0.4009209903:  error curl 0.000000000
+DEAL::-0.5231976677 -0.4156124105:  error curl 1.110223025e-16
+DEAL::-0.4752234518 -0.4263399913:  error curl 9.992007222e-16
+DEAL::
+DEAL::Cell with center -1.377727649e-16 -0.7500000000
+DEAL::-0.4960937500 -0.8592595803:  error curl -1.526556659e-16
+DEAL::-0.4247540009 -0.8623529487:  error curl -4.440892099e-16
+DEAL::-0.3568744593 -0.8602395646:  error curl 8.326672685e-16
+DEAL::-0.2924551254 -0.8532466611:  error curl -2.220446049e-16
+DEAL::-0.2314959991 -0.8417014709:  error curl -3.330669074e-16
+DEAL::-0.1739970804 -0.8259312270:  error curl -1.387778781e-16
+DEAL::-0.1199583694 -0.8062631623:  error curl -1.387778781e-16
+DEAL::-0.06937986592 -0.7830245096:  error curl -1.387778781e-17
+DEAL::-0.02226157007 -0.7565425018:  error curl 1.665334537e-16
+DEAL::0.02139651817 -0.7271443718:  error curl -9.020562075e-17
+DEAL::0.06159439879 -0.6951573525:  error curl -6.938893904e-17
+DEAL::0.09833207180 -0.6609086767:  error curl 1.387778781e-16
+DEAL::0.1316095372 -0.6247255774:  error curl 0.000000000
+DEAL::
+DEAL::Cell with center 0.6495190528 -0.3750000000
+DEAL::0.4960937500 -0.8592595803:  error curl 7.494005416e-16
+DEAL::0.5344425601 -0.7990242294:  error curl 4.718447855e-16
+DEAL::0.5665520866 -0.7391821301:  error curl -1.082467449e-15
+DEAL::0.5927057215 -0.6798968986:  error curl -3.885780586e-16
+DEAL::0.6131868567 -0.6213321516:  error curl 4.996003611e-16
+DEAL::0.6282788842 -0.5636515054:  error curl 2.220446049e-16
+DEAL::0.6382651960 -0.5070185764:  error curl 0.000000000
+DEAL::0.6434291841 -0.4515969812:  error curl 3.330669074e-16
+DEAL::0.6440542406 -0.3975503361:  error curl -1.196959198e-16
+DEAL::0.6404237573 -0.3450422576:  error curl 2.220446049e-16
+DEAL::0.6328211262 -0.2942363621:  error curl 0.000000000
+DEAL::0.6215297395 -0.2452962662:  error curl -8.326672685e-16
+DEAL::0.6068329890 -0.1983855861:  error curl 1.221245327e-15
+DEAL::
+DEAL::Mapping of degree 6
+DEAL::Cell with center 0.6495190528 0.3750000000
+DEAL::0.9921875000 0.000000000:  error curl 8.318619625e-15
+DEAL::0.9609496573 0.05926936446:  error curl -5.759281940e-15
+DEAL::0.9262895583 0.1146994816:  error curl -3.594347042e-15
+DEAL::0.8885607198 0.1661007475:  error curl 1.715294573e-14
+DEAL::0.8481266485 0.2133123969:  error curl 3.441691376e-15
+DEAL::0.8053586039 0.2562030323:  error curl 6.328271240e-15
+DEAL::0.7606333492 0.2946709691:  error curl 2.081668171e-15
+DEAL::0.7143309018 0.3286444014:  error curl 9.436895709e-16
+DEAL::0.6668322927 0.3580813945:  error curl -2.303712776e-15
+DEAL::0.6185173440 0.3829697098:  error curl 3.885780586e-16
+DEAL::0.5697624746 0.4033264664:  error curl -1.276756478e-15
+DEAL::0.5209385434 0.4191976473:  error curl -1.720845688e-15
+DEAL::0.4724087397 0.4306574530:  error curl -5.551115123e-16
+DEAL::
+DEAL::Cell with center 4.592425497e-17 0.7500000000
+DEAL::0.4960937500 0.8592595803:  error curl -2.622901896e-15
+DEAL::0.4291460533 0.8618414972:  error curl 1.915134717e-15
+DEAL::0.3638121143 0.8595400296:  error curl -2.400857291e-15
+DEAL::0.3004328929 0.8525665299:  error curl -2.359223927e-15
+DEAL::0.2393293696 0.8411554217:  error curl -8.465450563e-16
+DEAL::0.1808009674 0.8255625263:  error curl -2.498001805e-15
+DEAL::0.1251241296 0.8060632879:  error curl -9.992007222e-16
+DEAL::0.07255105048 0.7829509083:  error curl -3.580469254e-15
+DEAL::0.02330856204 0.7565344028:  error curl -9.020562075e-17
+DEAL::-0.02240282557 0.7271365875:  error curl -2.185751580e-16
+DEAL::-0.06440972863 0.6950920104:  error curl -1.568190022e-15
+DEAL::-0.1025665400 0.6607448360:  error curl 1.498801083e-15
+DEAL::-0.1367559248 0.6244466961:  error curl 1.720845688e-15
+DEAL::
+DEAL::Cell with center -0.6495190528 0.3750000000
+DEAL::-0.4960937500 0.8592595803:  error curl -9.520162436e-15
+DEAL::-0.5318036039 0.8025721327:  error curl -1.192795862e-14
+DEAL::-0.5624774440 0.7448405480:  error curl -3.454701802e-15
+DEAL::-0.5881278268 0.6864657824:  error curl -6.904199434e-16
+DEAL::-0.6087972789 0.6278430248:  error curl -2.484124018e-15
+DEAL::-0.6245576365 0.5693594940:  error curl 1.875236078e-15
+DEAL::-0.6355092196 0.5113923188:  error curl 3.741798538e-15
+DEAL::-0.6417798513 0.4543065070:  error curl -1.276756478e-15
+DEAL::-0.6435237306 0.3984530082:  error curl 2.969846591e-15
+DEAL::-0.6409201696 0.3441668777:  error curl -2.359223927e-15
+DEAL::-0.6341722033 0.2917655439:  error curl -2.775557562e-16
+DEAL::-0.6235050834 0.2415471888:  error curl 3.497202528e-15
+DEAL::-0.6091646645 0.1937892431:  error curl 5.828670879e-15
+DEAL::
+DEAL::Cell with center -0.6495190528 -0.3750000000
+DEAL::-0.9921875000 1.215079246e-16:  error curl 3.620703861e-15
+DEAL::-0.9609496573 -0.05926936446:  error curl -1.942890293e-16
+DEAL::-0.9262895583 -0.1146994816:  error curl -2.532696275e-15
+DEAL::-0.8885607198 -0.1661007475:  error curl -6.002143227e-16
+DEAL::-0.8481266485 -0.2133123969:  error curl 8.881784197e-16
+DEAL::-0.8053586039 -0.2562030323:  error curl -2.289834988e-16
+DEAL::-0.7606333492 -0.2946709691:  error curl -9.089951014e-16
+DEAL::-0.7143309018 -0.3286444014:  error curl -1.845745778e-15
+DEAL::-0.6668322927 -0.3580813945:  error curl -6.245004514e-16
+DEAL::-0.6185173440 -0.3829697098:  error curl 2.386979503e-15
+DEAL::-0.5697624746 -0.4033264664:  error curl -1.110223025e-15
+DEAL::-0.5209385434 -0.4191976473:  error curl -7.216449660e-16
+DEAL::-0.4724087397 -0.4306574530:  error curl -3.885780586e-16
+DEAL::
+DEAL::Cell with center -1.377727649e-16 -0.7500000000
+DEAL::-0.4960937500 -0.8592595803:  error curl -3.920475056e-15
+DEAL::-0.4291460533 -0.8618414972:  error curl -5.898059818e-17
+DEAL::-0.3638121143 -0.8595400296:  error curl 6.938893904e-18
+DEAL::-0.3004328929 -0.8525665299:  error curl -5.256212132e-16
+DEAL::-0.2393293696 -0.8411554217:  error curl 1.769417945e-16
+DEAL::-0.1808009674 -0.8255625263:  error curl 5.412337245e-16
+DEAL::-0.1251241296 -0.8060632879:  error curl -6.938893904e-17
+DEAL::-0.07255105048 -0.7829509083:  error curl -3.122502257e-16
+DEAL::-0.02330856204 -0.7565344028:  error curl -4.857225733e-17
+DEAL::0.02240282557 -0.7271365875:  error curl 7.285838599e-17
+DEAL::0.06440972863 -0.6950920104:  error curl 1.110223025e-16
+DEAL::0.1025665400 -0.6607448360:  error curl 3.053113318e-16
+DEAL::0.1367559248 -0.6244466961:  error curl 6.175615574e-16
+DEAL::
+DEAL::Cell with center 0.6495190528 -0.3750000000
+DEAL::0.4960937500 -0.8592595803:  error curl -1.533495553e-15
+DEAL::0.5318036039 -0.8025721327:  error curl -2.151057110e-16
+DEAL::0.5624774440 -0.7448405480:  error curl 3.605622745e-15
+DEAL::0.5881278268 -0.6864657824:  error curl 2.151057110e-15
+DEAL::0.6087972789 -0.6278430248:  error curl -6.591949209e-16
+DEAL::0.6245576365 -0.5693594940:  error curl 1.497066360e-15
+DEAL::0.6355092196 -0.5113923188:  error curl 1.417269080e-15
+DEAL::0.6417798513 -0.4543065070:  error curl 0.000000000
+DEAL::0.6435237306 -0.3984530082:  error curl -1.387778781e-17
+DEAL::0.6409201696 -0.3441668777:  error curl -7.494005416e-16
+DEAL::0.6341722033 -0.2917655439:  error curl -2.803313137e-15
+DEAL::0.6235050834 -0.2415471888:  error curl 3.275157923e-15
+DEAL::0.6091646645 -0.1937892431:  error curl 1.887379142e-15
+DEAL::
+DEAL::Mapping of degree 5
+DEAL::Cell with center 0.000000000 0.000000000 -0.7500000000
+DEAL::-0.5759520655 -0.5527659841 -0.5759520655:  error curl -2.396520482e-15 -2.867497906e-15 -2.870967353e-16
+DEAL::-0.5254898676 -0.5029040719 -0.6187865537:  error curl 6.656133977e-15 -9.063930162e-16 -6.678685383e-15
+DEAL::-0.4685426073 -0.4464599093 -0.6616551149:  error curl 1.026803832e-15 -2.722865336e-15 -1.922832553e-15
+DEAL::-0.4051542510 -0.3835292107 -0.7012978956:  error curl 7.540842950e-15 -1.366962099e-15 4.041905699e-16
+DEAL::-0.3359932318 -0.3148223149 -0.7343028440:  error curl -4.121702979e-15 2.662800536e-16 2.477185124e-15
+DEAL::-0.2623766962 -0.2416891159 -0.7575161861:  error curl 2.185751580e-16 -1.802377692e-15 -5.655198532e-16
+DEAL::-0.1861925434 -0.1660412603 -0.7684224851:  error curl 6.582625070e-15 -9.410874857e-16 8.001412033e-17
+DEAL::-0.1097387169 -0.09019150189 -0.7654365722:  error curl 1.040834086e-15 -2.234323837e-15 1.075528555e-15
+DEAL::-0.03550798632 -0.01663873419 -0.7480737719:  error curl 2.928213227e-15 -7.285838599e-16 9.488937414e-16
+DEAL::0.03404970570 0.05216906594 -0.7169856075:  error curl -1.543903894e-15 1.619364365e-15 -4.148157512e-16
+DEAL::0.09675558003 0.1140621810 -0.6738655685:  error curl 1.491862189e-15 -9.844555726e-16 -2.921924855e-16
+DEAL::0.1508930055 0.1673387801 -0.6212435337:  error curl -5.641320744e-15 1.179611964e-16 -8.604228441e-16
+DEAL::0.1953394317 0.2108952300 -0.5621980946:  error curl 3.053113318e-15 -9.992007222e-16 -1.401656569e-15
+DEAL::
+DEAL::Cell with center 0.7500000000 0.000000000 -4.163336342e-17
+DEAL::0.5879922187 -0.5879922187 -0.5412602679:  error curl -1.085069534e-15 4.478188653e-15 -4.575333168e-15
+DEAL::0.6317656581 -0.5359927398 -0.4904247681:  error curl -1.266955291e-14 1.231783772e-14 -5.610746243e-17
+DEAL::0.6752050943 -0.4774156808 -0.4328300536:  error curl -1.442639411e-15 -3.407864269e-15 -2.994457980e-15
+DEAL::0.7149693757 -0.4123477836 -0.3686646352:  error curl -2.730454751e-15 6.980527267e-15 -5.290906602e-16
+DEAL::0.7476112822 -0.3415220448 -0.2987437738:  error curl -1.492079030e-15 2.002738253e-15 1.501836849e-15
+DEAL::0.7699931685 -0.2663289165 -0.2245218407:  error curl 1.595945598e-16 4.466912951e-15 -3.756543687e-15
+DEAL::0.7796633598 -0.1887262238 -0.1480024603:  error curl -4.943961907e-16 -1.864394056e-15 3.255859124e-15
+DEAL::0.7751382896 -0.1110685560 -0.07156802410:  error curl 1.524821935e-15 8.187894807e-16 -1.321046137e-15
+DEAL::0.7560592058 -0.03588516627 0.002241853775:  error curl 9.142263769e-16 -3.044006019e-15 2.420733686e-15
+DEAL::0.7232127397 0.03436118254 0.07097021235:  error curl 2.561210792e-15 2.067790383e-15 -1.703579769e-15
+DEAL::0.6784217227 0.09750269006 0.1324691326:  error curl -2.021820211e-15 -3.330669074e-15 -2.950330952e-15
+DEAL::0.6243263560 0.1518546010 0.1850850191:  error curl 6.229392002e-15 1.171979180e-14 2.811119393e-15
+DEAL::0.5640861853 0.1963417794 0.2277818661:  error curl 4.010680676e-15 2.664535259e-15 4.842480583e-15
+DEAL::
+DEAL::Cell with center 0.000000000 -4.163336342e-17 0.7500000000
+DEAL::-0.5879922187 -0.5412602679 0.5879922187:  error curl -2.003345406e-14 9.534040224e-15 7.277164982e-16
+DEAL::-0.5359927398 -0.4904247681 0.6317656581:  error curl 1.081426615e-14 2.357402468e-14 -4.208872834e-15
+DEAL::-0.4774156808 -0.4328300536 0.6752050943:  error curl -1.174916013e-15 4.354589606e-15 3.986611388e-15
+DEAL::-0.4123477836 -0.3686646352 0.7149693757:  error curl -5.342948306e-16 -1.226449498e-15 -3.623837341e-15
+DEAL::-0.3415220448 -0.2987437738 0.7476112822:  error curl 3.663735981e-15 5.322131624e-15 -4.319461455e-16
+DEAL::-0.2663289165 -0.2245218407 0.7699931685:  error curl 7.136652380e-15 2.513614317e-15 1.905593738e-15
+DEAL::-0.1887262238 -0.1480024603 0.7796633598:  error curl -4.495427468e-15 -4.281080698e-15 -1.428544782e-15
+DEAL::-0.1110685560 -0.07156802410 0.7751382896:  error curl 2.539635169e-15 6.331740687e-17 2.541369892e-16
+DEAL::-0.03588516627 0.002241853775 0.7560592058:  error curl -8.779001831e-15 -1.442532168e-15 5.395396586e-16
+DEAL::0.03436118254 0.07097021235 0.7232127397:  error curl -4.121702979e-15 2.643962523e-16 1.373467312e-15
+DEAL::0.09750269006 0.1324691326 0.6784217227:  error curl -3.275157923e-15 3.342161617e-15 3.773023560e-16
+DEAL::0.1518546010 0.1850850191 0.6243263560:  error curl 2.485425060e-15 -4.878909776e-15 1.785030457e-15
+DEAL::0.1963417794 0.2277818661 0.5640861853:  error curl -9.436895709e-16 -2.107689023e-15 3.851086117e-16
+DEAL::
+DEAL::Cell with center -0.7500000000 0.000000000 0.000000000
+DEAL::-0.6012515503 -0.5766780283 -0.5531118850:  error curl 1.005966144e-14 -9.065664885e-15 -1.131039706e-15
+DEAL::-0.6458946236 -0.5236491363 -0.5006487306:  error curl 1.975416358e-15 -8.640223953e-15 9.182975561e-15
+DEAL::-0.6898189280 -0.4638689601 -0.4413495827:  error curl -4.548444954e-15 -2.452898995e-15 6.098420380e-15
+DEAL::-0.7296048967 -0.3975220470 -0.3754490693:  error curl -5.342948306e-16 4.052314040e-15 2.376571162e-16
+DEAL::-0.7617762381 -0.3254496235 -0.3038288032:  error curl -2.896554524e-15 1.283695372e-15 -8.838416110e-16
+DEAL::-0.7832188522 -0.2491480826 -0.2280163559:  error curl 1.894318036e-15 -2.554380318e-15 1.619798046e-15
+DEAL::-0.7915521155 -0.1706669598 -0.1500833118:  error curl 1.321588238e-15 6.893791094e-15 5.594483210e-17
+DEAL::-0.7854001879 -0.09242877894 -0.07246518026:  error curl -1.348747503e-15 1.540434447e-15 -2.815320676e-15
+DEAL::-0.7645345314 -0.01700078778 0.002266559507:  error curl -4.265776507e-16 -1.958936485e-15 1.352642001e-15
+DEAL::-0.7298790026 0.05314854421 0.07164727507:  error curl -9.647230931e-16 3.330669074e-15 -5.224499219e-16
+DEAL::-0.6833856700 0.1158765656 0.1335451305:  error curl -1.172673070e-15 3.996802889e-15 -1.229268423e-15
+DEAL::-0.6278029298 0.1695497377 0.1863429372:  error curl -9.818534874e-16 2.137179322e-15 3.504575102e-15
+DEAL::-0.5663675440 0.2131627953 0.2290556314:  error curl -1.675742878e-15 -4.829470157e-15 -5.734128450e-15
+DEAL::
+DEAL::Cell with center 0.000000000 -0.7500000000 -4.163336342e-17
+DEAL::-0.5879922187 -0.5879922187 -0.5412602679:  error curl -3.737461729e-15 7.163540594e-15 4.051880359e-15
+DEAL::-0.5359927398 -0.6317656581 -0.4904247681:  error curl 4.904171687e-15 7.656202061e-15 -3.730956516e-15
+DEAL::-0.4774156808 -0.6752050943 -0.4328300536:  error curl 6.158268340e-17 8.364836601e-15 -1.478147032e-15
+DEAL::-0.4123477836 -0.7149693757 -0.3686646352:  error curl 5.093148125e-15 1.863093013e-15 1.518750403e-15
+DEAL::-0.3415220448 -0.7476112822 -0.2987437738:  error curl -8.283304598e-16 -1.892800153e-15 1.602125550e-15
+DEAL::-0.2663289165 -0.7699931685 -0.2245218407:  error curl 1.010129480e-14 5.516420654e-16 1.932915633e-15
+DEAL::-0.1887262238 -0.7796633598 -0.1480024603:  error curl 3.438221929e-15 1.369130503e-15 -2.772304955e-16
+DEAL::-0.1110685560 -0.7751382896 -0.07156802410:  error curl 4.926614672e-15 -2.498001805e-15 -2.839525490e-16
+DEAL::-0.03588516627 -0.7560592058 0.002241853775:  error curl -2.796807924e-15 4.975573176e-16 3.128066866e-16
+DEAL::0.03436118254 -0.7232127397 0.07097021235:  error curl 2.872702076e-15 -6.374024572e-16 2.511012231e-16
+DEAL::0.09750269006 -0.6784217227 0.1324691326:  error curl -1.387778781e-15 1.662732452e-15 -1.154024792e-15
+DEAL::0.1518546010 -0.6243263560 0.1850850191:  error curl 1.711131237e-14 -8.621575676e-16 -3.157196726e-16
+DEAL::0.1963417794 -0.5640861853 0.2277818661:  error curl -2.775557562e-16 -2.289834988e-15 5.048045315e-16
+DEAL::
+DEAL::Cell with center 0.000000000 0.7500000000 0.000000000
+DEAL::-0.5766780283 0.6012515503 -0.5531118850:  error curl -4.690692279e-15 2.860385540e-14 -5.568462358e-16
+DEAL::-0.5236491363 0.6458946236 -0.5006487306:  error curl -1.453134488e-14 1.591175108e-15 -8.960280434e-15
+DEAL::-0.4638689601 0.6898189280 -0.4413495827:  error curl -1.539957398e-14 -1.934216676e-16 1.350265386e-15
+DEAL::-0.3975220470 0.7296048967 -0.3754490693:  error curl 6.685624276e-15 1.498801083e-15 3.979455654e-15
+DEAL::-0.3254496235 0.7617762381 -0.3038288032:  error curl 2.029626467e-16 -1.096778918e-15 2.881484114e-15
+DEAL::-0.2491480826 0.7832188522 -0.2280163559:  error curl -1.874368716e-15 -1.682681772e-16 1.754672796e-15
+DEAL::-0.1706669598 0.7915521155 -0.1500833118:  error curl -3.375771884e-15 -1.998130394e-15 -1.943432394e-17
+DEAL::-0.09242877894 0.7854001879 -0.07246518026:  error curl 1.845745778e-15 1.824495416e-15 2.299321757e-16
+DEAL::-0.01700078778 0.7645345314 0.002266559507:  error curl -3.115563363e-15 -1.075684409e-15 -2.299347486e-15
+DEAL::0.05314854421 0.7298790026 0.07164727507:  error curl 7.216449660e-16 1.392766111e-15 -2.396086801e-17
+DEAL::0.1158765656 0.6833856700 0.1335451305:  error curl 7.382983114e-15 1.173540431e-15 -1.338556002e-15
+DEAL::0.1695497377 0.6278029298 0.1863429372:  error curl 1.332267630e-15 -2.619432449e-16 -2.034830637e-15
+DEAL::0.2131627953 0.5663675440 0.2290556314:  error curl 2.275957200e-15 1.453698273e-15 -6.192962809e-16
+DEAL::