From 8d8fb9c8388ff35237ba183cd2af698e20ac2ec2 Mon Sep 17 00:00:00 2001 From: Sebastian Stark Date: Fri, 22 Feb 2019 15:09:35 +0200 Subject: [PATCH] interpolate_to_different_mesh for hp::DoFHandler Make interpolate_to_different_mesh work with hp::DoFHandler objects. Added also four tests for the implementation as well as the relevant instantiations. --- doc/news/changes/minor/20190316SebastianStark | 3 + .../deal.II/numerics/vector_tools.templates.h | 34 +- .../numerics/vector_tools_interpolate.inst.in | 23 + .../interpolate_to_different_mesh_01.cc | 152 ++ .../interpolate_to_different_mesh_01.output | 1069 ++++++++ .../interpolate_to_different_mesh_02.cc | 164 ++ .../interpolate_to_different_mesh_02.output | 2173 +++++++++++++++++ .../interpolate_to_different_mesh_03.cc | 109 + .../interpolate_to_different_mesh_03.output | 17 + .../interpolate_to_different_mesh_04.cc | 133 + .../interpolate_to_different_mesh_04.output | 25 + 11 files changed, 3891 insertions(+), 11 deletions(-) create mode 100644 doc/news/changes/minor/20190316SebastianStark create mode 100644 tests/numerics/interpolate_to_different_mesh_01.cc create mode 100644 tests/numerics/interpolate_to_different_mesh_01.output create mode 100644 tests/numerics/interpolate_to_different_mesh_02.cc create mode 100644 tests/numerics/interpolate_to_different_mesh_02.output create mode 100644 tests/numerics/interpolate_to_different_mesh_03.cc create mode 100644 tests/numerics/interpolate_to_different_mesh_03.output create mode 100644 tests/numerics/interpolate_to_different_mesh_04.cc create mode 100644 tests/numerics/interpolate_to_different_mesh_04.output diff --git a/doc/news/changes/minor/20190316SebastianStark b/doc/news/changes/minor/20190316SebastianStark new file mode 100644 index 0000000000..8ee817eb59 --- /dev/null +++ b/doc/news/changes/minor/20190316SebastianStark @@ -0,0 +1,3 @@ +New: interpolate_to_different_mesh works now with hp::DoFHandler objects +
+(Sebastian Stark, 2019/03/16) diff --git a/include/deal.II/numerics/vector_tools.templates.h b/include/deal.II/numerics/vector_tools.templates.h index 30e5fa5f1f..3b74b1b5cb 100644 --- a/include/deal.II/numerics/vector_tools.templates.h +++ b/include/deal.II/numerics/vector_tools.templates.h @@ -812,7 +812,6 @@ namespace VectorTools } } // namespace internal - template active() && !cell2->is_locally_owned()) continue; - Assert( - cell1->get_fe().get_name() == cell2->get_fe().get_name(), - ExcMessage( - "Source and destination cells need to use the same finite element")); - - cache.reinit(cell1->get_fe().dofs_per_cell); - // Get and set the corresponding // dof_values by interpolation. - cell1->get_interpolated_dof_values(u1, cache); - cell2->set_dof_values_by_interpolation(cache, u2); + if (cell1->active()) + { + cache.reinit(cell1->get_fe().dofs_per_cell); + cell1->get_interpolated_dof_values(u1, + cache, + cell1->active_fe_index()); + cell2->set_dof_values_by_interpolation(cache, + u2, + cell1->active_fe_index()); + } + else + { + cache.reinit(cell2->get_fe().dofs_per_cell); + cell1->get_interpolated_dof_values(u1, + cache, + cell2->active_fe_index()); + cell2->set_dof_values_by_interpolation(cache, + u2, + cell2->active_fe_index()); + } } // finish the work on parallel vectors @@ -893,7 +906,6 @@ namespace VectorTools constraints.distribute(u2); } - namespace internal { /** diff --git a/source/numerics/vector_tools_interpolate.inst.in b/source/numerics/vector_tools_interpolate.inst.in index 2fe1aba080..7d08e95194 100644 --- a/source/numerics/vector_tools_interpolate.inst.in +++ b/source/numerics/vector_tools_interpolate.inst.in @@ -132,6 +132,29 @@ for (VEC : VECTOR_TYPES; deal_II_dimension : DIMENSIONS; const VEC &, const AffineConstraints &, VEC &); + + template void + interpolate_to_different_mesh( + const hp::DoFHandler &, + const VEC &, + const hp::DoFHandler &, + VEC &); + + template void + interpolate_to_different_mesh( + const hp::DoFHandler &, + const VEC &, + const hp::DoFHandler &, + const AffineConstraints &, + VEC &); + + template void + interpolate_to_different_mesh( + const InterGridMap< + hp::DoFHandler> &, + const VEC &, + const AffineConstraints &, + VEC &); \} #endif } diff --git a/tests/numerics/interpolate_to_different_mesh_01.cc b/tests/numerics/interpolate_to_different_mesh_01.cc new file mode 100644 index 0000000000..08663a9041 --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_01.cc @@ -0,0 +1,152 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2000 - 2018 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE.md at +// the top level directory of deal.II. +// +// --------------------------------------------------------------------- + +// test for interpolate_to_different_mesh with dealii::DoFHandler + +#include + +#include + +#include +#include + +#include + +#include + +#include +#include + +#include "../tests.h" + +using namespace std; +using namespace dealii; + +template +class F1 : public Function +{ +public: + F1() + : Function() + {} + virtual double + value(const Point &p, const unsigned int component = 0) const + { + return p(0) * p(0); + } +}; + +template +void +check(const unsigned int refinement_1, const unsigned int refinement_2) +{ + MappingQGeneric mapping(1); + + Triangulation tria_1, tria_2; + GridGenerator::hyper_cube(tria_1); + GridGenerator::hyper_cube(tria_2); + + tria_1.refine_global(refinement_1); + tria_2.refine_global(refinement_2); + + DoFHandler dof_handler_1(tria_1); + DoFHandler dof_handler_2(tria_2); + + FE_Q fe(1); + dof_handler_1.distribute_dofs(fe); + dof_handler_2.distribute_dofs(fe); + + Vector u_1(dof_handler_1.n_dofs()), u_2(dof_handler_2.n_dofs()); + F1 f_1; + VectorTools::interpolate(mapping, dof_handler_1, f_1, u_1); + VectorTools::interpolate_to_different_mesh(dof_handler_1, + u_1, + dof_handler_2, + u_2); + + Point support_point, unit_support_point; + vector local_dof_indices(fe.n_dofs_per_cell()); + deallog << "dof values of dof_handler_1:" << endl; + for (const auto &cell : dof_handler_1.active_cell_iterators()) + { + cell->get_dof_indices(local_dof_indices); + for (unsigned int shapefun = 0; + shapefun < cell->get_fe().n_dofs_per_cell(); + ++shapefun) + { + unit_support_point = cell->get_fe().unit_support_point(shapefun); + support_point = + mapping.transform_unit_to_real_cell(cell, unit_support_point); + deallog << " " << support_point << ":\n " + << u_1[local_dof_indices[shapefun]] << endl + << endl; + } + } + deallog << endl << "dof values of dof_handler_2:" << endl; + for (const auto &cell : dof_handler_2.active_cell_iterators()) + + { + cell->get_dof_indices(local_dof_indices); + for (unsigned int shapefun = 0; + shapefun < cell->get_fe().n_dofs_per_cell(); + ++shapefun) + { + unit_support_point = cell->get_fe().unit_support_point(shapefun); + support_point = + mapping.transform_unit_to_real_cell(cell, unit_support_point); + deallog << " " << support_point << ":\n " + << u_2[local_dof_indices[shapefun]] << endl + << endl; + } + } +} + +int +main(int argc, char **argv) +{ + initlog(); + + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + + deallog + << "### 2D-Case, first cell refined once, second cell unrefined###\n\n"; + check<2>(1, 0); + deallog << endl; + + deallog + << "### 2D-Case, first cell unrefined, second cell refined once###\n\n"; + check<2>(0, 1); + deallog << endl; + + deallog + << "### 2D-Case, first cell refined once, second cell refined once###\n\n"; + check<2>(1, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell refined once, second cell unrefined###\n\n"; + check<3>(1, 0); + deallog << endl; + + deallog + << "### 3D-Case, first cell unrefined, second cell refined once###\n\n"; + check<3>(0, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell refined once, second cell refined once###\n\n"; + check<3>(1, 1); + deallog << endl; +} diff --git a/tests/numerics/interpolate_to_different_mesh_01.output b/tests/numerics/interpolate_to_different_mesh_01.output new file mode 100644 index 0000000000..e413654cda --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_01.output @@ -0,0 +1,1069 @@ + +DEAL::### 2D-Case, first cell refined once, second cell unrefined### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 2D-Case, first cell unrefined, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.500000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 2D-Case, first cell refined once, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 3D-Case, first cell refined once, second cell unrefined### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 3D-Case, first cell unrefined, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.500000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.500000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.500000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.500000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.500000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 3D-Case, first cell refined once, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: diff --git a/tests/numerics/interpolate_to_different_mesh_02.cc b/tests/numerics/interpolate_to_different_mesh_02.cc new file mode 100644 index 0000000000..1d58e9a318 --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_02.cc @@ -0,0 +1,164 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2000 - 2018 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE.md at +// the top level directory of deal.II. +// +// --------------------------------------------------------------------- + +// test for interpolate_to_different_mesh with dealii::hp::DoFHandler + +#include + +#include +#include + +#include + +#include +#include + +#include + +#include +#include + +#include "../tests.h" + +using namespace std; +using namespace dealii; + +template +class F1 : public Function +{ +public: + F1() + : Function() + {} + virtual double + value(const Point &p, const unsigned int component = 0) const + { + return p(0) * p(0); + } +}; + +template +void +check(const unsigned int refinement_1, const unsigned int refinement_2) +{ + MappingQGeneric mapping(1); + + Triangulation tria_1, tria_2; + GridGenerator::hyper_cube(tria_1); + GridGenerator::hyper_cube(tria_2); + + tria_1.refine_global(refinement_1); + tria_2.refine_global(refinement_2); + + hp::DoFHandler dof_handler_1(tria_1); + hp::DoFHandler dof_handler_2(tria_2); + + hp::FECollection fe_collection; + fe_collection.push_back(FE_Q(2)); + fe_collection.push_back(FE_Q(1)); + + for (const auto &cell : dof_handler_1.active_cell_iterators()) + cell->set_active_fe_index(0); + + for (const auto &cell : dof_handler_2.active_cell_iterators()) + cell->set_active_fe_index(1); + + dof_handler_1.distribute_dofs(fe_collection); + dof_handler_2.distribute_dofs(fe_collection); + + Vector u_1(dof_handler_1.n_dofs()), u_2(dof_handler_2.n_dofs()); + F1 f_1; + VectorTools::interpolate(mapping, dof_handler_1, f_1, u_1); + VectorTools::interpolate_to_different_mesh(dof_handler_1, + u_1, + dof_handler_2, + u_2); + + Point support_point, unit_support_point; + vector local_dof_indices_1( + fe_collection[0].n_dofs_per_cell()); + vector local_dof_indices_2( + fe_collection[1].n_dofs_per_cell()); + deallog << "dof values of dof_handler_1:" << endl; + for (const auto &cell : dof_handler_1.active_cell_iterators()) + { + cell->get_dof_indices(local_dof_indices_1); + for (unsigned int shapefun = 0; + shapefun < cell->get_fe().n_dofs_per_cell(); + ++shapefun) + { + unit_support_point = cell->get_fe().unit_support_point(shapefun); + support_point = + mapping.transform_unit_to_real_cell(cell, unit_support_point); + deallog << " " << support_point << ":\n " + << u_1[local_dof_indices_1[shapefun]] << endl + << endl; + } + } + deallog << endl << "dof values of dof_handler_2:" << endl; + for (const auto &cell : dof_handler_2.active_cell_iterators()) + { + cell->get_dof_indices(local_dof_indices_2); + for (unsigned int shapefun = 0; + shapefun < cell->get_fe().n_dofs_per_cell(); + ++shapefun) + { + unit_support_point = cell->get_fe().unit_support_point(shapefun); + support_point = + mapping.transform_unit_to_real_cell(cell, unit_support_point); + deallog << " " << support_point << ":\n " + << u_2[local_dof_indices_2[shapefun]] << endl + << endl; + } + } +} + +int +main(int argc, char **argv) +{ + initlog(); + + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + + deallog + << "### 2D-Case, first cell refined once, second cell unrefined###\n\n"; + check<2>(1, 0); + deallog << endl; + + deallog + << "### 2D-Case, first cell unrefined, second cell refined once###\n\n"; + check<2>(0, 1); + deallog << endl; + + deallog + << "### 2D-Case, first cell refined once, second cell refined once###\n\n"; + check<2>(1, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell refined once, second cell unrefined###\n\n"; + check<3>(1, 0); + deallog << endl; + + deallog + << "### 3D-Case, first cell unrefined, second cell refined once###\n\n"; + check<3>(0, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell refined once, second cell refined once###\n\n"; + check<3>(1, 1); + deallog << endl; +} diff --git a/tests/numerics/interpolate_to_different_mesh_02.output b/tests/numerics/interpolate_to_different_mesh_02.output new file mode 100644 index 0000000000..cc4e6eab72 --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_02.output @@ -0,0 +1,2173 @@ + +DEAL::### 2D-Case, first cell refined once, second cell unrefined### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000: + 0.562500 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 2D-Case, first cell unrefined, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 2D-Case, first cell refined once, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000: + 0.562500 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 3D-Case, first cell refined once, second cell unrefined### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.250000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.00000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.750000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.500000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.250000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.00000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.250000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.750000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.500000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.750000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.750000: + 0.562500 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 3D-Case, first cell unrefined, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: +DEAL::### 3D-Case, first cell refined once, second cell refined once### + +dof values of dof_handler_1: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.250000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.00000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.750000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.500000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.250000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.250000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.250000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.500000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.250000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.250000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.250000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.250000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.00000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.250000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.00000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.500000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.250000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.250000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.250000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.250000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.00000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.500000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.250000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.750000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.00000 0.500000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.750000 0.750000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.750000: + 0.250000 +DEAL:: +DEAL:: 0.250000 0.500000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 1.00000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.500000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 1.00000: + 0.0625000 +DEAL:: +DEAL:: 0.250000 0.750000 0.750000: + 0.0625000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.750000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.500000 0.500000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.750000 0.750000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.750000 0.750000: + 1.00000 +DEAL:: +DEAL:: 0.750000 0.500000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 1.00000 0.750000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.500000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 1.00000: + 0.562500 +DEAL:: +DEAL:: 0.750000 0.750000 0.750000: + 0.562500 +DEAL:: +DEAL:: +DEAL::dof values of dof_handler_2: +DEAL:: 0.00000 0.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.00000 0.500000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 0.500000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 0.00000 0.500000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.00000 1.00000 1.00000: + 0.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 0.500000 0.500000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 0.500000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 0.500000: + 1.00000 +DEAL:: +DEAL:: 0.500000 0.500000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 0.500000 1.00000: + 1.00000 +DEAL:: +DEAL:: 0.500000 1.00000 1.00000: + 0.250000 +DEAL:: +DEAL:: 1.00000 1.00000 1.00000: + 1.00000 +DEAL:: +DEAL:: diff --git a/tests/numerics/interpolate_to_different_mesh_03.cc b/tests/numerics/interpolate_to_different_mesh_03.cc new file mode 100644 index 0000000000..0d34358e6a --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_03.cc @@ -0,0 +1,109 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2000 - 2018 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE.md at +// the top level directory of deal.II. +// +// --------------------------------------------------------------------- + +// test for interpolate_to_different_mesh with dealii::DoFHandler +// checks that interpolation from coarse to fine and back to coarse +// is the identity operation + +#include + +#include + +#include +#include + +#include + +#include + +#include +#include +#include +#include + +#include "../tests.h" + +using namespace std; +using namespace dealii; + +template +void +check(const unsigned int refinement_1, const unsigned int refinement_2) +{ + MappingQGeneric mapping(1); + + Triangulation tria_1, tria_2; + GridGenerator::hyper_cube(tria_1); + GridGenerator::hyper_cube(tria_2); + + tria_1.refine_global(refinement_1); + tria_2.refine_global(refinement_2); + + DoFHandler dof_handler_1(tria_1); + DoFHandler dof_handler_2(tria_2); + + FE_Q fe(1); + dof_handler_1.distribute_dofs(fe); + dof_handler_2.distribute_dofs(fe); + + Vector u_1(dof_handler_1.n_dofs()), u_2(dof_handler_2.n_dofs()), + u_1_back(dof_handler_1.n_dofs()), d(dof_handler_1.n_dofs()); + for (unsigned int i1 = 0; i1 < u_1.size(); ++i1) + u_1[i1] = 2. * rand() / RAND_MAX - 1.; + + VectorTools::interpolate_to_different_mesh(dof_handler_1, + u_1, + dof_handler_2, + u_2); + + VectorTools::interpolate_to_different_mesh(dof_handler_2, + u_2, + dof_handler_1, + u_1_back); + + d = u_1; + d.sadd(-1.0, u_1_back); + deallog << "distance=" << d.l2_norm() << endl; +} + +int +main(int argc, char **argv) +{ + initlog(); + + srand(static_cast(clock())); + + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + + deallog + << "### 2D-Case, first cell unrefined, second cell refined once###\n\n"; + check<2>(0, 1); + deallog << endl; + + deallog + << "### 2D-Case, first cell refined once, second cell refined three times###\n\n"; + check<2>(1, 3); + deallog << endl; + + deallog + << "### 3D-Case, first cell unrefined, second cell refined once###\n\n"; + check<3>(0, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell refined once, second cell refined three times###\n\n"; + check<3>(1, 3); + deallog << endl; +} diff --git a/tests/numerics/interpolate_to_different_mesh_03.output b/tests/numerics/interpolate_to_different_mesh_03.output new file mode 100644 index 0000000000..9803fd2d68 --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_03.output @@ -0,0 +1,17 @@ + +DEAL::### 2D-Case, first cell unrefined, second cell refined once### + +distance=0.00000 +DEAL:: +DEAL::### 2D-Case, first cell refined once, second cell refined three times### + +distance=0.00000 +DEAL:: +DEAL::### 3D-Case, first cell unrefined, second cell refined once### + +distance=0.00000 +DEAL:: +DEAL::### 3D-Case, first cell refined once, second cell refined three times### + +distance=0.00000 +DEAL:: diff --git a/tests/numerics/interpolate_to_different_mesh_04.cc b/tests/numerics/interpolate_to_different_mesh_04.cc new file mode 100644 index 0000000000..f47e03980f --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_04.cc @@ -0,0 +1,133 @@ +// --------------------------------------------------------------------- +// +// Copyright (C) 2000 - 2018 by the deal.II authors +// +// This file is part of the deal.II library. +// +// The deal.II library is free software; you can use it, redistribute +// it, and/or modify it under the terms of the GNU Lesser General +// Public License as published by the Free Software Foundation; either +// version 2.1 of the License, or (at your option) any later version. +// The full text of the license can be found in the file LICENSE.md at +// the top level directory of deal.II. +// +// --------------------------------------------------------------------- + +// test for interpolate_to_different_mesh with dealii::hp::DoFHandler +// checks that interpolation from coarse to fine and back to coarse +// is the identity operation + +#include + +#include +#include + +#include + +#include +#include + +#include + +#include +#include +#include +#include + +#include "../tests.h" + +using namespace std; +using namespace dealii; + +template +void +check(const unsigned int refinement_1, const unsigned int refinement_2) +{ + MappingQGeneric mapping(1); + + Triangulation tria_1, tria_2; + GridGenerator::hyper_cube(tria_1); + GridGenerator::hyper_cube(tria_2); + + tria_1.refine_global(refinement_1); + tria_2.refine_global(refinement_2); + + hp::DoFHandler dof_handler_1(tria_1); + hp::DoFHandler dof_handler_2(tria_2); + + hp::FECollection fe_collection; + fe_collection.push_back(FE_Q(1)); + fe_collection.push_back(FE_Q(2)); + + for (const auto &cell : dof_handler_1.active_cell_iterators()) + cell->set_active_fe_index(0); + + for (const auto &cell : dof_handler_2.active_cell_iterators()) + cell->set_active_fe_index(1); + + dof_handler_1.distribute_dofs(fe_collection); + dof_handler_2.distribute_dofs(fe_collection); + + Vector u_1(dof_handler_1.n_dofs()), u_2(dof_handler_2.n_dofs()), + u_1_back(dof_handler_1.n_dofs()), d(dof_handler_1.n_dofs()); + for (unsigned int i1 = 0; i1 < u_1.size(); ++i1) + u_1[i1] = 2. * rand() / RAND_MAX - 1.; + + VectorTools::interpolate_to_different_mesh(dof_handler_1, + u_1, + dof_handler_2, + u_2); + + VectorTools::interpolate_to_different_mesh(dof_handler_2, + u_2, + dof_handler_1, + u_1_back); + + d = u_1; + d.sadd(-1.0, u_1_back); + deallog << "distance=" << d.l2_norm() << endl; +} + +int +main(int argc, char **argv) +{ + initlog(); + + srand(static_cast(clock())); + + Utilities::MPI::MPI_InitFinalize mpi_initialization(argc, argv, 1); + + deallog + << "### 2D-Case, first cell unrefined, second cell refined once###\n\n"; + check<2>(0, 1); + deallog << endl; + + deallog + << "### 2D-Case, first cell refined once, second cell refined three times###\n\n"; + check<2>(1, 3); + deallog << endl; + + // this should still be an identity operation although interpolation is from + // fine to coarse mesh because second mesh uses quadratic elements + deallog + << "### 2D-Case, first cell refined twice, second cell refined once###\n\n"; + check<2>(2, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell unrefined, second cell refined once###\n\n"; + check<3>(0, 1); + deallog << endl; + + deallog + << "### 3D-Case, first cell refined once, second cell refined three times###\n\n"; + check<3>(1, 3); + deallog << endl; + + // this should still be an identity operation although interpolation is from + // fine to coarse mesh because second mesh uses quadratic elements + deallog + << "### 3D-Case, first cell refined twice, second cell refined once###\n\n"; + check<3>(2, 1); + deallog << endl; +} diff --git a/tests/numerics/interpolate_to_different_mesh_04.output b/tests/numerics/interpolate_to_different_mesh_04.output new file mode 100644 index 0000000000..b199f69ae9 --- /dev/null +++ b/tests/numerics/interpolate_to_different_mesh_04.output @@ -0,0 +1,25 @@ + +DEAL::### 2D-Case, first cell unrefined, second cell refined once### + +distance=0.00000 +DEAL:: +DEAL::### 2D-Case, first cell refined once, second cell refined three times### + +distance=0.00000 +DEAL:: +DEAL::### 2D-Case, first cell refined twice, second cell refined once### + +distance=0.00000 +DEAL:: +DEAL::### 3D-Case, first cell unrefined, second cell refined once### + +distance=0.00000 +DEAL:: +DEAL::### 3D-Case, first cell refined once, second cell refined three times### + +distance=0.00000 +DEAL:: +DEAL::### 3D-Case, first cell refined twice, second cell refined once### + +distance=0.00000 +DEAL:: -- 2.39.5