From: 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-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=76f85d9d9e834cd8b64ad84dc9a9e20fac97248d;p=dealii-svn.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