From: kronbichler Date: Wed, 24 Aug 2011 12:14:21 +0000 (+0000) Subject: Modify shell radius to make test fail. X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=cf920dad87b520683ff02061cae8c530bf93e996;p=dealii-svn.git Modify shell radius to make test fail. git-svn-id: https://svn.dealii.org/trunk@24175 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/mpi/no_flux_constraints.cc b/tests/mpi/no_flux_constraints.cc index 671c2bf7b3..f66155fb7f 100644 --- a/tests/mpi/no_flux_constraints.cc +++ b/tests/mpi/no_flux_constraints.cc @@ -46,11 +46,14 @@ void test() parallel::distributed::Triangulation 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(), - 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 fe(FE_Q (2), dim); + static const FESystem fe(FE_Q (1), dim); DoFHandler 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 range = - vector.local_range(); - for (unsigned int i=range.first; i range = vector.local_range(); + for (unsigned int i=range.first; i