Fix tests.

Co-authored-by: David Wells <drwells@vt.edu>
Co-authored-by: Jean-Paul Pelteret <jppelteret@gmail.com>
Co-authored-by: Daniel Arndt <daniel.arndt@iwr.uni-heidelberg.de>

diff --git a/tests/aniso/mesh_3d_21.cc b/tests/aniso/mesh_3d_21.cc
index 7db7498f5f41d4e1a052e19f1ff30aefbbd8f8f3..982f197b09c5413f37829139cc4956f2a5ea1029 100644 (file)
--- a/tests/aniso/mesh_3d_21.cc
+++ b/tests/aniso/mesh_3d_21.cc
@@ -203,6 +203,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_manifold(0);
     check (coarse_grid);
   }
 

diff --git a/tests/aniso/refine_and_coarsen_anisotropic.cc b/tests/aniso/refine_and_coarsen_anisotropic.cc
index 350d20ce897ed3e287b0384ef68bde88d394d22c..d2b3057d17e8cb07cddaa25e8ab3a4bfbd747873 100644 (file)
--- a/tests/aniso/refine_and_coarsen_anisotropic.cc
+++ b/tests/aniso/refine_and_coarsen_anisotropic.cc
@@ -67,6 +67,7 @@ void test_isotropic (int type, std::ostream *logfile)
   else
     GridGenerator::hyper_ball (tria);
 
+  tria.reset_manifold(0);
   tria.refine_global(1);
 
   *logfile << "cycle: 0, number of cells: "<<tria.n_cells()

diff --git a/tests/bits/face_orientation_and_fe_nedelec_01.cc b/tests/bits/face_orientation_and_fe_nedelec_01.cc
index 08c0618a512967078f11a26fd1e31876387a1236..2e36d192408bfe6d832f33a0d7a571e2424362c7 100644 (file)
--- a/tests/bits/face_orientation_and_fe_nedelec_01.cc
+++ b/tests/bits/face_orientation_and_fe_nedelec_01.cc
@@ -178,6 +178,7 @@ void test_with_wrong_face_orientation (const FiniteElement<dim> &fe,
     {
       Triangulation<dim>     triangulation;
       GridGenerator::hyper_ball (triangulation);
+      triangulation.reset_manifold(0);
       typename Triangulation<dim>::active_cell_iterator
       cell = triangulation.begin_active();
       std::advance (cell, i);

diff --git a/tests/bits/face_orientation_and_fe_q_01.cc b/tests/bits/face_orientation_and_fe_q_01.cc
index 993b9dca50f6d92080f9f271511a1bb695ce174e..3fe5a2d37fa363ae00e6794a505273b595ba9e80 100644 (file)
--- a/tests/bits/face_orientation_and_fe_q_01.cc
+++ b/tests/bits/face_orientation_and_fe_q_01.cc
@@ -181,6 +181,7 @@ void test_with_wrong_face_orientation (const FiniteElement<dim> &fe,
     {
       Triangulation<dim>     triangulation;
       GridGenerator::hyper_ball (triangulation);
+      triangulation.reset_manifold(0);
       typename Triangulation<dim>::active_cell_iterator
       cell = triangulation.begin_active();
       std::advance (cell, i);

diff --git a/tests/bits/find_cell_3.cc b/tests/bits/find_cell_3.cc
index 90e750ed721e4ff1efd17669c00e392faf604eb9..ea35d2cb415e4dad5f32478b9340d4719c33bf5a 100644 (file)
--- a/tests/bits/find_cell_3.cc
+++ b/tests/bits/find_cell_3.cc
@@ -67,6 +67,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_all_manifolds();
     coarse_grid.refine_global (1);
     check (coarse_grid);
   }

diff --git a/tests/bits/find_cell_alt_3.cc b/tests/bits/find_cell_alt_3.cc
index 00538d6acb85556351b057415fb20375bb6ab6b1..a3dc77c5423a3eede4e14d484ab2a54fd47ef595 100644 (file)
--- a/tests/bits/find_cell_alt_3.cc
+++ b/tests/bits/find_cell_alt_3.cc
@@ -72,6 +72,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_all_manifolds();
     coarse_grid.refine_global (1);
     check (coarse_grid);
   }

diff --git a/tests/bits/step-2.cc b/tests/bits/step-2.cc
index 877aa351f9421f6ca57c6878e7764947d853485a..b8bc7317a6f327a6455bddd87511788850fb3a96 100644 (file)
--- a/tests/bits/step-2.cc
+++ b/tests/bits/step-2.cc
@@ -22,6 +22,7 @@
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/hp/dof_handler.h>
@@ -44,6 +45,8 @@ void make_grid (Triangulation<2> &triangulation)
   GridGenerator::hyper_shell (triangulation,
                               center, inner_radius, outer_radius,
                               10);
+  //triangulation.reset_all_manifolds();
+  GridTools::copy_boundary_to_manifold_id(triangulation);
 
   static const SphericalManifold<2> boundary_description(center);
   triangulation.set_manifold (0, boundary_description);

diff --git a/tests/bits/step-5.cc b/tests/bits/step-5.cc
index b88afc71813f9b626862ba1f7c2635f472b8d9eb..a6b777f29c538842e5ceab1f06e002d4fbe737dc 100644 (file)
--- a/tests/bits/step-5.cc
+++ b/tests/bits/step-5.cc
@@ -32,6 +32,7 @@
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/fe/fe_q.h>
@@ -298,6 +299,7 @@ void LaplaceProblem<dim>::run ()
           Assert (dim==2, ExcInternalError());
 
           grid_in.read_ucd (input_file);
+          GridTools::copy_boundary_to_manifold_id(triangulation);
 
           static const SphericalManifold<dim> boundary;
           triangulation.set_manifold (0, boundary);

diff --git a/tests/codim_one/bem.cc b/tests/codim_one/bem.cc
index f3d95e137d0eb0ff33b881f963418d8f3d908de5..fb281683b52a822aa1f95b0cc1d7906fd188c6c0 100644 (file)
--- a/tests/codim_one/bem.cc
+++ b/tests/codim_one/bem.cc
@@ -35,6 +35,7 @@
 #include <deal.II/fe/mapping_q1.h>
 #include <deal.II/grid/grid_in.h>
 #include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/manifold_lib.h>
@@ -128,7 +129,6 @@ BEM<spacedim>::run()
       // works up to cycle<9, but for testin purpose, we stop at 4
       for (unsigned int cycle=0; cycle<4; ++cycle)
         {
-
           tria.set_manifold(1, boundary);
           tria.refine_global(1);
 
@@ -465,6 +465,9 @@ BEM<spacedim>::read_grid(std::string filename)
   gi.attach_triangulation (tria);
   std::ifstream in (filename.c_str());
   gi.read_ucd (in);
+
+  for (auto cell = tria.begin_active(); cell != tria.end(); ++cell)
+    cell->set_all_manifold_ids(1);
 }
 
 

diff --git a/tests/codim_one/boundary_indicator_01.cc b/tests/codim_one/boundary_indicator_01.cc
index 5b2602ca2f637191a4efed98e3c18cea9d42ca3e..0fac459f16e8e01e488f97da448d2f6c010cb8e6 100644 (file)
--- a/tests/codim_one/boundary_indicator_01.cc
+++ b/tests/codim_one/boundary_indicator_01.cc
@@ -60,6 +60,7 @@ int main ()
     const SphericalManifold<dim> boundary_description;
     Triangulation<dim> volume_mesh;
     GridGenerator::hyper_ball(volume_mesh);
+    GridTools::copy_boundary_to_manifold_id(volume_mesh);
     volume_mesh.set_manifold (0, boundary_description);
 
     // exclude one of the 6 faces
@@ -107,6 +108,7 @@ int main ()
       for (unsigned int f=0; f<GeometryInfo<dim-1>::faces_per_cell; ++f)
         if (cell->at_boundary(f))
           cell->face(f)->set_boundary_id(1);
+    GridTools::copy_boundary_to_manifold_id(boundary_mesh);
 
     boundary_mesh.refine_global (2);
 

diff --git a/tests/codim_one/extract_boundary_mesh_01.cc b/tests/codim_one/extract_boundary_mesh_01.cc
index 18d89cbe64f5a62b1c5f74443ba0edade12bea8b..3a903ed6359de15b714bd1a77b9830f98fccded8 100644 (file)
--- a/tests/codim_one/extract_boundary_mesh_01.cc
+++ b/tests/codim_one/extract_boundary_mesh_01.cc
@@ -102,13 +102,14 @@ int main ()
     // Extract the boundary of a hyper-sphere
 
     const int dim = 3;
-    deallog << "Testing hyper_cube in dim: " << dim << "..."<< endl;
+    deallog << "Testing hyper_ball in dim: " << dim << "..."<< endl;
 
     map< Triangulation<dim-1,dim>::cell_iterator,
          Triangulation<dim,dim>::face_iterator>
          surface_to_volume_mapping;
     Triangulation<dim> volume_mesh;
     GridGenerator::hyper_ball(volume_mesh);
+    volume_mesh.reset_manifold (0);
 
     volume_mesh.refine_global (1);
 

diff --git a/tests/codim_one/extract_boundary_mesh_01.output b/tests/codim_one/extract_boundary_mesh_01.output
index 13519be2dce0a3ef2cde60acbe48eb04391a8597..5a040e28c02755d6a9d9f51c2f78eb87577a8219 100644 (file)
--- a/tests/codim_one/extract_boundary_mesh_01.output
+++ b/tests/codim_one/extract_boundary_mesh_01.output
@@ -1,5 +1,5 @@
 
-DEAL::Testing hyper_cube in dim: 3...
+DEAL::Testing hyper_ball in dim: 3...
 DEAL::The last column should be 0
 DEAL::Vol faces                Surf cell               Distance
 DEAL::

diff --git a/tests/codim_one/extract_boundary_mesh_03.cc b/tests/codim_one/extract_boundary_mesh_03.cc
index bb0a0e7f8eb0de8e1cae9f619e063b7ca6cd8ca6..4f2e8765240abaf9c8fdc45ae496599c1d08c3f5 100644 (file)
--- a/tests/codim_one/extract_boundary_mesh_03.cc
+++ b/tests/codim_one/extract_boundary_mesh_03.cc
@@ -97,7 +97,7 @@ int main ()
     // Extract the boundary of a hyper-sphere
 
     const int dim = 3;
-    deallog << "Testing hyper_cube in dim: " << dim << "..."<< endl;
+    deallog << "Testing hyper_ball in dim: " << dim << "..."<< endl;
 
     map< Triangulation<dim-1,dim>::cell_iterator,
          Triangulation<dim,dim>::face_iterator>
@@ -111,6 +111,7 @@ int main ()
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         if (cell->at_boundary(f))
           cell->face(f)->set_all_boundary_ids (1);
+    GridTools::copy_boundary_to_manifold_id(volume_mesh);
     volume_mesh.set_manifold (1, boundary_description);
     volume_mesh.refine_global (1);
 

diff --git a/tests/codim_one/extract_boundary_mesh_03.output b/tests/codim_one/extract_boundary_mesh_03.output
index 70d72f9f42d855ab2fdedf79a09ef4e70f769d9b..3f9e2c745f88e399dbe537f7419c0aa5a070a618 100644 (file)
--- a/tests/codim_one/extract_boundary_mesh_03.output
+++ b/tests/codim_one/extract_boundary_mesh_03.output
@@ -1,5 +1,5 @@
 
-DEAL::Testing hyper_cube in dim: 3...
+DEAL::Testing hyper_ball in dim: 3...
 DEAL::56
 DEAL::24
 -0.577350 -0.577350 -0.577350 1 1

diff --git a/tests/codim_one/extract_boundary_mesh_10.cc b/tests/codim_one/extract_boundary_mesh_10.cc
index 45623c267e0441aaefe7bae60c10a5878252f1da..acf880c24393a533664f25fe2ed9d436a3347e28 100644 (file)
--- a/tests/codim_one/extract_boundary_mesh_10.cc
+++ b/tests/codim_one/extract_boundary_mesh_10.cc
@@ -38,13 +38,18 @@
 #include <deal.II/grid/tria_boundary_lib.h>
 
 
-
 void test()
 {
   const int dim=3;
 
   Triangulation<dim>   triangulation;
   GridGenerator::cylinder(triangulation, 100, 200);
+  GridTools::copy_boundary_to_manifold_id(triangulation);
+  for (auto cell = triangulation.begin_active(); cell != triangulation.end(); ++cell)
+    for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
+      if (cell->face(f)->boundary_id() != numbers::invalid_boundary_id)
+        for (unsigned int e=0; e<GeometryInfo<dim>::lines_per_face; ++e)
+          cell->face(f)->line(e)->set_manifold_id(cell->face(f)->boundary_id());
 
   static const CylinderBoundary<dim> outer_cylinder (100,0);
   triangulation.set_manifold(0,outer_cylinder);
@@ -62,6 +67,9 @@ void test()
   triangulation_surface.refine_global (1);
   GridOut().write_gnuplot(triangulation_surface, deallog.get_file_stream());
 
+  std::ofstream output ("grid.vtk");
+  GridOut().write_vtk(triangulation_surface, output);
+
   deallog << triangulation_surface.n_used_vertices() << std::endl;
   deallog << triangulation_surface.n_active_cells() << std::endl;
 }

diff --git a/tests/codim_one/grid_refinement.cc b/tests/codim_one/grid_refinement.cc
index 7da831efa24900fa62ae5d18188365af7f20cdd0..eaaa2ac2eac9c8fed002c3740e1702dfb821ad24 100644 (file)
--- a/tests/codim_one/grid_refinement.cc
+++ b/tests/codim_one/grid_refinement.cc
@@ -24,6 +24,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/grid_in.h>
 #include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/manifold_lib.h>
 
 #include <string>
@@ -35,11 +36,12 @@ void test(std::string filename)
 {
   SphericalManifold<dim, spacedim> boundary;
   Triangulation<dim, spacedim> tria;
-  tria.set_manifold(1, boundary);
   GridIn<dim, spacedim> gi;
   gi.attach_triangulation (tria);
   std::ifstream in (filename.c_str());
   gi.read_ucd (in);
+  tria.set_all_manifold_ids(1);
+  tria.set_manifold(1, boundary);
 
   GridOut grid_out;
   grid_out.set_flags (GridOutFlags::Ucd(true));

diff --git a/tests/dofs/dof_renumbering.cc b/tests/dofs/dof_renumbering.cc
index 12417d7a09794ef4415976b881a8a72c1afd8daa..91c40b4e674932ce4de52146490542f7c1e7763a 100644 (file)
--- a/tests/dofs/dof_renumbering.cc
+++ b/tests/dofs/dof_renumbering.cc
@@ -153,6 +153,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_all_manifolds();
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/fe/cell_similarity_05.cc b/tests/fe/cell_similarity_05.cc
index 332515644edb13f23d2e0c7b735edf2cba0244d4..18d03e129f3d4ac8332b418f96121a2303d9c336 100644 (file)
--- a/tests/fe/cell_similarity_05.cc
+++ b/tests/fe/cell_similarity_05.cc
@@ -34,6 +34,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/manifold_lib.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -155,7 +156,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_06.cc b/tests/fe/cell_similarity_06.cc
index 6f992234349ebdb3a73169b756feeaf6850a8198..751785ff2999d0e24fc18e96eeb81ec923ca9a90 100644 (file)
--- a/tests/fe/cell_similarity_06.cc
+++ b/tests/fe/cell_similarity_06.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_07.cc b/tests/fe/cell_similarity_07.cc
index bcbc0670c8df2d11d9b833859fd9786cf950bdfe..5c05329aeb069f18dbc3ff0cd27605ee5def05b5 100644 (file)
--- a/tests/fe/cell_similarity_07.cc
+++ b/tests/fe/cell_similarity_07.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_08.cc b/tests/fe/cell_similarity_08.cc
index 1e9beadd34251e01a0b61306eef13d5a7e59c1ac..abfc2438e15ab1399c4e90abb7cf21d140e2527b 100644 (file)
--- a/tests/fe/cell_similarity_08.cc
+++ b/tests/fe/cell_similarity_08.cc
@@ -34,6 +34,7 @@
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/manifold_lib.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_q.h>
 #include <deal.II/fe/fe_dgq.h>
@@ -148,6 +149,7 @@ void test()
   std::vector<unsigned int> subdivisions (dim, 1);
   subdivisions[0] = 2;
   GridGenerator::subdivided_hyper_rectangle(tr, subdivisions, p1, p2);
+  GridTools::copy_boundary_to_manifold_id(tr);
 
   static const SphericalManifold<dim> boundary(tr.begin_active()->center());
   tr.set_manifold (1, boundary);
@@ -155,7 +157,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_11.cc b/tests/fe/cell_similarity_11.cc
index eb61db2120c8390399d866fe55dc008fafebe2fa..f25ef7fcbeca095d89d8b3f9842b3438afaeca0a 100644 (file)
--- a/tests/fe/cell_similarity_11.cc
+++ b/tests/fe/cell_similarity_11.cc
@@ -110,7 +110,7 @@ void test()
       if (cell->face(f)->at_boundary() &&
           std::abs(cell->face(f)->center()[0]+1.) < 1e-12 &&
           std::abs(cell->face(f)->center()[1]) < 1e-12)
-        cell->face(f)->set_boundary_id(1);
+        cell->face(f)->set_manifold_id(1);
 
   static const SphericalManifold<dim> boundary;
   tr.set_manifold (1, boundary);

diff --git a/tests/fe/cell_similarity_dgp_monomial_05.cc b/tests/fe/cell_similarity_dgp_monomial_05.cc
index 516fcdc0615e6e074269690051bf08fbee3c0216..06a43a006b30d2036166a0f3e40c5617c5d4ce4b 100644 (file)
--- a/tests/fe/cell_similarity_dgp_monomial_05.cc
+++ b/tests/fe/cell_similarity_dgp_monomial_05.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_monomial_06.cc b/tests/fe/cell_similarity_dgp_monomial_06.cc
index 0fc66c172d4df36e3d090a662da2bf4cfdd15eb6..92ecf3119a10ac11cf243f276c8ff7701d8ba1a8 100644 (file)
--- a/tests/fe/cell_similarity_dgp_monomial_06.cc
+++ b/tests/fe/cell_similarity_dgp_monomial_06.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_monomial_07.cc b/tests/fe/cell_similarity_dgp_monomial_07.cc
index 659f128946303fd1a475299938597e89f349d603..be640e7c517d96baa6b307bf86a1ebb031910816 100644 (file)
--- a/tests/fe/cell_similarity_dgp_monomial_07.cc
+++ b/tests/fe/cell_similarity_dgp_monomial_07.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_monomial_08.cc b/tests/fe/cell_similarity_dgp_monomial_08.cc
index f1ea5f93ca22116e324cf124e22626c94624fb96..4862fe312a7f676939fccbb12fc90934b66779fd 100644 (file)
--- a/tests/fe/cell_similarity_dgp_monomial_08.cc
+++ b/tests/fe/cell_similarity_dgp_monomial_08.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_nonparametric_05.cc b/tests/fe/cell_similarity_dgp_nonparametric_05.cc
index 71863c2a6666e3897502419c7722d31bb3ed6ebc..efc11763d1e9c73d667e4947cedbd0eb0e9b97e3 100644 (file)
--- a/tests/fe/cell_similarity_dgp_nonparametric_05.cc
+++ b/tests/fe/cell_similarity_dgp_nonparametric_05.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_nonparametric_06.cc b/tests/fe/cell_similarity_dgp_nonparametric_06.cc
index 4c327a865f4721445fd9e0b122688697155a49a3..a11b5de3ed9378fab8caba971b45dbd39ba7d269 100644 (file)
--- a/tests/fe/cell_similarity_dgp_nonparametric_06.cc
+++ b/tests/fe/cell_similarity_dgp_nonparametric_06.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_nonparametric_07.cc b/tests/fe/cell_similarity_dgp_nonparametric_07.cc
index 9c5c19aab8134a3e2506114de252ab93d7e21931..2c75b4f904af4d3f3e766af820bb1a5689259390 100644 (file)
--- a/tests/fe/cell_similarity_dgp_nonparametric_07.cc
+++ b/tests/fe/cell_similarity_dgp_nonparametric_07.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/cell_similarity_dgp_nonparametric_08.cc b/tests/fe/cell_similarity_dgp_nonparametric_08.cc
index ee19c1c38bd2be6fd224861f03688348e299f2f7..d7c1d29344b25bc248691ba506275d47f2ebe442 100644 (file)
--- a/tests/fe/cell_similarity_dgp_nonparametric_08.cc
+++ b/tests/fe/cell_similarity_dgp_nonparametric_08.cc
@@ -155,7 +155,7 @@ void test()
   // set boundary id on cell 1
   for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
     if (tr.begin_active()->at_boundary(f))
-      tr.begin_active()->face(f)->set_boundary_id (1);
+      tr.begin_active()->face(f)->set_manifold_id (1);
 
   test(tr);
 }

diff --git a/tests/fe/mass_matrix_arbitrary_nodes.cc b/tests/fe/mass_matrix_arbitrary_nodes.cc
index 44d669105234f0c06c7edaa7c7831e315f830eeb..39b95d9ca39b94d7f15c39e1ba10c5b08cc5a12e 100644 (file)
--- a/tests/fe/mass_matrix_arbitrary_nodes.cc
+++ b/tests/fe/mass_matrix_arbitrary_nodes.cc
@@ -44,6 +44,7 @@ void check ()
         GridGenerator::hyper_cube (tr);
       else
         GridGenerator::hyper_ball (tr);
+      tr.reset_manifold(0);
       tr.refine_global (3-dim);
 
       FE_DGQArbitraryNodes<dim> fe ((QGauss<1>(q_points)));

diff --git a/tests/grid/grid_hyper_shell.cc b/tests/grid/grid_hyper_shell.cc
index a767e6e3739370b5403f3ee38df0e7bdee8b5947..874b880627a30efc128e2300aad2ced74b84b789 100644 (file)
--- a/tests/grid/grid_hyper_shell.cc
+++ b/tests/grid/grid_hyper_shell.cc
@@ -45,12 +45,15 @@ void check (double r1, double r2, unsigned int n)
   // 3d. the test was written in a way that it tests with this bug present, so
   // undo the big fix here
   GridGenerator::hyper_shell (tria, center, r1, r2, n, true);
+  tria.reset_manifold(0);
+  tria.set_all_manifold_ids(numbers::flat_manifold_id);
+  GridTools::copy_boundary_to_manifold_id(tria);
   if (dim==3)
     for (typename Triangulation<dim>::active_cell_iterator c=tria.begin_active(); c!=tria.end(); ++c)
       for (unsigned int f=0; f<GeometryInfo<dim>::faces_per_cell; ++f)
         if (c->face(f)->at_boundary())
           for (unsigned int e=0; e<GeometryInfo<dim>::lines_per_face; ++e)
-            c->face(f)->line(e)->set_boundary_id(0);
+            c->face(f)->line(e)->set_manifold_id(0);
 
   static const SphericalManifold<dim> boundary(center);
   tria.set_manifold(0, boundary);

diff --git a/tests/grid/grid_out_svg_01.cc b/tests/grid/grid_out_svg_01.cc
index 8aea71c5eeed960a66b8d9d6a0e55d2ed2edb265..be939cb54c57edfbe66570f85f1d315663de1a17 100644 (file)
--- a/tests/grid/grid_out_svg_01.cc
+++ b/tests/grid/grid_out_svg_01.cc
@@ -35,6 +35,7 @@ Triangulation<2,2> create_grid()
   Point<2> center(0., 0.);
 
   GridGenerator::hyper_cube_with_cylindrical_hole(triangulation, inner_radius, outer_radius);
+  triangulation.reset_manifold(0);
   triangulation.refine_global(1);
 
   Triangulation<2>::active_cell_iterator

diff --git a/tests/grid/grid_out_svg_02.cc b/tests/grid/grid_out_svg_02.cc
index e85b91ea229713cedd8c57fc7f24bceeb537a8b8..ef82acab1e8073b97b2756661ae775358979ca0d 100644 (file)
--- a/tests/grid/grid_out_svg_02.cc
+++ b/tests/grid/grid_out_svg_02.cc
@@ -35,6 +35,7 @@ Triangulation<2,2> create_grid()
   Point<2> center(0., 0.);
 
   GridGenerator::hyper_cube_with_cylindrical_hole(triangulation, inner_radius, outer_radius);
+  triangulation.reset_manifold(0);
   triangulation.refine_global(1);
 
   for (unsigned int l=0; l<3; ++l)

diff --git a/tests/grid/grid_output_input.cc b/tests/grid/grid_output_input.cc
index d1aad878aa951184c7d19780ccf54f9bf8ed8a30..6622c08c198cafff52e34cd8aa9e4df6f79ea540 100644 (file)
--- a/tests/grid/grid_output_input.cc
+++ b/tests/grid/grid_output_input.cc
@@ -24,6 +24,7 @@
 #include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/grid_in.h>
+#include <deal.II/grid/grid_tools.h>
 
 
 
@@ -36,6 +37,8 @@ void test(std::ostream &out)
   Triangulation<dim> tr;
 
   GridGenerator::hyper_cube_with_cylindrical_hole(tr, .3, .4, 1, 1, false);
+  tr.reset_manifold(0);
+  GridTools::copy_boundary_to_manifold_id(tr);
   CylindricalManifold<dim> boundary(2);
   tr.set_manifold(1, boundary);
   {
@@ -46,8 +49,6 @@ void test(std::ostream &out)
   tr.refine_global(1);
   deallog << "Writing refined from constructor" << std::endl;
   go.write_ucd(tr, out);
-
-  tr.reset_manifold(1);
   tr.clear();
 
   GridIn<dim> gi;
@@ -59,7 +60,8 @@ void test(std::ostream &out)
     grid_file.close();
   }
 
-  tr.set_manifold(1, boundary);
+  GridTools::map_boundary_to_manifold_ids({1}, {0},tr);
+  tr.set_manifold(0, boundary);
   tr.refine_global(1);
   deallog << "Writing refined from file" << std::endl;
   go.write_ucd(tr, out);
@@ -69,8 +71,6 @@ void test(std::ostream &out)
     go.write_msh(tr, grid_file);
     grid_file.close();
   }
-
-  tr.reset_manifold(1);
   tr.clear();
 }
 

