]> https://gitweb.dealii.org/ - dealii.git/commitdiff
Switch to parallel tria
authorPeter Munch <peterrmuench@gmail.com>
Sat, 22 May 2021 16:20:57 +0000 (18:20 +0200)
committerPeter Munch <peterrmuench@gmail.com>
Sat, 22 May 2021 16:20:57 +0000 (18:20 +0200)
tests/multigrid-global-coarsening/multigrid_p_02.cc
tests/multigrid-global-coarsening/multigrid_p_02.mpirun=1.with_trilinos=true.output

index 2a3b197cff604004102249a9440588bb6d42b2c3..13111dc3782f95dc7f16daece43ce65a520cdf8e 100644 (file)
@@ -27,8 +27,10 @@ test(const unsigned int n_refinements, const unsigned int fe_degree_fine)
 {
   using VectorType = LinearAlgebra::distributed::Vector<Number>;
 
-  Triangulation<dim> tria(
-    Triangulation<dim>::MeshSmoothing::limit_level_difference_at_vertices);
+  parallel::distributed::Triangulation<dim> tria(
+    MPI_COMM_WORLD,
+    Triangulation<dim>::MeshSmoothing::none,
+    parallel::distributed::Triangulation<dim>::construct_multigrid_hierarchy);
   GridGenerator::subdivided_hyper_cube(
     tria, 2 * Utilities::pow<unsigned int>(2, n_refinements));
 
@@ -88,6 +90,8 @@ test(const unsigned int n_refinements, const unsigned int fe_degree_fine)
         mg_constrained_dofs.get_boundary_indices(0 /*level*/));
       constraint.close();
 
+      constraint.print(std::cout);
+
       // set up operator
       op.reinit(*mapping, dof_handler, *quad, constraint, 0 /*level*/);
     }
index 24dc959efecaf8269900fd2ffa588307dbdc2377..e460d32df9d9359be120c36b0c8b0f068f30ca70 100644 (file)
@@ -1,10 +1,10 @@
 
-DEAL:0::2 2 2 4
-DEAL:0::2 3 2 6
-DEAL:0::2 4 2 8
-DEAL:0::2 2 3 6
-DEAL:0::2 3 3 11
-DEAL:0::2 4 3 17
-DEAL:0::2 2 4 13
-DEAL:0::2 3 4 22
-DEAL:0::2 4 4 32
+DEAL:0::2 2 2 3
+DEAL:0::2 3 2 3
+DEAL:0::2 4 2 3
+DEAL:0::2 2 3 3
+DEAL:0::2 3 3 3
+DEAL:0::2 4 3 3
+DEAL:0::2 2 4 3
+DEAL:0::2 3 4 3
+DEAL:0::2 4 4 3

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.