Output all faces and edges to vtk

diff --git a/doc/news/changes/minor/20191011PeterMunch b/doc/news/changes/minor/20191011PeterMunch
index c2af4e67043ed0f84a4c54185298bd144c4c3001..263ecdd251da5d2e6232f06d9a9836939438eba6 100644 (file)
--- a/doc/news/changes/minor/20191011PeterMunch
+++ b/doc/news/changes/minor/20191011PeterMunch
@@ -1,4 +1,4 @@
-Improved: Make the output of cells, faces and co-faces optional in
-GridOut::write_vtk.
+Improved: Make the output of cells, (all non-internal) faces and (all non-internal) 
+co-faces optional in GridOut::write_vtk.
 <br>
 (Peter Munch, Niklas Fehn, 2019/10/11)

diff --git a/include/deal.II/grid/grid_out.h b/include/deal.II/grid/grid_out.h
index b05c7e57f73685f82e37b324c3577dc90888783c..08ce73cc2da4da13c11e0018145c60a7ba3ad9f6 100644 (file)
--- a/include/deal.II/grid/grid_out.h
+++ b/include/deal.II/grid/grid_out.h
@@ -869,12 +869,14 @@ namespace GridOutFlags
     /**
      * Default constructor.
      */
-    Vtk(const bool output_cells    = true,
-        const bool output_faces    = true,
-        const bool output_co_faces = true)
+    Vtk(const bool output_cells         = true,
+        const bool output_faces         = true,
+        const bool output_edges         = true,
+        const bool output_only_relevant = true)
       : output_cells(output_cells)
       , output_faces(output_faces)
-      , output_co_faces(output_co_faces)
+      , output_edges(output_edges)
+      , output_only_relevant(output_only_relevant)
     {}
 
     /**
@@ -890,7 +892,13 @@ namespace GridOutFlags
     /**
      * Output co-faces/edges.
      */
-    bool output_co_faces;
+    bool output_edges;
+
+    /**
+     * Output only faces/co-faces that differ from the default settings
+     * (e.g boundary_id).
+     */
+    bool output_only_relevant;
   };
 
 

diff --git a/source/grid/grid_out.cc b/source/grid/grid_out.cc
index 96c063561201fd8da738c68da343b81b99ea4d65..66874f19df3b370b4d2d290d67e52456c94cb63c 100644 (file)
--- a/source/grid/grid_out.cc
+++ b/source/grid/grid_out.cc
@@ -2951,13 +2951,53 @@ namespace
     return v;
   }
 
+  /**
+   * Return all boundary lines of non-internal faces in three dimension.
+   */
+  std::vector<typename Triangulation<3, 3>::active_line_iterator>
+  get_boundary_edge_iterators(const Triangulation<3, 3> &tria)
+  {
+    std::vector<typename Triangulation<3, 3>::active_line_iterator> res;
+
+    std::vector<bool> flags;
+    tria.save_user_flags_line(flags);
+    const_cast<Triangulation<3, 3> &>(tria).clear_user_flags_line();
+
+    for (auto face : tria.active_face_iterators())
+      for (unsigned int l = 0; l < GeometryInfo<3>::lines_per_face; ++l)
+        {
+          const auto line = face->line(l);
+          if (line->user_flag_set() || line->has_children())
+            continue;
+          else
+            line->set_user_flag();
+          if (line->at_boundary())
+            res.emplace_back(line);
+        }
+    const_cast<Triangulation<3, 3> &>(tria).load_user_flags_line(flags);
+    return res;
+  }
+
+
+
+  /**
+   * Same as above, for 1 and 2 dimensional grids. Does nothing.
+   */
+  template <int dim, int spacedim>
+  std::vector<typename Triangulation<dim, spacedim>::active_line_iterator>
+  get_boundary_edge_iterators(const Triangulation<dim, spacedim> &)
+  {
+    return {};
+  }
+
+
 
   /**
    * Return all lines of a face in three dimension that have a non-standard
    * boundary indicator (!=0), or a non-flat manifold indicator.
    */
   std::vector<typename Triangulation<3, 3>::active_line_iterator>