diff --git a/tests/grid/grid_transform.cc b/tests/grid/grid_transform.cc
index aa5589795564bfdec8edd8efdfa645680f57d974..ba3d89abf67d3267da44720e0ff6e5347a3dd0bc 100644 (file)
--- a/tests/grid/grid_transform.cc
+++ b/tests/grid/grid_transform.cc
@@ -33,13 +33,14 @@ int main ()
 {
   const unsigned int dim=2;
   Point<dim> origin;
-  SphericalManifold<dim> boundary(origin);
   MappingQ<dim> mapping(2);
   Triangulation<dim> tria;
-  tria.set_manifold(0, boundary);
   const double inner_radius=1.;
   const double outer_radius=5.;
   GridGenerator::hyper_shell(tria, origin, inner_radius, outer_radius, 8);
+  GridTools::copy_boundary_to_manifold_id(tria);
+  SphericalManifold<dim> boundary(origin);
+  tria.set_manifold(0, boundary);
   tria.refine_global(2);
 
   // build up a map of vertex indices
@@ -94,7 +95,7 @@ int main ()
                 }
           }
     }
-
+  GridTools::copy_boundary_to_manifold_id(tria);
   GridTools::laplace_transform (new_points, tria);
   SphericalManifold<dim> inner_ball(n_center);
   tria.set_manifold(1, inner_ball);

diff --git a/tests/grid/grid_transform_coefficient.cc b/tests/grid/grid_transform_coefficient.cc
index a30c11acd4cf522d3fa9c30c6405f0dca305d17e..00afc48fb8cee5205275af27b037b8c555d42b88 100644 (file)
--- a/tests/grid/grid_transform_coefficient.cc
+++ b/tests/grid/grid_transform_coefficient.cc
@@ -47,10 +47,12 @@ int main ()
   SphericalManifold<dim> boundary(origin);
   MappingQ<dim> mapping(2);
   Triangulation<dim> tria;
-  tria.set_manifold(0, boundary);
   const double inner_radius=1.;
   const double outer_radius=5.;
   GridGenerator::hyper_shell(tria, origin, inner_radius, outer_radius, 8);
+  tria.set_all_manifold_ids(numbers::flat_manifold_id);
+  GridTools::copy_boundary_to_manifold_id(tria);
+  tria.set_manifold(0, boundary);
   tria.refine_global(2);
 
   // build up a map of vertex indices
@@ -98,7 +100,7 @@ int main ()
                       // circle.
                       new_points.insert(std::pair<types::global_dof_index, Point<dim> > (
                                           face->vertex_index(vertex_no), n_radius/inner_radius*v+n_center));
-                      face->set_boundary_id(1);
+                      face->set_manifold_id(1);
                     }
                   else
                     Assert(false, ExcInternalError());

diff --git a/tests/grid/manifold_ids_on_boundary_01.cc b/tests/grid/manifold_ids_on_boundary_01.cc
index fcb22668c2396f5443e1787b6888b7a7415bf5a7..7e571304735b2b2c30c6bd10777ebedb08403226 100644 (file)
--- a/tests/grid/manifold_ids_on_boundary_01.cc
+++ b/tests/grid/manifold_ids_on_boundary_01.cc
@@ -21,6 +21,7 @@
 #include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
 
 #include <iostream>
 
@@ -34,8 +35,10 @@ void dim_2(std::ostream &os)
   Triangulation<d> tr;
 
   GridGenerator::hyper_cube_with_cylindrical_hole(tr,inner,outer,outer,true);
+  tr.set_all_manifold_ids_on_boundary(numbers::flat_manifold_id);
+  tr.reset_manifold(0);
   static SphericalManifold<d> boundary;
-  tr.set_all_manifold_ids_on_boundary(1,1);
+  GridTools::copy_boundary_to_manifold_id(tr);
   tr.set_manifold(1,boundary);
 
   tr.refine_global(2);
@@ -55,8 +58,10 @@ void dim_3(std::ostream &os)
   Triangulation<d> tr;
 
   GridGenerator::hyper_cube_with_cylindrical_hole(tr,inner,outer,outer,true);
+  tr.set_all_manifold_ids_on_boundary(numbers::flat_manifold_id);
+  tr.reset_manifold(0);
   static CylindricalManifold<d> boundary(2);
-  tr.set_all_manifold_ids_on_boundary(1,1);
+  GridTools::copy_boundary_to_manifold_id(tr);
   tr.set_manifold(1,boundary);
 
   tr.refine_global(1);

diff --git a/tests/grid/manifold_ids_on_boundary_01.output b/tests/grid/manifold_ids_on_boundary_01.output
index 193c94030efc286206d752b8efb219f39320d385..544274ac9a3e3dd4913ead9eb7448bb9b83a2e88 100644 (file)
--- a/tests/grid/manifold_ids_on_boundary_01.output
+++ b/tests/grid/manifold_ids_on_boundary_01.output
@@ -678,7 +678,7 @@ POINTS 512 double
 1.00000 0.500000 0.500000
 0.600000 0.00000 1.00000
 0.592388 0.288268 1.00000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.592388 0.288268 0.500000
 1.00000 0.500000 1.00000
 1.00000 1.00000 1.00000
