Merge pull request #15486 from kronbichler/test_robustness2

MatrixFree inverse mass matrix: Make tests more robust
dealii · Jun 26, 2023 · 1ffb60c · 1ffb60c
2 parents 27f5a85 + 9892dc5
commit 1ffb60c
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 14 deletions.
diff --git a/tests/matrix_free/inverse_mass_10.cc b/tests/matrix_free/inverse_mass_10.cc
@@ -150,7 +150,7 @@ class MatrixFreeTest
                   make_vectorized_array(random_value<Number>(1.0e-2, 1.0e-1));
                 // make diagonal dominant because inverse is needed
                 if (d1 == d2)
-                  dyadic_coefficients[q][d1][d2] *= 1.0e2;
+                  dyadic_coefficients[q][d1][d2] += 5.0;
               }
           }
         inverse_dyadic_coefficients[q] = invert(dyadic_coefficients[q]);
@@ -220,12 +220,6 @@ test()
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    if (cell->center().norm() < 1e-8)
-      cell->set_refine_flag();
-  tria.execute_coarsening_and_refinement();
 
   const unsigned int degree = fe_degree;
   FESystem<dim>      fe(FE_DGQ<dim>(degree), dim);

diff --git a/tests/matrix_free/inverse_mass_11.cc b/tests/matrix_free/inverse_mass_11.cc
@@ -151,7 +151,7 @@ class MatrixFreeTest
                   make_vectorized_array(random_value<Number>(1.0e-2, 1.0e-1));
                 // make diagonal dominant because inverse is needed
                 if (d1 == d2)
-                  dyadic_coefficients[q][d1][d2] *= 1.0e2;
+                  dyadic_coefficients[q][d1][d2] += 5.0;
               }
           }
         inverse_dyadic_coefficients[q] = invert(dyadic_coefficients[q]);
@@ -221,12 +221,6 @@ test()
   tria.begin(tria.n_levels() - 1)->set_refine_flag();
   tria.last()->set_refine_flag();
   tria.execute_coarsening_and_refinement();
-  typename Triangulation<dim>::active_cell_iterator cell = tria.begin_active(),
-                                                    endc = tria.end();
-  for (; cell != endc; ++cell)
-    if (cell->center().norm() < 1e-8)
-      cell->set_refine_flag();
-  tria.execute_coarsening_and_refinement();
 
   const unsigned int degree = fe_degree;
   FESystem<dim>      fe(FE_DGQ<dim>(degree), dim);