]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Modify shell radius to make test fail.
authorkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 24 Aug 2011 12:14:21 +0000 (12:14 +0000)
committerkronbichler <kronbichler@0785d39b-7218-0410-832d-ea1e28bc413d>
Wed, 24 Aug 2011 12:14:21 +0000 (12:14 +0000)
git-svn-id: https://svn.dealii.org/trunk@24175 0785d39b-7218-0410-832d-ea1e28bc413d

tests/mpi/no_flux_constraints.cc
tests/mpi/no_flux_constraints/ncpu_9/cmp/generic

index 671c2bf7b3d14e7809b0d8a491ef77e9215a35a8..f66155fb7fd7a558edd3d02eb838c92d85ac28d1 100644 (file)
@@ -46,11 +46,14 @@ void test()
   parallel::distributed::Triangulation<dim> triangulation(MPI_COMM_WORLD);
 
                                // create hypershell mesh and refine some
-                               // cells
+                               // cells. use some large numbers which make
+                               // round-off influence more pronounced
+  const double R0      = 6371000.-2890000.;     /* m          */
+  const double R1      = 6371000.-  35000.;     /* m          */
   GridGenerator::hyper_shell (triangulation,
                              Point<dim>(),
-                             0.5,
-                             1.,
+                             R0,
+                             R1,
                              96, true);
   triangulation.refine_global (2);
 
@@ -58,7 +61,7 @@ void test()
         cell = triangulation.begin_active();
        cell != triangulation.end(); ++cell)
     if (!cell->is_ghost() && !cell->is_artificial())
-      if (cell->center()[2] > 0.75)
+      if (cell->center()[2] > 0.75 * R1)
        {
          cell->set_refine_flag();
          for (unsigned int c=0; c<8; ++c)
@@ -90,7 +93,7 @@ void test()
 
                                // create FE_System and fill in no-normal flux
                                // conditions on boundary 1 (outer)
-  static const FESystem<dim> fe(FE_Q<dim> (2), dim);
+  static const FESystem<dim> fe(FE_Q<dim> (1), dim);
   DoFHandler<dim> dofh(triangulation);
   dofh.distribute_dofs (fe);
   if (myid == 0)
@@ -151,10 +154,8 @@ void test()
                                // now check that no entries were generated
                                // for constrained entries on the locally
                                // owned range.
-  const std::pair<unsigned int,unsigned int> range =
-    vector.local_range();
-  for (unsigned int i=range.first; i<range.second;
-       ++i)
+  const std::pair<unsigned int,unsigned int> range = vector.local_range();
+  for (unsigned int i=range.first; i<range.second; ++i)
     if (constraints.is_constrained(i))
       Assert (vector(i)==0, ExcInternalError());
 
index 62e305614465f2fb25b0ee3dc2767d7782cb4e6a..ba60b2ea3d4f90eeed32aee3f4435a64beb8e6cb 100644 (file)
@@ -1,14 +1,14 @@
 
 DEAL:0:3d::# flagged cells = 224
 DEAL:0:3d::#cells = 7712
-DEAL:0:3d::#dofs = 202422
-DEAL:0:3d::#constraints on 0: 607
-DEAL:0:3d::#constraints on 1: 4388
-DEAL:0:3d::#constraints on 2: 2518
-DEAL:0:3d::#constraints on 3: 4774
-DEAL:0:3d::#constraints on 4: 2557
-DEAL:0:3d::#constraints on 5: 4330
-DEAL:0:3d::#constraints on 6: 2922
-DEAL:0:3d::#constraints on 7: 3888
-DEAL:0:3d::#constraints on 8: 1147
+DEAL:0:3d::#dofs = 27018
+DEAL:0:3d::#constraints on 0: 168
+DEAL:0:3d::#constraints on 1: 1017
+DEAL:0:3d::#constraints on 2: 617
+DEAL:0:3d::#constraints on 3: 1096
+DEAL:0:3d::#constraints on 4: 626
+DEAL:0:3d::#constraints on 5: 1032
+DEAL:0:3d::#constraints on 6: 695
+DEAL:0:3d::#constraints on 7: 907
+DEAL:0:3d::#constraints on 8: 308
 DEAL:0:3d::OK

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.