From 1e97766e5977ef5610f03293992dd41e19604a87 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Wed, 20 Mar 2024 16:23:56 +0100 Subject: [PATCH] Add new test case --- tests/matrix_free/point_evaluation_30.cc | 152 +++++++ tests/matrix_free/point_evaluation_30.output | 277 ++++++++++++ tests/matrix_free/point_evaluation_31.cc | 160 +++++++ tests/matrix_free/point_evaluation_31.output | 180 ++++++++ tests/matrix_free/point_evaluation_32.cc | 167 +++++++ tests/matrix_free/point_evaluation_32.output | 439 +++++++++++++++++++ 6 files changed, 1375 insertions(+) create mode 100644 tests/matrix_free/point_evaluation_30.cc create mode 100644 tests/matrix_free/point_evaluation_30.output create mode 100644 tests/matrix_free/point_evaluation_31.cc create mode 100644 tests/matrix_free/point_evaluation_31.output create mode 100644 tests/matrix_free/point_evaluation_32.cc create mode 100644 tests/matrix_free/point_evaluation_32.output diff --git a/tests/matrix_free/point_evaluation_30.cc b/tests/matrix_free/point_evaluation_30.cc new file mode 100644 index 0000000000..36420fcd8c --- /dev/null +++ b/tests/matrix_free/point_evaluation_30.cc @@ -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 + +#include + +#include +#include +#include + +#include +#include + +#include + +#include + +#include + +#include + +#include "../tests.h" + + + +template +void +test(const unsigned int degree) +{ + Triangulation tria; + + if constexpr (dim + 1 == spacedim) + GridGenerator::hyper_sphere(tria, Point(), 1.); + else if constexpr (dim == 1) + GridGenerator::hyper_rectangle(tria, Point(), Point(1.2)); + else + static_assert("Case not implemented"); + + MappingQ mapping(degree); + deallog << "Mapping of degree " << degree << 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); + } + + FE_Q fe(degree); + FEValues fe_values(mapping, + fe, + Quadrature(unit_points), + update_values | update_gradients); + + DoFHandler dof_handler(tria); + dof_handler.distribute_dofs(fe); + Vector 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(exponents), + vector); + + std::vector solution_values(fe.dofs_per_cell); + std::vector function_values(unit_points.size()); + std::vector> 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 ? 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 index 0000000000..4246ff2f2f --- /dev/null +++ b/tests/matrix_free/point_evaluation_30.output @@ -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 index 0000000000..0a97d929b3 --- /dev/null +++ b/tests/matrix_free/point_evaluation_31.cc @@ -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 + +#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(const unsigned int degree) +{ + Triangulation tria; + + if constexpr (dim + 1 == spacedim) + GridGenerator::hyper_sphere(tria, Point(), 1.); + else if constexpr (dim == 1) + GridGenerator::hyper_rectangle(tria, Point(), Point(1.2)); + else + static_assert("Case not implemented"); + + MappingQ mapping(degree); + deallog << "Mapping of degree " << degree << 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); + } + + FESystem fe(FE_Q(degree), spacedim); + FEValues fe_values(mapping, + fe, + Quadrature(unit_points), + update_values | update_gradients); + + DoFHandler dof_handler(tria); + dof_handler.distribute_dofs(fe); + Vector vector(dof_handler.n_dofs()); + + FEPointEvaluation evaluator( + mapping, fe, update_values | update_gradients); + + VectorTools::interpolate(mapping, + dof_handler, + MyFunction(), + vector); + + std::vector solution_values(fe.dofs_per_cell); + std::vector> function_values(unit_points.size()); + std::vector> 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 index 0000000000..28850f770e --- /dev/null +++ b/tests/matrix_free/point_evaluation_31.output @@ -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 index 0000000000..b69106e0df --- /dev/null +++ b/tests/matrix_free/point_evaluation_32.cc @@ -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 + +#include + +#include +#include +#include + +#include +#include + +#include + +#include + +#include + +#include + +#include "../tests.h" + + + +template +class MyFunction : public Function +{ +public: + MyFunction() + : Function() + {} + + double + value(const Point &p, const unsigned int component) const override + { + return component + p[component]; + } +}; + + + +template +void +test(const unsigned int degree) +{ + Triangulation tria; + + + if constexpr (dim + 1 == spacedim) + GridGenerator::hyper_sphere(tria, Point(), 1.); + else if constexpr (dim == 1) + GridGenerator::hyper_rectangle(tria, Point(), Point(1.2)); + else + static_assert("Case not implemented"); + + MappingQ mapping(degree); + deallog << "Mapping of degree " << degree << std::endl; + + std::vector> unit_points; + for (unsigned int i = 0; i < (dim > 1 ? 13 : 1); ++i) + { + Point p; + for (unsigned int d = 0; d < dim - 1; ++d) + p[d] = static_cast(i) / 17. + 0.015625 * d; + unit_points.push_back(p); + } + + FE_Q fe(degree); + DoFHandler dof_handler(tria); + dof_handler.distribute_dofs(fe); + Vector vector(dof_handler.n_dofs()); + + NonMatching::MappingInfo mapping_info(mapping, + update_values | + update_gradients); + + std::vector>> 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(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(), + vector); + + std::vector solution_values_in(fe.dofs_per_cell); + std::vector 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 index 0000000000..ce5302dcb2 --- /dev/null +++ b/tests/matrix_free/point_evaluation_32.output @@ -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 -- 2.39.5