From 4ac9f5d70e60737f63e65d453cd842ca2c114eb5 Mon Sep 17 00:00:00 2001 From: Joerg Frohne 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