--- /dev/null
+// ---------------------------------------------------------------------
+//
+// Copyright (C) 2021 - 2022 by the deal.II authors
+//
+// This file is part of the deal.II library.
+//
+// The deal.II library is free software; you can use it, redistribute
+// it, and/or modify it under the terms of the GNU Lesser General
+// Public License as published by the Free Software Foundation; either
+// version 2.1 of the License, or (at your option) any later version.
+// The full text of the license can be found in the file LICENSE.md at
+// the top level directory of deal.II.
+//
+// ---------------------------------------------------------------------
+
+
+// Check FEPointEvaluation::evaluate()/integrate() with precomputed mapping
+// with different location of unit points provided my NonMatching::MappingInfo.
+
+#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_cartesian.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()
+{
+ using namespace dealii;
+ Triangulation<dim> tria;
+
+ GridGenerator::subdivided_hyper_cube(tria, 2, 0, 1);
+
+ MappingCartesian<dim> mapping;
+ deallog << "Cartesian linear mapping" << 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);
+ }
+
+ std::vector<Point<dim>> unit_points2;
+ 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) / 19. + 0.015625 * d;
+ unit_points2.push_back(p);
+ }
+
+ FESystem<dim> fe(FE_Q<dim>(2), dim);
+ DoFHandler<dim> dof_handler(tria);
+ dof_handler.distribute_dofs(fe);
+ Vector<double> vector(dof_handler.n_dofs());
+
+ // mapping info object for precomputed mapping
+ NonMatching::MappingInfo<dim, dim> mapping_info(mapping,
+ update_values |
+ update_gradients);
+ // First evaluator to use precomputed mapping
+ FEPointEvaluation<dim, dim> evaluator1(mapping_info, fe);
+
+ VectorTools::interpolate(mapping, dof_handler, MyFunction<dim>(), vector);
+
+ std::vector<double> solution_values(fe.dofs_per_cell);
+
+ for (const auto &cell : dof_handler.active_cell_iterators())
+ {
+ cell->get_dof_values(vector,
+ solution_values.begin(),
+ solution_values.end());
+
+ mapping_info.reinit(cell, unit_points);
+
+ evaluator1.evaluate(solution_values,
+ EvaluationFlags::values | EvaluationFlags::gradients);
+
+ deallog << "Cell with center " << cell->center(true) << std::endl;
+ for (unsigned int i = 0; i < unit_points.size(); ++i)
+ deallog << mapping.transform_unit_to_real_cell(cell, unit_points[i])
+ << ": " << evaluator1.get_value(i) << std::endl;
+ deallog << std::endl;
+
+ mapping_info.reinit(cell, unit_points2);
+
+ evaluator1.evaluate(solution_values,
+ EvaluationFlags::values | EvaluationFlags::gradients);
+
+ deallog << "Cell with center " << cell->center(true) << std::endl;
+ for (unsigned int i = 0; i < unit_points2.size(); ++i)
+ deallog << mapping.transform_unit_to_real_cell(cell, unit_points2[i])
+ << ": " << evaluator1.get_value(i) << std::endl;
+ deallog << std::endl;
+ }
+}
+
+
+
+int
+main()
+{
+ initlog();
+ deallog << std::setprecision(10);
+
+ test<2>();
+ test<3>();
+}
--- /dev/null
+
+DEAL::Cartesian linear mapping
+DEAL::Cell with center 0.2500000000 0.2500000000
+DEAL::0.000000000 0.007812500000: 0.000000000 1.007812500
+DEAL::0.02941176471 0.03722426471: 0.02941176471 1.037224265
+DEAL::0.05882352941 0.06663602941: 0.05882352941 1.066636029
+DEAL::0.08823529412 0.09604779412: 0.08823529412 1.096047794
+DEAL::0.1176470588 0.1254595588: 0.1176470588 1.125459559
+DEAL::0.1470588235 0.1548713235: 0.1470588235 1.154871324
+DEAL::0.1764705882 0.1842830882: 0.1764705882 1.184283088
+DEAL::0.2058823529 0.2136948529: 0.2058823529 1.213694853
+DEAL::0.2352941176 0.2431066176: 0.2352941176 1.243106618
+DEAL::0.2647058824 0.2725183824: 0.2647058824 1.272518382
+DEAL::0.2941176471 0.3019301471: 0.2941176471 1.301930147
+DEAL::0.3235294118 0.3313419118: 0.3235294118 1.331341912
+DEAL::0.3529411765 0.3607536765: 0.3529411765 1.360753676
+DEAL::
+DEAL::Cell with center 0.2500000000 0.2500000000
+DEAL::0.000000000 0.007812500000: 0.000000000 1.007812500
+DEAL::0.02631578947 0.03412828947: 0.02631578947 1.034128289
+DEAL::0.05263157895 0.06044407895: 0.05263157895 1.060444079
+DEAL::0.07894736842 0.08675986842: 0.07894736842 1.086759868
+DEAL::0.1052631579 0.1130756579: 0.1052631579 1.113075658
+DEAL::0.1315789474 0.1393914474: 0.1315789474 1.139391447
+DEAL::0.1578947368 0.1657072368: 0.1578947368 1.165707237
+DEAL::0.1842105263 0.1920230263: 0.1842105263 1.192023026
+DEAL::0.2105263158 0.2183388158: 0.2105263158 1.218338816
+DEAL::0.2368421053 0.2446546053: 0.2368421053 1.244654605
+DEAL::0.2631578947 0.2709703947: 0.2631578947 1.270970395
+DEAL::0.2894736842 0.2972861842: 0.2894736842 1.297286184
+DEAL::0.3157894737 0.3236019737: 0.3157894737 1.323601974
+DEAL::
+DEAL::Cell with center 0.7500000000 0.2500000000
+DEAL::0.5000000000 0.007812500000: 0.5000000000 1.007812500
+DEAL::0.5294117647 0.03722426471: 0.5294117647 1.037224265
+DEAL::0.5588235294 0.06663602941: 0.5588235294 1.066636029
+DEAL::0.5882352941 0.09604779412: 0.5882352941 1.096047794
+DEAL::0.6176470588 0.1254595588: 0.6176470588 1.125459559
+DEAL::0.6470588235 0.1548713235: 0.6470588235 1.154871324
+DEAL::0.6764705882 0.1842830882: 0.6764705882 1.184283088
+DEAL::0.7058823529 0.2136948529: 0.7058823529 1.213694853
+DEAL::0.7352941176 0.2431066176: 0.7352941176 1.243106618
+DEAL::0.7647058824 0.2725183824: 0.7647058824 1.272518382
+DEAL::0.7941176471 0.3019301471: 0.7941176471 1.301930147
+DEAL::0.8235294118 0.3313419118: 0.8235294118 1.331341912
+DEAL::0.8529411765 0.3607536765: 0.8529411765 1.360753676
+DEAL::
+DEAL::Cell with center 0.7500000000 0.2500000000
+DEAL::0.5000000000 0.007812500000: 0.5000000000 1.007812500
+DEAL::0.5263157895 0.03412828947: 0.5263157895 1.034128289
+DEAL::0.5526315789 0.06044407895: 0.5526315789 1.060444079
+DEAL::0.5789473684 0.08675986842: 0.5789473684 1.086759868
+DEAL::0.6052631579 0.1130756579: 0.6052631579 1.113075658
+DEAL::0.6315789474 0.1393914474: 0.6315789474 1.139391447
+DEAL::0.6578947368 0.1657072368: 0.6578947368 1.165707237
+DEAL::0.6842105263 0.1920230263: 0.6842105263 1.192023026
+DEAL::0.7105263158 0.2183388158: 0.7105263158 1.218338816
+DEAL::0.7368421053 0.2446546053: 0.7368421053 1.244654605
+DEAL::0.7631578947 0.2709703947: 0.7631578947 1.270970395
+DEAL::0.7894736842 0.2972861842: 0.7894736842 1.297286184
+DEAL::0.8157894737 0.3236019737: 0.8157894737 1.323601974
+DEAL::
+DEAL::Cell with center 0.2500000000 0.7500000000
+DEAL::0.000000000 0.5078125000: 0.000000000 1.507812500
+DEAL::0.02941176471 0.5372242647: 0.02941176471 1.537224265
+DEAL::0.05882352941 0.5666360294: 0.05882352941 1.566636029
+DEAL::0.08823529412 0.5960477941: 0.08823529412 1.596047794
+DEAL::0.1176470588 0.6254595588: 0.1176470588 1.625459559
+DEAL::0.1470588235 0.6548713235: 0.1470588235 1.654871324
+DEAL::0.1764705882 0.6842830882: 0.1764705882 1.684283088
+DEAL::0.2058823529 0.7136948529: 0.2058823529 1.713694853
+DEAL::0.2352941176 0.7431066176: 0.2352941176 1.743106618
+DEAL::0.2647058824 0.7725183824: 0.2647058824 1.772518382
+DEAL::0.2941176471 0.8019301471: 0.2941176471 1.801930147
+DEAL::0.3235294118 0.8313419118: 0.3235294118 1.831341912
+DEAL::0.3529411765 0.8607536765: 0.3529411765 1.860753676
+DEAL::
+DEAL::Cell with center 0.2500000000 0.7500000000
+DEAL::0.000000000 0.5078125000: 0.000000000 1.507812500
+DEAL::0.02631578947 0.5341282895: 0.02631578947 1.534128289
+DEAL::0.05263157895 0.5604440789: 0.05263157895 1.560444079
+DEAL::0.07894736842 0.5867598684: 0.07894736842 1.586759868
+DEAL::0.1052631579 0.6130756579: 0.1052631579 1.613075658
+DEAL::0.1315789474 0.6393914474: 0.1315789474 1.639391447
+DEAL::0.1578947368 0.6657072368: 0.1578947368 1.665707237
+DEAL::0.1842105263 0.6920230263: 0.1842105263 1.692023026
+DEAL::0.2105263158 0.7183388158: 0.2105263158 1.718338816
+DEAL::0.2368421053 0.7446546053: 0.2368421053 1.744654605
+DEAL::0.2631578947 0.7709703947: 0.2631578947 1.770970395
+DEAL::0.2894736842 0.7972861842: 0.2894736842 1.797286184
+DEAL::0.3157894737 0.8236019737: 0.3157894737 1.823601974
+DEAL::
+DEAL::Cell with center 0.7500000000 0.7500000000
+DEAL::0.5000000000 0.5078125000: 0.5000000000 1.507812500
+DEAL::0.5294117647 0.5372242647: 0.5294117647 1.537224265
+DEAL::0.5588235294 0.5666360294: 0.5588235294 1.566636029
+DEAL::0.5882352941 0.5960477941: 0.5882352941 1.596047794
+DEAL::0.6176470588 0.6254595588: 0.6176470588 1.625459559
+DEAL::0.6470588235 0.6548713235: 0.6470588235 1.654871324
+DEAL::0.6764705882 0.6842830882: 0.6764705882 1.684283088
+DEAL::0.7058823529 0.7136948529: 0.7058823529 1.713694853
+DEAL::0.7352941176 0.7431066176: 0.7352941176 1.743106618
+DEAL::0.7647058824 0.7725183824: 0.7647058824 1.772518382
+DEAL::0.7941176471 0.8019301471: 0.7941176471 1.801930147
+DEAL::0.8235294118 0.8313419118: 0.8235294118 1.831341912
+DEAL::0.8529411765 0.8607536765: 0.8529411765 1.860753676
+DEAL::
+DEAL::Cell with center 0.7500000000 0.7500000000
+DEAL::0.5000000000 0.5078125000: 0.5000000000 1.507812500
+DEAL::0.5263157895 0.5341282895: 0.5263157895 1.534128289
+DEAL::0.5526315789 0.5604440789: 0.5526315789 1.560444079
+DEAL::0.5789473684 0.5867598684: 0.5789473684 1.586759868
+DEAL::0.6052631579 0.6130756579: 0.6052631579 1.613075658
+DEAL::0.6315789474 0.6393914474: 0.6315789474 1.639391447
+DEAL::0.6578947368 0.6657072368: 0.6578947368 1.665707237
+DEAL::0.6842105263 0.6920230263: 0.6842105263 1.692023026
+DEAL::0.7105263158 0.7183388158: 0.7105263158 1.718338816
+DEAL::0.7368421053 0.7446546053: 0.7368421053 1.744654605
+DEAL::0.7631578947 0.7709703947: 0.7631578947 1.770970395
+DEAL::0.7894736842 0.7972861842: 0.7894736842 1.797286184
+DEAL::0.8157894737 0.8236019737: 0.8157894737 1.823601974
+DEAL::
+DEAL::Cartesian linear mapping
+DEAL::Cell with center 0.2500000000 0.2500000000 0.2500000000
+DEAL::0.000000000 0.007812500000 0.01562500000: 0.000000000 1.007812500 2.015625000
+DEAL::0.02941176471 0.03722426471 0.04503676471: 0.02941176471 1.037224265 2.045036765
+DEAL::0.05882352941 0.06663602941 0.07444852941: 0.05882352941 1.066636029 2.074448529
+DEAL::0.08823529412 0.09604779412 0.1038602941: 0.08823529412 1.096047794 2.103860294
+DEAL::0.1176470588 0.1254595588 0.1332720588: 0.1176470588 1.125459559 2.133272059
+DEAL::0.1470588235 0.1548713235 0.1626838235: 0.1470588235 1.154871324 2.162683824
+DEAL::0.1764705882 0.1842830882 0.1920955882: 0.1764705882 1.184283088 2.192095588
+DEAL::0.2058823529 0.2136948529 0.2215073529: 0.2058823529 1.213694853 2.221507353
+DEAL::0.2352941176 0.2431066176 0.2509191176: 0.2352941176 1.243106618 2.250919118
+DEAL::0.2647058824 0.2725183824 0.2803308824: 0.2647058824 1.272518382 2.280330882
+DEAL::0.2941176471 0.3019301471 0.3097426471: 0.2941176471 1.301930147 2.309742647
+DEAL::0.3235294118 0.3313419118 0.3391544118: 0.3235294118 1.331341912 2.339154412
+DEAL::0.3529411765 0.3607536765 0.3685661765: 0.3529411765 1.360753676 2.368566176
+DEAL::
+DEAL::Cell with center 0.2500000000 0.2500000000 0.2500000000
+DEAL::0.000000000 0.007812500000 0.01562500000: 0.000000000 1.007812500 2.015625000
+DEAL::0.02631578947 0.03412828947 0.04194078947: 0.02631578947 1.034128289 2.041940789
+DEAL::0.05263157895 0.06044407895 0.06825657895: 0.05263157895 1.060444079 2.068256579
+DEAL::0.07894736842 0.08675986842 0.09457236842: 0.07894736842 1.086759868 2.094572368
+DEAL::0.1052631579 0.1130756579 0.1208881579: 0.1052631579 1.113075658 2.120888158
+DEAL::0.1315789474 0.1393914474 0.1472039474: 0.1315789474 1.139391447 2.147203947
+DEAL::0.1578947368 0.1657072368 0.1735197368: 0.1578947368 1.165707237 2.173519737
+DEAL::0.1842105263 0.1920230263 0.1998355263: 0.1842105263 1.192023026 2.199835526
+DEAL::0.2105263158 0.2183388158 0.2261513158: 0.2105263158 1.218338816 2.226151316
+DEAL::0.2368421053 0.2446546053 0.2524671053: 0.2368421053 1.244654605 2.252467105
+DEAL::0.2631578947 0.2709703947 0.2787828947: 0.2631578947 1.270970395 2.278782895
+DEAL::0.2894736842 0.2972861842 0.3050986842: 0.2894736842 1.297286184 2.305098684
+DEAL::0.3157894737 0.3236019737 0.3314144737: 0.3157894737 1.323601974 2.331414474
+DEAL::
+DEAL::Cell with center 0.7500000000 0.2500000000 0.2500000000
+DEAL::0.5000000000 0.007812500000 0.01562500000: 0.5000000000 1.007812500 2.015625000
+DEAL::0.5294117647 0.03722426471 0.04503676471: 0.5294117647 1.037224265 2.045036765
+DEAL::0.5588235294 0.06663602941 0.07444852941: 0.5588235294 1.066636029 2.074448529
+DEAL::0.5882352941 0.09604779412 0.1038602941: 0.5882352941 1.096047794 2.103860294
+DEAL::0.6176470588 0.1254595588 0.1332720588: 0.6176470588 1.125459559 2.133272059
+DEAL::0.6470588235 0.1548713235 0.1626838235: 0.6470588235 1.154871324 2.162683824
+DEAL::0.6764705882 0.1842830882 0.1920955882: 0.6764705882 1.184283088 2.192095588
+DEAL::0.7058823529 0.2136948529 0.2215073529: 0.7058823529 1.213694853 2.221507353
+DEAL::0.7352941176 0.2431066176 0.2509191176: 0.7352941176 1.243106618 2.250919118
+DEAL::0.7647058824 0.2725183824 0.2803308824: 0.7647058824 1.272518382 2.280330882
+DEAL::0.7941176471 0.3019301471 0.3097426471: 0.7941176471 1.301930147 2.309742647
+DEAL::0.8235294118 0.3313419118 0.3391544118: 0.8235294118 1.331341912 2.339154412
+DEAL::0.8529411765 0.3607536765 0.3685661765: 0.8529411765 1.360753676 2.368566176
+DEAL::
+DEAL::Cell with center 0.7500000000 0.2500000000 0.2500000000
+DEAL::0.5000000000 0.007812500000 0.01562500000: 0.5000000000 1.007812500 2.015625000
+DEAL::0.5263157895 0.03412828947 0.04194078947: 0.5263157895 1.034128289 2.041940789
+DEAL::0.5526315789 0.06044407895 0.06825657895: 0.5526315789 1.060444079 2.068256579
+DEAL::0.5789473684 0.08675986842 0.09457236842: 0.5789473684 1.086759868 2.094572368
+DEAL::0.6052631579 0.1130756579 0.1208881579: 0.6052631579 1.113075658 2.120888158
+DEAL::0.6315789474 0.1393914474 0.1472039474: 0.6315789474 1.139391447 2.147203947
+DEAL::0.6578947368 0.1657072368 0.1735197368: 0.6578947368 1.165707237 2.173519737
+DEAL::0.6842105263 0.1920230263 0.1998355263: 0.6842105263 1.192023026 2.199835526
+DEAL::0.7105263158 0.2183388158 0.2261513158: 0.7105263158 1.218338816 2.226151316
+DEAL::0.7368421053 0.2446546053 0.2524671053: 0.7368421053 1.244654605 2.252467105
+DEAL::0.7631578947 0.2709703947 0.2787828947: 0.7631578947 1.270970395 2.278782895
+DEAL::0.7894736842 0.2972861842 0.3050986842: 0.7894736842 1.297286184 2.305098684
+DEAL::0.8157894737 0.3236019737 0.3314144737: 0.8157894737 1.323601974 2.331414474
+DEAL::
+DEAL::Cell with center 0.2500000000 0.7500000000 0.2500000000
+DEAL::0.000000000 0.5078125000 0.01562500000: 0.000000000 1.507812500 2.015625000
+DEAL::0.02941176471 0.5372242647 0.04503676471: 0.02941176471 1.537224265 2.045036765
+DEAL::0.05882352941 0.5666360294 0.07444852941: 0.05882352941 1.566636029 2.074448529
+DEAL::0.08823529412 0.5960477941 0.1038602941: 0.08823529412 1.596047794 2.103860294
+DEAL::0.1176470588 0.6254595588 0.1332720588: 0.1176470588 1.625459559 2.133272059
+DEAL::0.1470588235 0.6548713235 0.1626838235: 0.1470588235 1.654871324 2.162683824
+DEAL::0.1764705882 0.6842830882 0.1920955882: 0.1764705882 1.684283088 2.192095588
+DEAL::0.2058823529 0.7136948529 0.2215073529: 0.2058823529 1.713694853 2.221507353
+DEAL::0.2352941176 0.7431066176 0.2509191176: 0.2352941176 1.743106618 2.250919118
+DEAL::0.2647058824 0.7725183824 0.2803308824: 0.2647058824 1.772518382 2.280330882
+DEAL::0.2941176471 0.8019301471 0.3097426471: 0.2941176471 1.801930147 2.309742647
+DEAL::0.3235294118 0.8313419118 0.3391544118: 0.3235294118 1.831341912 2.339154412
+DEAL::0.3529411765 0.8607536765 0.3685661765: 0.3529411765 1.860753676 2.368566176
+DEAL::
+DEAL::Cell with center 0.2500000000 0.7500000000 0.2500000000
+DEAL::0.000000000 0.5078125000 0.01562500000: 0.000000000 1.507812500 2.015625000
+DEAL::0.02631578947 0.5341282895 0.04194078947: 0.02631578947 1.534128289 2.041940789
+DEAL::0.05263157895 0.5604440789 0.06825657895: 0.05263157895 1.560444079 2.068256579
+DEAL::0.07894736842 0.5867598684 0.09457236842: 0.07894736842 1.586759868 2.094572368
+DEAL::0.1052631579 0.6130756579 0.1208881579: 0.1052631579 1.613075658 2.120888158
+DEAL::0.1315789474 0.6393914474 0.1472039474: 0.1315789474 1.639391447 2.147203947
+DEAL::0.1578947368 0.6657072368 0.1735197368: 0.1578947368 1.665707237 2.173519737
+DEAL::0.1842105263 0.6920230263 0.1998355263: 0.1842105263 1.692023026 2.199835526
+DEAL::0.2105263158 0.7183388158 0.2261513158: 0.2105263158 1.718338816 2.226151316
+DEAL::0.2368421053 0.7446546053 0.2524671053: 0.2368421053 1.744654605 2.252467105
+DEAL::0.2631578947 0.7709703947 0.2787828947: 0.2631578947 1.770970395 2.278782895
+DEAL::0.2894736842 0.7972861842 0.3050986842: 0.2894736842 1.797286184 2.305098684
+DEAL::0.3157894737 0.8236019737 0.3314144737: 0.3157894737 1.823601974 2.331414474
+DEAL::
+DEAL::Cell with center 0.7500000000 0.7500000000 0.2500000000
+DEAL::0.5000000000 0.5078125000 0.01562500000: 0.5000000000 1.507812500 2.015625000
+DEAL::0.5294117647 0.5372242647 0.04503676471: 0.5294117647 1.537224265 2.045036765
+DEAL::0.5588235294 0.5666360294 0.07444852941: 0.5588235294 1.566636029 2.074448529
+DEAL::0.5882352941 0.5960477941 0.1038602941: 0.5882352941 1.596047794 2.103860294
+DEAL::0.6176470588 0.6254595588 0.1332720588: 0.6176470588 1.625459559 2.133272059
+DEAL::0.6470588235 0.6548713235 0.1626838235: 0.6470588235 1.654871324 2.162683824
+DEAL::0.6764705882 0.6842830882 0.1920955882: 0.6764705882 1.684283088 2.192095588
+DEAL::0.7058823529 0.7136948529 0.2215073529: 0.7058823529 1.713694853 2.221507353
+DEAL::0.7352941176 0.7431066176 0.2509191176: 0.7352941176 1.743106618 2.250919118
+DEAL::0.7647058824 0.7725183824 0.2803308824: 0.7647058824 1.772518382 2.280330882
+DEAL::0.7941176471 0.8019301471 0.3097426471: 0.7941176471 1.801930147 2.309742647
+DEAL::0.8235294118 0.8313419118 0.3391544118: 0.8235294118 1.831341912 2.339154412
+DEAL::0.8529411765 0.8607536765 0.3685661765: 0.8529411765 1.860753676 2.368566176
+DEAL::
+DEAL::Cell with center 0.7500000000 0.7500000000 0.2500000000
+DEAL::0.5000000000 0.5078125000 0.01562500000: 0.5000000000 1.507812500 2.015625000
+DEAL::0.5263157895 0.5341282895 0.04194078947: 0.5263157895 1.534128289 2.041940789
+DEAL::0.5526315789 0.5604440789 0.06825657895: 0.5526315789 1.560444079 2.068256579
+DEAL::0.5789473684 0.5867598684 0.09457236842: 0.5789473684 1.586759868 2.094572368
+DEAL::0.6052631579 0.6130756579 0.1208881579: 0.6052631579 1.613075658 2.120888158
+DEAL::0.6315789474 0.6393914474 0.1472039474: 0.6315789474 1.639391447 2.147203947
+DEAL::0.6578947368 0.6657072368 0.1735197368: 0.6578947368 1.665707237 2.173519737
+DEAL::0.6842105263 0.6920230263 0.1998355263: 0.6842105263 1.692023026 2.199835526
+DEAL::0.7105263158 0.7183388158 0.2261513158: 0.7105263158 1.718338816 2.226151316
+DEAL::0.7368421053 0.7446546053 0.2524671053: 0.7368421053 1.744654605 2.252467105
+DEAL::0.7631578947 0.7709703947 0.2787828947: 0.7631578947 1.770970395 2.278782895
+DEAL::0.7894736842 0.7972861842 0.3050986842: 0.7894736842 1.797286184 2.305098684
+DEAL::0.8157894737 0.8236019737 0.3314144737: 0.8157894737 1.823601974 2.331414474
+DEAL::
+DEAL::Cell with center 0.2500000000 0.2500000000 0.7500000000
+DEAL::0.000000000 0.007812500000 0.5156250000: 0.000000000 1.007812500 2.515625000
+DEAL::0.02941176471 0.03722426471 0.5450367647: 0.02941176471 1.037224265 2.545036765
+DEAL::0.05882352941 0.06663602941 0.5744485294: 0.05882352941 1.066636029 2.574448529
+DEAL::0.08823529412 0.09604779412 0.6038602941: 0.08823529412 1.096047794 2.603860294
+DEAL::0.1176470588 0.1254595588 0.6332720588: 0.1176470588 1.125459559 2.633272059
+DEAL::0.1470588235 0.1548713235 0.6626838235: 0.1470588235 1.154871324 2.662683824
+DEAL::0.1764705882 0.1842830882 0.6920955882: 0.1764705882 1.184283088 2.692095588
+DEAL::0.2058823529 0.2136948529 0.7215073529: 0.2058823529 1.213694853 2.721507353
+DEAL::0.2352941176 0.2431066176 0.7509191176: 0.2352941176 1.243106618 2.750919118
+DEAL::0.2647058824 0.2725183824 0.7803308824: 0.2647058824 1.272518382 2.780330882
+DEAL::0.2941176471 0.3019301471 0.8097426471: 0.2941176471 1.301930147 2.809742647
+DEAL::0.3235294118 0.3313419118 0.8391544118: 0.3235294118 1.331341912 2.839154412
+DEAL::0.3529411765 0.3607536765 0.8685661765: 0.3529411765 1.360753676 2.868566176
+DEAL::
+DEAL::Cell with center 0.2500000000 0.2500000000 0.7500000000
+DEAL::0.000000000 0.007812500000 0.5156250000: 0.000000000 1.007812500 2.515625000
+DEAL::0.02631578947 0.03412828947 0.5419407895: 0.02631578947 1.034128289 2.541940789
+DEAL::0.05263157895 0.06044407895 0.5682565789: 0.05263157895 1.060444079 2.568256579
+DEAL::0.07894736842 0.08675986842 0.5945723684: 0.07894736842 1.086759868 2.594572368
+DEAL::0.1052631579 0.1130756579 0.6208881579: 0.1052631579 1.113075658 2.620888158
+DEAL::0.1315789474 0.1393914474 0.6472039474: 0.1315789474 1.139391447 2.647203947
+DEAL::0.1578947368 0.1657072368 0.6735197368: 0.1578947368 1.165707237 2.673519737
+DEAL::0.1842105263 0.1920230263 0.6998355263: 0.1842105263 1.192023026 2.699835526
+DEAL::0.2105263158 0.2183388158 0.7261513158: 0.2105263158 1.218338816 2.726151316
+DEAL::0.2368421053 0.2446546053 0.7524671053: 0.2368421053 1.244654605 2.752467105
+DEAL::0.2631578947 0.2709703947 0.7787828947: 0.2631578947 1.270970395 2.778782895
+DEAL::0.2894736842 0.2972861842 0.8050986842: 0.2894736842 1.297286184 2.805098684
+DEAL::0.3157894737 0.3236019737 0.8314144737: 0.3157894737 1.323601974 2.831414474
+DEAL::
+DEAL::Cell with center 0.7500000000 0.2500000000 0.7500000000
+DEAL::0.5000000000 0.007812500000 0.5156250000: 0.5000000000 1.007812500 2.515625000
+DEAL::0.5294117647 0.03722426471 0.5450367647: 0.5294117647 1.037224265 2.545036765
+DEAL::0.5588235294 0.06663602941 0.5744485294: 0.5588235294 1.066636029 2.574448529
+DEAL::0.5882352941 0.09604779412 0.6038602941: 0.5882352941 1.096047794 2.603860294
+DEAL::0.6176470588 0.1254595588 0.6332720588: 0.6176470588 1.125459559 2.633272059
+DEAL::0.6470588235 0.1548713235 0.6626838235: 0.6470588235 1.154871324 2.662683824
+DEAL::0.6764705882 0.1842830882 0.6920955882: 0.6764705882 1.184283088 2.692095588
+DEAL::0.7058823529 0.2136948529 0.7215073529: 0.7058823529 1.213694853 2.721507353
+DEAL::0.7352941176 0.2431066176 0.7509191176: 0.7352941176 1.243106618 2.750919118
+DEAL::0.7647058824 0.2725183824 0.7803308824: 0.7647058824 1.272518382 2.780330882
+DEAL::0.7941176471 0.3019301471 0.8097426471: 0.7941176471 1.301930147 2.809742647
+DEAL::0.8235294118 0.3313419118 0.8391544118: 0.8235294118 1.331341912 2.839154412
+DEAL::0.8529411765 0.3607536765 0.8685661765: 0.8529411765 1.360753676 2.868566176
+DEAL::
+DEAL::Cell with center 0.7500000000 0.2500000000 0.7500000000
+DEAL::0.5000000000 0.007812500000 0.5156250000: 0.5000000000 1.007812500 2.515625000
+DEAL::0.5263157895 0.03412828947 0.5419407895: 0.5263157895 1.034128289 2.541940789
+DEAL::0.5526315789 0.06044407895 0.5682565789: 0.5526315789 1.060444079 2.568256579
+DEAL::0.5789473684 0.08675986842 0.5945723684: 0.5789473684 1.086759868 2.594572368
+DEAL::0.6052631579 0.1130756579 0.6208881579: 0.6052631579 1.113075658 2.620888158
+DEAL::0.6315789474 0.1393914474 0.6472039474: 0.6315789474 1.139391447 2.647203947
+DEAL::0.6578947368 0.1657072368 0.6735197368: 0.6578947368 1.165707237 2.673519737
+DEAL::0.6842105263 0.1920230263 0.6998355263: 0.6842105263 1.192023026 2.699835526
+DEAL::0.7105263158 0.2183388158 0.7261513158: 0.7105263158 1.218338816 2.726151316
+DEAL::0.7368421053 0.2446546053 0.7524671053: 0.7368421053 1.244654605 2.752467105
+DEAL::0.7631578947 0.2709703947 0.7787828947: 0.7631578947 1.270970395 2.778782895
+DEAL::0.7894736842 0.2972861842 0.8050986842: 0.7894736842 1.297286184 2.805098684
+DEAL::0.8157894737 0.3236019737 0.8314144737: 0.8157894737 1.323601974 2.831414474
+DEAL::
+DEAL::Cell with center 0.2500000000 0.7500000000 0.7500000000
+DEAL::0.000000000 0.5078125000 0.5156250000: 0.000000000 1.507812500 2.515625000
+DEAL::0.02941176471 0.5372242647 0.5450367647: 0.02941176471 1.537224265 2.545036765
+DEAL::0.05882352941 0.5666360294 0.5744485294: 0.05882352941 1.566636029 2.574448529
+DEAL::0.08823529412 0.5960477941 0.6038602941: 0.08823529412 1.596047794 2.603860294
+DEAL::0.1176470588 0.6254595588 0.6332720588: 0.1176470588 1.625459559 2.633272059
+DEAL::0.1470588235 0.6548713235 0.6626838235: 0.1470588235 1.654871324 2.662683824
+DEAL::0.1764705882 0.6842830882 0.6920955882: 0.1764705882 1.684283088 2.692095588
+DEAL::0.2058823529 0.7136948529 0.7215073529: 0.2058823529 1.713694853 2.721507353
+DEAL::0.2352941176 0.7431066176 0.7509191176: 0.2352941176 1.743106618 2.750919118
+DEAL::0.2647058824 0.7725183824 0.7803308824: 0.2647058824 1.772518382 2.780330882
+DEAL::0.2941176471 0.8019301471 0.8097426471: 0.2941176471 1.801930147 2.809742647
+DEAL::0.3235294118 0.8313419118 0.8391544118: 0.3235294118 1.831341912 2.839154412
+DEAL::0.3529411765 0.8607536765 0.8685661765: 0.3529411765 1.860753676 2.868566176
+DEAL::
+DEAL::Cell with center 0.2500000000 0.7500000000 0.7500000000
+DEAL::0.000000000 0.5078125000 0.5156250000: 0.000000000 1.507812500 2.515625000
+DEAL::0.02631578947 0.5341282895 0.5419407895: 0.02631578947 1.534128289 2.541940789
+DEAL::0.05263157895 0.5604440789 0.5682565789: 0.05263157895 1.560444079 2.568256579
+DEAL::0.07894736842 0.5867598684 0.5945723684: 0.07894736842 1.586759868 2.594572368
+DEAL::0.1052631579 0.6130756579 0.6208881579: 0.1052631579 1.613075658 2.620888158
+DEAL::0.1315789474 0.6393914474 0.6472039474: 0.1315789474 1.639391447 2.647203947
+DEAL::0.1578947368 0.6657072368 0.6735197368: 0.1578947368 1.665707237 2.673519737
+DEAL::0.1842105263 0.6920230263 0.6998355263: 0.1842105263 1.692023026 2.699835526
+DEAL::0.2105263158 0.7183388158 0.7261513158: 0.2105263158 1.718338816 2.726151316
+DEAL::0.2368421053 0.7446546053 0.7524671053: 0.2368421053 1.744654605 2.752467105
+DEAL::0.2631578947 0.7709703947 0.7787828947: 0.2631578947 1.770970395 2.778782895
+DEAL::0.2894736842 0.7972861842 0.8050986842: 0.2894736842 1.797286184 2.805098684
+DEAL::0.3157894737 0.8236019737 0.8314144737: 0.3157894737 1.823601974 2.831414474
+DEAL::
+DEAL::Cell with center 0.7500000000 0.7500000000 0.7500000000
+DEAL::0.5000000000 0.5078125000 0.5156250000: 0.5000000000 1.507812500 2.515625000
+DEAL::0.5294117647 0.5372242647 0.5450367647: 0.5294117647 1.537224265 2.545036765
+DEAL::0.5588235294 0.5666360294 0.5744485294: 0.5588235294 1.566636029 2.574448529
+DEAL::0.5882352941 0.5960477941 0.6038602941: 0.5882352941 1.596047794 2.603860294
+DEAL::0.6176470588 0.6254595588 0.6332720588: 0.6176470588 1.625459559 2.633272059
+DEAL::0.6470588235 0.6548713235 0.6626838235: 0.6470588235 1.654871324 2.662683824
+DEAL::0.6764705882 0.6842830882 0.6920955882: 0.6764705882 1.684283088 2.692095588
+DEAL::0.7058823529 0.7136948529 0.7215073529: 0.7058823529 1.713694853 2.721507353
+DEAL::0.7352941176 0.7431066176 0.7509191176: 0.7352941176 1.743106618 2.750919118
+DEAL::0.7647058824 0.7725183824 0.7803308824: 0.7647058824 1.772518382 2.780330882
+DEAL::0.7941176471 0.8019301471 0.8097426471: 0.7941176471 1.801930147 2.809742647
+DEAL::0.8235294118 0.8313419118 0.8391544118: 0.8235294118 1.831341912 2.839154412
+DEAL::0.8529411765 0.8607536765 0.8685661765: 0.8529411765 1.860753676 2.868566176
+DEAL::
+DEAL::Cell with center 0.7500000000 0.7500000000 0.7500000000
+DEAL::0.5000000000 0.5078125000 0.5156250000: 0.5000000000 1.507812500 2.515625000
+DEAL::0.5263157895 0.5341282895 0.5419407895: 0.5263157895 1.534128289 2.541940789
+DEAL::0.5526315789 0.5604440789 0.5682565789: 0.5526315789 1.560444079 2.568256579
+DEAL::0.5789473684 0.5867598684 0.5945723684: 0.5789473684 1.586759868 2.594572368
+DEAL::0.6052631579 0.6130756579 0.6208881579: 0.6052631579 1.613075658 2.620888158
+DEAL::0.6315789474 0.6393914474 0.6472039474: 0.6315789474 1.639391447 2.647203947
+DEAL::0.6578947368 0.6657072368 0.6735197368: 0.6578947368 1.665707237 2.673519737
+DEAL::0.6842105263 0.6920230263 0.6998355263: 0.6842105263 1.692023026 2.699835526
+DEAL::0.7105263158 0.7183388158 0.7261513158: 0.7105263158 1.718338816 2.726151316
+DEAL::0.7368421053 0.7446546053 0.7524671053: 0.7368421053 1.744654605 2.752467105
+DEAL::0.7631578947 0.7709703947 0.7787828947: 0.7631578947 1.770970395 2.778782895
+DEAL::0.7894736842 0.7972861842 0.8050986842: 0.7894736842 1.797286184 2.805098684
+DEAL::0.8157894737 0.8236019737 0.8314144737: 0.8157894737 1.823601974 2.831414474
+DEAL::