From: Rene Gassmoeller Date: Thu, 30 Mar 2023 12:58:32 +0000 (-0400) Subject: Add test X-Git-Tag: v9.5.0-rc1~385^2~2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=f121f893f95c8171ba463284da7f5a02faace640;p=dealii.git Add test --- diff --git a/tests/matrix_free/point_evaluation_20.cc b/tests/matrix_free/point_evaluation_20.cc new file mode 100644 index 0000000000..f7d65e7e2e --- /dev/null +++ b/tests/matrix_free/point_evaluation_20.cc @@ -0,0 +1,146 @@ +// --------------------------------------------------------------------- +// +// 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 + +#include + +#include +#include +#include + +#include +#include + +#include + +#include + +#include + +#include + +#include "../tests.h" + + + +template +class MyFunction : public Function +{ +public: + MyFunction() + : Function(dim) + {} + + double + value(const Point &p, const unsigned int component) const override + { + return component + p[component]; + } +}; + + + +template +void +test() +{ + using namespace dealii; + Triangulation tria; + + GridGenerator::subdivided_hyper_cube(tria, 2, 0, 1); + + MappingCartesian mapping; + deallog << "Cartesian linear mapping" << std::endl; + + std::vector> unit_points; + for (unsigned int i = 0; i < 13; ++i) + { + Point p; + for (unsigned int d = 0; d < dim; ++d) + p[d] = static_cast(i) / 17. + 0.015625 * d; + unit_points.push_back(p); + } + + std::vector> unit_points2; + for (unsigned int i = 0; i < 13; ++i) + { + Point p; + for (unsigned int d = 0; d < dim; ++d) + p[d] = static_cast(i) / 19. + 0.015625 * d; + unit_points2.push_back(p); + } + + FESystem fe(FE_Q(2), dim); + DoFHandler dof_handler(tria); + dof_handler.distribute_dofs(fe); + Vector vector(dof_handler.n_dofs()); + + // mapping info object for precomputed mapping + NonMatching::MappingInfo mapping_info(mapping, + update_values | + update_gradients); + // First evaluator to use precomputed mapping + FEPointEvaluation evaluator1(mapping_info, fe); + + VectorTools::interpolate(mapping, dof_handler, MyFunction(), vector); + + std::vector 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>(); +} diff --git a/tests/matrix_free/point_evaluation_20.output b/tests/matrix_free/point_evaluation_20.output new file mode 100644 index 0000000000..f740b951ea --- /dev/null +++ b/tests/matrix_free/point_evaluation_20.output @@ -0,0 +1,363 @@ + +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::