From: Guido Kanschat Date: Sun, 26 May 2013 22:04:50 +0000 (+0000) Subject: add data exchange to KElly X-Git-Tag: v8.0.0~414 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=5be6b2088096fe6003204cf0252db58188ad5ef2;p=dealii.git add data exchange to KElly git-svn-id: https://svn.dealii.org/trunk@29623 0785d39b-7218-0410-832d-ea1e28bc413d --- diff --git a/deal.II/examples/step-50/step-50.cc b/deal.II/examples/step-50/step-50.cc index bd9ff81261..7bdde8d208 100644 --- a/deal.II/examples/step-50/step-50.cc +++ b/deal.II/examples/step-50/step-50.cc @@ -950,10 +950,17 @@ namespace Step50 { Vector estimated_error_per_cell (triangulation.n_active_cells()); + TrilinosWrappers::MPI::Vector temp_solution; + IndexSet idx; + DoFTools::extract_locally_relevant_dofs (mg_dof_handler, + idx); + temp_solution.reinit(idx, MPI_COMM_WORLD); + temp_solution = solution; + KellyErrorEstimator::estimate (static_cast&>(mg_dof_handler), QGauss(3), typename FunctionMap::type(), - solution, + temp_solution, estimated_error_per_cell); GridRefinement::refine_and_coarsen_fixed_number (triangulation, estimated_error_per_cell, @@ -1017,7 +1024,7 @@ namespace Step50 } else triangulation.refine_global (1); - //refine_grid (); + //refine_grid (); deallog << " Number of active cells: "