From 8edb4a9e9ca0e8708cf3c03eebd4ac5ad5606dd3 Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Wed, 21 Oct 2020 11:30:49 +0200 Subject: [PATCH] Test case for transformation --- .../mapping_q_real_to_unit_internal.cc | 138 +++ .../mapping_q_real_to_unit_internal.output | 1002 +++++++++++++++++ 2 files changed, 1140 insertions(+) create mode 100644 tests/mappings/mapping_q_real_to_unit_internal.cc create mode 100644 tests/mappings/mapping_q_real_to_unit_internal.output diff --git a/tests/mappings/mapping_q_real_to_unit_internal.cc b/tests/mappings/mapping_q_real_to_unit_internal.cc new file mode 100644 index 0000000000..a02462d8ad --- /dev/null +++ b/tests/mappings/mapping_q_real_to_unit_internal.cc @@ -0,0 +1,138 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2020 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE at +// the top level of the deal.II distribution. +// +// --------------------------------------------------------------------- + + +// Check internal implementation of +// MappingQGeneric::transform_real_to_unit_point by printing Newton iteration +// information. This test is sensitive to roundoff errors by the nature of +// what gets tested, which can cause one more or one less iteration, +// especially due to FMA + +#include +#include +#include +#include + +#include +#include +#include +#include +#include +#include +#include + +#include +#include + +#include + +#include "../tests.h" + + + +template +void +print_result(const unsigned int mapping_degree, + const Triangulation &tria, + const Point p) +{ + deallog << "Testing " << dim << "D with point " << p << std::endl; + + FE_Q dummy(mapping_degree); + MappingQGeneric mapping(mapping_degree); + + FEValues fe_values(mapping, + dummy, + Quadrature(dummy.get_unit_support_points()), + update_quadrature_points); + + std::vector> polynomials = + Polynomials::generate_complete_Lagrange_basis( + QGaussLobatto<1>(mapping_degree + 1).get_points()); + std::vector renumber = + FETools::lexicographic_to_hierarchic_numbering(mapping_degree); + + for (const auto &cell : tria.active_cell_iterators()) + { + fe_values.reinit(cell); + deallog << "Testing on cell " << cell->id() << " with center " + << cell->center(true) << std::endl; + if (GeometryInfo::distance_to_unit_cell( + cell->real_to_unit_cell_affine_approximation(p)) < + (-0.6 + 1.3 * dim)) + internal::MappingQGenericImplementation:: + do_transform_real_to_unit_cell_internal( + p, + cell->real_to_unit_cell_affine_approximation(p), + fe_values.get_quadrature_points(), + polynomials, + renumber, + /* print_iterations = */ true); + deallog << std::endl; + } + deallog << std::endl; +} + + + +template +void +test(const unsigned mapping_degree, const unsigned int n_ref) +{ + Triangulation tria; + GridGenerator::hyper_shell(tria, Point(), 0.8, 1., dim == 2 ? 3 : 6); + tria.refine_global(n_ref); + { + const double phi = 0.56 * numbers::PI; + const double theta = 0.49 * numbers::PI; + for (unsigned int i = 0; i <= 4; ++i) + { + const double r = 0.8 + 0.2 * i / 4; + Point p; + p[0] = std::cos(phi) * std::sin(theta) * r; + p[1] = std::sin(phi) * std::sin(theta) * r; + if (dim == 3) + p[2] = std::cos(theta) * r; + print_result(mapping_degree, tria, p); + } + } + if (dim == 3) + { + const double phi = 0.3 * numbers::PI; + const double theta = 0.1; + for (unsigned int i = 0; i <= 6; ++i) + { + const double r = 0.8 + 0.2 * i / 6; + Point p; + p[0] = std::cos(phi) * std::sin(theta) * r; + p[1] = std::sin(phi) * std::sin(theta) * r; + if (dim == 3) + p[2] = std::cos(theta) * r; + print_result(mapping_degree, tria, p); + } + } +} + + + +int +main() +{ + initlog(); + deallog << std::setprecision(10); + + test<2, 2>(6, 1); + test<3, 3>(4, 0); +} diff --git a/tests/mappings/mapping_q_real_to_unit_internal.output b/tests/mappings/mapping_q_real_to_unit_internal.output new file mode 100644 index 0000000000..67c0a772fc --- /dev/null +++ b/tests/mappings/mapping_q_real_to_unit_internal.output @@ -0,0 +1,1002 @@ + +DEAL::Testing 2D with point -0.1498310826 0.7854420410 +DEAL::Testing on cell 0_1:0 with center 0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 0_1:1 with center 5.817072296e-17 0.9500000000 +DEAL::Newton iteration 0 for unit point guess 0.6577169290 0.9304965238 +DEAL:: delta=-0.01964393136 -1.075627827 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1091688320 +DEAL:: ||f*|| =0.002220292790 +DEAL:: ||f*||_A =0.002879917279 +DEAL::Newton iteration 1 for unit point guess 0.6773608604 2.006124351 +DEAL:: delta=-0.002639886633 0.002146340320 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.002220292790 +DEAL:: ||f*|| =3.111968626e-06 +DEAL:: ||f*||_A =3.056234667e-05 +DEAL::Newton iteration 2 for unit point guess 0.6800007470 2.003978011 +DEAL:: delta=7.153584652e-07 3.053776106e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.111968626e-06 +DEAL:: ||f*|| =2.475371550e-12 +DEAL:: ||f*||_A =5.596813561e-12 +DEAL::Iteration converged for p_unit = [ 0.6800000317 2.003947473 ] and iteration error 5.596813561e-12 +DEAL:: +DEAL::Testing on cell 0_1:2 with center 0.7361215932 0.4250000000 +DEAL:: +DEAL::Testing on cell 0_1:3 with center 5.204748896e-17 0.8500000000 +DEAL::Newton iteration 0 for unit point guess 0.6762718619 -0.06950347624 +DEAL:: delta=-0.003286901688 -1.073511888 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1073965686 +DEAL:: ||f*|| =0.0003695425182 +DEAL:: ||f*||_A =0.0003919325980 +DEAL::Newton iteration 1 for unit point guess 0.6795587635 1.004008412 +DEAL:: delta=-0.0004412714742 6.008501490e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0003695425182 +DEAL:: ||f*|| =8.541647302e-08 +DEAL:: ||f*||_A =8.537238397e-07 +DEAL::Newton iteration 2 for unit point guess 0.6800000350 1.003948327 +DEAL:: delta=3.347866320e-09 8.537045958e-07 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =8.541647302e-08 +DEAL:: ||f*|| =3.342213889e-16 +DEAL:: ||f*||_A =1.746985821e-15 +DEAL::Iteration converged for p_unit = [ 0.6800000317 1.003947473 ] and iteration error 1.746985821e-15 +DEAL:: +DEAL::Testing on cell 1_1:0 with center -0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:1 with center -0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:2 with center -0.7361215932 0.4250000000 +DEAL::Newton iteration 0 for unit point guess -0.2121143758 2.966937436 +DEAL:: delta=0.1426653570 1.911921754 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2113822431 +DEAL:: ||f*|| =0.02949335980 +DEAL:: ||f*||_A =0.09899499681 +DEAL::Newton iteration 1 for unit point guess -0.3547797328 1.055015682 +DEAL:: delta=-0.03501372833 0.04561257975 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.02949335980 +DEAL:: ||f*|| =0.0005592232088 +DEAL:: ||f*||_A =0.005362688549 +DEAL::Newton iteration 2 for unit point guess -0.3197660045 1.009403103 +DEAL:: delta=0.0002139315510 0.005297328002 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0005592232088 +DEAL:: ||f*|| =1.204134814e-07 +DEAL:: ||f*||_A =2.454543669e-07 +DEAL::Newton iteration 3 for unit point guess -0.3199799360 1.004105775 +DEAL:: delta=-1.417819883e-07 2.000356681e-07 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.204134814e-07 +DEAL:: ||f*|| =2.123533725e-14 +DEAL:: ||f*||_A =1.942358461e-13 +DEAL::Iteration converged for p_unit = [ -0.3199797942 1.004105575 ] and iteration error 1.942358461e-13 +DEAL:: +DEAL::Testing on cell 1_1:3 with center -0.7361215932 -0.4250000000 +DEAL:: +DEAL::Testing on cell 2_1:0 with center -1.018281668e-15 -0.9500000000 +DEAL:: +DEAL::Testing on cell 2_1:1 with center 0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 2_1:2 with center -9.110941236e-16 -0.8500000000 +DEAL:: +DEAL::Testing on cell 2_1:3 with center 0.7361215932 -0.4250000000 +DEAL:: +DEAL:: +DEAL::Testing 2D with point -0.1591955252 0.8345321686 +DEAL::Testing on cell 0_1:0 with center 0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 0_1:1 with center 5.817072296e-17 0.9500000000 +DEAL::Newton iteration 0 for unit point guess 0.6675742371 0.3636525565 +DEAL:: delta=-0.01095479005 -1.141260873 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1146602330 +DEAL:: ||f*|| =0.001310641447 +DEAL:: ||f*||_A =0.001412896727 +DEAL::Newton iteration 1 for unit point guess 0.6785290271 1.504913430 +DEAL:: delta=-0.001471128494 0.0007091597170 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.001310641447 +DEAL:: ||f*|| =1.014013368e-06 +DEAL:: ||f*||_A =1.008240256e-05 +DEAL::Newton iteration 2 for unit point guess 0.6800001556 1.504204270 +DEAL:: delta=1.239682307e-07 1.007997444e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.014013368e-06 +DEAL:: ||f*|| =1.543231669e-13 +DEAL:: ||f*||_A =2.394277731e-13 +DEAL::Iteration converged for p_unit = [ 0.6800000317 1.504194190 ] and iteration error 2.394277731e-13 +DEAL:: +DEAL::Testing on cell 0_1:2 with center 0.7361215932 0.4250000000 +DEAL:: +DEAL::Testing on cell 0_1:3 with center 5.204748896e-17 0.8500000000 +DEAL::Newton iteration 0 for unit point guess 0.6872888532 -0.6363474435 +DEAL:: delta=0.006426066544 -1.140789116 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1142630680 +DEAL:: ||f*|| =0.0007679627762 +DEAL:: ||f*||_A =0.0007847271756 +DEAL::Newton iteration 1 for unit point guess 0.6808627867 0.5044416722 +DEAL:: delta=0.0008627800205 0.0002440148031 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0007679627762 +DEAL:: ||f*|| =3.474553334e-07 +DEAL:: ||f*||_A =3.467706919e-06 +DEAL::Newton iteration 2 for unit point guess 0.6800000067 0.5041976574 +DEAL:: delta=-2.501310646e-08 3.467419577e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.474553334e-07 +DEAL:: ||f*|| =9.128842154e-15 +DEAL:: ||f*||_A =1.103190943e-14 +DEAL::Iteration converged for p_unit = [ 0.6800000317 0.5041941900 ] and iteration error 1.103190943e-14 +DEAL:: +DEAL::Testing on cell 1_1:0 with center -0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:1 with center -0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:2 with center -0.7361215932 0.4250000000 +DEAL::Newton iteration 0 for unit point guess -0.2566215243 2.589871026 +DEAL:: delta=0.08391429210 2.066838582 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2142321563 +DEAL:: ||f*|| =0.01836658887 +DEAL:: ||f*||_A =0.04233118966 +DEAL::Newton iteration 1 for unit point guess -0.3405358164 0.5230324440 +DEAL:: delta=-0.02059922788 0.01670588339 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.01836658887 +DEAL:: ||f*|| =0.0002001773451 +DEAL:: ||f*||_A =0.001972573880 +DEAL::Newton iteration 2 for unit point guess -0.3199365885 0.5063265606 +DEAL:: delta=4.321573340e-05 0.001964378821 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0002001773451 +DEAL:: ||f*|| =8.935712282e-09 +DEAL:: ||f*||_A =1.324222136e-08 +DEAL::Newton iteration 3 for unit point guess -0.3199798042 0.5043621818 +DEAL:: delta=-9.994412487e-09 8.680613328e-09 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =8.935712282e-09 +DEAL:: ||f*|| =7.850462293e-15 +DEAL:: ||f*||_A =9.547950136e-15 +DEAL::Iteration converged for p_unit = [ -0.3199797942 0.5043621731 ] and iteration error 9.547950136e-15 +DEAL:: +DEAL::Testing on cell 1_1:3 with center -0.7361215932 -0.4250000000 +DEAL:: +DEAL::Testing on cell 2_1:0 with center -1.018281668e-15 -0.9500000000 +DEAL:: +DEAL::Testing on cell 2_1:1 with center 0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 2_1:2 with center -9.110941236e-16 -0.8500000000 +DEAL:: +DEAL::Testing on cell 2_1:3 with center 0.7361215932 -0.4250000000 +DEAL:: +DEAL:: +DEAL::Testing 2D with point -0.1685599679 0.8836222961 +DEAL::Testing on cell 0_1:0 with center 0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 0_1:1 with center 5.817072296e-17 0.9500000000 +DEAL::Newton iteration 0 for unit point guess 0.6774315452 -0.2031914108 +DEAL:: delta=-0.002264482121 -1.207664857 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1207907202 +DEAL:: ||f*|| =0.0002863941104 +DEAL:: ||f*||_A =0.0002692200273 +DEAL::Newton iteration 1 for unit point guess 0.6796960273 1.004473446 +DEAL:: delta=-0.0003040054862 3.208354073e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0002863941104 +DEAL:: ||f*|| =4.559590695e-08 +DEAL:: ||f*||_A =4.558469573e-07 +DEAL::Newton iteration 2 for unit point guess 0.6800000328 1.004441363 +DEAL:: delta=1.094818045e-09 4.558424191e-07 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =4.559590695e-08 +DEAL:: ||f*|| =4.518280360e-16 +DEAL:: ||f*||_A =4.518258255e-15 +DEAL::Iteration converged for p_unit = [ 0.6800000317 1.004440907 ] and iteration error 4.518258255e-15 +DEAL:: +DEAL::Testing on cell 0_1:2 with center 0.7361215932 0.4250000000 +DEAL:: +DEAL::Testing on cell 0_1:3 with center 5.204748896e-17 0.8500000000 +DEAL::Newton iteration 0 for unit point guess 0.6983058446 -1.203191411 +DEAL:: delta=0.01613817896 -1.209285116 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1221516894 +DEAL:: ||f*|| =0.002048430507 +DEAL:: ||f*||_A =0.002304452941 +DEAL::Newton iteration 1 for unit point guess 0.6821676656 0.006093705129 +DEAL:: delta=0.002168030425 0.001629622689 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.002048430507 +DEAL:: ||f*|| =2.347437833e-06 +DEAL:: ||f*||_A =2.318708831e-05 +DEAL::Newton iteration 2 for unit point guess 0.6799996352 0.004464082440 +DEAL:: delta=-3.964683032e-07 2.317537267e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.347437833e-06 +DEAL:: ||f*|| =9.653137149e-13 +DEAL:: ||f*||_A =1.283764955e-12 +DEAL::Iteration converged for p_unit = [ 0.6800000317 0.004440907068 ] and iteration error 1.283764955e-12 +DEAL:: +DEAL::Testing on cell 1_1:0 with center -0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:1 with center -0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:2 with center -0.7361215932 0.4250000000 +DEAL::Newton iteration 0 for unit point guess -0.3011286728 2.212804616 +DEAL:: delta=0.02498309308 2.206435877 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2213639846 +DEAL:: ||f*|| =0.005779652832 +DEAL:: ||f*||_A =0.008680805469 +DEAL::Newton iteration 1 for unit point guess -0.3261117659 0.006368738899 +DEAL:: delta=-0.006133082543 0.001564845993 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.005779652832 +DEAL:: ||f*|| =1.854230994e-05 +DEAL:: ||f*||_A =0.0001851872864 +DEAL::Newton iteration 2 for unit point guess -0.3199786833 0.004803892906 +DEAL:: delta=1.110956583e-06 0.0001851213886 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.854230994e-05 +DEAL:: ||f*|| =2.156950753e-11 +DEAL:: ||f*||_A =2.371298577e-11 +DEAL::Newton iteration 3 for unit point guess -0.3199797943 0.004618771518 +DEAL:: delta=-2.288319216e-11 6.216154557e-12 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.156950753e-11 +DEAL:: ||f*|| =2.992994541e-14 +DEAL:: ||f*||_A =2.485467891e-13 +DEAL::Iteration converged for p_unit = [ -0.3199797942 0.004618771511 ] and iteration error 2.485467891e-13 +DEAL:: +DEAL::Testing on cell 1_1:3 with center -0.7361215932 -0.4250000000 +DEAL:: +DEAL::Testing on cell 2_1:0 with center -1.018281668e-15 -0.9500000000 +DEAL:: +DEAL::Testing on cell 2_1:1 with center 0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 2_1:2 with center -9.110941236e-16 -0.8500000000 +DEAL:: +DEAL::Testing on cell 2_1:3 with center 0.7361215932 -0.4250000000 +DEAL:: +DEAL:: +DEAL::Testing 2D with point -0.1779244106 0.9327124237 +DEAL::Testing on cell 0_1:0 with center 0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 0_1:1 with center 5.817072296e-17 0.9500000000 +DEAL::Newton iteration 0 for unit point guess 0.6872888532 -0.7700353780 +DEAL:: delta=0.006426066546 -1.274999600 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1277057819 +DEAL:: ||f*|| =0.0008583113365 +DEAL:: ||f*||_A =0.0007905789066 +DEAL::Newton iteration 1 for unit point guess 0.6808627867 0.5049642219 +DEAL:: delta=0.0008627800189 0.0002727224011 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0008583113365 +DEAL:: ||f*|| =3.883324301e-07 +DEAL:: ||f*||_A =3.875652879e-06 +DEAL::Newton iteration 2 for unit point guess 0.6800000067 0.5046914995 +DEAL:: delta=-2.501310404e-08 3.875351281e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.883324301e-07 +DEAL:: ||f*|| =1.018035770e-14 +DEAL:: ||f*||_A =1.282440299e-14 +DEAL::Iteration converged for p_unit = [ 0.6800000317 0.5046876241 ] and iteration error 1.282440299e-14 +DEAL:: +DEAL::Testing on cell 0_1:2 with center 0.7361215932 0.4250000000 +DEAL:: +DEAL::Testing on cell 0_1:3 with center 5.204748896e-17 0.8500000000 +DEAL::Newton iteration 0 for unit point guess 0.7093228359 -1.770035378 +DEAL:: delta=0.02584814312 -1.279199240 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1311997489 +DEAL:: ||f*|| =0.003483086009 +DEAL:: ||f*||_A =0.004640017241 +DEAL::Newton iteration 1 for unit point guess 0.6834746928 -0.4908361385 +DEAL:: delta=0.003476292242 0.004413379362 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.003483086009 +DEAL:: ||f*|| =6.491662569e-06 +DEAL:: ||f*||_A =6.293718017e-05 +DEAL::Newton iteration 2 for unit point guess 0.6799984006 -0.4952495179 +DEAL:: delta=-1.631085849e-06 6.285798976e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =6.491662569e-06 +DEAL:: ||f*|| =1.083964827e-11 +DEAL:: ||f*||_A =1.812776080e-11 +DEAL::Newton iteration 3 for unit point guess 0.6800000317 -0.4953123759 +DEAL:: delta=1.080251783e-11 1.455728296e-11 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.083964827e-11 +DEAL:: ||f*|| =5.432055468e-14 +DEAL:: ||f*||_A =5.431295206e-13 +DEAL::Iteration converged for p_unit = [ 0.6800000317 -0.4953123759 ] and iteration error 5.431295206e-13 +DEAL:: +DEAL::Testing on cell 1_1:0 with center -0.8227241336 0.4750000000 +DEAL::Newton iteration 0 for unit point guess -0.2566215243 2.835738205 +DEAL:: delta=0.08391430091 2.309996094 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2394359386 +DEAL:: ||f*|| =0.02052737244 +DEAL:: ||f*||_A =0.04536287509 +DEAL::Newton iteration 1 for unit point guess -0.3405358252 0.5257421110 +DEAL:: delta=-0.02059923664 0.01867124727 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.02052737244 +DEAL:: ||f*|| =0.0002237277955 +DEAL:: ||f*||_A =0.002204559446 +DEAL::Newton iteration 2 for unit point guess -0.3199365886 0.5070708637 +DEAL:: delta=4.321568008e-05 0.002195484093 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0002237277955 +DEAL:: ||f*|| =9.986962090e-09 +DEAL:: ||f*||_A =1.393354314e-08 +DEAL::Newton iteration 3 for unit point guess -0.3199798042 0.5048753796 +DEAL:: delta=-9.994403471e-09 9.701703420e-09 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =9.986962090e-09 +DEAL:: ||f*|| =1.175801663e-14 +DEAL:: ||f*||_A =7.866427655e-14 +DEAL::Iteration converged for p_unit = [ -0.3199797942 0.5048753699 ] and iteration error 7.866427655e-14 +DEAL:: +DEAL::Testing on cell 1_1:1 with center -0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:2 with center -0.7361215932 0.4250000000 +DEAL::Newton iteration 0 for unit point guess -0.3456358213 1.835738205 +DEAL:: delta=-0.03399801319 2.327443055 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2341373451 +DEAL:: ||f*|| =0.008302069701 +DEAL:: ||f*||_A =0.01256092621 +DEAL::Newton iteration 1 for unit point guess -0.3116378081 -0.4917048501 +DEAL:: delta=0.008344955586 0.003058087920 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.008302069701 +DEAL:: ||f*|| =3.628962144e-05 +DEAL:: ||f*||_A =0.0003619456002 +DEAL::Newton iteration 2 for unit point guess -0.3199827637 -0.4947629380 +DEAL:: delta=-2.969540907e-06 0.0003616920856 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.628962144e-05 +DEAL:: ||f*|| =1.129384400e-10 +DEAL:: ||f*||_A =1.384346268e-10 +DEAL::Newton iteration 3 for unit point guess -0.3199797941 -0.4951246301 +DEAL:: delta=1.132731643e-10 7.957265822e-11 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.129384400e-10 +DEAL:: ||f*|| =2.738251725e-12 +DEAL:: ||f*||_A =2.058225514e-11 +DEAL::Iteration converged for p_unit = [ -0.3199797942 -0.4951246301 ] and iteration error 2.058225514e-11 +DEAL:: +DEAL::Testing on cell 1_1:3 with center -0.7361215932 -0.4250000000 +DEAL:: +DEAL::Testing on cell 2_1:0 with center -1.018281668e-15 -0.9500000000 +DEAL:: +DEAL::Testing on cell 2_1:1 with center 0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 2_1:2 with center -9.110941236e-16 -0.8500000000 +DEAL:: +DEAL::Testing on cell 2_1:3 with center 0.7361215932 -0.4250000000 +DEAL:: +DEAL:: +DEAL::Testing 2D with point -0.1872888532 0.9818025513 +DEAL::Testing on cell 0_1:0 with center 0.8227241336 0.4750000000 +DEAL:: +DEAL::Testing on cell 0_1:1 with center 5.817072296e-17 0.9500000000 +DEAL::Newton iteration 0 for unit point guess 0.6971461613 -1.336879345 +DEAL:: delta=0.01511593008 -1.343424826 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.1355357751 +DEAL:: ||f*|| =0.002130899758 +DEAL:: ||f*||_A =0.002185874340 +DEAL::Newton iteration 1 for unit point guess 0.6820302312 0.006545480637 +DEAL:: delta=0.002030525298 0.001588550825 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.002130899758 +DEAL:: ||f*|| =2.284452314e-06 +DEAL:: ||f*||_A =2.259810758e-05 +DEAL::Newton iteration 2 for unit point guess 0.6799997059 0.004956929811 +DEAL:: delta=-3.257664409e-07 2.258862513e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.284452314e-06 +DEAL:: ||f*|| =7.727227667e-13 +DEAL:: ||f*||_A =9.381448999e-13 +DEAL::Iteration converged for p_unit = [ 0.6800000317 0.004934341186 ] and iteration error 9.381448999e-13 +DEAL:: +DEAL::Testing on cell 0_1:2 with center 0.7361215932 0.4250000000 +DEAL:: +DEAL::Testing on cell 0_1:3 with center 5.204748896e-17 0.8500000000 +DEAL:: +DEAL::Testing on cell 1_1:0 with center -0.8227241336 0.4750000000 +DEAL::Newton iteration 0 for unit point guess -0.2964437098 2.458671795 +DEAL:: delta=0.03119099912 2.450508799 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2462942128 +DEAL:: ||f*|| =0.008018379590 +DEAL:: ||f*||_A =0.01144192579 +DEAL::Newton iteration 1 for unit point guess -0.3276347089 0.008162996156 +DEAL:: delta=-0.007657091941 0.002710484511 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.008018379590 +DEAL:: ||f*|| =3.213641675e-05 +DEAL:: ||f*||_A =0.0003207212052 +DEAL::Newton iteration 2 for unit point guess -0.3199776170 0.005452511644 +DEAL:: delta=2.177360430e-06 0.0003205432723 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.213641675e-05 +DEAL:: ||f*|| =7.317141595e-11 +DEAL:: ||f*||_A =7.455483695e-11 +DEAL::Newton iteration 3 for unit point guess -0.3199797943 0.005131968372 +DEAL:: delta=-6.984991340e-11 2.605790417e-11 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =7.317141595e-11 +DEAL:: ||f*|| =1.867816199e-14 +DEAL:: ||f*||_A =1.781362949e-13 +DEAL::Iteration converged for p_unit = [ -0.3199797942 0.005131968346 ] and iteration error 1.781362949e-13 +DEAL:: +DEAL::Testing on cell 1_1:1 with center -0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 1_1:2 with center -0.7361215932 0.4250000000 +DEAL::Newton iteration 0 for unit point guess -0.3901429698 1.458671795 +DEAL:: delta=-0.09289797061 2.426643373 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2534993347 +DEAL:: ||f*|| =0.02392012964 +DEAL:: ||f*||_A =0.05560061155 +DEAL::Newton iteration 1 for unit point guess -0.2972449991 -0.9679715781 +DEAL:: delta=0.02279447828 0.02406831932 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.02392012964 +DEAL:: ||f*|| =0.0002896231357 +DEAL:: ||f*||_A =0.002842698490 +DEAL::Newton iteration 2 for unit point guess -0.3200394774 -0.9920398974 +DEAL:: delta=-5.970025909e-05 0.002828115558 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0002896231357 +DEAL:: ||f*|| =1.774864200e-08 +DEAL:: ||f*||_A =2.573172498e-08 +DEAL::Newton iteration 3 for unit point guess -0.3199797772 -0.9948680130 +DEAL:: delta=1.685125470e-08 1.943121730e-08 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.774864200e-08 +DEAL:: ||f*|| =2.007891506e-11 +DEAL:: ||f*||_A =1.928835934e-10 +DEAL::Newton iteration 4 for unit point guess -0.3199797940 -0.9948680324 +DEAL:: delta=-5.347814829e-12 -1.928094443e-10 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.007891506e-11 +DEAL:: ||f*|| =8.271995739e-11 +DEAL:: ||f*||_A =8.035405799e-11 +DEAL::Iteration converged for p_unit = [ -0.3199797940 -0.9948680322 ] and iteration error 8.035405799e-11 +DEAL:: +DEAL::Testing on cell 1_1:3 with center -0.7361215932 -0.4250000000 +DEAL:: +DEAL::Testing on cell 2_1:0 with center -1.018281668e-15 -0.9500000000 +DEAL:: +DEAL::Testing on cell 2_1:1 with center 0.8227241336 -0.4750000000 +DEAL:: +DEAL::Testing on cell 2_1:2 with center -9.110941236e-16 -0.8500000000 +DEAL:: +DEAL::Testing on cell 2_1:3 with center 0.7361215932 -0.4250000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point -0.1498310826 0.7854420410 0.02512860726 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL::Newton iteration 0 for unit point guess 3.702424762 1.255791956 0.5241800136 +DEAL:: delta=2.856388591 -0.7329178827 -0.07384132325 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5585835069 +DEAL:: ||f*|| =0.9496238329 +DEAL:: ||f*||_A =4.671562098 +DEAL:: step_length=0.5000000000 +DEAL:: ||f || =0.5585835069 +DEAL:: ||f*|| =0.3989025574 +DEAL:: ||f*||_A =2.107213631 +DEAL::Newton iteration 1 for unit point guess 2.274230467 1.622250898 0.5611006752 +DEAL::Newton iteration 1 for unit point guess 1.000000000 1.000000000 0.5611006752 +DEAL:: delta=-0.6775152184 -0.3855878006 -0.001349240320 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.4709858772 +DEAL:: ||f*|| =0.2020861084 +DEAL:: ||f*||_A =0.5469157809 +DEAL::Newton iteration 2 for unit point guess 1.677515218 1.385587801 0.5624499156 +DEAL:: delta=1.266079264 -0.2156126809 -0.008282391338 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2020861084 +DEAL:: ||f*|| =0.1049223850 +DEAL:: ||f*||_A =0.6628409621 +DEAL::Newton iteration 3 for unit point guess 0.4114359540 1.601200481 0.5707323069 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL::Newton iteration 0 for unit point guess -1.802127607 0.3558249736 0.5241800136 +DEAL:: delta=-2.805512532 -0.01580652674 0.002582563452 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5620799794 +DEAL:: ||f*|| =0.01436684397 +DEAL:: ||f*||_A =0.006730022603 +DEAL::Newton iteration 1 for unit point guess 1.003384925 0.3716315003 0.5215974502 +DEAL:: delta=0.003147124861 -0.01106989251 0.001788042630 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.01436684397 +DEAL:: ||f*|| =0.0001299736755 +DEAL:: ||f*||_A =0.0006449800586 +DEAL::Newton iteration 2 for unit point guess 1.000237800 0.3827013928 0.5198094075 +DEAL:: delta=0.0006433441903 1.326333806e-05 -5.908202149e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0001299736755 +DEAL:: ||f*|| =2.979104282e-09 +DEAL:: ||f*||_A =2.561313531e-09 +DEAL::Newton iteration 3 for unit point guess 0.9995944559 0.3826881295 0.5198153157 +DEAL:: delta=1.078284357e-09 -2.127435559e-09 9.323111453e-10 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.979104282e-09 +DEAL:: ||f*|| =1.163170575e-16 +DEAL:: ||f*||_A =5.681413772e-16 +DEAL::Iteration converged for p_unit = [ 0.9995944549 0.3826881316 0.5198153148 ] and iteration error 5.681413772e-16 +DEAL:: +DEAL:: +DEAL::Testing 3D with point -0.1591955252 0.8345321686 0.02669914522 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL::Newton iteration 0 for unit point guess 3.621326310 1.303028954 0.5256912645 +DEAL:: delta=3.194052890 -0.6600674790 -0.07966380171 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5920447676 +DEAL:: ||f*|| =0.9165285292 +DEAL:: ||f*||_A =4.978904026 +DEAL:: step_length=0.5000000000 +DEAL:: ||f || =0.5920447676 +DEAL:: ||f*|| =0.4201155126 +DEAL:: ||f*||_A =2.306653962 +DEAL::Newton iteration 1 for unit point guess 2.024299865 1.633062693 0.5655231653 +DEAL::Newton iteration 1 for unit point guess 1.000000000 1.000000000 0.5655231653 +DEAL:: delta=-0.4706798901 -0.4099576263 0.0004091119754 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.4882994847 +DEAL:: ||f*|| =0.2160200439 +DEAL:: ||f*||_A =0.6434558405 +DEAL::Newton iteration 2 for unit point guess 1.470679890 1.409957626 0.5651140533 +DEAL:: delta=1.453515393 -0.2064922549 -0.007935435237 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2160200439 +DEAL:: ||f*|| =0.1139084658 +DEAL:: ||f*||_A =0.7744783373 +DEAL::Newton iteration 3 for unit point guess 0.01716449706 1.616449881 0.5730494886 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL::Newton iteration 0 for unit point guess -2.227260583 0.3468140344 0.5256912645 +DEAL:: delta=-2.984004804 -0.02111069061 0.003484766185 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5987206642 +DEAL:: ||f*|| =0.02039700604 +DEAL:: ||f*||_A =0.009167287651 +DEAL::Newton iteration 1 for unit point guess 0.7567442217 0.3679247250 0.5222064983 +DEAL:: delta=0.005956738232 -0.01479312618 0.002402897797 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.02039700604 +DEAL:: ||f*|| =0.0002474545613 +DEAL:: ||f*||_A =0.001223428251 +DEAL::Newton iteration 2 for unit point guess 0.7507874835 0.3827178512 0.5198036005 +DEAL:: delta=0.001218369679 2.972808970e-05 -1.171759931e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0002474545613 +DEAL:: ||f*|| =1.243143907e-08 +DEAL:: ||f*||_A =1.066824344e-08 +DEAL::Newton iteration 3 for unit point guess 0.7495691138 0.3826881231 0.5198153181 +DEAL:: delta=5.551615407e-09 -8.494137003e-09 3.280125957e-09 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.243143907e-08 +DEAL:: ||f*|| =3.469446952e-18 +DEAL:: ||f*||_A =2.631025672e-18 +DEAL::Iteration converged for p_unit = [ 0.7495691083 0.3826881316 0.5198153148 ] and iteration error 2.631025672e-18 +DEAL:: +DEAL:: +DEAL::Testing 3D with point -0.1685599679 0.8836222961 0.02826968317 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL::Newton iteration 0 for unit point guess 3.540227857 1.350265951 0.5272025153 +DEAL:: delta=3.563108903 -0.5947761734 -0.08726992220 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6279729539 +DEAL:: ||f*|| =0.8983345290 +DEAL:: ||f*||_A =5.504053175 +DEAL:: step_length=0.5000000000 +DEAL:: ||f || =0.6279729539 +DEAL:: ||f*|| =0.4439967188 +DEAL:: ||f*||_A =2.555152222 +DEAL::Newton iteration 1 for unit point guess 1.758673406 1.647654038 0.5708374764 +DEAL::Newton iteration 1 for unit point guess 1.000000000 1.000000000 0.5708374764 +DEAL:: delta=-0.2642670141 -0.4344439561 0.002772053466 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5100712180 +DEAL:: ||f*|| =0.2331871882 +DEAL:: ||f*||_A =0.7513804837 +DEAL::Newton iteration 2 for unit point guess 1.264267014 1.434443956 0.5680654229 +DEAL:: delta=1.670051431 -0.2002124112 -0.007724121199 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2331871882 +DEAL:: ||f*|| =0.1274152335 +DEAL:: ||f*||_A =0.9518211610 +DEAL::Newton iteration 3 for unit point guess -0.4057844171 1.634656367 0.5757895441 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL::Newton iteration 0 for unit point guess -2.652393558 0.3378030953 0.5272025153 +DEAL:: delta=-3.163815489 -0.02641498558 0.004391680086 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6359961941 +DEAL:: ||f*|| =0.02704151679 +DEAL:: ||f*||_A =0.01179290428 +DEAL::Newton iteration 1 for unit point guess 0.5114219311 0.3642180808 0.5228108352 +DEAL:: delta=0.009859953041 -0.01852606978 0.003015648198 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.02704151679 +DEAL:: ||f*|| =0.0004125306161 +DEAL:: ||f*||_A =0.002030597423 +DEAL::Newton iteration 2 for unit point guess 0.5015619781 0.3827441506 0.5197951870 +DEAL:: delta=0.002018195959 5.604404835e-05 -2.013654738e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0004125306161 +DEAL:: ||f*|| =3.843012297e-08 +DEAL:: ||f*||_A =3.350019595e-08 +DEAL::Newton iteration 3 for unit point guess 0.4995437821 0.3826881066 0.5198153236 +DEAL:: delta=2.042052817e-08 -2.503750219e-08 8.776975331e-09 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.843012297e-08 +DEAL:: ||f*|| =7.098814035e-16 +DEAL:: ||f*||_A =3.132643602e-15 +DEAL::Iteration converged for p_unit = [ 0.4995437617 0.3826881316 0.5198153148 ] and iteration error 3.132643602e-15 +DEAL:: +DEAL:: +DEAL::Testing 3D with point -0.1779244106 0.9327124237 0.02984022112 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL::Newton iteration 0 for unit point guess 3.459129405 1.397502948 0.5287137662 +DEAL:: delta=3.979828397 -0.5423958605 -0.09739294423 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6668294077 +DEAL:: ||f*|| =0.9171795318 +DEAL:: ||f*||_A =6.576519582 +DEAL:: step_length=0.5000000000 +DEAL:: ||f || =0.6668294077 +DEAL:: ||f*|| =0.4737130779 +DEAL:: ||f*||_A =2.906475606 +DEAL::Newton iteration 1 for unit point guess 1.469215207 1.668700878 0.5774102383 +DEAL::Newton iteration 1 for unit point guess 1.000000000 1.000000000 0.5774102383 +DEAL:: delta=-0.05852563526 -0.4591073408 0.005969607497 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5358286371 +DEAL:: ||f*|| =0.2540023766 +DEAL:: ||f*||_A =0.8710470496 +DEAL::Newton iteration 2 for unit point guess 1.058525635 1.459107341 0.5714406308 +DEAL:: delta=1.928810109 -0.1985615361 -0.007803797499 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2540023766 +DEAL:: ||f*|| =0.1498020872 +DEAL:: ||f*||_A =1.270586809 +DEAL::Newton iteration 3 for unit point guess -0.8702844736 1.657668877 0.5792444283 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL::Newton iteration 0 for unit point guess -3.077526534 0.3287921561 0.5287137662 +DEAL:: delta=-3.345097096 -0.03171894519 0.005303684596 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6739766956 +DEAL:: ||f*|| =0.03430597763 +DEAL:: ||f*||_A =0.01468153345 +DEAL::Newton iteration 1 for unit point guess 0.2675705624 0.3605111013 0.5234100816 +DEAL:: delta=0.01498260400 -0.02227148874 0.003626316434 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.03430597763 +DEAL:: ||f*|| =0.0006321096333 +DEAL:: ||f*||_A =0.003095644194 +DEAL::Newton iteration 2 for unit point guess 0.2525879584 0.3827825900 0.5197837651 +DEAL:: delta=0.003069482877 9.451923784e-05 -3.156939906e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0006321096333 +DEAL:: ||f*|| =9.796087467e-08 +DEAL:: ||f*||_A =8.800136986e-08 +DEAL::Newton iteration 3 for unit point guess 0.2495184755 0.3826880708 0.5198153345 +DEAL:: delta=6.035640101e-08 -6.080731181e-08 1.973325093e-08 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =9.796087467e-08 +DEAL:: ||f*|| =5.689389486e-15 +DEAL:: ||f*||_A =2.742618229e-14 +DEAL::Iteration converged for p_unit = [ 0.2495184151 0.3826881316 0.5198153148 ] and iteration error 2.742618229e-14 +DEAL:: +DEAL:: +DEAL::Testing 3D with point -0.1872888532 0.9818025513 0.03141075908 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL::Newton iteration 0 for unit point guess 3.378030953 1.444739945 0.5302250170 +DEAL:: delta=4.481775076 -0.5164923720 -0.1114133430 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.7090545024 +DEAL:: ||f*|| =1.040376963 +DEAL:: ||f*||_A =9.431636223 +DEAL:: step_length=0.5000000000 +DEAL:: ||f || =0.7090545024 +DEAL:: ||f*|| =0.5189006404 +DEAL:: ||f*||_A =3.548765003 +DEAL::Newton iteration 1 for unit point guess 1.137143415 1.702986131 0.5859316885 +DEAL::Newton iteration 1 for unit point guess 1.000000000 1.000000000 0.5859316885 +DEAL:: delta=0.1460488594 -0.4840658966 0.01041830046 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5651606508 +DEAL:: ||f*|| =0.2788717609 +DEAL:: ||f*||_A =1.003004700 +DEAL::Newton iteration 2 for unit point guess 0.8539511406 1.484065897 0.5755133880 +DEAL:: delta=2.256510264 -0.2053129307 -0.008548455215 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.2788717609 +DEAL:: ||f*|| =0.1920633747 +DEAL:: ||f*||_A =1.953837842 +DEAL::Newton iteration 3 for unit point guess -1.402559123 1.689378827 0.5840618433 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.04694448799 0.06461354454 0.7960033322 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL::Newton iteration 0 for unit point guess 0.5451723546 -1.893591072 0.5621744122 +DEAL:: delta=0.004969656478 -2.894555770 0.006855209102 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.5791702826 +DEAL:: ||f*|| =0.007834938277 +DEAL:: ||f*||_A =0.003571810880 +DEAL::Newton iteration 1 for unit point guess 0.5402026982 1.000964698 0.5553192031 +DEAL:: delta=0.003587465407 0.0009021921021 0.004953715025 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.007834938277 +DEAL:: ||f*|| =3.884066112e-05 +DEAL:: ||f*||_A =0.0001929674524 +DEAL::Newton iteration 2 for unit point guess 0.5366152327 1.000062506 0.5503654881 +DEAL:: delta=-2.660056699e-06 0.0001927096002 -2.683324683e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =3.884066112e-05 +DEAL:: ||f*|| =2.326926478e-10 +DEAL:: ||f*||_A =1.956074025e-10 +DEAL::Newton iteration 3 for unit point guess 0.5366178928 0.9998697960 0.5503681714 +DEAL:: delta=1.274583800e-10 7.366186838e-11 1.287881821e-10 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.326926478e-10 +DEAL:: ||f*|| =3.379462172e-16 +DEAL:: ||f*||_A =1.639763864e-15 +DEAL::Iteration converged for p_unit = [ 0.5366178927 0.9998697959 0.5503681713 ] and iteration error 1.639763864e-15 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.04890050833 0.06730577556 0.8291701377 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL::Newton iteration 0 for unit point guess 0.5470545361 -2.180824033 0.5647650127 +DEAL:: delta=0.006064776973 -3.015717358 0.008360635756 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6035598825 +DEAL:: ||f*|| =0.009951748567 +DEAL:: ||f*||_A =0.004371343605 +DEAL::Newton iteration 1 for unit point guess 0.5409897591 0.8348933245 0.5564043770 +DEAL:: delta=0.004376158946 0.001397238017 0.006040640167 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.009951748567 +DEAL:: ||f*|| =6.023554465e-05 +DEAL:: ||f*||_A =0.0002989039364 +DEAL::Newton iteration 2 for unit point guess 0.5366136002 0.8334960865 0.5503637368 +DEAL:: delta=-4.292825751e-06 0.0002983822338 -4.434798063e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =6.023554465e-05 +DEAL:: ||f*|| =5.885477122e-10 +DEAL:: ||f*||_A =4.851415704e-10 +DEAL::Newton iteration 3 for unit point guess 0.5366178930 0.8331977043 0.5503681716 +DEAL:: delta=3.055682214e-10 2.047476826e-10 3.162746966e-10 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =5.885477122e-10 +DEAL:: ||f*|| =5.003707553e-17 +DEAL:: ||f*||_A =3.836115168e-17 +DEAL::Iteration converged for p_unit = [ 0.5366178927 0.8331977041 0.5503681713 ] and iteration error 3.836115168e-17 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.05085652866 0.06999800658 0.8623369432 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL::Newton iteration 0 for unit point guess 0.5489367175 -2.468056995 0.5673556132 +DEAL:: delta=0.007160560813 -3.137037000 0.009866511490 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6280236718 +DEAL:: ||f*|| =0.01221120457 +DEAL:: ||f*||_A =0.005181873655 +DEAL::Newton iteration 1 for unit point guess 0.5417761567 0.6689800054 0.5574891018 +DEAL:: delta=0.005164708567 0.002022568282 0.007127722446 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.01221120457 +DEAL:: ||f*|| =8.732694438e-05 +DEAL:: ||f*||_A =0.0004327777066 +DEAL::Newton iteration 2 for unit point guess 0.5366114481 0.6669574371 0.5503613793 +DEAL:: delta=-6.445176340e-06 0.0004318243265 -6.792659995e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =8.732694438e-05 +DEAL:: ||f*|| =1.292239418e-09 +DEAL:: ||f*||_A =1.049589754e-09 +DEAL::Newton iteration 3 for unit point guess 0.5366178933 0.6665256127 0.5503681720 +DEAL:: delta=6.380390351e-10 4.900768629e-10 6.738742064e-10 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =1.292239418e-09 +DEAL:: ||f*|| =1.187747315e-16 +DEAL:: ||f*||_A =5.386767756e-16 +DEAL::Iteration converged for p_unit = [ 0.5366178927 0.6665256123 0.5503681713 ] and iteration error 5.386767756e-16 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.05281254899 0.07269023761 0.8955037488 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL::Newton iteration 0 for unit point guess 0.5508188990 -2.755289956 0.5699462137 +DEAL:: delta=0.008257034834 -3.258527571 0.01137284273 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6525676654 +DEAL:: ||f*|| =0.01461334118 +DEAL:: ||f*||_A =0.006006826699 +DEAL::Newton iteration 1 for unit point guess 0.5425618641 0.5032376147 0.5585733710 +DEAL:: delta=0.005953152018 0.002788806072 0.008215041053 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.01461334118 +DEAL:: ||f*|| =0.0001206127230 +DEAL:: ||f*||_A =0.0005969040839 +DEAL::Newton iteration 2 for unit point guess 0.5366087121 0.5004488087 0.5503583300 +DEAL:: delta=-9.181772873e-06 0.0005952871909 -9.842617085e-06 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0001206127230 +DEAL:: ||f*|| =2.560943991e-09 +DEAL:: ||f*||_A =2.059237066e-09 +DEAL::Newton iteration 3 for unit point guess 0.5366178939 0.4998535215 0.5503681726 +DEAL:: delta=1.205909696e-09 1.051589847e-09 1.295765321e-09 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =2.560943991e-09 +DEAL:: ||f*|| =2.330476870e-16 +DEAL:: ||f*||_A =1.082227180e-15 +DEAL::Iteration converged for p_unit = [ 0.5366178927 0.4998535204 0.5503681713 ] and iteration error 1.082227180e-15 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.05476856932 0.07538246863 0.9286705543 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL::Newton iteration 0 for unit point guess 0.5527010804 -3.042522917 0.5725368143 +DEAL:: delta=0.009354226180 -3.380201895 0.01287963629 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.6771978331 +DEAL:: ||f*|| =0.01715820711 +DEAL:: ||f*||_A =0.006850101626 +DEAL::Newton iteration 1 for unit point guess 0.5433468542 0.3376789780 0.5596571780 +DEAL:: delta=0.006741526547 0.003706534631 0.009302674610 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.01715820711 +DEAL:: ||f*|| =0.0001605975923 +DEAL:: ||f*||_A =0.0007936021325 +DEAL::Newton iteration 2 for unit point guess 0.5366053277 0.3339724433 0.5503545034 +DEAL:: delta=-1.256711300e-05 0.0007910126895 -1.367023795e-05 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =0.0001605975923 +DEAL:: ||f*|| =4.695015766e-09 +DEAL:: ||f*||_A =3.754227887e-09 +DEAL::Newton iteration 3 for unit point guess 0.5366178948 0.3331814307 0.5503681736 +DEAL:: delta=2.113650368e-09 2.075279878e-09 2.305167360e-09 +DEAL:: step_length=1.000000000 +DEAL:: ||f || =4.695015766e-09 +DEAL:: ||f*|| =1.316562551e-16 +DEAL:: ||f*||_A =5.860138287e-16 +DEAL::Iteration converged for p_unit = [ 0.5366178927 0.3331814286 0.5503681713 ] and iteration error 5.860138287e-16 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.05672458966 0.07807469965 0.9618373598 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: +DEAL::Testing 3D with point 0.05868060999 0.08076693067 0.9950041653 +DEAL::Testing on cell 0_0: with center -9.723960037e-18 -9.723960037e-18 -0.9000000000 +DEAL:: +DEAL::Testing on cell 1_0: with center 0.9000000000 -1.696561603e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 2_0: with center -1.696561603e-17 0.000000000 0.9000000000 +DEAL:: +DEAL::Testing on cell 3_0: with center -0.9000000000 -5.358086410e-17 0.000000000 +DEAL:: +DEAL::Testing on cell 4_0: with center -1.334478803e-17 -0.9000000000 0.000000000 +DEAL:: +DEAL::Testing on cell 5_0: with center -5.358086410e-17 0.9000000000 0.000000000 +DEAL:: +DEAL:: -- 2.39.5