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.

diff --git a/doc/news/changes/minor/20190316SebastianStark b/doc/news/changes/minor/20190316SebastianStark
new file mode 100644 (file)
index 0000000..8ee817e
--- /dev/null
+++ b/doc/news/changes/minor/20190316SebastianStark
@@ -0,0 +1,3 @@
+New: interpolate_to_different_mesh works now with hp::DoFHandler objects
+<br>
+(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 30e5fa5f1fe5e531e2d97c05a731c8a973ac0d1b..3b74b1b5cb307c13b60a44694988d681fa77683a 100644 (file)
--- 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 <int dim,
             int spacedim,
             typename VectorType,
@@ -829,6 +828,9 @@ namespace VectorTools
       intergridmap.get_destination_grid();
     (void)dof2;
 
+    Assert(dof1.get_fe_collection() == dof2.get_fe_collection(),
+           ExcMessage(
+             "The FECollections of both DoFHandler objects must match"));
     Assert(u1.size() == dof1.n_dofs(),
            ExcDimensionMismatch(u1.size(), dof1.n_dofs()));
     Assert(u2.size() == dof2.n_dofs(),
@@ -874,17 +876,28 @@ namespace VectorTools
         if (cell2->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 2fe1aba080878d56732e6e31bab747f864cebe96..7d08e95194c8d5313a7e9315db288c332a67d36e 100644 (file)
--- 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::value_type> &,
         VEC &);
+
+      template void
+      interpolate_to_different_mesh(
+        const hp::DoFHandler<deal_II_dimension, deal_II_space_dimension> &,
+        const VEC &,
+        const hp::DoFHandler<deal_II_dimension, deal_II_space_dimension> &,
+        VEC &);
+
+      template void
+      interpolate_to_different_mesh(
+        const hp::DoFHandler<deal_II_dimension, deal_II_space_dimension> &,
+        const VEC &,
+        const hp::DoFHandler<deal_II_dimension, deal_II_space_dimension> &,
+        const AffineConstraints<VEC::value_type> &,
+        VEC &);
+
+      template void
+      interpolate_to_different_mesh(
+        const InterGridMap<
+          hp::DoFHandler<deal_II_dimension, deal_II_space_dimension>> &,
+        const VEC &,
+        const AffineConstraints<VEC::value_type> &,
+        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 (file)
index 0000000..08663a9
--- /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 <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/dofs/dof_handler.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/mapping_q_generic.h>
+
+#include <deal.II/grid/grid_generator.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <iostream>
+#include <vector>
+
+#include "../tests.h"
+
+using namespace std;
+using namespace dealii;
+
+template <int spacedim>
+class F1 : public Function<spacedim>
+{
+public:
+  F1()
+    : Function<spacedim>()
+  {}
+  virtual double
+  value(const Point<spacedim> &p, const unsigned int component = 0) const
+  {
+    return p(0) * p(0);
+  }
+};
+
+template <unsigned int spacedim>
+void
+check(const unsigned int refinement_1, const unsigned int refinement_2)
+{
+  MappingQGeneric<spacedim> mapping(1);
+
+  Triangulation<spacedim> 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<spacedim> dof_handler_1(tria_1);
+  DoFHandler<spacedim> dof_handler_2(tria_2);
+
+  FE_Q<spacedim> fe(1);
+  dof_handler_1.distribute_dofs(fe);
+  dof_handler_2.distribute_dofs(fe);
+
+  Vector<double> u_1(dof_handler_1.n_dofs()), u_2(dof_handler_2.n_dofs());
+  F1<spacedim>   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<spacedim>                 support_point, unit_support_point;
+  vector<types::global_dof_index> 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 (file)
index 0000000..e413654
--- /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 (file)
index 0000000..1d58e9a
--- /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 <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/mapping_q_generic.h>
+
+#include <deal.II/grid/grid_generator.h>
+
+#include <deal.II/hp/dof_handler.h>
+#include <deal.II/hp/fe_collection.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <iostream>
+#include <vector>
+
+#include "../tests.h"
+
+using namespace std;
+using namespace dealii;
+
+template <int spacedim>
+class F1 : public Function<spacedim>
+{
+public:
+  F1()
+    : Function<spacedim>()
+  {}
+  virtual double
+  value(const Point<spacedim> &p, const unsigned int component = 0) const
+  {
+    return p(0) * p(0);
+  }
+};
+
+template <unsigned int spacedim>
+void
+check(const unsigned int refinement_1, const unsigned int refinement_2)
+{
+  MappingQGeneric<spacedim> mapping(1);
+
+  Triangulation<spacedim> 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<spacedim> dof_handler_1(tria_1);
+  hp::DoFHandler<spacedim> dof_handler_2(tria_2);
+
+  hp::FECollection<spacedim> fe_collection;
+  fe_collection.push_back(FE_Q<spacedim>(2));
+  fe_collection.push_back(FE_Q<spacedim>(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<double> u_1(dof_handler_1.n_dofs()), u_2(dof_handler_2.n_dofs());
+  F1<spacedim>   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<spacedim>                 support_point, unit_support_point;
+  vector<types::global_dof_index> local_dof_indices_1(
+    fe_collection[0].n_dofs_per_cell());
+  vector<types::global_dof_index> 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 (file)
index 0000000..cc4e6ea
--- /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 (file)
index 0000000..0d34358
--- /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 <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/dofs/dof_handler.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/mapping_q_generic.h>
+
+#include <deal.II/grid/grid_generator.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <cstdlib>
+#include <ctime>
+#include <iostream>
+#include <vector>
+
+#include "../tests.h"
+
+using namespace std;
+using namespace dealii;
+
+template <unsigned int spacedim>
+void
+check(const unsigned int refinement_1, const unsigned int refinement_2)
+{
+  MappingQGeneric<spacedim> mapping(1);
+
+  Triangulation<spacedim> 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<spacedim> dof_handler_1(tria_1);
+  DoFHandler<spacedim> dof_handler_2(tria_2);
+
+  FE_Q<spacedim> fe(1);
+  dof_handler_1.distribute_dofs(fe);
+  dof_handler_2.distribute_dofs(fe);
+
+  Vector<double> 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<unsigned int>(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 (file)
index 0000000..9803fd2
--- /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 (file)
index 0000000..f47e039
--- /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 <deal.II/base/quadrature_lib.h>
+
+#include <deal.II/fe/fe_q.h>
+#include <deal.II/fe/mapping_q_generic.h>
+
+#include <deal.II/grid/grid_generator.h>
+
+#include <deal.II/hp/dof_handler.h>
+#include <deal.II/hp/fe_collection.h>
+
+#include <deal.II/numerics/vector_tools.h>
+
+#include <cstdlib>
+#include <ctime>
+#include <iostream>
+#include <vector>
+
+#include "../tests.h"
+
+using namespace std;
+using namespace dealii;
+
+template <unsigned int spacedim>
+void
+check(const unsigned int refinement_1, const unsigned int refinement_2)
+{
+  MappingQGeneric<spacedim> mapping(1);
+
+  Triangulation<spacedim> 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<spacedim> dof_handler_1(tria_1);
+  hp::DoFHandler<spacedim> dof_handler_2(tria_2);
+
+  hp::FECollection<spacedim> fe_collection;
+  fe_collection.push_back(FE_Q<spacedim>(1));
+  fe_collection.push_back(FE_Q<spacedim>(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<double> 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<unsigned int>(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 (file)
index 0000000..b199f69
--- /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::