From: Wolfgang Bangerth Date: Wed, 21 Nov 2012 20:04:42 +0000 (+0000) Subject: We can't use PETSc with threads, so can't use the Kelly error estimator here. Use... X-Git-Tag: v8.0.0~1777 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=d0a0733f4d97f83098e91145140552f768d8415d;p=dealii.git We can't use PETSc with threads, so can't use the Kelly error estimator here. Use globally refined meshes instead. git-svn-id: https://svn.dealii.org/trunk@27664 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/tests/mpi/step-40.cc b/tests/mpi/step-40.cc index caef596fea..5a2a2b6715 100644 --- a/tests/mpi/step-40.cc +++ b/tests/mpi/step-40.cc @@ -270,17 +270,7 @@ namespace Step40 template void LaplaceProblem::refine_grid () { - Vector estimated_error_per_cell (triangulation.n_active_cells()); - KellyErrorEstimator::estimate (dof_handler, - QGauss(3), - typename FunctionMap::type(), - locally_relevant_solution, - estimated_error_per_cell); - parallel::distributed::GridRefinement:: - refine_and_coarsen_fixed_number (triangulation, - estimated_error_per_cell, - 0.3, 0.03); - triangulation.execute_coarsening_and_refinement (); + triangulation.refine_global (1); } @@ -289,7 +279,7 @@ namespace Step40 template void LaplaceProblem::run () { - const unsigned int n_cycles = 3; + const unsigned int n_cycles = 2; for (unsigned int cycle=0; cycle