@@ -692,7 +692,7 @@ POINTS 512 double
 1.00000 0.500000 0.500000
 1.00000 0.00000 0.00000
 1.00000 0.500000 0.00000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.592388 0.288268 0.500000
 0.600000 0.00000 0.00000
 0.592388 0.288268 0.00000
@@ -706,11 +706,11 @@ POINTS 512 double
 0.570711 0.570711 0.00000
 0.600000 0.00000 1.00000
 0.592388 0.288268 1.00000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.592388 0.288268 0.500000
 0.200000 0.00000 1.00000
 0.184776 0.0765367 1.00000
-0.200000 0.00000 0.500000
+0.200000 -1.38778e-17 0.500000
 0.184776 0.0765367 0.500000
 0.592388 0.288268 1.00000
 0.570711 0.570711 1.00000
@@ -720,11 +720,11 @@ POINTS 512 double
 0.141421 0.141421 1.00000
 0.184776 0.0765367 0.500000
 0.141421 0.141421 0.500000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.592388 0.288268 0.500000
 0.600000 0.00000 0.00000
 0.592388 0.288268 0.00000
-0.200000 0.00000 0.500000
+0.200000 -1.38778e-17 0.500000
 0.184776 0.0765367 0.500000
 0.200000 0.00000 0.00000
 0.184776 0.0765367 0.00000
@@ -751,7 +751,7 @@ POINTS 512 double
 0.288268 0.592388 1.00000
 3.67394e-17 0.600000 1.00000
 0.288268 0.592388 0.500000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 1.00000 1.00000 0.500000
 0.500000 1.00000 0.500000
 1.00000 1.00000 0.00000
@@ -765,7 +765,7 @@ POINTS 512 double
 0.500000 1.00000 0.00000
 6.12323e-17 1.00000 0.00000
 0.288268 0.592388 0.500000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 0.288268 0.592388 0.00000
 3.67394e-17 0.600000 0.00000
 0.570711 0.570711 1.00000
@@ -779,11 +779,11 @@ POINTS 512 double
 0.288268 0.592388 1.00000
 3.67394e-17 0.600000 1.00000
 0.288268 0.592388 0.500000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 0.0765367 0.184776 1.00000
 1.22465e-17 0.200000 1.00000
 0.0765367 0.184776 0.500000
-1.22465e-17 0.200000 0.500000
+1.38778e-17 0.200000 0.500000
 0.570711 0.570711 0.500000
 0.288268 0.592388 0.500000
 0.570711 0.570711 0.00000
@@ -793,11 +793,11 @@ POINTS 512 double
 0.141421 0.141421 0.00000
 0.0765367 0.184776 0.00000
 0.288268 0.592388 0.500000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 0.288268 0.592388 0.00000
 3.67394e-17 0.600000 0.00000
 0.0765367 0.184776 0.500000
-1.22465e-17 0.200000 0.500000
+1.38778e-17 0.200000 0.500000
 0.0765367 0.184776 0.00000
 1.22465e-17 0.200000 0.00000
 6.12323e-17 1.00000 1.00000
@@ -806,7 +806,7 @@ POINTS 512 double
 -0.500000 1.00000 0.500000
 3.67394e-17 0.600000 1.00000
 -0.288268 0.592388 1.00000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 -0.288268 0.592388 0.500000
 -0.500000 1.00000 1.00000
 -1.00000 1.00000 1.00000
@@ -820,7 +820,7 @@ POINTS 512 double
 -0.500000 1.00000 0.500000
 6.12323e-17 1.00000 0.00000
 -0.500000 1.00000 0.00000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 -0.288268 0.592388 0.500000
 3.67394e-17 0.600000 0.00000
 -0.288268 0.592388 0.00000
@@ -834,11 +834,11 @@ POINTS 512 double
 -0.570711 0.570711 0.00000
 3.67394e-17 0.600000 1.00000
 -0.288268 0.592388 1.00000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 -0.288268 0.592388 0.500000
 1.22465e-17 0.200000 1.00000
 -0.0765367 0.184776 1.00000
-1.22465e-17 0.200000 0.500000
+1.38778e-17 0.200000 0.500000
 -0.0765367 0.184776 0.500000
 -0.288268 0.592388 1.00000
 -0.570711 0.570711 1.00000
@@ -848,11 +848,11 @@ POINTS 512 double
 -0.141421 0.141421 1.00000
 -0.0765367 0.184776 0.500000
 -0.141421 0.141421 0.500000
-3.67394e-17 0.600000 0.500000
+3.75551e-17 0.600000 0.500000
 -0.288268 0.592388 0.500000
 3.67394e-17 0.600000 0.00000
 -0.288268 0.592388 0.00000
-1.22465e-17 0.200000 0.500000
+1.38778e-17 0.200000 0.500000
 -0.0765367 0.184776 0.500000
 1.22465e-17 0.200000 0.00000
 -0.0765367 0.184776 0.00000
@@ -879,7 +879,7 @@ POINTS 512 double
 -0.592388 0.288268 1.00000
 -0.600000 7.34788e-17 1.00000
 -0.592388 0.288268 0.500000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -1.00000 1.00000 0.500000
 -1.00000 0.500000 0.500000
 -1.00000 1.00000 0.00000
@@ -893,7 +893,7 @@ POINTS 512 double
 -1.00000 0.500000 0.00000
 -1.00000 1.22465e-16 0.00000
 -0.592388 0.288268 0.500000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.592388 0.288268 0.00000
 -0.600000 7.34788e-17 0.00000
 -0.570711 0.570711 1.00000
@@ -907,11 +907,11 @@ POINTS 512 double
 -0.592388 0.288268 1.00000
 -0.600000 7.34788e-17 1.00000
 -0.592388 0.288268 0.500000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.184776 0.0765367 1.00000
 -0.200000 2.44929e-17 1.00000
 -0.184776 0.0765367 0.500000
--0.200000 2.44929e-17 0.500000
+-0.200000 4.16334e-17 0.500000
 -0.570711 0.570711 0.500000
 -0.592388 0.288268 0.500000
 -0.570711 0.570711 0.00000
@@ -921,11 +921,11 @@ POINTS 512 double
 -0.141421 0.141421 0.00000
 -0.184776 0.0765367 0.00000
 -0.592388 0.288268 0.500000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.592388 0.288268 0.00000
 -0.600000 7.34788e-17 0.00000
 -0.184776 0.0765367 0.500000
--0.200000 2.44929e-17 0.500000
+-0.200000 4.16334e-17 0.500000
 -0.184776 0.0765367 0.00000
 -0.200000 2.44929e-17 0.00000
 -1.00000 1.22465e-16 1.00000
@@ -934,7 +934,7 @@ POINTS 512 double
 -1.00000 -0.500000 0.500000
 -0.600000 7.34788e-17 1.00000
 -0.592388 -0.288268 1.00000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.592388 -0.288268 0.500000
 -1.00000 -0.500000 1.00000
 -1.00000 -1.00000 1.00000
@@ -948,7 +948,7 @@ POINTS 512 double
 -1.00000 -0.500000 0.500000
 -1.00000 1.22465e-16 0.00000
 -1.00000 -0.500000 0.00000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.592388 -0.288268 0.500000
 -0.600000 7.34788e-17 0.00000
 -0.592388 -0.288268 0.00000
@@ -962,11 +962,11 @@ POINTS 512 double
 -0.570711 -0.570711 0.00000
 -0.600000 7.34788e-17 1.00000
 -0.592388 -0.288268 1.00000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.592388 -0.288268 0.500000
 -0.200000 2.44929e-17 1.00000
 -0.184776 -0.0765367 1.00000
--0.200000 2.44929e-17 0.500000
+-0.200000 4.16334e-17 0.500000
 -0.184776 -0.0765367 0.500000
 -0.592388 -0.288268 1.00000
 -0.570711 -0.570711 1.00000
@@ -976,11 +976,11 @@ POINTS 512 double
 -0.141421 -0.141421 1.00000
 -0.184776 -0.0765367 0.500000
 -0.141421 -0.141421 0.500000
--0.600000 7.34788e-17 0.500000
+-0.600000 8.20490e-17 0.500000
 -0.592388 -0.288268 0.500000
 -0.600000 7.34788e-17 0.00000
 -0.592388 -0.288268 0.00000
--0.200000 2.44929e-17 0.500000
+-0.200000 4.16334e-17 0.500000
 -0.184776 -0.0765367 0.500000
 -0.200000 2.44929e-17 0.00000
 -0.184776 -0.0765367 0.00000
@@ -1007,7 +1007,7 @@ POINTS 512 double
 -0.288268 -0.592388 1.00000
 -1.10218e-16 -0.600000 1.00000
 -0.288268 -0.592388 0.500000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 -1.00000 -1.00000 0.500000
 -0.500000 -1.00000 0.500000
 -1.00000 -1.00000 0.00000
@@ -1021,7 +1021,7 @@ POINTS 512 double
 -0.500000 -1.00000 0.00000
 -1.83697e-16 -1.00000 0.00000
 -0.288268 -0.592388 0.500000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 -0.288268 -0.592388 0.00000
 -1.10218e-16 -0.600000 0.00000
 -0.570711 -0.570711 1.00000
@@ -1035,11 +1035,11 @@ POINTS 512 double
 -0.288268 -0.592388 1.00000
 -1.10218e-16 -0.600000 1.00000
 -0.288268 -0.592388 0.500000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 -0.0765367 -0.184776 1.00000
 -3.67394e-17 -0.200000 1.00000
 -0.0765367 -0.184776 0.500000
--3.67394e-17 -0.200000 0.500000
+-1.38778e-17 -0.200000 0.500000
 -0.570711 -0.570711 0.500000
 -0.288268 -0.592388 0.500000
 -0.570711 -0.570711 0.00000
@@ -1049,11 +1049,11 @@ POINTS 512 double
 -0.141421 -0.141421 0.00000
 -0.0765367 -0.184776 0.00000
 -0.288268 -0.592388 0.500000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 -0.288268 -0.592388 0.00000
 -1.10218e-16 -0.600000 0.00000
 -0.0765367 -0.184776 0.500000
--3.67394e-17 -0.200000 0.500000
+-1.38778e-17 -0.200000 0.500000
 -0.0765367 -0.184776 0.00000
 -3.67394e-17 -0.200000 0.00000
 -1.83697e-16 -1.00000 1.00000
@@ -1062,7 +1062,7 @@ POINTS 512 double
 0.500000 -1.00000 0.500000
 -1.10218e-16 -0.600000 1.00000
 0.288268 -0.592388 1.00000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 0.288268 -0.592388 0.500000
 0.500000 -1.00000 1.00000
 1.00000 -1.00000 1.00000
@@ -1076,7 +1076,7 @@ POINTS 512 double
 0.500000 -1.00000 0.500000
 -1.83697e-16 -1.00000 0.00000
 0.500000 -1.00000 0.00000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 0.288268 -0.592388 0.500000
 -1.10218e-16 -0.600000 0.00000
 0.288268 -0.592388 0.00000
@@ -1090,11 +1090,11 @@ POINTS 512 double
 0.570711 -0.570711 0.00000
 -1.10218e-16 -0.600000 1.00000
 0.288268 -0.592388 1.00000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 0.288268 -0.592388 0.500000
 -3.67394e-17 -0.200000 1.00000
 0.0765367 -0.184776 1.00000
--3.67394e-17 -0.200000 0.500000
+-1.38778e-17 -0.200000 0.500000
 0.0765367 -0.184776 0.500000
 0.288268 -0.592388 1.00000
 0.570711 -0.570711 1.00000
@@ -1104,11 +1104,11 @@ POINTS 512 double
 0.141421 -0.141421 1.00000
 0.0765367 -0.184776 0.500000
 0.141421 -0.141421 0.500000
--1.10218e-16 -0.600000 0.500000
+-9.87874e-17 -0.600000 0.500000
 0.288268 -0.592388 0.500000
 -1.10218e-16 -0.600000 0.00000
 0.288268 -0.592388 0.00000
--3.67394e-17 -0.200000 0.500000
+-1.38778e-17 -0.200000 0.500000
 0.0765367 -0.184776 0.500000
 -3.67394e-17 -0.200000 0.00000
 0.0765367 -0.184776 0.00000
@@ -1135,7 +1135,7 @@ POINTS 512 double
 0.592388 -0.288268 1.00000
 0.600000 0.00000 1.00000
 0.592388 -0.288268 0.500000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 1.00000 -1.00000 0.500000
 1.00000 -0.500000 0.500000
 1.00000 -1.00000 0.00000
@@ -1149,7 +1149,7 @@ POINTS 512 double
 1.00000 -0.500000 0.00000
 1.00000 0.00000 0.00000
 0.592388 -0.288268 0.500000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.592388 -0.288268 0.00000
 0.600000 0.00000 0.00000
 0.570711 -0.570711 1.00000
@@ -1163,11 +1163,11 @@ POINTS 512 double
 0.592388 -0.288268 1.00000
 0.600000 0.00000 1.00000
 0.592388 -0.288268 0.500000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.184776 -0.0765367 1.00000
 0.200000 0.00000 1.00000
 0.184776 -0.0765367 0.500000
-0.200000 0.00000 0.500000
+0.200000 -1.38778e-17 0.500000
 0.570711 -0.570711 0.500000
 0.592388 -0.288268 0.500000
 0.570711 -0.570711 0.00000
@@ -1177,11 +1177,11 @@ POINTS 512 double
 0.141421 -0.141421 0.00000
 0.184776 -0.0765367 0.00000
 0.592388 -0.288268 0.500000
-0.600000 0.00000 0.500000
+0.600000 -6.93889e-18 0.500000
 0.592388 -0.288268 0.00000
 0.600000 0.00000 0.00000
 0.184776 -0.0765367 0.500000
-0.200000 0.00000 0.500000
+0.200000 -1.38778e-17 0.500000
 0.184776 -0.0765367 0.00000
 0.200000 0.00000 0.00000
 
@@ -1259,7 +1259,7 @@ LOOKUP_TABLE default
 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 1.00000 
 SCALARS manifold double 1
 LOOKUP_TABLE default
--1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 -1.00000 
+0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 
 SCALARS material double 1
 LOOKUP_TABLE default
 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 

