From 6275894b3f9b1f15fb88a656d7f99e4a40196fcc Mon Sep 17 00:00:00 2001
From: frohne <frohne@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Wed, 20 Feb 2013 22:28:57 +0000
Subject: [PATCH] fix the computation of cell_constitution for the parallel
 case

git-svn-id: https://svn.dealii.org/trunk@28496 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/examples/step-42/step-42.cc | 9 ++++++++-
 1 file changed, 8 insertions(+), 1 deletion(-)

diff --git a/deal.II/examples/step-42/step-42.cc b/deal.II/examples/step-42/step-42.cc
index f2506e2b5d..a24556b2c2 100644
--- a/deal.II/examples/step-42/step-42.cc
+++ b/deal.II/examples/step-42/step-42.cc
@@ -896,7 +896,7 @@ namespace Step42
     unsigned int elast_points = 0;
     unsigned int plast_points = 0;
     double       yield = 0;
-    unsigned int cell_number = 0;
+    double       cell_number = 0;
     cell_constitution = 0;
 
     for (; cell!=endc; ++cell)
@@ -918,6 +918,7 @@ namespace Step42
 
               plast_lin_hard->plast_linear_hardening (stress_strain_tensor, strain_tensor[q_point],
                                                       elast_points, plast_points, yield);
+
               cell_constitution (cell_number) += yield;
               for (unsigned int i=0; i<dofs_per_cell; ++i)
                 {
@@ -961,9 +962,15 @@ namespace Step42
               system_rhs_newton);
 
           cell_number += 1;
+        }
+      else
+        {
+          cell_constitution (cell_number) = 0;
+          cell_number += 1;
         };
 
     cell_constitution /= n_q_points;
+    cell_constitution.compress (VectorOperation::add);
     system_rhs_newton.compress (VectorOperation::add);
 
     unsigned int sum_elast_points = Utilities::MPI::sum(elast_points, mpi_communicator);
-- 
2.39.5