-  relevant_co_faces(const Triangulation<3, 3> &tria)
+  get_relevant_edge_iterators(const Triangulation<3, 3> &tria)
   {
     std::vector<typename Triangulation<3, 3>::active_line_iterator> res;
 
@@ -2965,10 +3005,10 @@ namespace
     tria.save_user_flags_line(flags);
     const_cast<Triangulation<3, 3> &>(tria).clear_user_flags_line();
 
-    for (auto face = tria.begin_active_face(); face != tria.end_face(); ++face)
+    for (auto face : tria.active_face_iterators())
       for (unsigned int l = 0; l < GeometryInfo<3>::lines_per_face; ++l)
         {
-          auto line = face->line(l);
+          const auto line = face->line(l);
           if (line->user_flag_set() || line->has_children())
             continue;
           else
@@ -2988,26 +3028,47 @@ namespace
    */
   template <int dim, int spacedim>
   std::vector<typename Triangulation<dim, spacedim>::active_line_iterator>
-  relevant_co_faces(const Triangulation<dim, spacedim> &)
+  get_relevant_edge_iterators(const Triangulation<dim, spacedim> &)
   {
     return {};
   }
 
 
 
+  /**
+   * Return all boundary faces of a triangulation.
+   */
+  template <int dim, int spacedim>
+  std::vector<typename Triangulation<dim, spacedim>::active_face_iterator>
+  get_boundary_face_iterators(const Triangulation<dim, spacedim> &tria)
+  {
+    std::vector<typename Triangulation<dim, spacedim>::active_face_iterator>
+      res;
+    if (dim == 1)
+      return res;
+    for (auto face : tria.active_face_iterators())
+      {
+        if (face->boundary_id() != numbers::invalid_boundary_id)
+          res.push_back(face);
+      }
+    return res;
+  }
+
+
+
   /**
    * Return all faces of a triangulation that have a non-standard
    * boundary indicator (!=0), or a non-flat manifold indicator.
    */
   template <int dim, int spacedim>
   std::vector<typename Triangulation<dim, spacedim>::active_face_iterator>
-  relevant_faces(const Triangulation<dim, spacedim> &tria)
+  get_relevant_face_iterators(const Triangulation<dim, spacedim> &tria)
   {
     std::vector<typename Triangulation<dim, spacedim>::active_face_iterator>
       res;
     if (dim == 1)
       return res;
-    for (auto face = tria.begin_active_face(); face != tria.end_face(); ++face)
+    for (auto face : tria.active_face_iterators())
       {
         if (face->manifold_id() != numbers::flat_manifold_id ||
             (face->boundary_id() != 0 &&
@@ -3047,19 +3108,23 @@ GridOut::write_vtk(const Triangulation<dim, spacedim> &tria,
       out << '\n';
     }
 
-  const auto faces    = relevant_faces(tria);
-  const auto co_faces = relevant_co_faces(tria);
+  const auto faces = vtk_flags.output_only_relevant ?
+                       get_relevant_face_iterators(tria) :
+                       get_boundary_face_iterators(tria);
+  const auto edges = vtk_flags.output_only_relevant ?
+                       get_relevant_edge_iterators(tria) :
+                       get_boundary_edge_iterators(tria);
 
   AssertThrow(
-    vtk_flags.output_cells || (dim > 2 && vtk_flags.output_faces) ||
-      (dim > 3 && vtk_flags.output_co_faces),
+    vtk_flags.output_cells || (dim >= 2 && vtk_flags.output_faces) ||
+      (dim >= 3 && vtk_flags.output_edges),
     ExcMessage(
-      "At least one of the flags (output_cells, output_faces, output_co_faces) has to be enabled!"));
+      "At least one of the flags (output_cells, output_faces, output_edges) has to be enabled!"));
 
   // Write cells preamble
   const int n_cells = (vtk_flags.output_cells ? tria.n_active_cells() : 0) +
                       (vtk_flags.output_faces ? faces.size() : 0) +
-                      (vtk_flags.output_co_faces ? co_faces.size() : 0);
+                      (vtk_flags.output_edges ? edges.size() : 0);
 
   // VTK now expects a number telling the total storage requirement to read all
   // cell connectivity information. The connectivity information is read cell by
@@ -3074,8 +3139,8 @@ GridOut::write_vtk(const Triangulation<dim, spacedim> &tria,
     (vtk_flags.output_faces ?
        faces.size() * (GeometryInfo<dim>::vertices_per_face + 1) :
        0) +
-    (vtk_flags.output_co_faces ? co_faces.size() * (3) :
-                                 0); // only in 3d, otherwise it is always zero.
+    (vtk_flags.output_edges ? edges.size() * (3) :
+                              0); // only in 3d, otherwise it is always zero.
 
   AssertThrow(cells_size > 0, ExcMessage("No cells given to be output!"));
 
@@ -3117,12 +3182,12 @@ GridOut::write_vtk(const Triangulation<dim, spacedim> &tria,
           }
         out << '\n';
       }
-  if (vtk_flags.output_co_faces)
-    for (const auto &co_face : co_faces)
+  if (vtk_flags.output_edges)
+    for (const auto &edge : edges)
       {
         out << 2;
         for (unsigned int i = 0; i < 2; ++i)
-          out << ' ' << co_face->vertex_index(i);
+          out << ' ' << edge->vertex_index(i);
         out << '\n';
       }
 
@@ -3144,9 +3209,9 @@ GridOut::write_vtk(const Triangulation<dim, spacedim> &tria,
         }
       out << '\n';
     }
-  if (vtk_flags.output_co_faces)
+  if (vtk_flags.output_edges)
     {
-      for (unsigned int i = 0; i < co_faces.size(); ++i)
+      for (unsigned int i = 0; i < edges.size(); ++i)
         {
           out << co_face_type << ' ';
         }
@@ -3176,12 +3241,12 @@ GridOut::write_vtk(const Triangulation<dim, spacedim> &tria,
         }
       out << '\n';
     }
-  if (vtk_flags.output_co_faces)
+  if (vtk_flags.output_edges)
     {
-      for (const auto &co_face : co_faces)
+      for (const auto &edge : edges)
         {
           out << static_cast<std::make_signed<types::boundary_id>::type>(
-                   co_face->boundary_id())
+                   edge->boundary_id())
               << ' ';
         }
     }
@@ -3210,12 +3275,12 @@ GridOut::write_vtk(const Triangulation<dim, spacedim> &tria,
         }
       out << '\n';
     }
-  if (vtk_flags.output_co_faces)
+  if (vtk_flags.output_edges)
     {
-      for (const auto &co_face : co_faces)
+      for (const auto &edge : edges)
         {
           out << static_cast<std::make_signed<types::boundary_id>::type>(
-                   co_face->manifold_id())
+                   edge->manifold_id())
               << ' ';
         }
       out << '\n';

diff --git a/tests/grid/grid_out_vtk_03.cc b/tests/grid/grid_out_vtk_03.cc
index f5c9d6b524cebb2651a8496fa01cf16c15a40f48..b70244e654c843a332d2f8840de287b1e17033a6 100644 (file)
--- a/tests/grid/grid_out_vtk_03.cc
+++ b/tests/grid/grid_out_vtk_03.cc
@@ -37,7 +37,8 @@ void
 do_test(std::ostream &logfile,
         const bool    output_cells,
         const bool    output_faces,
-        const bool    output_co_faces)
+        const bool    output_edges,
+        const bool    output_only_relevant)
 {
   Triangulation<dim, spacedim> tria_1;
   GridGenerator::hyper_cube(tria_1, 0, 1, true);
@@ -47,15 +48,15 @@ do_test(std::ostream &logfile,
   tria_1.execute_coarsening_and_refinement();
 
   GridOutFlags::Vtk flags;
-  flags.output_cells    = output_cells;
-  flags.output_faces    = output_faces;
-  flags.output_co_faces = output_co_faces;
+  flags.output_cells         = output_cells;
+  flags.output_faces         = output_faces;
+  flags.output_edges         = output_edges;
+  flags.output_only_relevant = output_only_relevant;
 
   GridOut grid_out_1;
   grid_out_1.set_flags(flags);
   grid_out_1.write_vtk(tria_1, logfile);
 
-
   // write to buffer
   std::ostringstream buf;
   grid_out_1.write_vtk(tria_1, buf);
@@ -82,10 +83,15 @@ template <int dim, int spacedim>
 void
 test(std::ostream &logfile)
 {
-  do_test<dim, spacedim>(logfile, true, false, false);
-  do_test<dim, spacedim>(logfile, true, false, true);
-  do_test<dim, spacedim>(logfile, true, true, false);
-  do_test<dim, spacedim>(logfile, true, true, true);
+  do_test<dim, spacedim>(logfile, true, false, false, true);
+  do_test<dim, spacedim>(logfile, true, false, true, true);
+  do_test<dim, spacedim>(logfile, true, true, false, true);
+  do_test<dim, spacedim>(logfile, true, true, true, true);
+
+  do_test<dim, spacedim>(logfile, true, false, false, false);
+  do_test<dim, spacedim>(logfile, true, false, true, false);
+  do_test<dim, spacedim>(logfile, true, true, false, false);
+  do_test<dim, spacedim>(logfile, true, true, true, false);
 }
 
 

diff --git a/tests/grid/grid_out_vtk_03.output b/tests/grid/grid_out_vtk_03.output
index aa95b5534fe44fa5d4d53480078200fb325bd726..92140c5435abdb8743c9ffa52674c081f700588d 100644 (file)
--- a/tests/grid/grid_out_vtk_03.output
+++ b/tests/grid/grid_out_vtk_03.output
@@ -252,10 +252,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -280,10 +280,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -312,10 +312,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -341,10 +341,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -373,10 +373,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -404,10 +404,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -436,10 +436,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -468,10 +468,10 @@ Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
 POINTS 4 double
-0.00000 0.00000 0
-1.00000 0.00000 0
-0.500000 0.00000 0
-0.250000 0.00000 0
+0.00000 0 0
+1.00000 0 0
+0.500000 0 0
+0.250000 0 0
 
 CELLS 3 9
 2 2 1
@@ -499,778 +499,2974 @@ LOOKUP_TABLE default
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
 
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
 
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 4 1
-2 1 6
-2 6 3
-2 2 7
-2 7 3
-2 0 9
-2 9 4
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 14
+
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 1 3 3 2 2 
+0 0 0 
+
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
+
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 0 9
-2 1 6
-2 2 7
-2 4 1
-2 6 3
-2 7 3
-2 9 4
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 14
+
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 3 2 1 3 2 
+0 0 0 
+
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
+
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 4 1
-2 1 6
-2 6 3
-2 2 7
-2 7 3
-2 0 9
-2 9 4
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 14
+
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 1 3 3 2 2 
+0 0 0 
+
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
+
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 4 double
 0.00000 0.00000 0
 1.00000 0.00000 0
-0.00000 1.00000 0
-1.00000 1.00000 0
 0.500000 0.00000 0
-0.00000 0.500000 0
-1.00000 0.500000 0
-0.500000 1.00000 0
-0.500000 0.500000 0
 0.250000 0.00000 0
-0.00000 0.250000 0
-0.500000 0.250000 0
-0.250000 0.500000 0
-0.250000 0.250000 0
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 0 9
-2 1 6
-2 2 7
-2 4 1
-2 6 3
-2 7 3
-2 9 4
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 14
+
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 3 2 1 3 2 
+0 0 0 
+
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
+
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
-0.00000 0.00000 0.00000
-1.00000 0.00000 0.00000
-0.00000 1.00000 0.00000
-1.00000 1.00000 0.00000
-0.500000 0.00000 0.00000
-0.00000 0.500000 0.00000
-1.00000 0.500000 0.00000
-0.500000 1.00000 0.00000
-0.500000 0.500000 0.00000
-0.250000 0.00000 0.00000
-0.00000 0.250000 0.00000
-0.500000 0.250000 0.00000
-0.250000 0.500000 0.00000
-0.250000 0.250000 0.00000
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
-0.00000 0.00000 0.00000
-1.00000 0.00000 0.00000
-0.00000 1.00000 0.00000
-1.00000 1.00000 0.00000
-0.500000 0.00000 0.00000
-0.00000 0.500000 0.00000
-1.00000 0.500000 0.00000
-0.500000 1.00000 0.00000
-0.500000 0.500000 0.00000
-0.250000 0.00000 0.00000
-0.00000 0.250000 0.00000
-0.500000 0.250000 0.00000
-0.250000 0.500000 0.00000
-0.250000 0.250000 0.00000
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 
 
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
-0.00000 0.00000 0.00000
-1.00000 0.00000 0.00000
-0.00000 1.00000 0.00000
-1.00000 1.00000 0.00000
-0.500000 0.00000 0.00000
-0.00000 0.500000 0.00000
-1.00000 0.500000 0.00000
-0.500000 1.00000 0.00000
-0.500000 0.500000 0.00000
-0.250000 0.00000 0.00000
-0.00000 0.250000 0.00000
-0.500000 0.250000 0.00000
-0.250000 0.500000 0.00000
-0.250000 0.250000 0.00000
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 3
+3 3 3 
 
 
-CELL_DATA 7
+CELL_DATA 3
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
+0 0 0 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
+-1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
+
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
+
+CELL_TYPES 3
+3 3 3 
+
+
+
+CELL_DATA 3
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
+
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
+
+CELL_TYPES 3
+3 3 3 
+
+
+
+CELL_DATA 3
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
+
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
+
+CELL_TYPES 3
+3 3 3 
+
+
+
+CELL_DATA 3
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
+
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
+
+CELL_TYPES 3
+3 3 3 
+
+
+
+CELL_DATA 3
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 4 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.500000 0.00000 0
+0.250000 0.00000 0
+
+CELLS 3 9
+2 2 1
+2 0 3
+2 3 2
+
+CELL_TYPES 3
+3 3 3 
+
+
+
+CELL_DATA 3
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 1 3 3 2 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 1 3 3 2 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 17 65
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 5 2
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+2 0 10
+2 10 5
+
+CELL_TYPES 17
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 
+
+
+CELL_DATA 17
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 0 1 1 3 3 2 2 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 17 65
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 5 2
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+2 0 10
+2 10 5
+
+CELL_TYPES 17
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 
+
+
+CELL_DATA 17
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 0 1 1 3 3 2 2 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0
+1.00000 0.00000 0
+0.00000 1.00000 0
+1.00000 1.00000 0
+0.500000 0.00000 0
+0.00000 0.500000 0
+1.00000 0.500000 0
+0.500000 1.00000 0
+0.500000 0.500000 0
+0.250000 0.00000 0
+0.00000 0.250000 0
+0.500000 0.250000 0
+0.250000 0.500000 0
+0.250000 0.250000 0
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 1 3 3 2 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 1 3 3 2 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 7 35
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+
+CELL_TYPES 7
+9 9 9 9 9 9 9 
+
+
+
+CELL_DATA 7
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 17 65
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 5 2
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+2 0 10
+2 10 5
+
+CELL_TYPES 17
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 
+
+
+CELL_DATA 17
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 0 1 1 3 3 2 2 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 17 65
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 4 1
+2 5 2
+2 1 6
+2 6 3
+2 2 7
+2 7 3
+2 0 9
+2 9 4
+2 0 10
+2 10 5
+
+CELL_TYPES 17
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 
+
+
+CELL_DATA 17
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 0 1 1 3 3 2 2 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 14 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+1.00000 0.500000 0.00000
+0.500000 1.00000 0.00000
+0.500000 0.500000 0.00000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.250000 0.250000 0.00000
+
+CELLS 14 56
+4 4 1 6 8
+4 5 8 7 2
+4 8 6 3 7
+4 0 9 13 10
+4 9 4 11 13
+4 10 13 12 5
+4 13 11 8 12
+2 0 9
+2 1 6
+2 2 7
+2 4 1
+2 6 3
+2 7 3
+2 9 4
+
+CELL_TYPES 14
+9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 
+
+
+CELL_DATA 14
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 
+2 1 3 2 1 3 2 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 
 -1 -1 -1 -1 -1 -1 -1 
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 46 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
+1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
+0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
+0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
+0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
+
+CELLS 15 135
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+
+CELL_TYPES 15
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+
+
+CELL_DATA 15
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 46 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
+1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
+0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
+0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
+0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
+
+CELLS 15 135
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+
+CELL_TYPES 15
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+
+
+
+CELL_DATA 15
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+
+
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+
+
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 46 double
 0.00000 0.00000 0.00000
 1.00000 0.00000 0.00000
 0.00000 1.00000 0.00000
 1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
 0.500000 0.00000 0.00000
 0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
 1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
 0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
 0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
 0.250000 0.00000 0.00000
 0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
 0.500000 0.250000 0.00000
 0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
 0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
 
-CELLS 7 35
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
+CELLS 47 231
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+2 20 16
+2 20 12
+2 21 11
+2 21 13
+2 11 23
+2 12 23
+2 23 15
+2 23 18
+2 14 24
+2 13 24
+2 24 19
+2 24 15
+2 16 25
+2 17 25
+2 25 18
+2 25 19
+2 29 39
+2 27 39
+2 10 30
+2 39 30
+2 39 31
+2 8 31
+2 30 20
+2 31 20
+2 27 40
+2 28 40
+2 8 32
+2 40 32
+2 40 33
+2 9 33
+2 32 21
+2 33 21
 
-CELL_TYPES 7
-9 9 9 9 9 9 9 
+CELL_TYPES 47
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
+
+CELL_DATA 47
+SCALARS MaterialID int 1
+LOOKUP_TABLE default
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+2 2 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 2 2 2 2 4 4 4 4 4 4 4 4 
+
+SCALARS ManifoldID int 1
+LOOKUP_TABLE default
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+# vtk DataFile Version 3.0
+Triangulation generated with deal.II
+ASCII
+DATASET UNSTRUCTURED_GRID
+POINTS 46 double
+0.00000 0.00000 0.00000
+1.00000 0.00000 0.00000
+0.00000 1.00000 0.00000
+1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
+0.500000 0.00000 0.00000
+0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
+1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
+0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
+0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
+0.250000 0.00000 0.00000
+0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
+0.500000 0.250000 0.00000
+0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
+0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
+
+CELLS 47 231
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+2 29 39
+2 27 39
+2 27 40
+2 28 40
+2 11 23
+2 12 23
+2 14 24
+2 13 24
+2 16 25
+2 17 25
+2 20 12
+2 21 11
+2 8 31
+2 8 32
+2 21 13
+2 9 33
+2 20 16
+2 10 30
+2 23 15
+2 23 18
+2 24 15
+2 24 19
+2 25 18
+2 25 19
+2 39 31
+2 40 32
+2 40 33
+2 39 30
+2 30 20
+2 31 20
+2 32 21
+2 33 21
 
+CELL_TYPES 47
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
 
+3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
 
-CELL_DATA 7
+CELL_DATA 47
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
 
+2 2 4 4 1 1 3 3 5 5 2 4 2 4 4 4 2 2 1 1 3 3 5 5 2 4 4 2 2 2 4 4 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
-
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 46 double
 0.00000 0.00000 0.00000
 1.00000 0.00000 0.00000
 0.00000 1.00000 0.00000
 1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
 0.500000 0.00000 0.00000
 0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
 1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
 0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
 0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
 0.250000 0.00000 0.00000
 0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
 0.500000 0.250000 0.00000
 0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
 0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 4 1
-2 1 6
-2 6 3
-2 2 7
-2 7 3
-2 0 9
-2 9 4
+CELLS 41 265
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+4 10 4 16 20
+4 8 20 12 1
+4 20 16 5 12
+4 8 1 11 21
+4 9 21 13 2
+4 21 11 3 13
+4 1 11 23 12
+4 11 3 15 23
+4 12 23 18 5
+4 23 15 7 18
+4 2 14 24 13
+4 14 6 19 24
+4 13 24 15 3
+4 24 19 7 15
+4 4 16 25 17
+4 16 5 18 25
+4 17 25 19 6
+4 25 18 7 19
+4 0 29 39 27
+4 29 10 30 39
+4 27 39 31 8
+4 39 30 20 31
+4 0 27 40 28
+4 27 8 32 40
+4 28 40 33 9
+4 40 32 21 33
 
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELL_TYPES 41
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
 
 
-CELL_DATA 14
+CELL_DATA 41
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 1 3 3 2 2 
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+2 2 2 4 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 4 4 4 4 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 46 double
 0.00000 0.00000 0.00000
 1.00000 0.00000 0.00000
 0.00000 1.00000 0.00000
 1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
 0.500000 0.00000 0.00000
 0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
 1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
 0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
 0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
 0.250000 0.00000 0.00000
 0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
 0.500000 0.250000 0.00000
 0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
 0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 0 9
-2 1 6
-2 2 7
-2 4 1
-2 6 3
-2 7 3
-2 9 4
+CELLS 41 265
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+4 0 29 39 27
+4 0 27 40 28
+4 1 11 23 12
+4 2 14 24 13
+4 4 16 25 17
+4 8 20 12 1
+4 8 1 11 21
+4 9 21 13 2
+4 10 4 16 20
+4 11 3 15 23
+4 12 23 18 5
+4 13 24 15 3
+4 14 6 19 24
+4 16 5 18 25
+4 17 25 19 6
+4 20 16 5 12
+4 21 11 3 13
+4 23 15 7 18
+4 24 19 7 15
+4 25 18 7 19
+4 27 39 31 8
+4 27 8 32 40
+4 28 40 33 9
+4 29 10 30 39
+4 39 30 20 31
+4 40 32 21 33
 
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELL_TYPES 41
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
 
 
-CELL_DATA 14
+CELL_DATA 41
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 3 2 1 3 2 
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+2 4 1 3 5 2 4 4 2 1 1 3 3 5 5 2 4 1 3 5 2 4 4 2 2 4 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 46 double
 0.00000 0.00000 0.00000
 1.00000 0.00000 0.00000
 0.00000 1.00000 0.00000
 1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
 0.500000 0.00000 0.00000
 0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
 1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
 0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
 0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
 0.250000 0.00000 0.00000
 0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
 0.500000 0.250000 0.00000
 0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
 0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 4 1
-2 1 6
-2 6 3
-2 2 7
-2 7 3
-2 0 9
-2 9 4
-
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELLS 73 361
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+4 10 4 16 20
+4 8 20 12 1
+4 20 16 5 12
+4 8 1 11 21
+4 9 21 13 2
+4 21 11 3 13
+4 1 11 23 12
+4 11 3 15 23
+4 12 23 18 5
+4 23 15 7 18
+4 2 14 24 13
+4 14 6 19 24
+4 13 24 15 3
+4 24 19 7 15
+4 4 16 25 17
+4 16 5 18 25
+4 17 25 19 6
+4 25 18 7 19
+4 0 29 39 27
+4 29 10 30 39
+4 27 39 31 8
+4 39 30 20 31
+4 0 27 40 28
+4 27 8 32 40
+4 28 40 33 9
+4 40 32 21 33
+2 20 16
+2 20 12
+2 21 11
+2 21 13
+2 11 23
+2 12 23
+2 23 15
+2 23 18
+2 14 24
+2 13 24
+2 24 19
+2 24 15
+2 16 25
+2 17 25
+2 25 18
+2 25 19
+2 29 39
+2 27 39
+2 10 30
+2 39 30
+2 39 31
+2 8 31
+2 30 20
+2 31 20
+2 27 40
+2 28 40
+2 8 32
+2 40 32
+2 40 33
+2 9 33
+2 32 21
+2 33 21
 
+CELL_TYPES 73
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
 
-CELL_DATA 14
+CELL_DATA 73
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 1 3 3 2 2 
-
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+2 2 2 4 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 4 4 4 4 
+2 2 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 2 2 2 2 4 4 4 4 4 4 4 4 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
-
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
 DATASET UNSTRUCTURED_GRID
-POINTS 14 double
+POINTS 46 double
 0.00000 0.00000 0.00000
 1.00000 0.00000 0.00000
 0.00000 1.00000 0.00000
 1.00000 1.00000 0.00000
+0.00000 0.00000 1.00000
+1.00000 0.00000 1.00000
+0.00000 1.00000 1.00000
+1.00000 1.00000 1.00000
 0.500000 0.00000 0.00000
 0.00000 0.500000 0.00000
+0.00000 0.00000 0.500000
 1.00000 0.500000 0.00000
+1.00000 0.00000 0.500000
 0.500000 1.00000 0.00000
+0.00000 1.00000 0.500000
+1.00000 1.00000 0.500000
+0.500000 0.00000 1.00000
+0.00000 0.500000 1.00000
+1.00000 0.500000 1.00000
+0.500000 1.00000 1.00000
+0.500000 0.00000 0.500000
 0.500000 0.500000 0.00000
+0.00000 0.500000 0.500000
+1.00000 0.500000 0.500000
+0.500000 1.00000 0.500000
+0.500000 0.500000 1.00000
+0.500000 0.500000 0.500000
 0.250000 0.00000 0.00000
 0.00000 0.250000 0.00000
+0.00000 0.00000 0.250000
+0.250000 0.00000 0.500000
+0.500000 0.00000 0.250000
 0.500000 0.250000 0.00000
 0.250000 0.500000 0.00000
+0.00000 0.500000 0.250000
+0.00000 0.250000 0.500000
+0.500000 0.500000 0.250000
+0.250000 0.500000 0.500000
+0.500000 0.250000 0.500000
+0.250000 0.00000 0.250000
 0.250000 0.250000 0.00000
+0.00000 0.250000 0.250000
+0.250000 0.250000 0.500000
+0.250000 0.500000 0.250000
+0.500000 0.250000 0.250000
+0.250000 0.250000 0.250000
 
-CELLS 14 56
-4 4 1 6 8
-4 5 8 7 2
-4 8 6 3 7
-4 0 9 13 10
-4 9 4 11 13
-4 10 13 12 5
-4 13 11 8 12
-2 0 9
-2 1 6
-2 2 7
-2 4 1
-2 6 3
-2 7 3
-2 9 4
-
-CELL_TYPES 14
-9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 
+CELLS 73 361
+8 8 1 12 20 21 11 23 26
+8 9 21 26 22 2 13 24 14
+8 21 11 23 26 13 3 15 24
+8 10 20 16 4 22 26 25 17
+8 20 12 5 16 26 23 18 25
+8 22 26 25 17 14 24 19 6
+8 26 23 18 25 24 15 7 19
+8 0 27 39 29 28 40 45 41
+8 27 8 31 39 40 32 44 45
+8 28 40 45 41 9 33 43 34
+8 40 32 44 45 33 21 36 43
+8 29 39 30 10 41 45 42 35
+8 39 31 20 30 45 44 38 42
+8 41 45 42 35 34 43 37 22
+8 45 44 38 42 43 36 26 37
+4 0 29 39 27
+4 0 27 40 28
+4 1 11 23 12
+4 2 14 24 13
+4 4 16 25 17
+4 8 20 12 1
+4 8 1 11 21
+4 9 21 13 2
+4 10 4 16 20
+4 11 3 15 23
+4 12 23 18 5
+4 13 24 15 3
+4 14 6 19 24
+4 16 5 18 25
+4 17 25 19 6
+4 20 16 5 12
+4 21 11 3 13
+4 23 15 7 18
+4 24 19 7 15
+4 25 18 7 19
+4 27 39 31 8
+4 27 8 32 40
+4 28 40 33 9
+4 29 10 30 39
+4 39 30 20 31
+4 40 32 21 33
+2 29 39
+2 27 39
+2 27 40
+2 28 40
+2 11 23
+2 12 23
+2 14 24
+2 13 24
+2 16 25
+2 17 25
+2 20 12
+2 21 11
+2 8 31
+2 8 32
+2 21 13
+2 9 33
+2 20 16
+2 10 30
+2 23 15
+2 23 18
+2 24 15
+2 24 19
+2 25 18
+2 25 19
+2 39 31
+2 40 32
+2 40 33
+2 39 30
+2 30 20
+2 31 20
+2 32 21
+2 33 21
 
+CELL_TYPES 73
+12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
+9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
 
-CELL_DATA 14
+CELL_DATA 73
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
-0 0 0 0 0 0 0 
-2 1 3 2 1 3 2 
-
+0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
+2 4 1 3 5 2 4 4 2 1 1 3 3 5 5 2 4 1 3 5 2 4 4 2 2 4 
+2 2 4 4 1 1 3 3 5 5 2 4 2 4 4 4 2 2 1 1 3 3 5 5 2 4 4 2 2 2 4 4 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
--1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 
-
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
@@ -1491,7 +3687,7 @@ POINTS 46 double
 0.500000 0.250000 0.250000
 0.250000 0.250000 0.250000
 
-CELLS 47 231
+CELLS 84 342
 8 8 1 12 20 21 11 23 26
 8 9 21 26 22 2 13 24 14
 8 21 11 23 26 13 3 15 24
@@ -1507,53 +3703,90 @@ CELLS 47 231
 8 39 31 20 30 45 44 38 42
 8 41 45 42 35 34 43 37 22
 8 45 44 38 42 43 36 26 37
+2 4 16
+2 10 4
 2 20 16
+2 8 1
 2 20 12
+2 1 12
+2 16 5
+2 12 5
+2 1 11
 2 21 11
+2 9 2
 2 21 13
+2 2 13
+2 11 3
+2 13 3
+2 2 14
+2 22 14
+2 22 17
+2 4 17
+2 14 6
+2 17 6
 2 11 23
 2 12 23
+2 3 15
 2 23 15
 2 23 18
+2 5 18
+2 15 7
+2 18 7
 2 14 24
 2 13 24
+2 6 19
 2 24 19
 2 24 15
+2 19 7
 2 16 25
 2 17 25
 2 25 18
 2 25 19
+2 0 27
 2 29 39
+2 0 29
 2 27 39
 2 10 30
+2 29 10
 2 39 30
+2 27 8
 2 39 31
 2 8 31
 2 30 20
 2 31 20
+2 0 28
 2 27 40
 2 28 40
 2 8 32
 2 40 32
+2 28 9
 2 40 33
 2 9 33
 2 32 21
 2 33 21
-
-CELL_TYPES 47
+2 28 41
+2 29 41
+2 9 34
+2 41 34
+2 41 35
+2 10 35
+2 34 22
+2 35 22
+
+CELL_TYPES 84
 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
-3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
+3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
 
-CELL_DATA 47
+CELL_DATA 84
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
-2 2 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 2 2 2 2 4 4 4 4 4 4 4 4 
+0 0 2 0 2 0 0 0 0 4 0 4 0 0 0 0 0 0 0 0 0 1 1 0 1 1 0 0 0 3 3 0 3 3 0 5 5 5 5 0 2 0 2 2 0 2 0 2 2 2 2 0 4 4 4 4 0 4 4 4 4 0 0 0 0 0 0 0 0 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
@@ -1724,7 +3957,7 @@ POINTS 46 double
 0.500000 0.250000 0.250000
 0.250000 0.250000 0.250000
 
-CELLS 41 265
+CELLS 48 300
 8 8 1 12 20 21 11 23 26
 8 9 21 26 22 2 13 24 14
 8 21 11 23 26 13 3 15 24
@@ -1746,6 +3979,9 @@ CELLS 41 265
 4 8 1 11 21
 4 9 21 13 2
 4 21 11 3 13
+4 9 2 14 22
+4 10 22 17 4
+4 22 14 6 17
 4 1 11 23 12
 4 11 3 15 23
 4 12 23 18 5
@@ -1766,23 +4002,27 @@ CELLS 41 265
 4 27 8 32 40
 4 28 40 33 9
 4 40 32 21 33
+4 0 28 41 29
+4 28 9 34 41
+4 29 41 35 10
+4 41 34 22 35
 
-CELL_TYPES 41
+CELL_TYPES 48
 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
-9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
+9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
 
 
-CELL_DATA 41
+CELL_DATA 48
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
-2 2 2 4 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 4 4 4 4 
+2 2 2 4 4 4 0 0 0 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 4 4 4 4 0 0 0 0 
 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII
@@ -1947,7 +4187,7 @@ POINTS 46 double
 0.500000 0.250000 0.250000
 0.250000 0.250000 0.250000
 
-CELLS 73 361
+CELLS 117 507
 8 8 1 12 20 21 11 23 26
 8 9 21 26 22 2 13 24 14
 8 21 11 23 26 13 3 15 24
@@ -1969,6 +4209,9 @@ CELLS 73 361
 4 8 1 11 21
 4 9 21 13 2
 4 21 11 3 13
+4 9 2 14 22
+4 10 22 17 4
+4 22 14 6 17
 4 1 11 23 12
 4 11 3 15 23
 4 12 23 18 5
@@ -1989,56 +4232,97 @@ CELLS 73 361
 4 27 8 32 40
 4 28 40 33 9
 4 40 32 21 33
+4 0 28 41 29
+4 28 9 34 41
+4 29 41 35 10
+4 41 34 22 35
+2 4 16
+2 10 4
 2 20 16
+2 8 1
 2 20 12
+2 1 12
+2 16 5
+2 12 5
+2 1 11
 2 21 11
+2 9 2
 2 21 13
+2 2 13
+2 11 3
+2 13 3
+2 2 14
+2 22 14
+2 22 17
+2 4 17
+2 14 6
+2 17 6
 2 11 23
 2 12 23
+2 3 15
 2 23 15
 2 23 18
+2 5 18
+2 15 7
+2 18 7
 2 14 24
 2 13 24
+2 6 19
 2 24 19
 2 24 15
+2 19 7
 2 16 25
 2 17 25
 2 25 18
 2 25 19
+2 0 27
 2 29 39
+2 0 29
 2 27 39
 2 10 30
+2 29 10
 2 39 30
+2 27 8
 2 39 31
 2 8 31
 2 30 20
 2 31 20
+2 0 28
 2 27 40
 2 28 40
 2 8 32
 2 40 32
+2 28 9
 2 40 33
 2 9 33
 2 32 21
 2 33 21
-
-CELL_TYPES 73
+2 28 41
+2 29 41
+2 9 34
+2 41 34
+2 41 35
+2 10 35
+2 34 22
+2 35 22
+
+CELL_TYPES 117
 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 
-9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
-3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
+9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 9 
+3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 
 
-CELL_DATA 73
+CELL_DATA 117
 SCALARS MaterialID int 1
 LOOKUP_TABLE default
 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
-2 2 2 4 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 4 4 4 4 
-2 2 4 4 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 2 2 2 2 4 4 4 4 4 4 4 4 
+2 2 2 4 4 4 0 0 0 1 1 1 1 3 3 3 3 5 5 5 5 2 2 2 2 4 4 4 4 0 0 0 0 
+0 0 2 0 2 0 0 0 0 4 0 4 0 0 0 0 0 0 0 0 0 1 1 0 1 1 0 0 0 3 3 0 3 3 0 5 5 5 5 0 2 0 2 2 0 2 0 2 2 2 2 0 4 4 4 4 0 4 4 4 4 0 0 0 0 0 0 0 0 
 
 SCALARS ManifoldID int 1
 LOOKUP_TABLE default
 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
--1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
+-1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 
 # vtk DataFile Version 3.0
 Triangulation generated with deal.II
 ASCII