diff --git a/tests/grid/maximal_cell_diameter.cc b/tests/grid/maximal_cell_diameter.cc
index 91246c5c261397328c9f144a17f09ba1bb02b35c..76fc4760a45c61ed20fadc87e83367637b726074 100644 (file)
--- a/tests/grid/maximal_cell_diameter.cc
+++ b/tests/grid/maximal_cell_diameter.cc
@@ -55,6 +55,7 @@ void test1 ()
     {
       Triangulation<dim> tria;
       GridGenerator::hyper_ball(tria, Point<dim>(), 1);
+      tria.reset_manifold(0);
 
       for (unsigned int i=0; i<2; ++i)
         {

diff --git a/tests/grid/maximal_cell_diameter_mpi.cc b/tests/grid/maximal_cell_diameter_mpi.cc
index 36ea24eb55bec1120730c413b308d212251eea6e..4878a4e9e0fa3ad9746bca53211795ddef44030e 100644 (file)
--- a/tests/grid/maximal_cell_diameter_mpi.cc
+++ b/tests/grid/maximal_cell_diameter_mpi.cc
@@ -50,6 +50,7 @@ void test1 ()
     {
       parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
       GridGenerator::hyper_ball(tria, Point<dim>(), 1);
+      tria.reset_manifold(0);
 
       for (unsigned int i=0; i<2; ++i)
         {

diff --git a/tests/grid/mesh_3d_12.cc b/tests/grid/mesh_3d_12.cc
index 780597f0c0df34ae6526ea7581b66f06abb81747..d2b24665bf847aecf4bb9469de8d7815fee8676c 100644 (file)
--- a/tests/grid/mesh_3d_12.cc
+++ b/tests/grid/mesh_3d_12.cc
@@ -103,6 +103,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_manifold(0);
     check (coarse_grid);
   }
 

diff --git a/tests/grid/mesh_3d_14.cc b/tests/grid/mesh_3d_14.cc
index 32168e58e2dfdab9be306c2a854025586752d927..207b39a7a3b6089932a2eba642094f18bcd89d9e 100644 (file)
--- a/tests/grid/mesh_3d_14.cc
+++ b/tests/grid/mesh_3d_14.cc
@@ -137,6 +137,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_manifold(0);
     check (coarse_grid);
   }
 

diff --git a/tests/grid/mesh_3d_26.cc b/tests/grid/mesh_3d_26.cc
index 325c2377f3549beec5fec95aeb96fbca2b35d47d..0f8e61b64bb86571ef359da10fe8f269748aafd1 100644 (file)
--- a/tests/grid/mesh_3d_26.cc
+++ b/tests/grid/mesh_3d_26.cc
@@ -114,6 +114,7 @@ int main ()
 
   Triangulation<3> triangulation;
   GridGenerator::hyper_ball(triangulation);
+  triangulation.reset_manifold(0);
   triangulation.refine_global(1);
   test<3>(triangulation);
 }

diff --git a/tests/grid/mesh_3d_6.cc b/tests/grid/mesh_3d_6.cc
index 01c8ecca117a1e18d024a17382e22eeb1ccb95ba..a279f1f72a4a78bdc539e66fa78f87b72185baae 100644 (file)
--- a/tests/grid/mesh_3d_6.cc
+++ b/tests/grid/mesh_3d_6.cc
@@ -126,6 +126,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_manifold(0);
     check (coarse_grid);
   }
 

diff --git a/tests/grid/mesh_3d_7.cc b/tests/grid/mesh_3d_7.cc
index 3dfa2e5b9d023750e29de74748beda0fc85a8234..56bb6941cb886585a6f9ddd1384e4085a3d0e92f 100644 (file)
--- a/tests/grid/mesh_3d_7.cc
+++ b/tests/grid/mesh_3d_7.cc
@@ -135,6 +135,7 @@ int main ()
   {
     Triangulation<3> coarse_grid;
     GridGenerator::hyper_ball (coarse_grid);
+    coarse_grid.reset_manifold(0);
     check (coarse_grid);
   }
 

