From: Maien Hamed Date: Tue, 8 Sep 2015 10:53:09 +0000 (+0200) Subject: added tests for other elements X-Git-Tag: v8.4.0-rc2~461^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d6b8a0b5b19a81746c454715f77c86fa2b09eb6e;p=dealii.git added tests for other elements --- diff --git a/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc new file mode 100644 index 0000000000..0e0c35a146 --- /dev/null +++ b/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.cc @@ -0,0 +1,177 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2003 - 2015 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 the correctness of fe_values.shape_3rd_derivative for FE_DGP by comparing +// the integral of all shape third derivative components with the flux of the +// hessian over the boundary by the divergence theorem + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include + +template +Tensor<1,dim> ones () +{ + Tensor<1,dim> result; + for (unsigned int i=0; i +void test (const Triangulation &tr, + const FiniteElement &fe, + const double tolerance) +{ + DoFHandler dof(tr); + dof.distribute_dofs(fe); + + std::stringstream ss; + + deallog << "FE=" << fe.get_name() << std::endl; + + const QGauss quadrature(6); + FEValues fe_values (fe, quadrature, update_3rd_derivatives + | update_quadrature_points + | update_JxW_values); + + const QGauss face_quadrature(6); + FEFaceValues fe_face_values (fe, face_quadrature, + update_hessians + | update_quadrature_points + | update_normal_vectors + | update_JxW_values); + + for (typename DoFHandler::active_cell_iterator cell = dof.begin_active(); + cell != dof.end(); + ++cell) + { + fe_values.reinit (cell); + + deallog << "Cell nodes:" << std::endl; + for (unsigned int i=0; i::vertices_per_cell; ++i) + { + deallog << i << ": ( "; + for (unsigned int d=0; dvertex(i)[d] << " "; + deallog << ")" << std::endl; + } + + bool cell_ok = true; + + for (unsigned int c=0; c +void test_hyper_ball(const double tolerance) +{ + Triangulation tr; + GridGenerator::hyper_ball(tr); + + static const HyperBallBoundary boundary; + tr.set_boundary (0, boundary); + + tr.refine_global(1); + + FE_DGP fe(2); + test(tr, fe, tolerance); +} + + +int main() +{ + std::ofstream logfile ("output"); + deallog << std::setprecision (3); + + deallog.attach(logfile); + deallog.depth_console (0); + deallog.threshold_double(1.e-7); + + test_hyper_ball<2>(1e-6); + test_hyper_ball<3>(1e-6); + + deallog << "done..." << std::endl; +} + diff --git a/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.output b/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.output new file mode 100644 index 0000000000..e758ca449b --- /dev/null +++ b/tests/fe/fe_dgp_3rd_derivative_divergence_theorem.output @@ -0,0 +1,760 @@ + +DEAL::FE=FE_DGP<2>(2) +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( 0 -1.00 ) +DEAL::2: ( -0.500 -0.500 ) +DEAL::3: ( 0 -0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 ) +DEAL::1: ( 0.707 -0.707 ) +DEAL::2: ( 0 -0.646 ) +DEAL::3: ( 0.500 -0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( 0 -0.646 ) +DEAL::2: ( -0.293 -0.293 ) +DEAL::3: ( 0 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.646 ) +DEAL::1: ( 0.500 -0.500 ) +DEAL::2: ( 0 -0.293 ) +DEAL::3: ( 0.293 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( -0.500 -0.500 ) +DEAL::2: ( -1.00 0 ) +DEAL::3: ( -0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( -0.293 -0.293 ) +DEAL::2: ( -0.646 0 ) +DEAL::3: ( -0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 ) +DEAL::1: ( -0.646 0 ) +DEAL::2: ( -0.707 0.707 ) +DEAL::3: ( -0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.646 0 ) +DEAL::1: ( -0.293 0 ) +DEAL::2: ( -0.500 0.500 ) +DEAL::3: ( -0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 -0.293 ) +DEAL::1: ( 0 -0.293 ) +DEAL::2: ( -0.293 0 ) +DEAL::3: ( 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.293 ) +DEAL::1: ( 0.293 -0.293 ) +DEAL::2: ( 0 0 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 0 ) +DEAL::1: ( 0 0 ) +DEAL::2: ( -0.293 0.293 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 ) +DEAL::1: ( 0.293 0 ) +DEAL::2: ( 0 0.293 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 ) +DEAL::1: ( 1.00 0 ) +DEAL::2: ( 0.500 -0.500 ) +DEAL::3: ( 0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 ) +DEAL::1: ( 0.707 0.707 ) +DEAL::2: ( 0.646 0 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.500 -0.500 ) +DEAL::1: ( 0.646 0 ) +DEAL::2: ( 0.293 -0.293 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.646 0 ) +DEAL::1: ( 0.500 0.500 ) +DEAL::2: ( 0.293 0 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 ) +DEAL::1: ( -0.500 0.500 ) +DEAL::2: ( 0 1.00 ) +DEAL::3: ( 0 0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 0.500 ) +DEAL::1: ( -0.293 0.293 ) +DEAL::2: ( 0 0.646 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 ) +DEAL::1: ( 0 0.646 ) +DEAL::2: ( 0.707 0.707 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.646 ) +DEAL::1: ( 0 0.293 ) +DEAL::2: ( 0.500 0.500 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::FE=FE_DGP<3>(2) +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 -0.211 ) +DEAL::1: ( 0 -0.211 -0.211 ) +DEAL::2: ( -0.211 0 -0.211 ) +DEAL::3: ( 0 0 -0.211 ) +DEAL::4: ( -0.211 -0.211 0 ) +DEAL::5: ( 0 -0.211 0 ) +DEAL::6: ( -0.211 0 0 ) +DEAL::7: ( 0 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 -0.211 ) +DEAL::1: ( 0.211 -0.211 -0.211 ) +DEAL::2: ( 0 0 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0 -0.211 0 ) +DEAL::5: ( 0.211 -0.211 0 ) +DEAL::6: ( 0 0 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 -0.211 ) +DEAL::1: ( 0 0 -0.211 ) +DEAL::2: ( -0.211 0.211 -0.211 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.211 0 0 ) +DEAL::5: ( 0 0 0 ) +DEAL::6: ( -0.211 0.211 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.211 ) +DEAL::1: ( 0.211 0 -0.211 ) +DEAL::2: ( 0 0.211 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0 0 ) +DEAL::5: ( 0.211 0 0 ) +DEAL::6: ( 0 0.211 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 0 ) +DEAL::1: ( 0 -0.211 0 ) +DEAL::2: ( -0.211 0 0 ) +DEAL::3: ( 0 0 0 ) +DEAL::4: ( -0.211 -0.211 0.211 ) +DEAL::5: ( 0 -0.211 0.211 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 0 ) +DEAL::1: ( 0.211 -0.211 0 ) +DEAL::2: ( 0 0 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0 -0.211 0.211 ) +DEAL::5: ( 0.211 -0.211 0.211 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 0 ) +DEAL::1: ( 0 0 0 ) +DEAL::2: ( -0.211 0.211 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.211 0 0.211 ) +DEAL::5: ( 0 0 0.211 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0 ) +DEAL::1: ( 0.211 0 0 ) +DEAL::2: ( 0 0.211 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0 0.211 ) +DEAL::5: ( 0.211 0 0.211 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( 0 0 -1.00 ) +DEAL::4: ( -0.394 -0.394 -0.394 ) +DEAL::5: ( 0 -0.419 -0.419 ) +DEAL::6: ( -0.419 0 -0.419 ) +DEAL::7: ( 0 0 -0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 0 -1.00 ) +DEAL::3: ( 0.707 0 -0.707 ) +DEAL::4: ( 0 -0.419 -0.419 ) +DEAL::5: ( 0.394 -0.394 -0.394 ) +DEAL::6: ( 0 0 -0.457 ) +DEAL::7: ( 0.419 0 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( 0 0 -1.00 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( 0 0.707 -0.707 ) +DEAL::4: ( -0.419 0 -0.419 ) +DEAL::5: ( 0 0 -0.457 ) +DEAL::6: ( -0.394 0.394 -0.394 ) +DEAL::7: ( 0 0.419 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -1.00 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0.577 0.577 -0.577 ) +DEAL::4: ( 0 0 -0.457 ) +DEAL::5: ( 0.419 0 -0.419 ) +DEAL::6: ( 0 0.419 -0.419 ) +DEAL::7: ( 0.394 0.394 -0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( 0 0 -0.457 ) +DEAL::4: ( -0.211 -0.211 -0.211 ) +DEAL::5: ( 0 -0.211 -0.211 ) +DEAL::6: ( -0.211 0 -0.211 ) +DEAL::7: ( 0 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 0 -0.457 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0 -0.211 -0.211 ) +DEAL::5: ( 0.211 -0.211 -0.211 ) +DEAL::6: ( 0 0 -0.211 ) +DEAL::7: ( 0.211 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( 0 0 -0.457 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.211 0 -0.211 ) +DEAL::5: ( 0 0 -0.211 ) +DEAL::6: ( -0.211 0.211 -0.211 ) +DEAL::7: ( 0 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.457 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 0 -0.211 ) +DEAL::5: ( 0.211 0 -0.211 ) +DEAL::6: ( 0 0.211 -0.211 ) +DEAL::7: ( 0.211 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.577 -0.577 -0.577 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0.394 -0.394 -0.394 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0.707 -0.707 0 ) +DEAL::5: ( 1.00 0 0 ) +DEAL::6: ( 0.419 -0.419 0 ) +DEAL::7: ( 0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 0 -0.707 ) +DEAL::1: ( 0.577 0.577 -0.577 ) +DEAL::2: ( 0.419 0 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 1.00 0 0 ) +DEAL::5: ( 0.707 0.707 0 ) +DEAL::6: ( 0.457 0 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.394 -0.394 -0.394 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0.211 -0.211 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0.419 -0.419 0 ) +DEAL::5: ( 0.457 0 0 ) +DEAL::6: ( 0.211 -0.211 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 0 -0.419 ) +DEAL::1: ( 0.394 0.394 -0.394 ) +DEAL::2: ( 0.211 0 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0.457 0 0 ) +DEAL::5: ( 0.419 0.419 0 ) +DEAL::6: ( 0.211 0 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 0 ) +DEAL::1: ( 1.00 0 0 ) +DEAL::2: ( 0.419 -0.419 0 ) +DEAL::3: ( 0.457 0 0 ) +DEAL::4: ( 0.577 -0.577 0.577 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0.394 -0.394 0.394 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 0 ) +DEAL::1: ( 0.707 0.707 0 ) +DEAL::2: ( 0.457 0 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0.707 0 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0.419 0 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 -0.419 0 ) +DEAL::1: ( 0.457 0 0 ) +DEAL::2: ( 0.211 -0.211 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0.394 -0.394 0.394 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0.211 -0.211 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.457 0 0 ) +DEAL::1: ( 0.419 0.419 0 ) +DEAL::2: ( 0.211 0 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0.419 0 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0.211 0 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 0.577 ) +DEAL::1: ( 0 -0.707 0.707 ) +DEAL::2: ( -0.394 -0.394 0.394 ) +DEAL::3: ( 0 -0.419 0.419 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( 0 0 1.00 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( 0 0 0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 0.707 ) +DEAL::1: ( 0.577 -0.577 0.577 ) +DEAL::2: ( 0 -0.419 0.419 ) +DEAL::3: ( 0.394 -0.394 0.394 ) +DEAL::4: ( 0 0 1.00 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0 0 0.457 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 0.394 ) +DEAL::1: ( 0 -0.419 0.419 ) +DEAL::2: ( -0.211 -0.211 0.211 ) +DEAL::3: ( 0 -0.211 0.211 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( 0 0 0.457 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 0.419 ) +DEAL::1: ( 0.394 -0.394 0.394 ) +DEAL::2: ( 0 -0.211 0.211 ) +DEAL::3: ( 0.211 -0.211 0.211 ) +DEAL::4: ( 0 0 0.457 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 0.707 ) +DEAL::1: ( 0 0 1.00 ) +DEAL::2: ( -0.419 0 0.419 ) +DEAL::3: ( 0 0 0.457 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( 0 0.707 0.707 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 1.00 ) +DEAL::1: ( 0.707 0 0.707 ) +DEAL::2: ( 0 0 0.457 ) +DEAL::3: ( 0.419 0 0.419 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 0.419 ) +DEAL::1: ( 0 0 0.457 ) +DEAL::2: ( -0.211 0 0.211 ) +DEAL::3: ( 0 0 0.211 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0.457 ) +DEAL::1: ( 0.419 0 0.419 ) +DEAL::2: ( 0 0 0.211 ) +DEAL::3: ( 0.211 0 0.211 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( -0.394 -0.394 -0.394 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( -0.419 0 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( -0.419 -0.419 0 ) +DEAL::6: ( -1.00 0 0 ) +DEAL::7: ( -0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( -0.211 -0.211 -0.211 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( -0.211 0 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( -0.211 -0.211 0 ) +DEAL::6: ( -0.457 0 0 ) +DEAL::7: ( -0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( -0.419 0 -0.419 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( -0.394 0.394 -0.394 ) +DEAL::4: ( -1.00 0 0 ) +DEAL::5: ( -0.457 0 0 ) +DEAL::6: ( -0.707 0.707 0 ) +DEAL::7: ( -0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( -0.211 0 -0.211 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( -0.211 0.211 -0.211 ) +DEAL::4: ( -0.457 0 0 ) +DEAL::5: ( -0.211 0 0 ) +DEAL::6: ( -0.419 0.419 0 ) +DEAL::7: ( -0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( -0.419 -0.419 0 ) +DEAL::2: ( -1.00 0 0 ) +DEAL::3: ( -0.457 0 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( -0.394 -0.394 0.394 ) +DEAL::6: ( -0.707 0 0.707 ) +DEAL::7: ( -0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( -0.211 -0.211 0 ) +DEAL::2: ( -0.457 0 0 ) +DEAL::3: ( -0.211 0 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( -0.211 -0.211 0.211 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( -0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 0 ) +DEAL::1: ( -0.457 0 0 ) +DEAL::2: ( -0.707 0.707 0 ) +DEAL::3: ( -0.419 0.419 0 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( -0.419 0 0.419 ) +DEAL::6: ( -0.577 0.577 0.577 ) +DEAL::7: ( -0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.457 0 0 ) +DEAL::1: ( -0.211 0 0 ) +DEAL::2: ( -0.419 0.419 0 ) +DEAL::3: ( -0.211 0.211 0 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( -0.211 0 0.211 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( -0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.394 -0.394 -0.394 ) +DEAL::3: ( 0 -0.419 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( 0 -1.00 0 ) +DEAL::6: ( -0.419 -0.419 0 ) +DEAL::7: ( 0 -0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 -0.419 -0.419 ) +DEAL::3: ( 0.394 -0.394 -0.394 ) +DEAL::4: ( 0 -1.00 0 ) +DEAL::5: ( 0.707 -0.707 0 ) +DEAL::6: ( 0 -0.457 0 ) +DEAL::7: ( 0.419 -0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.211 -0.211 -0.211 ) +DEAL::3: ( 0 -0.211 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( 0 -0.457 0 ) +DEAL::6: ( -0.211 -0.211 0 ) +DEAL::7: ( 0 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 -0.211 -0.211 ) +DEAL::3: ( 0.211 -0.211 -0.211 ) +DEAL::4: ( 0 -0.457 0 ) +DEAL::5: ( 0.419 -0.419 0 ) +DEAL::6: ( 0 -0.211 0 ) +DEAL::7: ( 0.211 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( 0 -1.00 0 ) +DEAL::2: ( -0.419 -0.419 0 ) +DEAL::3: ( 0 -0.457 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( 0 -0.707 0.707 ) +DEAL::6: ( -0.394 -0.394 0.394 ) +DEAL::7: ( 0 -0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 0 ) +DEAL::1: ( 0.707 -0.707 0 ) +DEAL::2: ( 0 -0.457 0 ) +DEAL::3: ( 0.419 -0.419 0 ) +DEAL::4: ( 0 -0.707 0.707 ) +DEAL::5: ( 0.577 -0.577 0.577 ) +DEAL::6: ( 0 -0.419 0.419 ) +DEAL::7: ( 0.394 -0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( 0 -0.457 0 ) +DEAL::2: ( -0.211 -0.211 0 ) +DEAL::3: ( 0 -0.211 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( 0 -0.419 0.419 ) +DEAL::6: ( -0.211 -0.211 0.211 ) +DEAL::7: ( 0 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.457 0 ) +DEAL::1: ( 0.419 -0.419 0 ) +DEAL::2: ( 0 -0.211 0 ) +DEAL::3: ( 0.211 -0.211 0 ) +DEAL::4: ( 0 -0.419 0.419 ) +DEAL::5: ( 0.394 -0.394 0.394 ) +DEAL::6: ( 0 -0.211 0.211 ) +DEAL::7: ( 0.211 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 0.577 -0.577 ) +DEAL::1: ( -0.394 0.394 -0.394 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.707 0.707 0 ) +DEAL::5: ( -0.419 0.419 0 ) +DEAL::6: ( 0 1.00 0 ) +DEAL::7: ( 0 0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 0.394 -0.394 ) +DEAL::1: ( -0.211 0.211 -0.211 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.419 0.419 0 ) +DEAL::5: ( -0.211 0.211 0 ) +DEAL::6: ( 0 0.457 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.707 -0.707 ) +DEAL::1: ( 0 0.419 -0.419 ) +DEAL::2: ( 0.577 0.577 -0.577 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 1.00 0 ) +DEAL::5: ( 0 0.457 0 ) +DEAL::6: ( 0.707 0.707 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.419 -0.419 ) +DEAL::1: ( 0 0.211 -0.211 ) +DEAL::2: ( 0.394 0.394 -0.394 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0.457 0 ) +DEAL::5: ( 0 0.211 0 ) +DEAL::6: ( 0.419 0.419 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 0 ) +DEAL::1: ( -0.419 0.419 0 ) +DEAL::2: ( 0 1.00 0 ) +DEAL::3: ( 0 0.457 0 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( -0.394 0.394 0.394 ) +DEAL::6: ( 0 0.707 0.707 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0.419 0 ) +DEAL::1: ( -0.211 0.211 0 ) +DEAL::2: ( 0 0.457 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( -0.211 0.211 0.211 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 0 ) +DEAL::1: ( 0 0.457 0 ) +DEAL::2: ( 0.707 0.707 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( 0.577 0.577 0.577 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.457 0 ) +DEAL::1: ( 0 0.211 0 ) +DEAL::2: ( 0.419 0.419 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0 0.211 0.211 ) +DEAL::6: ( 0.394 0.394 0.394 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::done... diff --git a/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc new file mode 100644 index 0000000000..4d47c20d64 --- /dev/null +++ b/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.cc @@ -0,0 +1,176 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2003 - 2015 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 the correctness of fe_values.shape_3rd_derivative for FE_DGQ by comparing +// the integral of all shape third derivative components with the flux of the +// hessian over the boundary by the divergence theorem + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include + +template +Tensor<1,dim> ones () +{ + Tensor<1,dim> result; + for (unsigned int i=0; i +void test (const Triangulation &tr, + const FiniteElement &fe, + const double tolerance) +{ + DoFHandler dof(tr); + dof.distribute_dofs(fe); + + std::stringstream ss; + + deallog << "FE=" << fe.get_name() << std::endl; + + const QGauss quadrature(6); + FEValues fe_values (fe, quadrature, update_3rd_derivatives + | update_quadrature_points + | update_JxW_values); + + const QGauss face_quadrature(6); + FEFaceValues fe_face_values (fe, face_quadrature, + update_hessians + | update_quadrature_points + | update_normal_vectors + | update_JxW_values); + + for (typename DoFHandler::active_cell_iterator cell = dof.begin_active(); + cell != dof.end(); + ++cell) + { + fe_values.reinit (cell); + + deallog << "Cell nodes:" << std::endl; + for (unsigned int i=0; i::vertices_per_cell; ++i) + { + deallog << i << ": ( "; + for (unsigned int d=0; dvertex(i)[d] << " "; + deallog << ")" << std::endl; + } + + bool cell_ok = true; + + for (unsigned int c=0; c +void test_hyper_ball(const double tolerance) +{ + Triangulation tr; + GridGenerator::hyper_ball(tr); + + static const HyperBallBoundary boundary; + tr.set_boundary (0, boundary); + + tr.refine_global(1); + + FE_DGQ fe(2); + test(tr, fe, tolerance); +} + + +int main() +{ + std::ofstream logfile ("output"); + deallog << std::setprecision (3); + + deallog.attach(logfile); + deallog.depth_console (0); + deallog.threshold_double(1.e-7); + + test_hyper_ball<2>(1e-6); + test_hyper_ball<3>(1e-6); + + deallog << "done..." << std::endl; +} + diff --git a/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.output b/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.output new file mode 100644 index 0000000000..1e607eded7 --- /dev/null +++ b/tests/fe/fe_dgq_3rd_derivative_divergence_theorem.output @@ -0,0 +1,760 @@ + +DEAL::FE=FE_DGQ<2>(2) +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( 0 -1.00 ) +DEAL::2: ( -0.500 -0.500 ) +DEAL::3: ( 0 -0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 ) +DEAL::1: ( 0.707 -0.707 ) +DEAL::2: ( 0 -0.646 ) +DEAL::3: ( 0.500 -0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( 0 -0.646 ) +DEAL::2: ( -0.293 -0.293 ) +DEAL::3: ( 0 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.646 ) +DEAL::1: ( 0.500 -0.500 ) +DEAL::2: ( 0 -0.293 ) +DEAL::3: ( 0.293 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( -0.500 -0.500 ) +DEAL::2: ( -1.00 0 ) +DEAL::3: ( -0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( -0.293 -0.293 ) +DEAL::2: ( -0.646 0 ) +DEAL::3: ( -0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 ) +DEAL::1: ( -0.646 0 ) +DEAL::2: ( -0.707 0.707 ) +DEAL::3: ( -0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.646 0 ) +DEAL::1: ( -0.293 0 ) +DEAL::2: ( -0.500 0.500 ) +DEAL::3: ( -0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 -0.293 ) +DEAL::1: ( 0 -0.293 ) +DEAL::2: ( -0.293 0 ) +DEAL::3: ( 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.293 ) +DEAL::1: ( 0.293 -0.293 ) +DEAL::2: ( 0 0 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 0 ) +DEAL::1: ( 0 0 ) +DEAL::2: ( -0.293 0.293 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 ) +DEAL::1: ( 0.293 0 ) +DEAL::2: ( 0 0.293 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 ) +DEAL::1: ( 1.00 0 ) +DEAL::2: ( 0.500 -0.500 ) +DEAL::3: ( 0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 ) +DEAL::1: ( 0.707 0.707 ) +DEAL::2: ( 0.646 0 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.500 -0.500 ) +DEAL::1: ( 0.646 0 ) +DEAL::2: ( 0.293 -0.293 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.646 0 ) +DEAL::1: ( 0.500 0.500 ) +DEAL::2: ( 0.293 0 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 ) +DEAL::1: ( -0.500 0.500 ) +DEAL::2: ( 0 1.00 ) +DEAL::3: ( 0 0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 0.500 ) +DEAL::1: ( -0.293 0.293 ) +DEAL::2: ( 0 0.646 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 ) +DEAL::1: ( 0 0.646 ) +DEAL::2: ( 0.707 0.707 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.646 ) +DEAL::1: ( 0 0.293 ) +DEAL::2: ( 0.500 0.500 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::FE=FE_DGQ<3>(2) +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 -0.211 ) +DEAL::1: ( 0 -0.211 -0.211 ) +DEAL::2: ( -0.211 0 -0.211 ) +DEAL::3: ( 0 0 -0.211 ) +DEAL::4: ( -0.211 -0.211 0 ) +DEAL::5: ( 0 -0.211 0 ) +DEAL::6: ( -0.211 0 0 ) +DEAL::7: ( 0 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 -0.211 ) +DEAL::1: ( 0.211 -0.211 -0.211 ) +DEAL::2: ( 0 0 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0 -0.211 0 ) +DEAL::5: ( 0.211 -0.211 0 ) +DEAL::6: ( 0 0 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 -0.211 ) +DEAL::1: ( 0 0 -0.211 ) +DEAL::2: ( -0.211 0.211 -0.211 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.211 0 0 ) +DEAL::5: ( 0 0 0 ) +DEAL::6: ( -0.211 0.211 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.211 ) +DEAL::1: ( 0.211 0 -0.211 ) +DEAL::2: ( 0 0.211 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0 0 ) +DEAL::5: ( 0.211 0 0 ) +DEAL::6: ( 0 0.211 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 0 ) +DEAL::1: ( 0 -0.211 0 ) +DEAL::2: ( -0.211 0 0 ) +DEAL::3: ( 0 0 0 ) +DEAL::4: ( -0.211 -0.211 0.211 ) +DEAL::5: ( 0 -0.211 0.211 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 0 ) +DEAL::1: ( 0.211 -0.211 0 ) +DEAL::2: ( 0 0 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0 -0.211 0.211 ) +DEAL::5: ( 0.211 -0.211 0.211 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 0 ) +DEAL::1: ( 0 0 0 ) +DEAL::2: ( -0.211 0.211 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.211 0 0.211 ) +DEAL::5: ( 0 0 0.211 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0 ) +DEAL::1: ( 0.211 0 0 ) +DEAL::2: ( 0 0.211 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0 0.211 ) +DEAL::5: ( 0.211 0 0.211 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( 0 0 -1.00 ) +DEAL::4: ( -0.394 -0.394 -0.394 ) +DEAL::5: ( 0 -0.419 -0.419 ) +DEAL::6: ( -0.419 0 -0.419 ) +DEAL::7: ( 0 0 -0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 0 -1.00 ) +DEAL::3: ( 0.707 0 -0.707 ) +DEAL::4: ( 0 -0.419 -0.419 ) +DEAL::5: ( 0.394 -0.394 -0.394 ) +DEAL::6: ( 0 0 -0.457 ) +DEAL::7: ( 0.419 0 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( 0 0 -1.00 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( 0 0.707 -0.707 ) +DEAL::4: ( -0.419 0 -0.419 ) +DEAL::5: ( 0 0 -0.457 ) +DEAL::6: ( -0.394 0.394 -0.394 ) +DEAL::7: ( 0 0.419 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -1.00 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0.577 0.577 -0.577 ) +DEAL::4: ( 0 0 -0.457 ) +DEAL::5: ( 0.419 0 -0.419 ) +DEAL::6: ( 0 0.419 -0.419 ) +DEAL::7: ( 0.394 0.394 -0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( 0 0 -0.457 ) +DEAL::4: ( -0.211 -0.211 -0.211 ) +DEAL::5: ( 0 -0.211 -0.211 ) +DEAL::6: ( -0.211 0 -0.211 ) +DEAL::7: ( 0 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 0 -0.457 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0 -0.211 -0.211 ) +DEAL::5: ( 0.211 -0.211 -0.211 ) +DEAL::6: ( 0 0 -0.211 ) +DEAL::7: ( 0.211 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( 0 0 -0.457 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.211 0 -0.211 ) +DEAL::5: ( 0 0 -0.211 ) +DEAL::6: ( -0.211 0.211 -0.211 ) +DEAL::7: ( 0 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.457 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 0 -0.211 ) +DEAL::5: ( 0.211 0 -0.211 ) +DEAL::6: ( 0 0.211 -0.211 ) +DEAL::7: ( 0.211 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.577 -0.577 -0.577 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0.394 -0.394 -0.394 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0.707 -0.707 0 ) +DEAL::5: ( 1.00 0 0 ) +DEAL::6: ( 0.419 -0.419 0 ) +DEAL::7: ( 0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 0 -0.707 ) +DEAL::1: ( 0.577 0.577 -0.577 ) +DEAL::2: ( 0.419 0 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 1.00 0 0 ) +DEAL::5: ( 0.707 0.707 0 ) +DEAL::6: ( 0.457 0 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.394 -0.394 -0.394 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0.211 -0.211 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0.419 -0.419 0 ) +DEAL::5: ( 0.457 0 0 ) +DEAL::6: ( 0.211 -0.211 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 0 -0.419 ) +DEAL::1: ( 0.394 0.394 -0.394 ) +DEAL::2: ( 0.211 0 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0.457 0 0 ) +DEAL::5: ( 0.419 0.419 0 ) +DEAL::6: ( 0.211 0 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 0 ) +DEAL::1: ( 1.00 0 0 ) +DEAL::2: ( 0.419 -0.419 0 ) +DEAL::3: ( 0.457 0 0 ) +DEAL::4: ( 0.577 -0.577 0.577 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0.394 -0.394 0.394 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 0 ) +DEAL::1: ( 0.707 0.707 0 ) +DEAL::2: ( 0.457 0 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0.707 0 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0.419 0 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 -0.419 0 ) +DEAL::1: ( 0.457 0 0 ) +DEAL::2: ( 0.211 -0.211 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0.394 -0.394 0.394 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0.211 -0.211 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.457 0 0 ) +DEAL::1: ( 0.419 0.419 0 ) +DEAL::2: ( 0.211 0 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0.419 0 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0.211 0 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 0.577 ) +DEAL::1: ( 0 -0.707 0.707 ) +DEAL::2: ( -0.394 -0.394 0.394 ) +DEAL::3: ( 0 -0.419 0.419 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( 0 0 1.00 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( 0 0 0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 0.707 ) +DEAL::1: ( 0.577 -0.577 0.577 ) +DEAL::2: ( 0 -0.419 0.419 ) +DEAL::3: ( 0.394 -0.394 0.394 ) +DEAL::4: ( 0 0 1.00 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0 0 0.457 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 0.394 ) +DEAL::1: ( 0 -0.419 0.419 ) +DEAL::2: ( -0.211 -0.211 0.211 ) +DEAL::3: ( 0 -0.211 0.211 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( 0 0 0.457 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 0.419 ) +DEAL::1: ( 0.394 -0.394 0.394 ) +DEAL::2: ( 0 -0.211 0.211 ) +DEAL::3: ( 0.211 -0.211 0.211 ) +DEAL::4: ( 0 0 0.457 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 0.707 ) +DEAL::1: ( 0 0 1.00 ) +DEAL::2: ( -0.419 0 0.419 ) +DEAL::3: ( 0 0 0.457 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( 0 0.707 0.707 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 1.00 ) +DEAL::1: ( 0.707 0 0.707 ) +DEAL::2: ( 0 0 0.457 ) +DEAL::3: ( 0.419 0 0.419 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 0.419 ) +DEAL::1: ( 0 0 0.457 ) +DEAL::2: ( -0.211 0 0.211 ) +DEAL::3: ( 0 0 0.211 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0.457 ) +DEAL::1: ( 0.419 0 0.419 ) +DEAL::2: ( 0 0 0.211 ) +DEAL::3: ( 0.211 0 0.211 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( -0.394 -0.394 -0.394 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( -0.419 0 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( -0.419 -0.419 0 ) +DEAL::6: ( -1.00 0 0 ) +DEAL::7: ( -0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( -0.211 -0.211 -0.211 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( -0.211 0 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( -0.211 -0.211 0 ) +DEAL::6: ( -0.457 0 0 ) +DEAL::7: ( -0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( -0.419 0 -0.419 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( -0.394 0.394 -0.394 ) +DEAL::4: ( -1.00 0 0 ) +DEAL::5: ( -0.457 0 0 ) +DEAL::6: ( -0.707 0.707 0 ) +DEAL::7: ( -0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( -0.211 0 -0.211 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( -0.211 0.211 -0.211 ) +DEAL::4: ( -0.457 0 0 ) +DEAL::5: ( -0.211 0 0 ) +DEAL::6: ( -0.419 0.419 0 ) +DEAL::7: ( -0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( -0.419 -0.419 0 ) +DEAL::2: ( -1.00 0 0 ) +DEAL::3: ( -0.457 0 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( -0.394 -0.394 0.394 ) +DEAL::6: ( -0.707 0 0.707 ) +DEAL::7: ( -0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( -0.211 -0.211 0 ) +DEAL::2: ( -0.457 0 0 ) +DEAL::3: ( -0.211 0 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( -0.211 -0.211 0.211 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( -0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 0 ) +DEAL::1: ( -0.457 0 0 ) +DEAL::2: ( -0.707 0.707 0 ) +DEAL::3: ( -0.419 0.419 0 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( -0.419 0 0.419 ) +DEAL::6: ( -0.577 0.577 0.577 ) +DEAL::7: ( -0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.457 0 0 ) +DEAL::1: ( -0.211 0 0 ) +DEAL::2: ( -0.419 0.419 0 ) +DEAL::3: ( -0.211 0.211 0 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( -0.211 0 0.211 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( -0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.394 -0.394 -0.394 ) +DEAL::3: ( 0 -0.419 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( 0 -1.00 0 ) +DEAL::6: ( -0.419 -0.419 0 ) +DEAL::7: ( 0 -0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 -0.419 -0.419 ) +DEAL::3: ( 0.394 -0.394 -0.394 ) +DEAL::4: ( 0 -1.00 0 ) +DEAL::5: ( 0.707 -0.707 0 ) +DEAL::6: ( 0 -0.457 0 ) +DEAL::7: ( 0.419 -0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.211 -0.211 -0.211 ) +DEAL::3: ( 0 -0.211 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( 0 -0.457 0 ) +DEAL::6: ( -0.211 -0.211 0 ) +DEAL::7: ( 0 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 -0.211 -0.211 ) +DEAL::3: ( 0.211 -0.211 -0.211 ) +DEAL::4: ( 0 -0.457 0 ) +DEAL::5: ( 0.419 -0.419 0 ) +DEAL::6: ( 0 -0.211 0 ) +DEAL::7: ( 0.211 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( 0 -1.00 0 ) +DEAL::2: ( -0.419 -0.419 0 ) +DEAL::3: ( 0 -0.457 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( 0 -0.707 0.707 ) +DEAL::6: ( -0.394 -0.394 0.394 ) +DEAL::7: ( 0 -0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 0 ) +DEAL::1: ( 0.707 -0.707 0 ) +DEAL::2: ( 0 -0.457 0 ) +DEAL::3: ( 0.419 -0.419 0 ) +DEAL::4: ( 0 -0.707 0.707 ) +DEAL::5: ( 0.577 -0.577 0.577 ) +DEAL::6: ( 0 -0.419 0.419 ) +DEAL::7: ( 0.394 -0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( 0 -0.457 0 ) +DEAL::2: ( -0.211 -0.211 0 ) +DEAL::3: ( 0 -0.211 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( 0 -0.419 0.419 ) +DEAL::6: ( -0.211 -0.211 0.211 ) +DEAL::7: ( 0 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.457 0 ) +DEAL::1: ( 0.419 -0.419 0 ) +DEAL::2: ( 0 -0.211 0 ) +DEAL::3: ( 0.211 -0.211 0 ) +DEAL::4: ( 0 -0.419 0.419 ) +DEAL::5: ( 0.394 -0.394 0.394 ) +DEAL::6: ( 0 -0.211 0.211 ) +DEAL::7: ( 0.211 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 0.577 -0.577 ) +DEAL::1: ( -0.394 0.394 -0.394 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.707 0.707 0 ) +DEAL::5: ( -0.419 0.419 0 ) +DEAL::6: ( 0 1.00 0 ) +DEAL::7: ( 0 0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 0.394 -0.394 ) +DEAL::1: ( -0.211 0.211 -0.211 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.419 0.419 0 ) +DEAL::5: ( -0.211 0.211 0 ) +DEAL::6: ( 0 0.457 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.707 -0.707 ) +DEAL::1: ( 0 0.419 -0.419 ) +DEAL::2: ( 0.577 0.577 -0.577 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 1.00 0 ) +DEAL::5: ( 0 0.457 0 ) +DEAL::6: ( 0.707 0.707 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.419 -0.419 ) +DEAL::1: ( 0 0.211 -0.211 ) +DEAL::2: ( 0.394 0.394 -0.394 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0.457 0 ) +DEAL::5: ( 0 0.211 0 ) +DEAL::6: ( 0.419 0.419 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 0 ) +DEAL::1: ( -0.419 0.419 0 ) +DEAL::2: ( 0 1.00 0 ) +DEAL::3: ( 0 0.457 0 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( -0.394 0.394 0.394 ) +DEAL::6: ( 0 0.707 0.707 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0.419 0 ) +DEAL::1: ( -0.211 0.211 0 ) +DEAL::2: ( 0 0.457 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( -0.211 0.211 0.211 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 0 ) +DEAL::1: ( 0 0.457 0 ) +DEAL::2: ( 0.707 0.707 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( 0.577 0.577 0.577 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.457 0 ) +DEAL::1: ( 0 0.211 0 ) +DEAL::2: ( 0.419 0.419 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0 0.211 0.211 ) +DEAL::6: ( 0.394 0.394 0.394 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::done... diff --git a/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc new file mode 100644 index 0000000000..2ca982f372 --- /dev/null +++ b/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.cc @@ -0,0 +1,177 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2003 - 2015 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 the correctness of fe_values.shape_3rd_derivative for FE_Q_Hierarchical by comparing +// the integral of all shape third derivative components with the flux of the +// hessian over the boundary by the divergence theorem + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include + +template +Tensor<1,dim> ones () +{ + Tensor<1,dim> result; + for (unsigned int i=0; i +void test (const Triangulation &tr, + const FiniteElement &fe, + const double tolerance) +{ + DoFHandler dof(tr); + dof.distribute_dofs(fe); + + std::stringstream ss; + + deallog << "FE=" << fe.get_name() << std::endl; + + const QGauss quadrature(6); + FEValues fe_values (fe, quadrature, update_3rd_derivatives + | update_quadrature_points + | update_JxW_values); + + const QGauss face_quadrature(6); + FEFaceValues fe_face_values (fe, face_quadrature, + update_hessians + | update_quadrature_points + | update_normal_vectors + | update_JxW_values); + + for (typename DoFHandler::active_cell_iterator cell = dof.begin_active(); + cell != dof.end(); + ++cell) + { + fe_values.reinit (cell); + + deallog << "Cell nodes:" << std::endl; + for (unsigned int i=0; i::vertices_per_cell; ++i) + { + deallog << i << ": ( "; + for (unsigned int d=0; dvertex(i)[d] << " "; + deallog << ")" << std::endl; + } + + bool cell_ok = true; + + for (unsigned int c=0; c +void test_hyper_ball(const double tolerance) +{ + Triangulation tr; + GridGenerator::hyper_ball(tr); + + static const HyperBallBoundary boundary; + tr.set_boundary (0, boundary); + + tr.refine_global(1); + + FE_Q_Hierarchical fe(2); + test(tr, fe, tolerance); +} + + +int main() +{ + std::ofstream logfile ("output"); + deallog << std::setprecision (3); + + deallog.attach(logfile); + deallog.depth_console (0); + deallog.threshold_double(1.e-7); + + test_hyper_ball<2>(1e-6); + test_hyper_ball<3>(1e-6); + + deallog << "done..." << std::endl; +} + diff --git a/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.output b/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.output new file mode 100644 index 0000000000..f787888d12 --- /dev/null +++ b/tests/fe/fe_q_hierarchical_3rd_derivative_divergence_theorem.output @@ -0,0 +1,760 @@ + +DEAL::FE=FE_Q_Hierarchical<2>(2) +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( 0 -1.00 ) +DEAL::2: ( -0.500 -0.500 ) +DEAL::3: ( 0 -0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 ) +DEAL::1: ( 0.707 -0.707 ) +DEAL::2: ( 0 -0.646 ) +DEAL::3: ( 0.500 -0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( 0 -0.646 ) +DEAL::2: ( -0.293 -0.293 ) +DEAL::3: ( 0 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.646 ) +DEAL::1: ( 0.500 -0.500 ) +DEAL::2: ( 0 -0.293 ) +DEAL::3: ( 0.293 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( -0.500 -0.500 ) +DEAL::2: ( -1.00 0 ) +DEAL::3: ( -0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( -0.293 -0.293 ) +DEAL::2: ( -0.646 0 ) +DEAL::3: ( -0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 ) +DEAL::1: ( -0.646 0 ) +DEAL::2: ( -0.707 0.707 ) +DEAL::3: ( -0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.646 0 ) +DEAL::1: ( -0.293 0 ) +DEAL::2: ( -0.500 0.500 ) +DEAL::3: ( -0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 -0.293 ) +DEAL::1: ( 0 -0.293 ) +DEAL::2: ( -0.293 0 ) +DEAL::3: ( 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.293 ) +DEAL::1: ( 0.293 -0.293 ) +DEAL::2: ( 0 0 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 0 ) +DEAL::1: ( 0 0 ) +DEAL::2: ( -0.293 0.293 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 ) +DEAL::1: ( 0.293 0 ) +DEAL::2: ( 0 0.293 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 ) +DEAL::1: ( 1.00 0 ) +DEAL::2: ( 0.500 -0.500 ) +DEAL::3: ( 0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 ) +DEAL::1: ( 0.707 0.707 ) +DEAL::2: ( 0.646 0 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.500 -0.500 ) +DEAL::1: ( 0.646 0 ) +DEAL::2: ( 0.293 -0.293 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.646 0 ) +DEAL::1: ( 0.500 0.500 ) +DEAL::2: ( 0.293 0 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 ) +DEAL::1: ( -0.500 0.500 ) +DEAL::2: ( 0 1.00 ) +DEAL::3: ( 0 0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 0.500 ) +DEAL::1: ( -0.293 0.293 ) +DEAL::2: ( 0 0.646 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 ) +DEAL::1: ( 0 0.646 ) +DEAL::2: ( 0.707 0.707 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.646 ) +DEAL::1: ( 0 0.293 ) +DEAL::2: ( 0.500 0.500 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::FE=FE_Q_Hierarchical<3>(2) +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 -0.211 ) +DEAL::1: ( 0 -0.211 -0.211 ) +DEAL::2: ( -0.211 0 -0.211 ) +DEAL::3: ( 0 0 -0.211 ) +DEAL::4: ( -0.211 -0.211 0 ) +DEAL::5: ( 0 -0.211 0 ) +DEAL::6: ( -0.211 0 0 ) +DEAL::7: ( 0 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 -0.211 ) +DEAL::1: ( 0.211 -0.211 -0.211 ) +DEAL::2: ( 0 0 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0 -0.211 0 ) +DEAL::5: ( 0.211 -0.211 0 ) +DEAL::6: ( 0 0 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 -0.211 ) +DEAL::1: ( 0 0 -0.211 ) +DEAL::2: ( -0.211 0.211 -0.211 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.211 0 0 ) +DEAL::5: ( 0 0 0 ) +DEAL::6: ( -0.211 0.211 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.211 ) +DEAL::1: ( 0.211 0 -0.211 ) +DEAL::2: ( 0 0.211 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0 0 ) +DEAL::5: ( 0.211 0 0 ) +DEAL::6: ( 0 0.211 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 0 ) +DEAL::1: ( 0 -0.211 0 ) +DEAL::2: ( -0.211 0 0 ) +DEAL::3: ( 0 0 0 ) +DEAL::4: ( -0.211 -0.211 0.211 ) +DEAL::5: ( 0 -0.211 0.211 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 0 ) +DEAL::1: ( 0.211 -0.211 0 ) +DEAL::2: ( 0 0 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0 -0.211 0.211 ) +DEAL::5: ( 0.211 -0.211 0.211 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 0 ) +DEAL::1: ( 0 0 0 ) +DEAL::2: ( -0.211 0.211 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.211 0 0.211 ) +DEAL::5: ( 0 0 0.211 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0 ) +DEAL::1: ( 0.211 0 0 ) +DEAL::2: ( 0 0.211 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0 0.211 ) +DEAL::5: ( 0.211 0 0.211 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( 0 0 -1.00 ) +DEAL::4: ( -0.394 -0.394 -0.394 ) +DEAL::5: ( 0 -0.419 -0.419 ) +DEAL::6: ( -0.419 0 -0.419 ) +DEAL::7: ( 0 0 -0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 0 -1.00 ) +DEAL::3: ( 0.707 0 -0.707 ) +DEAL::4: ( 0 -0.419 -0.419 ) +DEAL::5: ( 0.394 -0.394 -0.394 ) +DEAL::6: ( 0 0 -0.457 ) +DEAL::7: ( 0.419 0 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( 0 0 -1.00 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( 0 0.707 -0.707 ) +DEAL::4: ( -0.419 0 -0.419 ) +DEAL::5: ( 0 0 -0.457 ) +DEAL::6: ( -0.394 0.394 -0.394 ) +DEAL::7: ( 0 0.419 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -1.00 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0.577 0.577 -0.577 ) +DEAL::4: ( 0 0 -0.457 ) +DEAL::5: ( 0.419 0 -0.419 ) +DEAL::6: ( 0 0.419 -0.419 ) +DEAL::7: ( 0.394 0.394 -0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( 0 0 -0.457 ) +DEAL::4: ( -0.211 -0.211 -0.211 ) +DEAL::5: ( 0 -0.211 -0.211 ) +DEAL::6: ( -0.211 0 -0.211 ) +DEAL::7: ( 0 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 0 -0.457 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0 -0.211 -0.211 ) +DEAL::5: ( 0.211 -0.211 -0.211 ) +DEAL::6: ( 0 0 -0.211 ) +DEAL::7: ( 0.211 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( 0 0 -0.457 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.211 0 -0.211 ) +DEAL::5: ( 0 0 -0.211 ) +DEAL::6: ( -0.211 0.211 -0.211 ) +DEAL::7: ( 0 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.457 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 0 -0.211 ) +DEAL::5: ( 0.211 0 -0.211 ) +DEAL::6: ( 0 0.211 -0.211 ) +DEAL::7: ( 0.211 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.577 -0.577 -0.577 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0.394 -0.394 -0.394 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0.707 -0.707 0 ) +DEAL::5: ( 1.00 0 0 ) +DEAL::6: ( 0.419 -0.419 0 ) +DEAL::7: ( 0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 0 -0.707 ) +DEAL::1: ( 0.577 0.577 -0.577 ) +DEAL::2: ( 0.419 0 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 1.00 0 0 ) +DEAL::5: ( 0.707 0.707 0 ) +DEAL::6: ( 0.457 0 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.394 -0.394 -0.394 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0.211 -0.211 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0.419 -0.419 0 ) +DEAL::5: ( 0.457 0 0 ) +DEAL::6: ( 0.211 -0.211 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 0 -0.419 ) +DEAL::1: ( 0.394 0.394 -0.394 ) +DEAL::2: ( 0.211 0 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0.457 0 0 ) +DEAL::5: ( 0.419 0.419 0 ) +DEAL::6: ( 0.211 0 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 0 ) +DEAL::1: ( 1.00 0 0 ) +DEAL::2: ( 0.419 -0.419 0 ) +DEAL::3: ( 0.457 0 0 ) +DEAL::4: ( 0.577 -0.577 0.577 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0.394 -0.394 0.394 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 0 ) +DEAL::1: ( 0.707 0.707 0 ) +DEAL::2: ( 0.457 0 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0.707 0 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0.419 0 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 -0.419 0 ) +DEAL::1: ( 0.457 0 0 ) +DEAL::2: ( 0.211 -0.211 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0.394 -0.394 0.394 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0.211 -0.211 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.457 0 0 ) +DEAL::1: ( 0.419 0.419 0 ) +DEAL::2: ( 0.211 0 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0.419 0 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0.211 0 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 0.577 ) +DEAL::1: ( 0 -0.707 0.707 ) +DEAL::2: ( -0.394 -0.394 0.394 ) +DEAL::3: ( 0 -0.419 0.419 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( 0 0 1.00 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( 0 0 0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 0.707 ) +DEAL::1: ( 0.577 -0.577 0.577 ) +DEAL::2: ( 0 -0.419 0.419 ) +DEAL::3: ( 0.394 -0.394 0.394 ) +DEAL::4: ( 0 0 1.00 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0 0 0.457 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 0.394 ) +DEAL::1: ( 0 -0.419 0.419 ) +DEAL::2: ( -0.211 -0.211 0.211 ) +DEAL::3: ( 0 -0.211 0.211 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( 0 0 0.457 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 0.419 ) +DEAL::1: ( 0.394 -0.394 0.394 ) +DEAL::2: ( 0 -0.211 0.211 ) +DEAL::3: ( 0.211 -0.211 0.211 ) +DEAL::4: ( 0 0 0.457 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 0.707 ) +DEAL::1: ( 0 0 1.00 ) +DEAL::2: ( -0.419 0 0.419 ) +DEAL::3: ( 0 0 0.457 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( 0 0.707 0.707 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 1.00 ) +DEAL::1: ( 0.707 0 0.707 ) +DEAL::2: ( 0 0 0.457 ) +DEAL::3: ( 0.419 0 0.419 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 0.419 ) +DEAL::1: ( 0 0 0.457 ) +DEAL::2: ( -0.211 0 0.211 ) +DEAL::3: ( 0 0 0.211 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0.457 ) +DEAL::1: ( 0.419 0 0.419 ) +DEAL::2: ( 0 0 0.211 ) +DEAL::3: ( 0.211 0 0.211 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( -0.394 -0.394 -0.394 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( -0.419 0 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( -0.419 -0.419 0 ) +DEAL::6: ( -1.00 0 0 ) +DEAL::7: ( -0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( -0.211 -0.211 -0.211 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( -0.211 0 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( -0.211 -0.211 0 ) +DEAL::6: ( -0.457 0 0 ) +DEAL::7: ( -0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( -0.419 0 -0.419 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( -0.394 0.394 -0.394 ) +DEAL::4: ( -1.00 0 0 ) +DEAL::5: ( -0.457 0 0 ) +DEAL::6: ( -0.707 0.707 0 ) +DEAL::7: ( -0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( -0.211 0 -0.211 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( -0.211 0.211 -0.211 ) +DEAL::4: ( -0.457 0 0 ) +DEAL::5: ( -0.211 0 0 ) +DEAL::6: ( -0.419 0.419 0 ) +DEAL::7: ( -0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( -0.419 -0.419 0 ) +DEAL::2: ( -1.00 0 0 ) +DEAL::3: ( -0.457 0 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( -0.394 -0.394 0.394 ) +DEAL::6: ( -0.707 0 0.707 ) +DEAL::7: ( -0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( -0.211 -0.211 0 ) +DEAL::2: ( -0.457 0 0 ) +DEAL::3: ( -0.211 0 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( -0.211 -0.211 0.211 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( -0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 0 ) +DEAL::1: ( -0.457 0 0 ) +DEAL::2: ( -0.707 0.707 0 ) +DEAL::3: ( -0.419 0.419 0 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( -0.419 0 0.419 ) +DEAL::6: ( -0.577 0.577 0.577 ) +DEAL::7: ( -0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.457 0 0 ) +DEAL::1: ( -0.211 0 0 ) +DEAL::2: ( -0.419 0.419 0 ) +DEAL::3: ( -0.211 0.211 0 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( -0.211 0 0.211 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( -0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.394 -0.394 -0.394 ) +DEAL::3: ( 0 -0.419 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( 0 -1.00 0 ) +DEAL::6: ( -0.419 -0.419 0 ) +DEAL::7: ( 0 -0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 -0.419 -0.419 ) +DEAL::3: ( 0.394 -0.394 -0.394 ) +DEAL::4: ( 0 -1.00 0 ) +DEAL::5: ( 0.707 -0.707 0 ) +DEAL::6: ( 0 -0.457 0 ) +DEAL::7: ( 0.419 -0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.211 -0.211 -0.211 ) +DEAL::3: ( 0 -0.211 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( 0 -0.457 0 ) +DEAL::6: ( -0.211 -0.211 0 ) +DEAL::7: ( 0 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 -0.211 -0.211 ) +DEAL::3: ( 0.211 -0.211 -0.211 ) +DEAL::4: ( 0 -0.457 0 ) +DEAL::5: ( 0.419 -0.419 0 ) +DEAL::6: ( 0 -0.211 0 ) +DEAL::7: ( 0.211 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( 0 -1.00 0 ) +DEAL::2: ( -0.419 -0.419 0 ) +DEAL::3: ( 0 -0.457 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( 0 -0.707 0.707 ) +DEAL::6: ( -0.394 -0.394 0.394 ) +DEAL::7: ( 0 -0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 0 ) +DEAL::1: ( 0.707 -0.707 0 ) +DEAL::2: ( 0 -0.457 0 ) +DEAL::3: ( 0.419 -0.419 0 ) +DEAL::4: ( 0 -0.707 0.707 ) +DEAL::5: ( 0.577 -0.577 0.577 ) +DEAL::6: ( 0 -0.419 0.419 ) +DEAL::7: ( 0.394 -0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( 0 -0.457 0 ) +DEAL::2: ( -0.211 -0.211 0 ) +DEAL::3: ( 0 -0.211 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( 0 -0.419 0.419 ) +DEAL::6: ( -0.211 -0.211 0.211 ) +DEAL::7: ( 0 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.457 0 ) +DEAL::1: ( 0.419 -0.419 0 ) +DEAL::2: ( 0 -0.211 0 ) +DEAL::3: ( 0.211 -0.211 0 ) +DEAL::4: ( 0 -0.419 0.419 ) +DEAL::5: ( 0.394 -0.394 0.394 ) +DEAL::6: ( 0 -0.211 0.211 ) +DEAL::7: ( 0.211 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 0.577 -0.577 ) +DEAL::1: ( -0.394 0.394 -0.394 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.707 0.707 0 ) +DEAL::5: ( -0.419 0.419 0 ) +DEAL::6: ( 0 1.00 0 ) +DEAL::7: ( 0 0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 0.394 -0.394 ) +DEAL::1: ( -0.211 0.211 -0.211 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.419 0.419 0 ) +DEAL::5: ( -0.211 0.211 0 ) +DEAL::6: ( 0 0.457 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.707 -0.707 ) +DEAL::1: ( 0 0.419 -0.419 ) +DEAL::2: ( 0.577 0.577 -0.577 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 1.00 0 ) +DEAL::5: ( 0 0.457 0 ) +DEAL::6: ( 0.707 0.707 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.419 -0.419 ) +DEAL::1: ( 0 0.211 -0.211 ) +DEAL::2: ( 0.394 0.394 -0.394 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0.457 0 ) +DEAL::5: ( 0 0.211 0 ) +DEAL::6: ( 0.419 0.419 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 0 ) +DEAL::1: ( -0.419 0.419 0 ) +DEAL::2: ( 0 1.00 0 ) +DEAL::3: ( 0 0.457 0 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( -0.394 0.394 0.394 ) +DEAL::6: ( 0 0.707 0.707 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0.419 0 ) +DEAL::1: ( -0.211 0.211 0 ) +DEAL::2: ( 0 0.457 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( -0.211 0.211 0.211 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 0 ) +DEAL::1: ( 0 0.457 0 ) +DEAL::2: ( 0.707 0.707 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( 0.577 0.577 0.577 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.457 0 ) +DEAL::1: ( 0 0.211 0 ) +DEAL::2: ( 0.419 0.419 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0 0.211 0.211 ) +DEAL::6: ( 0.394 0.394 0.394 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::done... diff --git a/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc b/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc new file mode 100644 index 0000000000..22a6a02d29 --- /dev/null +++ b/tests/fe/fe_system_3rd_derivative_divergence_theorem.cc @@ -0,0 +1,176 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2003 - 2015 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 the correctness of fe_values.shape_3rd_derivative for FE_System by comparing +// the integral of all shape third derivative components with the flux of the +// hessian over the boundary by the divergence theorem + +#include "../tests.h" +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include +#include + +#include +#include + +template +Tensor<1,dim> ones () +{ + Tensor<1,dim> result; + for (unsigned int i=0; i +void test (const Triangulation &tr, + const FiniteElement &fe, + const double tolerance) +{ + DoFHandler dof(tr); + dof.distribute_dofs(fe); + + std::stringstream ss; + + deallog << "FE=" << fe.get_name() << std::endl; + + const QGauss quadrature(6); + FEValues fe_values (fe, quadrature, update_3rd_derivatives + | update_quadrature_points + | update_JxW_values); + + const QGauss face_quadrature(6); + FEFaceValues fe_face_values (fe, face_quadrature, + update_hessians + | update_quadrature_points + | update_normal_vectors + | update_JxW_values); + + for (typename DoFHandler::active_cell_iterator cell = dof.begin_active(); + cell != dof.end(); + ++cell) + { + fe_values.reinit (cell); + + deallog << "Cell nodes:" << std::endl; + for (unsigned int i=0; i::vertices_per_cell; ++i) + { + deallog << i << ": ( "; + for (unsigned int d=0; dvertex(i)[d] << " "; + deallog << ")" << std::endl; + } + + bool cell_ok = true; + + for (unsigned int c=0; c +void test_hyper_ball(const double tolerance) +{ + Triangulation tr; + GridGenerator::hyper_ball(tr); + + static const HyperBallBoundary boundary; + tr.set_boundary (0, boundary); + + tr.refine_global(1); + + FESystem fe(FE_Q(2),1); + test(tr, fe, tolerance); +} + + +int main() +{ + std::ofstream logfile ("output"); + deallog << std::setprecision (3); + + deallog.attach(logfile); + deallog.depth_console (0); + deallog.threshold_double(1.e-7); + + test_hyper_ball<2>(1e-6); + test_hyper_ball<3>(1e-6); + + deallog << "done..." << std::endl; +} + diff --git a/tests/fe/fe_system_3rd_derivative_divergence_theorem.output b/tests/fe/fe_system_3rd_derivative_divergence_theorem.output new file mode 100644 index 0000000000..0b728545b5 --- /dev/null +++ b/tests/fe/fe_system_3rd_derivative_divergence_theorem.output @@ -0,0 +1,760 @@ + +DEAL::FE=FESystem<2>[FE_Q<2>(2)] +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( 0 -1.00 ) +DEAL::2: ( -0.500 -0.500 ) +DEAL::3: ( 0 -0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 ) +DEAL::1: ( 0.707 -0.707 ) +DEAL::2: ( 0 -0.646 ) +DEAL::3: ( 0.500 -0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( 0 -0.646 ) +DEAL::2: ( -0.293 -0.293 ) +DEAL::3: ( 0 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.646 ) +DEAL::1: ( 0.500 -0.500 ) +DEAL::2: ( 0 -0.293 ) +DEAL::3: ( 0.293 -0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 ) +DEAL::1: ( -0.500 -0.500 ) +DEAL::2: ( -1.00 0 ) +DEAL::3: ( -0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 -0.500 ) +DEAL::1: ( -0.293 -0.293 ) +DEAL::2: ( -0.646 0 ) +DEAL::3: ( -0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 ) +DEAL::1: ( -0.646 0 ) +DEAL::2: ( -0.707 0.707 ) +DEAL::3: ( -0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.646 0 ) +DEAL::1: ( -0.293 0 ) +DEAL::2: ( -0.500 0.500 ) +DEAL::3: ( -0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 -0.293 ) +DEAL::1: ( 0 -0.293 ) +DEAL::2: ( -0.293 0 ) +DEAL::3: ( 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.293 ) +DEAL::1: ( 0.293 -0.293 ) +DEAL::2: ( 0 0 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.293 0 ) +DEAL::1: ( 0 0 ) +DEAL::2: ( -0.293 0.293 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 ) +DEAL::1: ( 0.293 0 ) +DEAL::2: ( 0 0.293 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 ) +DEAL::1: ( 1.00 0 ) +DEAL::2: ( 0.500 -0.500 ) +DEAL::3: ( 0.646 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 ) +DEAL::1: ( 0.707 0.707 ) +DEAL::2: ( 0.646 0 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.500 -0.500 ) +DEAL::1: ( 0.646 0 ) +DEAL::2: ( 0.293 -0.293 ) +DEAL::3: ( 0.293 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.646 0 ) +DEAL::1: ( 0.500 0.500 ) +DEAL::2: ( 0.293 0 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 ) +DEAL::1: ( -0.500 0.500 ) +DEAL::2: ( 0 1.00 ) +DEAL::3: ( 0 0.646 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.500 0.500 ) +DEAL::1: ( -0.293 0.293 ) +DEAL::2: ( 0 0.646 ) +DEAL::3: ( 0 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 ) +DEAL::1: ( 0 0.646 ) +DEAL::2: ( 0.707 0.707 ) +DEAL::3: ( 0.500 0.500 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.646 ) +DEAL::1: ( 0 0.293 ) +DEAL::2: ( 0.500 0.500 ) +DEAL::3: ( 0.293 0.293 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::FE=FESystem<3>[FE_Q<3>(2)] +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 -0.211 ) +DEAL::1: ( 0 -0.211 -0.211 ) +DEAL::2: ( -0.211 0 -0.211 ) +DEAL::3: ( 0 0 -0.211 ) +DEAL::4: ( -0.211 -0.211 0 ) +DEAL::5: ( 0 -0.211 0 ) +DEAL::6: ( -0.211 0 0 ) +DEAL::7: ( 0 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 -0.211 ) +DEAL::1: ( 0.211 -0.211 -0.211 ) +DEAL::2: ( 0 0 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0 -0.211 0 ) +DEAL::5: ( 0.211 -0.211 0 ) +DEAL::6: ( 0 0 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 -0.211 ) +DEAL::1: ( 0 0 -0.211 ) +DEAL::2: ( -0.211 0.211 -0.211 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.211 0 0 ) +DEAL::5: ( 0 0 0 ) +DEAL::6: ( -0.211 0.211 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.211 ) +DEAL::1: ( 0.211 0 -0.211 ) +DEAL::2: ( 0 0.211 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0 0 ) +DEAL::5: ( 0.211 0 0 ) +DEAL::6: ( 0 0.211 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 -0.211 0 ) +DEAL::1: ( 0 -0.211 0 ) +DEAL::2: ( -0.211 0 0 ) +DEAL::3: ( 0 0 0 ) +DEAL::4: ( -0.211 -0.211 0.211 ) +DEAL::5: ( 0 -0.211 0.211 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.211 0 ) +DEAL::1: ( 0.211 -0.211 0 ) +DEAL::2: ( 0 0 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0 -0.211 0.211 ) +DEAL::5: ( 0.211 -0.211 0.211 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.211 0 0 ) +DEAL::1: ( 0 0 0 ) +DEAL::2: ( -0.211 0.211 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.211 0 0.211 ) +DEAL::5: ( 0 0 0.211 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0 ) +DEAL::1: ( 0.211 0 0 ) +DEAL::2: ( 0 0.211 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0 0.211 ) +DEAL::5: ( 0.211 0 0.211 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( 0 0 -1.00 ) +DEAL::4: ( -0.394 -0.394 -0.394 ) +DEAL::5: ( 0 -0.419 -0.419 ) +DEAL::6: ( -0.419 0 -0.419 ) +DEAL::7: ( 0 0 -0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 0 -1.00 ) +DEAL::3: ( 0.707 0 -0.707 ) +DEAL::4: ( 0 -0.419 -0.419 ) +DEAL::5: ( 0.394 -0.394 -0.394 ) +DEAL::6: ( 0 0 -0.457 ) +DEAL::7: ( 0.419 0 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( 0 0 -1.00 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( 0 0.707 -0.707 ) +DEAL::4: ( -0.419 0 -0.419 ) +DEAL::5: ( 0 0 -0.457 ) +DEAL::6: ( -0.394 0.394 -0.394 ) +DEAL::7: ( 0 0.419 -0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -1.00 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0.577 0.577 -0.577 ) +DEAL::4: ( 0 0 -0.457 ) +DEAL::5: ( 0.419 0 -0.419 ) +DEAL::6: ( 0 0.419 -0.419 ) +DEAL::7: ( 0.394 0.394 -0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( 0 0 -0.457 ) +DEAL::4: ( -0.211 -0.211 -0.211 ) +DEAL::5: ( 0 -0.211 -0.211 ) +DEAL::6: ( -0.211 0 -0.211 ) +DEAL::7: ( 0 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 0 -0.457 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0 -0.211 -0.211 ) +DEAL::5: ( 0.211 -0.211 -0.211 ) +DEAL::6: ( 0 0 -0.211 ) +DEAL::7: ( 0.211 0 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( 0 0 -0.457 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.211 0 -0.211 ) +DEAL::5: ( 0 0 -0.211 ) +DEAL::6: ( -0.211 0.211 -0.211 ) +DEAL::7: ( 0 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 -0.457 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 0 -0.211 ) +DEAL::5: ( 0.211 0 -0.211 ) +DEAL::6: ( 0 0.211 -0.211 ) +DEAL::7: ( 0.211 0.211 -0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.577 -0.577 -0.577 ) +DEAL::1: ( 0.707 0 -0.707 ) +DEAL::2: ( 0.394 -0.394 -0.394 ) +DEAL::3: ( 0.419 0 -0.419 ) +DEAL::4: ( 0.707 -0.707 0 ) +DEAL::5: ( 1.00 0 0 ) +DEAL::6: ( 0.419 -0.419 0 ) +DEAL::7: ( 0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 0 -0.707 ) +DEAL::1: ( 0.577 0.577 -0.577 ) +DEAL::2: ( 0.419 0 -0.419 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 1.00 0 0 ) +DEAL::5: ( 0.707 0.707 0 ) +DEAL::6: ( 0.457 0 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.394 -0.394 -0.394 ) +DEAL::1: ( 0.419 0 -0.419 ) +DEAL::2: ( 0.211 -0.211 -0.211 ) +DEAL::3: ( 0.211 0 -0.211 ) +DEAL::4: ( 0.419 -0.419 0 ) +DEAL::5: ( 0.457 0 0 ) +DEAL::6: ( 0.211 -0.211 0 ) +DEAL::7: ( 0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 0 -0.419 ) +DEAL::1: ( 0.394 0.394 -0.394 ) +DEAL::2: ( 0.211 0 -0.211 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0.457 0 0 ) +DEAL::5: ( 0.419 0.419 0 ) +DEAL::6: ( 0.211 0 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.707 -0.707 0 ) +DEAL::1: ( 1.00 0 0 ) +DEAL::2: ( 0.419 -0.419 0 ) +DEAL::3: ( 0.457 0 0 ) +DEAL::4: ( 0.577 -0.577 0.577 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0.394 -0.394 0.394 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 1.00 0 0 ) +DEAL::1: ( 0.707 0.707 0 ) +DEAL::2: ( 0.457 0 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0.707 0 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0.419 0 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.419 -0.419 0 ) +DEAL::1: ( 0.457 0 0 ) +DEAL::2: ( 0.211 -0.211 0 ) +DEAL::3: ( 0.211 0 0 ) +DEAL::4: ( 0.394 -0.394 0.394 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0.211 -0.211 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0.457 0 0 ) +DEAL::1: ( 0.419 0.419 0 ) +DEAL::2: ( 0.211 0 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0.419 0 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0.211 0 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 0.577 ) +DEAL::1: ( 0 -0.707 0.707 ) +DEAL::2: ( -0.394 -0.394 0.394 ) +DEAL::3: ( 0 -0.419 0.419 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( 0 0 1.00 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( 0 0 0.457 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 0.707 ) +DEAL::1: ( 0.577 -0.577 0.577 ) +DEAL::2: ( 0 -0.419 0.419 ) +DEAL::3: ( 0.394 -0.394 0.394 ) +DEAL::4: ( 0 0 1.00 ) +DEAL::5: ( 0.707 0 0.707 ) +DEAL::6: ( 0 0 0.457 ) +DEAL::7: ( 0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 0.394 ) +DEAL::1: ( 0 -0.419 0.419 ) +DEAL::2: ( -0.211 -0.211 0.211 ) +DEAL::3: ( 0 -0.211 0.211 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( 0 0 0.457 ) +DEAL::6: ( -0.211 0 0.211 ) +DEAL::7: ( 0 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 0.419 ) +DEAL::1: ( 0.394 -0.394 0.394 ) +DEAL::2: ( 0 -0.211 0.211 ) +DEAL::3: ( 0.211 -0.211 0.211 ) +DEAL::4: ( 0 0 0.457 ) +DEAL::5: ( 0.419 0 0.419 ) +DEAL::6: ( 0 0 0.211 ) +DEAL::7: ( 0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 0.707 ) +DEAL::1: ( 0 0 1.00 ) +DEAL::2: ( -0.419 0 0.419 ) +DEAL::3: ( 0 0 0.457 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( 0 0.707 0.707 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 1.00 ) +DEAL::1: ( 0.707 0 0.707 ) +DEAL::2: ( 0 0 0.457 ) +DEAL::3: ( 0.419 0 0.419 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0.577 0.577 0.577 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 0.419 ) +DEAL::1: ( 0 0 0.457 ) +DEAL::2: ( -0.211 0 0.211 ) +DEAL::3: ( 0 0 0.211 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( -0.211 0.211 0.211 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0 0.457 ) +DEAL::1: ( 0.419 0 0.419 ) +DEAL::2: ( 0 0 0.211 ) +DEAL::3: ( 0.211 0 0.211 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0.394 0.394 0.394 ) +DEAL::6: ( 0 0.211 0.211 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( -0.394 -0.394 -0.394 ) +DEAL::2: ( -0.707 0 -0.707 ) +DEAL::3: ( -0.419 0 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( -0.419 -0.419 0 ) +DEAL::6: ( -1.00 0 0 ) +DEAL::7: ( -0.457 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( -0.211 -0.211 -0.211 ) +DEAL::2: ( -0.419 0 -0.419 ) +DEAL::3: ( -0.211 0 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( -0.211 -0.211 0 ) +DEAL::6: ( -0.457 0 0 ) +DEAL::7: ( -0.211 0 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0 -0.707 ) +DEAL::1: ( -0.419 0 -0.419 ) +DEAL::2: ( -0.577 0.577 -0.577 ) +DEAL::3: ( -0.394 0.394 -0.394 ) +DEAL::4: ( -1.00 0 0 ) +DEAL::5: ( -0.457 0 0 ) +DEAL::6: ( -0.707 0.707 0 ) +DEAL::7: ( -0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0 -0.419 ) +DEAL::1: ( -0.211 0 -0.211 ) +DEAL::2: ( -0.394 0.394 -0.394 ) +DEAL::3: ( -0.211 0.211 -0.211 ) +DEAL::4: ( -0.457 0 0 ) +DEAL::5: ( -0.211 0 0 ) +DEAL::6: ( -0.419 0.419 0 ) +DEAL::7: ( -0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( -0.419 -0.419 0 ) +DEAL::2: ( -1.00 0 0 ) +DEAL::3: ( -0.457 0 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( -0.394 -0.394 0.394 ) +DEAL::6: ( -0.707 0 0.707 ) +DEAL::7: ( -0.419 0 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( -0.211 -0.211 0 ) +DEAL::2: ( -0.457 0 0 ) +DEAL::3: ( -0.211 0 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( -0.211 -0.211 0.211 ) +DEAL::6: ( -0.419 0 0.419 ) +DEAL::7: ( -0.211 0 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -1.00 0 0 ) +DEAL::1: ( -0.457 0 0 ) +DEAL::2: ( -0.707 0.707 0 ) +DEAL::3: ( -0.419 0.419 0 ) +DEAL::4: ( -0.707 0 0.707 ) +DEAL::5: ( -0.419 0 0.419 ) +DEAL::6: ( -0.577 0.577 0.577 ) +DEAL::7: ( -0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.457 0 0 ) +DEAL::1: ( -0.211 0 0 ) +DEAL::2: ( -0.419 0.419 0 ) +DEAL::3: ( -0.211 0.211 0 ) +DEAL::4: ( -0.419 0 0.419 ) +DEAL::5: ( -0.211 0 0.211 ) +DEAL::6: ( -0.394 0.394 0.394 ) +DEAL::7: ( -0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 -0.577 -0.577 ) +DEAL::1: ( 0 -0.707 -0.707 ) +DEAL::2: ( -0.394 -0.394 -0.394 ) +DEAL::3: ( 0 -0.419 -0.419 ) +DEAL::4: ( -0.707 -0.707 0 ) +DEAL::5: ( 0 -1.00 0 ) +DEAL::6: ( -0.419 -0.419 0 ) +DEAL::7: ( 0 -0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.707 -0.707 ) +DEAL::1: ( 0.577 -0.577 -0.577 ) +DEAL::2: ( 0 -0.419 -0.419 ) +DEAL::3: ( 0.394 -0.394 -0.394 ) +DEAL::4: ( 0 -1.00 0 ) +DEAL::5: ( 0.707 -0.707 0 ) +DEAL::6: ( 0 -0.457 0 ) +DEAL::7: ( 0.419 -0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 -0.394 -0.394 ) +DEAL::1: ( 0 -0.419 -0.419 ) +DEAL::2: ( -0.211 -0.211 -0.211 ) +DEAL::3: ( 0 -0.211 -0.211 ) +DEAL::4: ( -0.419 -0.419 0 ) +DEAL::5: ( 0 -0.457 0 ) +DEAL::6: ( -0.211 -0.211 0 ) +DEAL::7: ( 0 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.419 -0.419 ) +DEAL::1: ( 0.394 -0.394 -0.394 ) +DEAL::2: ( 0 -0.211 -0.211 ) +DEAL::3: ( 0.211 -0.211 -0.211 ) +DEAL::4: ( 0 -0.457 0 ) +DEAL::5: ( 0.419 -0.419 0 ) +DEAL::6: ( 0 -0.211 0 ) +DEAL::7: ( 0.211 -0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 -0.707 0 ) +DEAL::1: ( 0 -1.00 0 ) +DEAL::2: ( -0.419 -0.419 0 ) +DEAL::3: ( 0 -0.457 0 ) +DEAL::4: ( -0.577 -0.577 0.577 ) +DEAL::5: ( 0 -0.707 0.707 ) +DEAL::6: ( -0.394 -0.394 0.394 ) +DEAL::7: ( 0 -0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -1.00 0 ) +DEAL::1: ( 0.707 -0.707 0 ) +DEAL::2: ( 0 -0.457 0 ) +DEAL::3: ( 0.419 -0.419 0 ) +DEAL::4: ( 0 -0.707 0.707 ) +DEAL::5: ( 0.577 -0.577 0.577 ) +DEAL::6: ( 0 -0.419 0.419 ) +DEAL::7: ( 0.394 -0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 -0.419 0 ) +DEAL::1: ( 0 -0.457 0 ) +DEAL::2: ( -0.211 -0.211 0 ) +DEAL::3: ( 0 -0.211 0 ) +DEAL::4: ( -0.394 -0.394 0.394 ) +DEAL::5: ( 0 -0.419 0.419 ) +DEAL::6: ( -0.211 -0.211 0.211 ) +DEAL::7: ( 0 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 -0.457 0 ) +DEAL::1: ( 0.419 -0.419 0 ) +DEAL::2: ( 0 -0.211 0 ) +DEAL::3: ( 0.211 -0.211 0 ) +DEAL::4: ( 0 -0.419 0.419 ) +DEAL::5: ( 0.394 -0.394 0.394 ) +DEAL::6: ( 0 -0.211 0.211 ) +DEAL::7: ( 0.211 -0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.577 0.577 -0.577 ) +DEAL::1: ( -0.394 0.394 -0.394 ) +DEAL::2: ( 0 0.707 -0.707 ) +DEAL::3: ( 0 0.419 -0.419 ) +DEAL::4: ( -0.707 0.707 0 ) +DEAL::5: ( -0.419 0.419 0 ) +DEAL::6: ( 0 1.00 0 ) +DEAL::7: ( 0 0.457 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.394 0.394 -0.394 ) +DEAL::1: ( -0.211 0.211 -0.211 ) +DEAL::2: ( 0 0.419 -0.419 ) +DEAL::3: ( 0 0.211 -0.211 ) +DEAL::4: ( -0.419 0.419 0 ) +DEAL::5: ( -0.211 0.211 0 ) +DEAL::6: ( 0 0.457 0 ) +DEAL::7: ( 0 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.707 -0.707 ) +DEAL::1: ( 0 0.419 -0.419 ) +DEAL::2: ( 0.577 0.577 -0.577 ) +DEAL::3: ( 0.394 0.394 -0.394 ) +DEAL::4: ( 0 1.00 0 ) +DEAL::5: ( 0 0.457 0 ) +DEAL::6: ( 0.707 0.707 0 ) +DEAL::7: ( 0.419 0.419 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.419 -0.419 ) +DEAL::1: ( 0 0.211 -0.211 ) +DEAL::2: ( 0.394 0.394 -0.394 ) +DEAL::3: ( 0.211 0.211 -0.211 ) +DEAL::4: ( 0 0.457 0 ) +DEAL::5: ( 0 0.211 0 ) +DEAL::6: ( 0.419 0.419 0 ) +DEAL::7: ( 0.211 0.211 0 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.707 0.707 0 ) +DEAL::1: ( -0.419 0.419 0 ) +DEAL::2: ( 0 1.00 0 ) +DEAL::3: ( 0 0.457 0 ) +DEAL::4: ( -0.577 0.577 0.577 ) +DEAL::5: ( -0.394 0.394 0.394 ) +DEAL::6: ( 0 0.707 0.707 ) +DEAL::7: ( 0 0.419 0.419 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( -0.419 0.419 0 ) +DEAL::1: ( -0.211 0.211 0 ) +DEAL::2: ( 0 0.457 0 ) +DEAL::3: ( 0 0.211 0 ) +DEAL::4: ( -0.394 0.394 0.394 ) +DEAL::5: ( -0.211 0.211 0.211 ) +DEAL::6: ( 0 0.419 0.419 ) +DEAL::7: ( 0 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 1.00 0 ) +DEAL::1: ( 0 0.457 0 ) +DEAL::2: ( 0.707 0.707 0 ) +DEAL::3: ( 0.419 0.419 0 ) +DEAL::4: ( 0 0.707 0.707 ) +DEAL::5: ( 0 0.419 0.419 ) +DEAL::6: ( 0.577 0.577 0.577 ) +DEAL::7: ( 0.394 0.394 0.394 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::Cell nodes: +DEAL::0: ( 0 0.457 0 ) +DEAL::1: ( 0 0.211 0 ) +DEAL::2: ( 0.419 0.419 0 ) +DEAL::3: ( 0.211 0.211 0 ) +DEAL::4: ( 0 0.419 0.419 ) +DEAL::5: ( 0 0.211 0.211 ) +DEAL::6: ( 0.394 0.394 0.394 ) +DEAL::7: ( 0.211 0.211 0.211 ) +DEAL::OK: cell bulk and boundary integrals match... + +DEAL::done...