From 5be6b2088096fe6003204cf0252db58188ad5ef2 Mon Sep 17 00:00:00 2001 From: Guido Kanschat Date: Sun, 26 May 2013 22:04:50 +0000 Subject: [PATCH] add data exchange to KElly git-svn-id: https://svn.dealii.org/trunk@29623 0785d39b-7218-0410-832d-ea1e28bc413d --- deal.II/examples/step-50/step-50.cc | 11 +++++++++-- 1 file changed, 9 insertions(+), 2 deletions(-) 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: " -- 2.39.5