diff --git a/tests/grid/minimal_cell_diameter.cc b/tests/grid/minimal_cell_diameter.cc
index 4a120b8a03280a24e8c2d28d0042ac36623adcfe..381ad4418b7620625d77a88b8ee2a82b622620b2 100644 (file)
--- a/tests/grid/minimal_cell_diameter.cc
+++ b/tests/grid/minimal_cell_diameter.cc
@@ -51,6 +51,7 @@ void test1 ()
     {
       Triangulation<dim> tria;
       GridGenerator::hyper_ball(tria, Point<dim>(), 1);
+      tria.reset_manifold(0);
 
       for (unsigned int i=0; i<2; ++i)
         {

diff --git a/tests/grid/minimal_cell_diameter_mpi.cc b/tests/grid/minimal_cell_diameter_mpi.cc
index 1d3afdeafc9682c5771559c2e9aa983c2537959f..041b14c539dfc2502bb79c0ef2e688df981e63c5 100644 (file)
--- a/tests/grid/minimal_cell_diameter_mpi.cc
+++ b/tests/grid/minimal_cell_diameter_mpi.cc
@@ -46,6 +46,7 @@ void test1 ()
     {
       parallel::distributed::Triangulation<dim> tria(MPI_COMM_WORLD);
       GridGenerator::hyper_ball(tria, Point<dim>(), 1);
+      tria.reset_manifold(0);
 
       for (unsigned int i=0; i<2; ++i)
         {

diff --git a/tests/grid/recursively_set_material_id.cc b/tests/grid/recursively_set_material_id.cc
index ac95f1d390f0c953f6030dd5e0150a0ad648ccdf..0966a0f69e461f1d0baafdb361e21250cbe2fda9 100644 (file)
--- a/tests/grid/recursively_set_material_id.cc
+++ b/tests/grid/recursively_set_material_id.cc
@@ -35,7 +35,10 @@ void test()
 
       Triangulation<dim> tr;
       if (dim > 1)
-        GridGenerator::hyper_ball(tr);
+        {
+          GridGenerator::hyper_ball(tr);
+          tr.reset_manifold(0);
+        }
       else
         GridGenerator::hyper_cube(tr);
 

diff --git a/tests/hp/boundary_matrices.cc b/tests/hp/boundary_matrices.cc
index 779eba87b53fbc7f540b96b2c329396e8f1671f7..74321056041914c0f4566669e97a26dcfa0f7cfb 100644 (file)
--- a/tests/hp/boundary_matrices.cc
+++ b/tests/hp/boundary_matrices.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/boundary_matrices_hp.cc b/tests/hp/boundary_matrices_hp.cc
index 7723cb304b84902cdbb3aba4e7cf84e5ed0cb8ee..2ecccddc083f0b9661fe0d6af8a7649c2b8d6299 100644 (file)
--- a/tests/hp/boundary_matrices_hp.cc
+++ b/tests/hp/boundary_matrices_hp.cc
@@ -76,6 +76,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_01.cc b/tests/hp/create_laplace_matrix_01.cc
index 058691c08b28bd8c5da3fd31e1cfafbd66c17a8e..c89640139fd7fa6718b47d221e8c12aa142c8599 100644 (file)
--- a/tests/hp/create_laplace_matrix_01.cc
+++ b/tests/hp/create_laplace_matrix_01.cc
@@ -61,6 +61,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_01b.cc b/tests/hp/create_laplace_matrix_01b.cc
index fc7a45739233ba6227fef3df8e209eb6af26b577..8a2b04af0495d19f5721641e46cfec459de5ffa7 100644 (file)
--- a/tests/hp/create_laplace_matrix_01b.cc
+++ b/tests/hp/create_laplace_matrix_01b.cc
@@ -62,6 +62,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_02.cc b/tests/hp/create_laplace_matrix_02.cc
index 852e451a4176f24fd0b0d8720cf657a5490cd418..26f70f5b03ff2fd59247c998b9db87a46f09ccbd 100644 (file)
--- a/tests/hp/create_laplace_matrix_02.cc
+++ b/tests/hp/create_laplace_matrix_02.cc
@@ -60,6 +60,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_02b.cc b/tests/hp/create_laplace_matrix_02b.cc
index 0286f64dbb932dec45ca1ea31d5156e4422ac003..4791a29ea277d774e44688c978f3c4ac086fb6ea 100644 (file)
--- a/tests/hp/create_laplace_matrix_02b.cc
+++ b/tests/hp/create_laplace_matrix_02b.cc
@@ -61,6 +61,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_03.cc b/tests/hp/create_laplace_matrix_03.cc
index 51ec100ac72afd3f63e85186f6fa692158d35c8e..c1e5608dcdffd77f8f7537bb478f76f675310330 100644 (file)
--- a/tests/hp/create_laplace_matrix_03.cc
+++ b/tests/hp/create_laplace_matrix_03.cc
@@ -60,6 +60,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_03b.cc b/tests/hp/create_laplace_matrix_03b.cc
index 9340e051bcfcdb661616dea91cc120ea83dce87f..814e614e453b1ec982617a9bb561038551396c62 100644 (file)
--- a/tests/hp/create_laplace_matrix_03b.cc
+++ b/tests/hp/create_laplace_matrix_03b.cc
@@ -61,6 +61,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_04.cc b/tests/hp/create_laplace_matrix_04.cc
index 938fa314be01cba8b9713d79c8ec2c1f5489a257..da477709fc1e951308b9be2babcbfcc02e54a125 100644 (file)
--- a/tests/hp/create_laplace_matrix_04.cc
+++ b/tests/hp/create_laplace_matrix_04.cc
@@ -59,6 +59,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_laplace_matrix_04b.cc b/tests/hp/create_laplace_matrix_04b.cc
index a2fe2af6f71132fe7acfa9d8aa108b96c56d8367..9288b58ef11252350eea9a8e605a93bbdce9dc94 100644 (file)
--- a/tests/hp/create_laplace_matrix_04b.cc
+++ b/tests/hp/create_laplace_matrix_04b.cc
@@ -60,6 +60,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_01.cc b/tests/hp/create_mass_matrix_01.cc
index 0b8a8406792a2b4ed811e809f221560a607fb5b2..3bb93575c0bbbb62541db90be12b61410fa2912a 100644 (file)
--- a/tests/hp/create_mass_matrix_01.cc
+++ b/tests/hp/create_mass_matrix_01.cc
@@ -61,6 +61,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_01b.cc b/tests/hp/create_mass_matrix_01b.cc
index c85e6bb8edd0d6b66fc86dbe1636cdece81e3151..9e5e309245cd24090f1e23c4592f8394b36391aa 100644 (file)
--- a/tests/hp/create_mass_matrix_01b.cc
+++ b/tests/hp/create_mass_matrix_01b.cc
@@ -62,6 +62,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_02.cc b/tests/hp/create_mass_matrix_02.cc
index 1a61e741424f845a011f7568148c4e9f47d8f8e5..537c90fbbb93b2d12aa27cf08c8d3fd33bef6a1f 100644 (file)
--- a/tests/hp/create_mass_matrix_02.cc
+++ b/tests/hp/create_mass_matrix_02.cc
@@ -60,6 +60,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_02b.cc b/tests/hp/create_mass_matrix_02b.cc
index 3fe28037c119d1571ea48db76db0a077313f176e..77f96e9c800c8a326eb295e66ced542299c7f6c3 100644 (file)
--- a/tests/hp/create_mass_matrix_02b.cc
+++ b/tests/hp/create_mass_matrix_02b.cc
@@ -61,6 +61,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_03.cc b/tests/hp/create_mass_matrix_03.cc
index 32970c15900abe1eea21dc75a41526e1d13e65b4..8bc602ba7e2ae6e96c74e9fef9e53b398d5f4e31 100644 (file)
--- a/tests/hp/create_mass_matrix_03.cc
+++ b/tests/hp/create_mass_matrix_03.cc
@@ -60,6 +60,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_03b.cc b/tests/hp/create_mass_matrix_03b.cc
index 33d008dbeb8fb6fb52d58ca39a3ead47b941df33..44f16ec2c7a28a32d3910586ab8908fd96f3a23d 100644 (file)
--- a/tests/hp/create_mass_matrix_03b.cc
+++ b/tests/hp/create_mass_matrix_03b.cc
@@ -61,6 +61,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_04.cc b/tests/hp/create_mass_matrix_04.cc
index 9b390e1058f838744ff133dd5ab4031a1fd8f97c..1b7c55482f60212e5b68c6f3d6e722917d5e9243 100644 (file)
--- a/tests/hp/create_mass_matrix_04.cc
+++ b/tests/hp/create_mass_matrix_04.cc
@@ -59,6 +59,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_04b.cc b/tests/hp/create_mass_matrix_04b.cc
index 5e7e050252f52cc3b997e4afa7031afa4915a87d..4e4c5a27d34f1b1e97da630f6a0105f6c6ab63f7 100644 (file)
--- a/tests/hp/create_mass_matrix_04b.cc
+++ b/tests/hp/create_mass_matrix_04b.cc
@@ -60,6 +60,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/create_mass_matrix_05.cc b/tests/hp/create_mass_matrix_05.cc
index 01c92d38e28a5bf1b2fa76fcabafae8e420b97cc..3cf600b1b6a96cf00af7b01753b4d9dc3e2b4db6 100644 (file)
--- a/tests/hp/create_mass_matrix_05.cc
+++ b/tests/hp/create_mass_matrix_05.cc
@@ -52,6 +52,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/matrices.cc b/tests/hp/matrices.cc
index c038ee43b464449d98ca7f559e7f7271a618b837..4a04029f6aa14a40e07c6292d2cbd09a4bbaf44b 100644 (file)
--- a/tests/hp/matrices.cc
+++ b/tests/hp/matrices.cc
@@ -132,6 +132,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/matrices_hp.cc b/tests/hp/matrices_hp.cc
index 96feec31ffb3f1ace57a6d46d8c52bfe06f85287..1ac803fa998016ca0aa637b19e1a07bffc2f026c 100644 (file)
--- a/tests/hp/matrices_hp.cc
+++ b/tests/hp/matrices_hp.cc
@@ -131,6 +131,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/step-2.cc b/tests/hp/step-2.cc
index 6e2adbe4e7cfb50d9cad3aebac970e61f287f4ee..4647b36b86dc8519011939a0e0d5438cf740ca04 100644 (file)
--- a/tests/hp/step-2.cc
+++ b/tests/hp/step-2.cc
@@ -22,6 +22,7 @@
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/grid_generator.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/manifold_lib.h>
 
 #include <deal.II/hp/dof_handler.h>
@@ -44,6 +45,7 @@ void make_grid (Triangulation<2> &triangulation)
   GridGenerator::hyper_shell (triangulation,
                               center, inner_radius, outer_radius,
                               10);
+  GridTools::copy_boundary_to_manifold_id(triangulation);
 
   static const SphericalManifold<2> boundary_description(center);
   triangulation.set_manifold (0, boundary_description);

diff --git a/tests/hp/step-5.cc b/tests/hp/step-5.cc
index 8f3eaf755043f3671c8e72514f35a251005167c6..aa050a5eb8ea3b5be002557b14573fcb16e53beb 100644 (file)
--- a/tests/hp/step-5.cc
+++ b/tests/hp/step-5.cc
@@ -31,6 +31,7 @@
 #include <deal.II/hp/dof_handler.h>
 #include <deal.II/grid/tria_accessor.h>
 #include <deal.II/grid/tria_iterator.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/dofs/dof_tools.h>
 #include <deal.II/fe/fe_q.h>
@@ -295,7 +296,8 @@ void LaplaceProblem<dim>::run ()
           std::ifstream input_file(SOURCE_DIR "/grids/circle-grid.inp");
           Assert (dim==2, ExcInternalError());
 
-          grid_in.read_ucd (input_file);
+          grid_in.read_ucd (input_file, true);
+          GridTools::copy_boundary_to_manifold_id(triangulation);
 
           static const SphericalManifold<dim> boundary;
           triangulation.set_manifold (0, boundary);

diff --git a/tests/hp/vectors_boundary_rhs_01.cc b/tests/hp/vectors_boundary_rhs_01.cc
index de56955fea93e92a27b8a63132a037b3bcc7a1c3..adc44b4a09680045bf9409b374a4b610d554776a 100644 (file)
--- a/tests/hp/vectors_boundary_rhs_01.cc
+++ b/tests/hp/vectors_boundary_rhs_01.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_boundary_rhs_02.cc b/tests/hp/vectors_boundary_rhs_02.cc
index 53d8d0f12c7cdc545158767c5e74d56f8776f019..c6b7989459a003a11ac6cf2cb8d8e98a70c632ee 100644 (file)
--- a/tests/hp/vectors_boundary_rhs_02.cc
+++ b/tests/hp/vectors_boundary_rhs_02.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_boundary_rhs_03.cc b/tests/hp/vectors_boundary_rhs_03.cc
index d7b7e6eb89c9d933ed29911089b465d2abaf08fe..82003503ccff0931bf7da93e7ca4d0ba2d907e3f 100644 (file)
--- a/tests/hp/vectors_boundary_rhs_03.cc
+++ b/tests/hp/vectors_boundary_rhs_03.cc
@@ -73,6 +73,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_boundary_rhs_hp_01.cc b/tests/hp/vectors_boundary_rhs_hp_01.cc
index 1f8f86de574f1f164423dbe58953b0ac2d793a22..8e6ea100506b0d6727908c40dad53529041b7c6d 100644 (file)
--- a/tests/hp/vectors_boundary_rhs_hp_01.cc
+++ b/tests/hp/vectors_boundary_rhs_hp_01.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_boundary_rhs_hp_03.cc b/tests/hp/vectors_boundary_rhs_hp_03.cc
index b0cdda8f8547b96304711987e0150eff2e9f7e12..9c3d5de45695a2a4cac6632da1e5019d9118c6d2 100644 (file)
--- a/tests/hp/vectors_boundary_rhs_hp_03.cc
+++ b/tests/hp/vectors_boundary_rhs_hp_03.cc
@@ -73,6 +73,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_point_source_01.cc b/tests/hp/vectors_point_source_01.cc
index fd5e5d10cb11a82ab7d79e885cddd4bf64edf8a4..ada1b8a7187f50017dec02898944826f90f7180c 100644 (file)
--- a/tests/hp/vectors_point_source_01.cc
+++ b/tests/hp/vectors_point_source_01.cc
@@ -52,6 +52,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_point_source_hp_01.cc b/tests/hp/vectors_point_source_hp_01.cc
index b5d6576b57df6b60ab62c1a917203bab5c168cbe..f26eca9b514350819133ea546c7621265a34cfef 100644 (file)
--- a/tests/hp/vectors_point_source_hp_01.cc
+++ b/tests/hp/vectors_point_source_hp_01.cc
@@ -51,6 +51,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_rhs_01.cc b/tests/hp/vectors_rhs_01.cc
index a9fcd74ec35db8b1f3affe3fd30f4acc3b2588ba..8a274ba9e4510a497ab28dec22d53c492cac2d14 100644 (file)
--- a/tests/hp/vectors_rhs_01.cc
+++ b/tests/hp/vectors_rhs_01.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_rhs_02.cc b/tests/hp/vectors_rhs_02.cc
index 1713d4a95bb5929029ab02aa1cee876dbd605ef2..9e78259f922aaaf37c63c86b4ca4471f15f5ee9c 100644 (file)
--- a/tests/hp/vectors_rhs_02.cc
+++ b/tests/hp/vectors_rhs_02.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_rhs_03.cc b/tests/hp/vectors_rhs_03.cc
index 38c9ebf1fed260686b4f5386d106be5a9c41a856..735d2cd715460fc3b62262be6d6ab1a281841ea7 100644 (file)
--- a/tests/hp/vectors_rhs_03.cc
+++ b/tests/hp/vectors_rhs_03.cc
@@ -73,6 +73,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_rhs_hp_01.cc b/tests/hp/vectors_rhs_hp_01.cc
index fe6843f9838ba6b7402e687028e8613149d84886..f2cb7c760b69258164fccf95d67a1da2f27a44c5 100644 (file)
--- a/tests/hp/vectors_rhs_hp_01.cc
+++ b/tests/hp/vectors_rhs_hp_01.cc
@@ -74,6 +74,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/hp/vectors_rhs_hp_03.cc b/tests/hp/vectors_rhs_hp_03.cc
index a60d7fed59fe48bdd7d2838888f7d9257c9906af..241169d4b7bf267816c3546b5fd016f87b89ae4d 100644 (file)
--- a/tests/hp/vectors_rhs_hp_03.cc
+++ b/tests/hp/vectors_rhs_hp_03.cc
@@ -73,6 +73,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/lac/filtered_matrix_06.cc b/tests/lac/filtered_matrix_06.cc
index b1a75f4c0d3156e81ae1426f3b4b4475ab1757f2..3605109bfdef292b2b6ffbe33b7aca89620a0694 100644 (file)
--- a/tests/lac/filtered_matrix_06.cc
+++ b/tests/lac/filtered_matrix_06.cc
@@ -104,6 +104,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
 
   tr.refine_global (5-dim);
 

diff --git a/tests/lac/inhomogeneous_constraints_nonsymmetric.cc b/tests/lac/inhomogeneous_constraints_nonsymmetric.cc
index bb0cdfbff3d2b56330e76351bb8a2e9c304a9090..e9989d059135b44517698351c8c444801ffb7e2c 100644 (file)
--- a/tests/lac/inhomogeneous_constraints_nonsymmetric.cc
+++ b/tests/lac/inhomogeneous_constraints_nonsymmetric.cc
@@ -428,6 +428,7 @@ template <int dim>
 void AdvectionProblem<dim>::run ()
 {
   GridGenerator::hyper_ball (triangulation);
+  triangulation.reset_manifold(0);
   triangulation.refine_global (3-dim);
 
   // manually refine the first two cells

diff --git a/tests/lac/matrices_02.cc b/tests/lac/matrices_02.cc
index 8e9c3cbb5bd85b3165c8ae9349563c1eba8be68a..3fb7670b02e4d4a0020b74a38d80481e59ba4fd1 100644 (file)
--- a/tests/lac/matrices_02.cc
+++ b/tests/lac/matrices_02.cc
@@ -127,6 +127,7 @@ check ()
     GridGenerator::hyper_ball(tr, Point<dim>(), 1);
   else
     GridGenerator::hyper_cube(tr, -1,1);
+  tr.reset_manifold(0);
   tr.refine_global (1);
   tr.begin_active()->set_refine_flag ();
   tr.execute_coarsening_and_refinement ();

diff --git a/tests/manifold/polar_manifold_02.cc b/tests/manifold/polar_manifold_02.cc
index 25711da93fe8965b122f56754719eb6f1a828b3c..db36463986c7e73d0d4f2ab44e81111edc89c11d 100644 (file)
--- a/tests/manifold/polar_manifold_02.cc
+++ b/tests/manifold/polar_manifold_02.cc
@@ -37,6 +37,7 @@ void test(unsigned int ref=1)
 
   Triangulation<dim,spacedim> tria;
   GridGenerator::hyper_ball (tria);
+  tria.reset_manifold(0);
 
   typename Triangulation<dim,spacedim>::active_cell_iterator cell;
 

diff --git a/tests/manifold/transfinite_manifold_02.cc b/tests/manifold/transfinite_manifold_02.cc
index 3a7c6df9b870e0241894975f317e1a20fedbc465..5eb6c86a021a05fff6181a5fca78032260628707 100644 (file)
--- a/tests/manifold/transfinite_manifold_02.cc
+++ b/tests/manifold/transfinite_manifold_02.cc
@@ -18,6 +18,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/manifold_lib.h>
 
 
@@ -29,11 +30,19 @@ void do_test(const Triangulation<dim,spacedim> &tria)
     {
       deallog << "Lines on cell with center: " << cell->center() << std::endl;
       for (unsigned int line=0; line<GeometryInfo<dim>::lines_per_cell; ++line)
-        deallog << cell->line(line)->center(/*respect_manifold=*/true) << std::endl;
+        {
+          deallog << cell->line(line)->center(/*respect_manifold=*/true) << std::endl;
+          deallog << "Manifold id: " << cell->line(line)->manifold_id() << std::endl;
+        }
       deallog << "Faces on cell with center: " << cell->center() << std::endl;
       for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
-        deallog << cell->face(face)->center(/*respect_manifold=*/true) << std::endl;
+        {
+          deallog << cell->face(face)->center(/*respect_manifold=*/true) << std::endl;
+          deallog << "Manifold id: " << cell->face(face)->manifold_id() << std::endl;
+        }
+
       deallog << "Center with manifold: " << cell->center(true) << std::endl;
+      deallog << "Manifold id: " << cell->manifold_id() << std::endl;
       for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
         if (cell->at_boundary(face))
           {
@@ -49,6 +58,7 @@ void do_test(const Triangulation<dim,spacedim> &tria)
                                    make_array_view(weights));
             deallog << "Distance between cell manifold and face manifold: "
                     << (pref-p) << std::endl;
+            Assert ((pref-p).norm() < 1e-6, ExcInternalError());
             weights[0] = 0.55;
             weights[1] = 0.45;
             p = cell->get_manifold().get_new_point(make_array_view(points),
@@ -57,6 +67,7 @@ void do_test(const Triangulation<dim,spacedim> &tria)
                                                                   make_array_view(weights));
             deallog << "Distance between cell manifold and face manifold: "
                     << (pref-p) << std::endl;
+            Assert ((pref-p).norm() < 1e-6, ExcInternalError());
           }
     }
   deallog << std::endl;
@@ -73,6 +84,7 @@ void test_polar()
 
   Triangulation<dim,spacedim> tria;
   GridGenerator::hyper_ball (tria);
+  tria.set_all_manifold_ids(numbers::flat_manifold_id);
 
   // set all cells at the boundary to the transfinite manifold except for the
   // boundaries where we put the polar manifold
@@ -100,6 +112,7 @@ void test_spherical()
 
   Triangulation<dim,spacedim> tria;
   GridGenerator::hyper_ball (tria);
+  tria.set_all_manifold_ids(numbers::flat_manifold_id);
 
   // set all cells at the boundary to the transfinite manifold except for the
   // boundaries where we put the polar manifold
@@ -125,12 +138,14 @@ void test_cylinder(unsigned int ref=1)
   TransfiniteInterpolationManifold<dim,spacedim> manifold;
   Triangulation<dim,spacedim> tria;
   GridGenerator::cylinder (tria);
+  tria.reset_all_manifolds();
+  tria.set_all_manifold_ids(numbers::flat_manifold_id);
 
   for (typename Triangulation<dim,spacedim>::active_cell_iterator
        cell=tria.begin_active(); cell!=tria.end(); ++cell)
     {
       if (cell->at_boundary())
-        cell->set_manifold_id(1);
+        cell->set_all_manifold_ids(1);
       for (unsigned int face=0; face<GeometryInfo<dim>::faces_per_cell; ++face)
         if (cell->at_boundary(face))
           {

diff --git a/tests/manifold/transfinite_manifold_02.output b/tests/manifold/transfinite_manifold_02.output
index a5b642cdc5b9d9bbcc306992e330bf1828d9d390..16ceb78898e53f7dbeba2cbe2cd9fd7b24f8f12f 100644 (file)
--- a/tests/manifold/transfinite_manifold_02.output
+++ b/tests/manifold/transfinite_manifold_02.output
@@ -2,539 +2,952 @@
 DEAL::Testing with PolarManifold dim=2, spacedim=2
 DEAL::Lines on cell with center: 0.00000 -0.500000
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-1.83697e-16 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 -0.500000
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-1.83697e-16 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: -9.71445e-17 -0.646447
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Lines on cell with center: -0.500000 -1.38778e-17
 DEAL::-1.00000 1.22465e-16
+DEAL::Manifold id: 0
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: -0.500000 -1.38778e-17
 DEAL::-1.00000 1.22465e-16
+DEAL::Manifold id: 0
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
-DEAL::Center with manifold: -0.646447 4.16334e-17
+DEAL::Manifold id: 4294967295
+DEAL::Center with manifold: -0.646447 5.55112e-17
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Lines on cell with center: 0.00000 0.00000
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 0.00000
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 0.00000 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Lines on cell with center: 0.500000 0.00000
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::1.00000 0.00000
+DEAL::Manifold id: 0
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.500000 0.00000
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::1.00000 0.00000
+DEAL::Manifold id: 0
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 0.646447 0.00000
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Lines on cell with center: -1.38778e-17 0.500000
 DEAL::6.12323e-17 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: -1.38778e-17 0.500000
 DEAL::6.12323e-17 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 2.77556e-17 0.646447
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::
 DEAL::Testing with SphericalManifold dim=2, spacedim=2
 DEAL::Lines on cell with center: 0.00000 -0.500000
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-1.11022e-16 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 -0.500000
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-1.11022e-16 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: -5.55112e-17 -0.646447
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Lines on cell with center: -0.500000 -1.38778e-17
 DEAL::-1.00000 -1.11022e-16
+DEAL::Manifold id: 0
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: -0.500000 -1.38778e-17
 DEAL::-1.00000 -1.11022e-16
+DEAL::Manifold id: 0
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
-DEAL::Center with manifold: -0.646447 -1.29480e-14
+DEAL::Manifold id: 4294967295
+DEAL::Center with manifold: -0.646447 -5.55112e-17
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Lines on cell with center: 0.00000 0.00000
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 0.00000
 DEAL::-0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.292893
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 0.00000 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Lines on cell with center: 0.500000 0.00000
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::1.00000 -1.11022e-16
+DEAL::Manifold id: 0
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.500000 0.00000
 DEAL::0.500000 -0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::1.00000 -1.11022e-16
+DEAL::Manifold id: 0
 DEAL::0.292893 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 0.646447 -5.55112e-17
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Lines on cell with center: -1.38778e-17 0.500000
 DEAL::-1.11022e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: -1.38778e-17 0.500000
 DEAL::-1.11022e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.292893
+DEAL::Manifold id: 4294967295
 DEAL::-0.500000 0.500000
+DEAL::Manifold id: 4294967295
 DEAL::0.500000 0.500000
-DEAL::Center with manifold: -1.29480e-14 0.646447
+DEAL::Manifold id: 4294967295
+DEAL::Center with manifold: -5.55112e-17 0.646447
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000
 DEAL::
 DEAL::Testing with SphericalManifold dim=3, spacedim=3
 DEAL::Lines on cell with center: 0.00000 0.00000 -6.93889e-18
 DEAL::-0.211325 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.211325 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.211325 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.211325 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.211325 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.211325 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.211325 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.211325 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 0.00000 -6.93889e-18
 DEAL::-0.211325 0.00000 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.00000 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 0.00000 0.00000 -1.38778e-17
+DEAL::Manifold id: 4294967295
 DEAL::Lines on cell with center: 0.00000 0.00000 -0.394338
 DEAL::-0.707107 0.00000 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.707107 0.00000 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.211325 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.211325 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 -0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 -0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 0.00000 -0.394338
 DEAL::-0.410557 0.00000 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 0.00000 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::6.93889e-18 -0.410557 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::6.93889e-18 0.410557 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.00000 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: -1.94289e-16 2.37914e-13 -0.450853
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: -5.55112e-17 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: -6.93889e-18 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 1.38778e-17 0.00000 0.00000
 DEAL::Lines on cell with center: 0.394338 0.00000 2.08167e-17
 DEAL::0.394338 -0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.707107 0.00000 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.211325 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 -0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.707107 0.00000 0.707107
+DEAL::Manifold id: 0
 DEAL::0.211325 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.707107 -0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::0.707107 0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::0.211325 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.394338 0.00000 2.08167e-17
 DEAL::0.410557 -0.410557 6.93889e-18
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 0.410557 6.93889e-18
+DEAL::Manifold id: 4294967295
 DEAL::1.00000 0.00000 0.00000
+DEAL::Manifold id: 0
 DEAL::0.211325 0.00000 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 0.00000 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 0.00000 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 0.450853 1.19127e-13 0.00000
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 1.11022e-16 -1.11022e-16 5.55112e-17
-DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 -6.93889e-18
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 6.24500e-17
 DEAL::Lines on cell with center: 0.00000 2.08167e-17 0.394338
 DEAL::-0.394338 -0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 -0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.211325 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 0.211325 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.707107 0.00000 0.707107
+DEAL::Manifold id: 0
 DEAL::0.707107 0.00000 0.707107
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 2.08167e-17 0.394338
 DEAL::-0.410557 6.93889e-18 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 0.00000 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.00000 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.410557 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.410557 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 1.19127e-13 0.00000 0.450853
-DEAL::Distance between cell manifold and face manifold: -1.11022e-16 0.00000 1.11022e-16
-DEAL::Distance between cell manifold and face manifold: 0.00000 -6.93889e-18 0.00000
+DEAL::Manifold id: 1
+DEAL::Distance between cell manifold and face manifold: -1.11022e-16 5.55112e-17 1.11022e-16
+DEAL::Distance between cell manifold and face manifold: 0.00000 6.93889e-17 0.00000
 DEAL::Lines on cell with center: -0.394338 3.46945e-18 1.73472e-17
 DEAL::-0.707107 0.00000 -0.707107
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.00000 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 -0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::-0.707107 0.00000 0.707107
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.00000 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 -0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::-0.707107 -0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::-0.211325 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.707107 0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: -0.394338 3.46945e-18 1.73472e-17
 DEAL::-1.00000 0.00000 0.00000
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.00000 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.410557 -0.410557 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.410557 0.410557 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.410557 0.00000 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::-0.410557 6.93889e-18 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: -0.450853 1.19113e-13 -6.93889e-18
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 -5.55112e-17 0.00000
-DEAL::Distance between cell manifold and face manifold: 0.00000 -6.93889e-18 0.00000
+DEAL::Distance between cell manifold and face manifold: 0.00000 1.38778e-17 0.00000
 DEAL::Lines on cell with center: 0.00000 -0.394338 2.08167e-17
 DEAL::-0.394338 -0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 -0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.211325 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 -0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 -0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.211325 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.707107 -0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::0.707107 -0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::-0.211325 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.211325 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 0.00000 -0.394338 2.08167e-17
 DEAL::-0.410557 -0.410557 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 -0.410557 6.93889e-18
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -1.00000 0.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 -0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::6.93889e-18 -0.410557 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 -0.410557 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 1.19127e-13 -0.450853 0.00000
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: -1.11022e-16 -1.11022e-16 5.55112e-17
-DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 -6.93889e-18
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 6.24500e-17
 DEAL::Lines on cell with center: 3.46945e-18 0.394338 1.73472e-17
 DEAL::0.00000 0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 0.211325 -0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 0.394338 -0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.00000 0.211325 0.211325
+DEAL::Manifold id: 4294967295
 DEAL::-0.394338 0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::0.394338 0.394338 0.394338
+DEAL::Manifold id: 4294967295
 DEAL::-0.707107 0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::-0.211325 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.707107 0.707107 0.00000
+DEAL::Manifold id: 0
 DEAL::0.211325 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::Faces on cell with center: 3.46945e-18 0.394338 1.73472e-17
 DEAL::0.00000 1.00000 0.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.211325 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::-0.410557 0.410557 0.00000
+DEAL::Manifold id: 4294967295
 DEAL::0.410557 0.410557 6.93889e-18
+DEAL::Manifold id: 4294967295
 DEAL::6.93889e-18 0.410557 -0.410557
+DEAL::Manifold id: 4294967295
 DEAL::0.00000 0.410557 0.410557
+DEAL::Manifold id: 4294967295
 DEAL::Center with manifold: 1.19113e-13 0.450853 -6.93889e-18
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: -5.55112e-17 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: -6.93889e-18 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 1.38778e-17 0.00000 0.00000
 DEAL::
 DEAL::Testing with CylindricalManifold in 3d
 DEAL::Lines on cell with center: -0.500000 0.00000 -0.500000
 DEAL::-0.500000 0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::-0.500000 -0.707107 -0.707107
-DEAL::-1.00000 -7.85046e-17 -1.00000
-DEAL::0.00000 -7.85046e-17 -1.00000
+DEAL::Manifold id: 0
+DEAL::-1.00000 -5.55112e-17 -1.00000
+DEAL::Manifold id: 0
+DEAL::0.00000 -5.55112e-17 -1.00000
+DEAL::Manifold id: 0
 DEAL::-0.500000 0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::-1.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: -0.500000 0.00000 -0.500000
 DEAL::-0.500000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.500000 -0.500000
-DEAL::-1.00000 3.33067e-15 -0.646447
-DEAL::0.00000 0.00000 -0.536612
-DEAL::-0.500000 -7.85046e-17 -1.00000
+DEAL::Manifold id: 1
+DEAL::-1.00000 3.28904e-15 -0.646447
+DEAL::Manifold id: 1
+DEAL::0.00000 3.28904e-15 -0.646447
+DEAL::Manifold id: 1
+DEAL::-0.500000 -5.55112e-17 -1.00000
+DEAL::Manifold id: 0
 DEAL::-0.500000 0.00000 -0.292893
-DEAL::Center with manifold: -0.500000 1.54737e-15 -0.591529
+DEAL::Manifold id: 1
+DEAL::Center with manifold: -0.500000 3.26822e-15 -0.646447
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 -2.22045e-16
+DEAL::Distance between cell manifold and face manifold: 0.00000 1.11022e-16 -1.11022e-16
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: -0.500000 0.500000 -6.93889e-18
 DEAL::-0.500000 0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::-0.500000 0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::-0.500000 0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 0.500000
-DEAL::-1.00000 1.00000 -7.85046e-17
+DEAL::Manifold id: 1
+DEAL::-1.00000 1.00000 -5.55112e-17
+DEAL::Manifold id: 0
 DEAL::-1.00000 0.292893 6.93889e-18
-DEAL::0.00000 1.00000 -7.85046e-17
-DEAL::0.00000 0.292893 0.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 1.00000 -5.55112e-17
+DEAL::Manifold id: 0
+DEAL::0.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: -0.500000 0.500000 -6.93889e-18
-DEAL::-0.500000 1.00000 -7.85046e-17
-DEAL::-0.500000 0.292893 0.00000
-DEAL::-1.00000 0.646447 -4.16334e-17
-DEAL::0.00000 0.536612 -9.81308e-18
+DEAL::-0.500000 1.00000 -5.55112e-17
+DEAL::Manifold id: 0
+DEAL::-0.500000 0.292893 3.46945e-18
+DEAL::Manifold id: 1
+DEAL::-1.00000 0.646447 -2.77556e-17
+DEAL::Manifold id: 1
+DEAL::0.00000 0.646447 -2.77556e-17
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.500000 0.500000
-DEAL::Center with manifold: -0.500000 0.591529 -2.77556e-17
+DEAL::Manifold id: 1
+DEAL::Center with manifold: -0.500000 0.646447 -2.77556e-17
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: -2.77556e-17 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: -0.500000 0.00000 1.38778e-17
 DEAL::-0.500000 0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
 DEAL::-1.00000 -0.292893 0.00000
-DEAL::0.00000 0.292893 0.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: -0.500000 0.00000 1.38778e-17
-DEAL::-0.500000 0.292893 0.00000
+DEAL::-0.500000 0.292893 3.46945e-18
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.00000 0.00000
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 0.00000
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::Center with manifold: -0.500000 0.00000 1.38778e-17
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: -0.500000 -0.500000 0.00000
 DEAL::-0.500000 -0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::-0.500000 -0.707107 0.707107
-DEAL::-1.00000 -1.00000 2.35514e-16
-DEAL::0.00000 -1.00000 2.35514e-16
+DEAL::Manifold id: 0
+DEAL::-1.00000 -1.00000 2.22045e-16
+DEAL::Manifold id: 0
+DEAL::0.00000 -1.00000 2.22045e-16
+DEAL::Manifold id: 0
 DEAL::-0.500000 -0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::-1.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::-1.00000 -0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: -0.500000 -0.500000 0.00000
 DEAL::-0.500000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.500000 0.500000
-DEAL::-1.00000 -0.646447 -2.99760e-15
-DEAL::0.00000 -0.536612 2.08167e-17
-DEAL::-0.500000 -1.00000 2.35514e-16
+DEAL::Manifold id: 1
+DEAL::-1.00000 -0.646447 -3.02536e-15
+DEAL::Manifold id: 1
+DEAL::0.00000 -0.646447 -3.02536e-15
+DEAL::Manifold id: 1
+DEAL::-0.500000 -1.00000 2.22045e-16
+DEAL::Manifold id: 0
 DEAL::-0.500000 -0.292893 0.00000
-DEAL::Center with manifold: -0.500000 -0.591529 -1.51268e-15
+DEAL::Manifold id: 1
+DEAL::Center with manifold: -0.500000 -0.646447 -2.92127e-15
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: -0.500000 6.93889e-18 0.500000
 DEAL::-0.500000 0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::-0.500000 0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::-0.500000 -0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 -0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 0.500000
-DEAL::-1.00000 1.57009e-16 1.00000
+DEAL::Manifold id: 1
+DEAL::-1.00000 1.66533e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::-1.00000 0.00000 0.292893
-DEAL::0.00000 1.57009e-16 1.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 1.66533e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: -0.500000 6.93889e-18 0.500000
-DEAL::-0.500000 1.57009e-16 1.00000
+DEAL::-0.500000 4.99600e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::-0.500000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::-1.00000 8.32667e-17 0.646447
-DEAL::0.00000 1.96262e-17 0.536612
+DEAL::Manifold id: 1
+DEAL::0.00000 8.32667e-17 0.646447
+DEAL::Manifold id: 1
 DEAL::-0.500000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::-0.500000 -0.500000 0.500000
-DEAL::Center with manifold: -0.500000 -1.24900e-16 0.591529
+DEAL::Manifold id: 1
+DEAL::Center with manifold: -0.500000 2.84495e-16 0.646447
+DEAL::Manifold id: 1
 DEAL::Distance between cell manifold and face manifold: -2.77556e-17 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: 0.500000 0.00000 -0.500000
 DEAL::0.500000 0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.500000 -0.707107 -0.707107
-DEAL::0.00000 -7.85046e-17 -1.00000
-DEAL::1.00000 -7.85046e-17 -1.00000
+DEAL::Manifold id: 0
+DEAL::0.00000 -5.55112e-17 -1.00000
+DEAL::Manifold id: 0
+DEAL::1.00000 -5.55112e-17 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.500000 0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::1.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: 0.500000 0.00000 -0.500000
 DEAL::0.500000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.500000 -0.500000
-DEAL::0.00000 0.00000 -0.536612
-DEAL::1.00000 0.00000 -0.536612
-DEAL::0.500000 -7.85046e-17 -1.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 3.28904e-15 -0.646447
+DEAL::Manifold id: 1
+DEAL::1.00000 3.28904e-15 -0.646447
+DEAL::Manifold id: 1
+DEAL::0.500000 -5.55112e-17 -1.00000
+DEAL::Manifold id: 0
 DEAL::0.500000 0.00000 -0.292893
-DEAL::Center with manifold: 0.500000 -2.77556e-17 -0.536612
-DEAL::Distance between cell manifold and face manifold: -2.22045e-16 -5.55112e-17 1.11022e-16
-DEAL::Distance between cell manifold and face manifold: -2.22045e-16 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 -2.22045e-16
+DEAL::Manifold id: 1
+DEAL::Center with manifold: 0.500000 3.26822e-15 -0.646447
+DEAL::Manifold id: 1
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 0.00000 1.11022e-16 -1.11022e-16
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: 0.500000 0.500000 -6.93889e-18
 DEAL::0.500000 0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.500000 0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.500000 0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.500000 0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 0.500000 0.500000
-DEAL::0.00000 1.00000 -7.85046e-17
-DEAL::0.00000 0.292893 0.00000
-DEAL::1.00000 1.00000 -7.85046e-17
-DEAL::1.00000 0.292893 0.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 1.00000 -5.55112e-17
+DEAL::Manifold id: 0
+DEAL::0.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
+DEAL::1.00000 1.00000 -5.55112e-17
+DEAL::Manifold id: 0
+DEAL::1.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: 0.500000 0.500000 -6.93889e-18
-DEAL::0.500000 1.00000 -7.85046e-17
-DEAL::0.500000 0.292893 0.00000
-DEAL::0.00000 0.536612 -9.81308e-18
-DEAL::1.00000 0.536612 -9.81308e-18
+DEAL::0.500000 1.00000 -5.55112e-17
+DEAL::Manifold id: 0
+DEAL::0.500000 0.292893 3.46945e-18
+DEAL::Manifold id: 1
+DEAL::0.00000 0.646447 -2.77556e-17
+DEAL::Manifold id: 1
+DEAL::1.00000 0.646447 -2.77556e-17
+DEAL::Manifold id: 1
 DEAL::0.500000 0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.500000 0.500000 0.500000
-DEAL::Center with manifold: 0.500000 0.536612 -2.77556e-17
+DEAL::Manifold id: 1
+DEAL::Center with manifold: 0.500000 0.646447 -2.77556e-17
+DEAL::Manifold id: 1
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: 0.00000 -0.101910 -0.0220998
-DEAL::Distance between cell manifold and face manifold: 0.00000 -0.289811 0.00774842
 DEAL::Lines on cell with center: 0.500000 0.00000 1.38778e-17
 DEAL::0.500000 0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::1.00000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.500000 0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::1.00000 0.00000 0.292893
-DEAL::0.00000 0.292893 0.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.292893 0.00000
-DEAL::1.00000 0.292893 0.00000
+DEAL::Manifold id: 1
+DEAL::1.00000 0.292893 6.93889e-18
+DEAL::Manifold id: 1
 DEAL::1.00000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: 0.500000 0.00000 1.38778e-17
-DEAL::0.500000 0.292893 0.00000
+DEAL::0.500000 0.292893 3.46945e-18
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::0.00000 0.00000 0.00000
+DEAL::Manifold id: 1
 DEAL::1.00000 0.00000 0.00000
+DEAL::Manifold id: 1
 DEAL::0.500000 0.00000 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.500000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::Center with manifold: 0.500000 0.00000 1.38778e-17
+DEAL::Manifold id: 1
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 -6.93889e-18
 DEAL::Lines on cell with center: 0.500000 -0.500000 0.00000
 DEAL::0.500000 -0.707107 -0.707107
+DEAL::Manifold id: 0
 DEAL::0.500000 -0.707107 0.707107
-DEAL::0.00000 -1.00000 2.35514e-16
-DEAL::1.00000 -1.00000 2.35514e-16
+DEAL::Manifold id: 0
+DEAL::0.00000 -1.00000 2.22045e-16
+DEAL::Manifold id: 0
+DEAL::1.00000 -1.00000 2.22045e-16
+DEAL::Manifold id: 0
 DEAL::0.500000 -0.292893 -0.292893
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::1.00000 -0.292893 0.00000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 -0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: 0.500000 -0.500000 0.00000
 DEAL::0.500000 -0.500000 -0.500000
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.500000 0.500000
-DEAL::0.00000 -0.536612 2.08167e-17
-DEAL::1.00000 -0.536612 2.08167e-17
-DEAL::0.500000 -1.00000 2.35514e-16
+DEAL::Manifold id: 1
+DEAL::0.00000 -0.646447 -3.02536e-15
+DEAL::Manifold id: 1
+DEAL::1.00000 -0.646447 -3.02536e-15
+DEAL::Manifold id: 1
+DEAL::0.500000 -1.00000 2.22045e-16
+DEAL::Manifold id: 0
 DEAL::0.500000 -0.292893 0.00000
-DEAL::Center with manifold: 0.500000 -0.536612 5.55112e-17
-DEAL::Distance between cell manifold and face manifold: 0.00000 -5.55112e-17 -5.55112e-17
-DEAL::Distance between cell manifold and face manifold: 0.00000 -1.11022e-16 -1.11022e-16
+DEAL::Manifold id: 1
+DEAL::Center with manifold: 0.500000 -0.646447 -2.92127e-15
+DEAL::Manifold id: 1
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Lines on cell with center: 0.500000 6.93889e-18 0.500000
 DEAL::0.500000 0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.500000 0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.707107 0.707107
+DEAL::Manifold id: 0
 DEAL::0.500000 -0.292893 0.292893
+DEAL::Manifold id: 1
 DEAL::0.00000 -0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::1.00000 -0.500000 0.500000
-DEAL::0.00000 1.57009e-16 1.00000
+DEAL::Manifold id: 1
+DEAL::0.00000 1.66533e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::0.00000 0.00000 0.292893
-DEAL::1.00000 1.57009e-16 1.00000
+DEAL::Manifold id: 1
+DEAL::1.00000 1.66533e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::1.00000 0.00000 0.292893
+DEAL::Manifold id: 1
 DEAL::Faces on cell with center: 0.500000 6.93889e-18 0.500000
-DEAL::0.500000 1.57009e-16 1.00000
+DEAL::0.500000 4.99600e-16 1.00000
+DEAL::Manifold id: 0
 DEAL::0.500000 0.00000 0.292893
-DEAL::0.00000 1.96262e-17 0.536612
-DEAL::1.00000 1.96262e-17 0.536612
+DEAL::Manifold id: 1
+DEAL::0.00000 8.32667e-17 0.646447
+DEAL::Manifold id: 1
+DEAL::1.00000 8.32667e-17 0.646447
+DEAL::Manifold id: 1
 DEAL::0.500000 0.500000 0.500000
+DEAL::Manifold id: 1
 DEAL::0.500000 -0.500000 0.500000
-DEAL::Center with manifold: 0.500000 5.55112e-17 0.536612
+DEAL::Manifold id: 1
+DEAL::Center with manifold: 0.500000 2.77556e-16 0.646447
+DEAL::Manifold id: 1
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
+DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
 DEAL::Distance between cell manifold and face manifold: 0.00000 0.00000 0.00000
-DEAL::Distance between cell manifold and face manifold: 0.00000 0.0220998 -0.101910
-DEAL::Distance between cell manifold and face manifold: 0.00000 -0.00774842 -0.289811
 DEAL::

diff --git a/tests/mappings/mapping.cc b/tests/mappings/mapping.cc
index fb2963b6d99cec01ea04792515a321725ef4d40a..80f7c2e17c1bf52c87165e5f2a85b58d5bd83d6b 100644 (file)
--- a/tests/mappings/mapping.cc
+++ b/tests/mappings/mapping.cc
@@ -25,6 +25,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
 #include <deal.II/grid/manifold_lib.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/fe/mapping_cartesian.h>
 #include <deal.II/fe/mapping_q1.h>
 #include <deal.II/fe/mapping_q.h>
@@ -304,8 +305,8 @@ void create_triangulations(std::vector<Triangulation<2> *> &tria_ptr,
       v3(1) = 3.;
       tria->set_manifold(1,*boundary1);
       tria->set_manifold(2,*boundary2);
-      tria->begin_active()->face(1)->set_boundary_id(1);
-      tria->begin_active()->face(3)->set_boundary_id(2);
+      tria->begin_active()->face(1)->set_manifold_id(1);
+      tria->begin_active()->face(3)->set_manifold_id(2);
       double pi=std::acos(-1.);
       double alpha=2*std::atan(0.5);
       exact_areas.push_back(4+pi-2.5*(alpha-std::sin(alpha)));
@@ -342,7 +343,7 @@ void create_triangulations(std::vector<Triangulation<2> *> &tria_ptr,
       v3(0) = 0.5;
       v3(1) = 1.5;
       tria->set_manifold(1,*boundary1);
-      tria->begin_active()->face(1)->set_boundary_id(1);
+      tria->begin_active()->face(1)->set_manifold_id(1);
       exact_areas.push_back(0.);
       for (unsigned int i=0; i<=4; ++i)
         show[4][i]=1;
@@ -397,7 +398,7 @@ void create_triangulations(std::vector<Triangulation<3> *> &tria_ptr,
       tria_ptr.push_back(tria);
       GridGenerator::hyper_cube(*tria, 1., 3.);
       tria->set_manifold(1,*boundary1);
-      tria->begin_active()->face(1)->set_boundary_id(1);
+      tria->begin_active()->face(1)->set_manifold_id(1);
       exact_areas.push_back(8.+pi/3*h*h*(3*r-h));
     }
 

diff --git a/tests/mappings/mapping_c1.cc b/tests/mappings/mapping_c1.cc
index 69f5dcaa5b015f99fb7fc04d1c21b78be3674824..5b9712c34ab6b5ad379cc8a39fc887ec8ae48ac0 100644 (file)
--- a/tests/mappings/mapping_c1.cc
+++ b/tests/mappings/mapping_c1.cc
@@ -43,6 +43,7 @@ int main ()
   SphericalManifold<2> circle;
   Triangulation<2> tria;
   GridGenerator::hyper_cube(tria, -1, 1);
+  tria.set_all_manifold_ids_on_boundary(0);
   tria.set_manifold (0, circle);
 
 

diff --git a/tests/mappings/mapping_manifold_04.cc b/tests/mappings/mapping_manifold_04.cc
index fed66a09a650093b123e75f09887f43d2112837e..ac0833849f8e2ab409515ebbe5ef4d6476f5775e 100644 (file)
--- a/tests/mappings/mapping_manifold_04.cc
+++ b/tests/mappings/mapping_manifold_04.cc
@@ -45,7 +45,7 @@ void test()
                outer_radius = 1.0;
   GridGenerator::hyper_shell (triangulation,
                               center, inner_radius, outer_radius);
-
+  triangulation.reset_manifold(0);
   const QGauss<spacedim> quad(3);
 
   FE_Q<dim,spacedim> fe(1);

diff --git a/tests/mappings/mapping_manifold_05.cc b/tests/mappings/mapping_manifold_05.cc
index ad9efbb8715a45bc39f12fa509f6c076b8a927a1..44dd54f1c325183ffc8f0399a4bbe77436c43df8 100644 (file)
--- a/tests/mappings/mapping_manifold_05.cc
+++ b/tests/mappings/mapping_manifold_05.cc
@@ -46,6 +46,7 @@ void test()
                outer_radius = 1.0;
   GridGenerator::hyper_shell (triangulation,
                               center, inner_radius, outer_radius);
+  triangulation.reset_manifold(0);
 
   const QGauss<spacedim-1> quad(3);
 

diff --git a/tests/mpi/map_dofs_to_support_points.cc b/tests/mpi/map_dofs_to_support_points.cc
index d83d7b9757fef795c4abac8bcd40c27f6a63924a..3fb8d5c9e25d351c5e236db6c8e43bb3f6ce0f4a 100644 (file)
--- a/tests/mpi/map_dofs_to_support_points.cc
+++ b/tests/mpi/map_dofs_to_support_points.cc
@@ -33,6 +33,7 @@ void test()
   parallel::distributed::Triangulation<dim> tr(MPI_COMM_WORLD);
 
   GridGenerator::hyper_ball(tr);
+  tr.reset_manifold(0);
   tr.refine_global (1);
 
   const FE_Q<dim> fe(1);

diff --git a/tests/mpi/no_flux_constraints.cc b/tests/mpi/no_flux_constraints.cc
index f98d61929d1a3784c5024b39208dd60354a7fd98..5e5c588f1f5e55d51ce9693da2cd4bd9040d7ce4 100644 (file)
--- a/tests/mpi/no_flux_constraints.cc
+++ b/tests/mpi/no_flux_constraints.cc
@@ -56,6 +56,7 @@ void test()
                               R0,
                               R1,
                               96, true);
+  triangulation.reset_manifold(0);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator

diff --git a/tests/mpi/no_flux_constraints_02.cc b/tests/mpi/no_flux_constraints_02.cc
index 1e205602a973cdff33f6b1ec877fa60b62e29985..3a784526cc375692d6b51f67f94ea00714970900 100644 (file)
--- a/tests/mpi/no_flux_constraints_02.cc
+++ b/tests/mpi/no_flux_constraints_02.cc
@@ -57,6 +57,7 @@ void test()
                               R0,
                               R1,
                               96, true);
+  triangulation.reset_manifold(0);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator

diff --git a/tests/mpi/normal_flux_constraints.cc b/tests/mpi/normal_flux_constraints.cc
index f0c4020759393ff5ae0dac25f15d66824cb7d066..66323744264c49819456a43b4c0d9bae9b7b2f79 100644 (file)
--- a/tests/mpi/normal_flux_constraints.cc
+++ b/tests/mpi/normal_flux_constraints.cc
@@ -55,6 +55,7 @@ void test()
                               R0,
                               R1,
                               96, true);
+  triangulation.reset_manifold(0);
   triangulation.refine_global (2);
 
   for (typename Triangulation<dim>::active_cell_iterator

diff --git a/tests/mpi/p4est_2d_constraintmatrix_03.cc b/tests/mpi/p4est_2d_constraintmatrix_03.cc
index 697e929481bb7ecc33bf14ec235c9038799abeeb..15fd58c6a8f519f436f2cc128b6fe0c7f02f5480 100644 (file)
--- a/tests/mpi/p4est_2d_constraintmatrix_03.cc
+++ b/tests/mpi/p4est_2d_constraintmatrix_03.cc
@@ -29,6 +29,7 @@
 #include <deal.II/distributed/tria.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/dofs/dof_tools.h>
@@ -62,6 +63,7 @@ void test()
                               R1,
                               12,
                               true);
+  GridTools::copy_boundary_to_manifold_id(tr);
 
   static SphericalManifold<dim> boundary;
   tr.set_manifold (0, boundary);

diff --git a/tests/mpi/p4est_2d_constraintmatrix_04.cc b/tests/mpi/p4est_2d_constraintmatrix_04.cc
index c5f3838fdbecf174325ffdf0483a3bc83fba21b3..7db642ccf16827d0ba2ecece971df1745ae30712 100644 (file)
--- a/tests/mpi/p4est_2d_constraintmatrix_04.cc
+++ b/tests/mpi/p4est_2d_constraintmatrix_04.cc
@@ -143,6 +143,7 @@ void test()
                               R1,
                               12,
                               true);
+  tr.reset_all_manifolds();
   static SphericalManifold<dim> boundary;
 //  tr.set_manifold (0, boundary);
 // tr.set_manifold (1, boundary);

diff --git a/tests/mpi/p4est_3d_constraintmatrix_03.cc b/tests/mpi/p4est_3d_constraintmatrix_03.cc
index ccb4c8a3f68b52514755f688b724772ebac022bb..9db1cb99a94a51a30ff23434d6950a28001bee0a 100644 (file)
--- a/tests/mpi/p4est_3d_constraintmatrix_03.cc
+++ b/tests/mpi/p4est_3d_constraintmatrix_03.cc
@@ -142,7 +142,9 @@ void test()
                               R1,
                               12,
                               true);
-  static SphericalManifold<dim> boundary;
+  tr.reset_manifold(0);
+
+//  static SphericalManifold<dim> boundary;
 //  tr.set_manifold (0, boundary);
 //tr.set_manifold (1, boundary);
 

diff --git a/tests/mpi/p4est_3d_refine_04.cc b/tests/mpi/p4est_3d_refine_04.cc
index 12d6f0825144f13052ba4b82771c412c9176d388..a78d222a8aad10373dba898bf7a798cbb4a3a99f 100644 (file)
--- a/tests/mpi/p4est_3d_refine_04.cc
+++ b/tests/mpi/p4est_3d_refine_04.cc
@@ -56,6 +56,7 @@ void test()
                                   R1,
                                   (dim==3) ? 96 : 12,
                                   true);
+      triangulation.reset_manifold(0);
       triangulation.refine_global (2);
 
       // then flag all cells for

diff --git a/tests/mpi/periodicity_02.cc b/tests/mpi/periodicity_02.cc
index 5aa9308d701f3fd988bda99b1e4dc9cf3e4deb8f..1185f82ed0f49ccf6eeb7bf4abec7702c315f350 100644 (file)
--- a/tests/mpi/periodicity_02.cc
+++ b/tests/mpi/periodicity_02.cc
@@ -719,6 +719,8 @@ namespace Step22
                                         outer_radius,
                                         0,
                                         true);
+    triangulation.set_all_manifold_ids(numbers::flat_manifold_id);
+    GridTools::copy_boundary_to_manifold_id(triangulation);
 
 #ifdef PERIODIC
     std::vector<GridTools::PeriodicFacePair<typename parallel::distributed::Triangulation<dim>::cell_iterator> >

diff --git a/tests/mpi/renumber_z_order_02.cc b/tests/mpi/renumber_z_order_02.cc
index e1e9b68594b369c436bf4366dea0819ddf98f5ca..b06ddeaca24eaf67048815d1901ea8106f13cee4 100644 (file)
--- a/tests/mpi/renumber_z_order_02.cc
+++ b/tests/mpi/renumber_z_order_02.cc
@@ -53,6 +53,7 @@ void test()
                               0.3,
                               1.0,
                               12, true);
