From de237fbd7c75aae389a5e83ff9475f31952c2c4a Mon Sep 17 00:00:00 2001
From: maier <maier@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Mon, 8 Apr 2013 18:49:48 +0000
Subject: [PATCH] Well, do it correctly.

git-svn-id: https://svn.dealii.org/trunk@29221 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/include/deal.II/numerics/vector_tools.templates.h | 8 ++------
 1 file changed, 2 insertions(+), 6 deletions(-)

diff --git a/deal.II/include/deal.II/numerics/vector_tools.templates.h b/deal.II/include/deal.II/numerics/vector_tools.templates.h
index 33bde2ad2f..01b12b9a27 100644
--- a/deal.II/include/deal.II/numerics/vector_tools.templates.h
+++ b/deal.II/include/deal.II/numerics/vector_tools.templates.h
@@ -5645,18 +5645,14 @@ namespace VectorTools
     Assert(p_select.size() == 0 || p_select.size() == n,
            ExcDimensionMismatch(p_select.size(), n));
 
-    typename VECTOR::value_type s = 0.;
-
     if(p_select.size() == 0)
       {
         // In case of an empty boolean mask operate on the whole vector:
-        for (unsigned int i=0; i<n; ++i)
-          s += v(i);
-
-        v.add(-s/n);
+        v.add( - v.mean_value() );
       }
     else
       {
+        typename VECTOR::value_type s = 0.;
         unsigned int counter = 0;
         for (unsigned int i=0; i<n; ++i)
           if (p_select[i])
-- 
2.39.5