From 5be6b2088096fe6003204cf0252db58188ad5ef2 Mon Sep 17 00:00:00 2001
From: Guido Kanschat <dr.guido.kanschat@gmail.com>
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<float> 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<dim>::estimate (static_cast<DoFHandler<dim>&>(mg_dof_handler),
                                         QGauss<dim-1>(3),
                                         typename FunctionMap<dim>::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