+  tr.reset_all_manifolds();
 //  GridGenerator::hyper_cube (tr, -1.0, 1.0);
   tr.refine_global (1);
 

diff --git a/tests/numerics/create_laplace_matrix_01.cc b/tests/numerics/create_laplace_matrix_01.cc
index a154599250234c643319437594303ea9337e94b6..f508f2d7391a3768653770ed0a368125f1657b0d 100644 (file)
--- a/tests/numerics/create_laplace_matrix_01.cc
+++ b/tests/numerics/create_laplace_matrix_01.cc
@@ -55,7 +55,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_01b.cc b/tests/numerics/create_laplace_matrix_01b.cc
index aa5916b99b72b0f2243a310595f60ce412d38dd3..1dba5e56197aa27b781e61fca0182a0f0be382f8 100644 (file)
--- a/tests/numerics/create_laplace_matrix_01b.cc
+++ b/tests/numerics/create_laplace_matrix_01b.cc
@@ -56,7 +56,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_02.cc b/tests/numerics/create_laplace_matrix_02.cc
index 41dce03a2383df1d9c11e17dce5de120ef883151..dc3d08d785ee8aa768504c31179956c6a8581792 100644 (file)
--- a/tests/numerics/create_laplace_matrix_02.cc
+++ b/tests/numerics/create_laplace_matrix_02.cc
@@ -54,7 +54,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_02b.cc b/tests/numerics/create_laplace_matrix_02b.cc
index cc4209b0753ab70ec99b09fc6d8c9bf634db35e9..7959609cb7230c4dbb69ae397a807c63703c57be 100644 (file)
--- a/tests/numerics/create_laplace_matrix_02b.cc
+++ b/tests/numerics/create_laplace_matrix_02b.cc
@@ -55,7 +55,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_03.cc b/tests/numerics/create_laplace_matrix_03.cc
index 4a20d1c1180c361e6e9e20c06bcc8983d2d87783..769074088479d1752c2321f66eda72f25f4c8b46 100644 (file)
--- a/tests/numerics/create_laplace_matrix_03.cc
+++ b/tests/numerics/create_laplace_matrix_03.cc
@@ -54,7 +54,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_03b.cc b/tests/numerics/create_laplace_matrix_03b.cc
index 383df4bd1bec64a56402ab7a9f9e1e1b71835e64..6ad024fce545b9eb62986d73a8ad91bc6cbae037 100644 (file)
--- a/tests/numerics/create_laplace_matrix_03b.cc
+++ b/tests/numerics/create_laplace_matrix_03b.cc
@@ -55,7 +55,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_04.cc b/tests/numerics/create_laplace_matrix_04.cc
index 0eec67b2370d2a1834e8f05394496917a95c551b..ab390afd7d972ea2287b082a61c347db70dd4b15 100644 (file)
--- a/tests/numerics/create_laplace_matrix_04.cc
+++ b/tests/numerics/create_laplace_matrix_04.cc
@@ -53,7 +53,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_laplace_matrix_04b.cc b/tests/numerics/create_laplace_matrix_04b.cc
index 874b801876bf82537b9c670248ff2c758c66dead..d63c7231ef7b741fe6094cf47a9597ceed0621e8 100644 (file)
--- a/tests/numerics/create_laplace_matrix_04b.cc
+++ b/tests/numerics/create_laplace_matrix_04b.cc
@@ -54,7 +54,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_01.cc b/tests/numerics/create_mass_matrix_01.cc
index 3ccdb1701ec138fd58daf163ed21ec9041031589..35588c7850b8db9a4786f84cbd725abccb061be2 100644 (file)
--- a/tests/numerics/create_mass_matrix_01.cc
+++ b/tests/numerics/create_mass_matrix_01.cc
@@ -55,7 +55,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_01b.cc b/tests/numerics/create_mass_matrix_01b.cc
index d670686bdbb9d012925be8593b46544e76b8fc23..ee586539f977ee95753b763d1108101410148d57 100644 (file)
--- a/tests/numerics/create_mass_matrix_01b.cc
+++ b/tests/numerics/create_mass_matrix_01b.cc
@@ -56,7 +56,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_02.cc b/tests/numerics/create_mass_matrix_02.cc
index b3f3229aa8e2b88320294e4417d9dff9f99bfb70..808c5f6aa2fa95a20c2e31a28c8b5f2213511bfa 100644 (file)
--- a/tests/numerics/create_mass_matrix_02.cc
+++ b/tests/numerics/create_mass_matrix_02.cc
@@ -54,7 +54,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_02b.cc b/tests/numerics/create_mass_matrix_02b.cc
index afdf5cd675a55187c2b8a8781e6e3698dc7a8b67..9fec860c8ece3f9c51de9c0d29587462173a3dd1 100644 (file)
--- a/tests/numerics/create_mass_matrix_02b.cc
+++ b/tests/numerics/create_mass_matrix_02b.cc
@@ -55,7 +55,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_03.cc b/tests/numerics/create_mass_matrix_03.cc
index 5e0e71a2ee0f519f4d5d8cfba01c9286c13fb971..33ab57ed188cd8a73e137cc37c2181339f108328 100644 (file)
--- a/tests/numerics/create_mass_matrix_03.cc
+++ b/tests/numerics/create_mass_matrix_03.cc
@@ -54,7 +54,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_03b.cc b/tests/numerics/create_mass_matrix_03b.cc
index 08a2e0004cd19e82389d8c2cac95e30443c08f09..7c431071853ed57d46cadfc215cce6906fc84ae1 100644 (file)
--- a/tests/numerics/create_mass_matrix_03b.cc
+++ b/tests/numerics/create_mass_matrix_03b.cc
@@ -55,7 +55,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_04.cc b/tests/numerics/create_mass_matrix_04.cc
index 6bf851ae1c19d49775c8f8bb9db2c1a7328aa9b8..e35d1db6b428226faac8e9548a1c032527380d08 100644 (file)
--- a/tests/numerics/create_mass_matrix_04.cc
+++ b/tests/numerics/create_mass_matrix_04.cc
@@ -53,7 +53,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_04b.cc b/tests/numerics/create_mass_matrix_04b.cc
index 65d8f09cd90e3f8f1637fe47bc6b82f95e116a9f..7289f4db3c5a6ef52db0779a357ad6595d90f136 100644 (file)
--- a/tests/numerics/create_mass_matrix_04b.cc
+++ b/tests/numerics/create_mass_matrix_04b.cc
@@ -54,7 +54,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/create_mass_matrix_05.cc b/tests/numerics/create_mass_matrix_05.cc
index 0ca32353642a50414e25d73282d8df49474765c5..2f4660a586dd9d3466442610b5f55ecff0f8c1b3 100644 (file)
--- a/tests/numerics/create_mass_matrix_05.cc
+++ b/tests/numerics/create_mass_matrix_05.cc
@@ -46,7 +46,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/derivative_approximation_01.cc b/tests/numerics/derivative_approximation_01.cc
index 4f656d5513237c5e1b1e60ce85beab04a8f974b6..8a8b76a479543ed392fd13d1fcb0ea3037046355 100644 (file)
--- a/tests/numerics/derivative_approximation_01.cc
+++ b/tests/numerics/derivative_approximation_01.cc
@@ -42,7 +42,10 @@ check ()
 
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/derivatives.cc b/tests/numerics/derivatives.cc
index 51009385d8c7070d47081586226bbf8a907ab7e1..833e05a74503c5832b0c20f457d042a76d59938a 100644 (file)
--- a/tests/numerics/derivatives.cc
+++ b/tests/numerics/derivatives.cc
@@ -54,7 +54,10 @@ check (const unsigned int level,
 
   DoFHandler<dim> dof(tr);
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
 

diff --git a/tests/numerics/error_estimator.cc b/tests/numerics/error_estimator.cc
index 02bb273266423b29bc45ab1fc90776f9643f7737..915ffe09e88ccea395876d7b9052250c5919bdfc 100644 (file)
--- a/tests/numerics/error_estimator.cc
+++ b/tests/numerics/error_estimator.cc
@@ -74,7 +74,10 @@ check ()
 
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/no_flux_08.cc b/tests/numerics/no_flux_08.cc
index 8f4b134ea961b79894925962ea388cd9db66ea15..2c0357bf29c9c9189b705b208ead0bd9293eaf4a 100644 (file)
--- a/tests/numerics/no_flux_08.cc
+++ b/tests/numerics/no_flux_08.cc
@@ -49,6 +49,7 @@ check ()
                                    Point<dim>(),
                                    0.5,
                                    1,0);
+  tr.reset_manifold(0);
 
   ConstraintMatrix cm;
   MappingQ<dim> mapping(1);

diff --git a/tests/numerics/no_flux_09.cc b/tests/numerics/no_flux_09.cc
index bea8c7e501b30bd4418694160ab0a53e4d7f8b72..4e26587df88bee1cd904855f5f5a5ccde96fe93b 100644 (file)
--- a/tests/numerics/no_flux_09.cc
+++ b/tests/numerics/no_flux_09.cc
@@ -47,6 +47,7 @@ check ()
                                       Point<dim>(),
                                       0.5, 1.0,
                                       3, true);
+  tr.reset_manifold(0);
 
   ConstraintMatrix cm;
   MappingQ<dim> mapping(1);

diff --git a/tests/numerics/no_flux_10.cc b/tests/numerics/no_flux_10.cc
index 49ba02516a2dbe7f4b01d2d43a1ec59c3d222ea1..962dbd3498e2871049bcf8935077aaf17ce5727d 100644 (file)
--- a/tests/numerics/no_flux_10.cc
+++ b/tests/numerics/no_flux_10.cc
@@ -35,6 +35,7 @@
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/grid/manifold_lib.h>
 #include <deal.II/grid/grid_out.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_handler.h>
 #include <deal.II/fe/fe_system.h>
 #include <deal.II/fe/mapping_q.h>
@@ -169,6 +170,7 @@ void run()
                          Point<dim>(),
                          0.5,
                          1.0);
+  GridTools::copy_boundary_to_manifold_id(triangulation);
 
   static SphericalManifold<dim> boundary((Point<dim>()));
   triangulation.set_manifold (0, boundary);

diff --git a/tests/numerics/no_flux_12.cc b/tests/numerics/no_flux_12.cc
index db43580e347b2f7195b846e224baba172f034d05..5648ec6e727142e6b6715945796fbb9837743f27 100644 (file)
--- a/tests/numerics/no_flux_12.cc
+++ b/tests/numerics/no_flux_12.cc
@@ -73,6 +73,7 @@ void run()
 
   // indicator 6 = cylinder
   GridGenerator::hyper_cube_with_cylindrical_hole (tria, 0.25, 0.5, 0.5, 1, true);
+  tria.reset_manifold(0);
 
   /*  std::string filename = "Mesh.eps";
   std::ofstream output (filename.c_str());

diff --git a/tests/numerics/no_flux_14.cc b/tests/numerics/no_flux_14.cc
index fe1690edddb5a2c37dc50454d3e4160f6023cb5a..b0f34aa754a375fa768d9b6dfc4d11e6d072f989 100644 (file)
--- a/tests/numerics/no_flux_14.cc
+++ b/tests/numerics/no_flux_14.cc
@@ -40,6 +40,7 @@ check ()
                                       Point<dim>(),
                                       0.5, 1.0,
                                       3, true);
+  tr.reset_manifold(0);
 
   ConstraintMatrix cm;
   MappingQ<dim> mapping(1);

diff --git a/tests/numerics/project_boundary_01.cc b/tests/numerics/project_boundary_01.cc
index a9fdb476b72ae3565754a4434169734aae7dc700..36ab6999d24d5a50d3601a2f5401ff5ca1aeddaa 100644 (file)
--- a/tests/numerics/project_boundary_01.cc
+++ b/tests/numerics/project_boundary_01.cc
@@ -25,6 +25,7 @@
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/lac/vector.h>
 #include <deal.II/grid/tria.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/manifold_lib.h>
 #include <deal.II/dofs/dof_accessor.h>
@@ -110,6 +111,8 @@ check ()
     }
   else
     GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
+  GridTools::copy_boundary_to_manifold_id(tr);
+
   static const SphericalManifold<dim> boundary;
   if (dim != 1)
     {

diff --git a/tests/numerics/project_common.h b/tests/numerics/project_common.h
index 4abde25cd555887cba6c0c1f39606e381d1c2905..592aa831378ea2f74cbb8e573945ff98430a2493 100644 (file)
--- a/tests/numerics/project_common.h
+++ b/tests/numerics/project_common.h
@@ -204,6 +204,7 @@ void test_with_wrong_face_orientation (const FiniteElement<dim> &fe,
     {
       Triangulation<dim>     triangulation;
       GridGenerator::hyper_ball (triangulation);
+      triangulation.set_manifold(0);
       typename Triangulation<dim>::active_cell_iterator
       cell = triangulation.begin_active();
       std::advance (cell, i);

diff --git a/tests/numerics/project_parallel_common.h b/tests/numerics/project_parallel_common.h
index 39219ca9e4cfac161a422139e58af793f6bc3135..561e17abfe8e0bc7ae19b951e063bcb311528971 100644 (file)
--- a/tests/numerics/project_parallel_common.h
+++ b/tests/numerics/project_parallel_common.h
@@ -216,6 +216,7 @@ void test_with_wrong_face_orientation (const FiniteElement<dim> &fe,
     {
       parallel::distributed::Triangulation<dim> triangulation(MPI_COMM_WORLD);
       GridGenerator::hyper_ball (triangulation);
+      triangulation.reset_manifold(0);
       typename parallel::distributed::Triangulation<dim>::active_cell_iterator
       cell = triangulation.begin_active();
       std::advance (cell, i);

diff --git a/tests/numerics/vectors_boundary_rhs_01.cc b/tests/numerics/vectors_boundary_rhs_01.cc
index 45a7758c8afbde27c7023bbc423cb49db4e32bc8..8977b6824faa114cf6b57723a6e4ef25c4be85ba 100644 (file)
--- a/tests/numerics/vectors_boundary_rhs_01.cc
+++ b/tests/numerics/vectors_boundary_rhs_01.cc
@@ -67,7 +67,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/vectors_boundary_rhs_02.cc b/tests/numerics/vectors_boundary_rhs_02.cc
index d559f3c0350e4d8129627dfe86b4a85e5d35d3b3..f35b14608bc2fd47422cf032fabe7bb95c63b4b5 100644 (file)
--- a/tests/numerics/vectors_boundary_rhs_02.cc
+++ b/tests/numerics/vectors_boundary_rhs_02.cc
@@ -70,7 +70,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/vectors_boundary_rhs_03.cc b/tests/numerics/vectors_boundary_rhs_03.cc
index 6b6c53ae0bb29a4d6781b8dcb397ced72d070801..97249bc9afcdb95daa8d549e83e7fa03db0d688a 100644 (file)
--- a/tests/numerics/vectors_boundary_rhs_03.cc
+++ b/tests/numerics/vectors_boundary_rhs_03.cc
@@ -68,7 +68,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/vectors_point_source_01.cc b/tests/numerics/vectors_point_source_01.cc
index 032ea0eca211a72cf8aa2c7a6564c2f074eed160..e2692288d9089ff175c24f557222369150eb2c34 100644 (file)
--- a/tests/numerics/vectors_point_source_01.cc
+++ b/tests/numerics/vectors_point_source_01.cc
@@ -69,7 +69,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/vectors_rhs_01.cc b/tests/numerics/vectors_rhs_01.cc
index a1f023c7f599c1848fb8e79dbddc9c4e4cefd01f..fa6d996cd1ab49c577c48a4cb01ada7d5fc71ca4 100644 (file)
--- a/tests/numerics/vectors_rhs_01.cc
+++ b/tests/numerics/vectors_rhs_01.cc
@@ -69,7 +69,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/numerics/vectors_rhs_03.cc b/tests/numerics/vectors_rhs_03.cc
index f359b2de8eae603b3d21502ae0e064bcffbcd1f7..45ccadb25d254853ee88e1841aa9e0412be5b183 100644 (file)
--- a/tests/numerics/vectors_rhs_03.cc
+++ b/tests/numerics/vectors_rhs_03.cc
@@ -69,7 +69,10 @@ check ()
 {
   Triangulation<dim> tr;
   if (dim==2)
-    GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+    {
+      GridGenerator::hyper_ball(tr, Point<dim>(), 1);
+      tr.reset_manifold(0);
+    }
   else
     GridGenerator::hyper_cube(tr, -1,1);
   tr.refine_global (1);

diff --git a/tests/vector_tools/boundaries.cc b/tests/vector_tools/boundaries.cc
index 237dd13de6aa42ab3c49140c70d3b49d1d450d22..ee837ff881ac101b1930ae7ca12693f862560124 100644 (file)
--- a/tests/vector_tools/boundaries.cc
+++ b/tests/vector_tools/boundaries.cc
@@ -26,6 +26,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/manifold_lib.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -87,6 +88,7 @@ check ()
     }
   else
     GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
+  GridTools::copy_boundary_to_manifold_id(tr);
   static const SphericalManifold<dim> boundary;
   if (dim != 1)
     {

diff --git a/tests/vector_tools/boundaries_complex.cc b/tests/vector_tools/boundaries_complex.cc
index 037f83ea3f8e01c53fb94333b7d062703a914784..e181aa07dda82e440beab0a9a78cb9fd1f636bd4 100644 (file)
--- a/tests/vector_tools/boundaries_complex.cc
+++ b/tests/vector_tools/boundaries_complex.cc
@@ -28,6 +28,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/manifold_lib.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -90,6 +91,7 @@ check ()
     }
   else
     GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
+  GridTools::copy_boundary_to_manifold_id(tr);
   static const SphericalManifold<dim> boundary;
   if (dim != 1)
     {

diff --git a/tests/vector_tools/boundaries_complex_hp.cc b/tests/vector_tools/boundaries_complex_hp.cc
index d0c199113c58216d65bef80bdd8c43a21b297d2a..5491a37df413f568609fb09fd21a1e6a2a98df31 100644 (file)
--- a/tests/vector_tools/boundaries_complex_hp.cc
+++ b/tests/vector_tools/boundaries_complex_hp.cc
@@ -26,6 +26,7 @@
 #include <deal.II/grid/tria.h>
 #include <deal.II/grid/tria_iterator.h>
 #include <deal.II/grid/manifold_lib.h>
+#include <deal.II/grid/grid_tools.h>
 #include <deal.II/dofs/dof_accessor.h>
 #include <deal.II/grid/grid_generator.h>
 #include <deal.II/dofs/dof_handler.h>
@@ -88,6 +89,7 @@ check ()
     }
   else
     GridGenerator::hyper_cube(tr, -1./std::sqrt(static_cast<double>(dim)),1./std::sqrt(static_cast<double>(dim)));
+  GridTools::copy_boundary_to_manifold_id(tr);
   static const SphericalManifold<dim> boundary;
   if (dim != 1)
     {