From: Matthias Maier <tamiko@43-1.org>
Date: Thu, 3 May 2018 18:40:15 +0000 (-0500)
Subject: Bugfix: Use the general vector interface
X-Git-Tag: v9.0.0-rc1~28^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F6431%2Fhead;p=dealii.git

Bugfix: Use the general vector interface

We should really try hard to avoid using implementation specific low
level access via operator[].
---

diff --git a/include/deal.II/lac/linear_operator.h b/include/deal.II/lac/linear_operator.h
index 874d795d5a..97a5a20946 100644
--- a/include/deal.II/lac/linear_operator.h
+++ b/include/deal.II/lac/linear_operator.h
@@ -866,15 +866,17 @@ mean_value_filter(const std::function<void(Range &, bool)> &reinit_vector)
   return_op.vmult = [](Range &v, const Range &u)
   {
     const auto mean = u.mean_value();
-    for (types::global_dof_index i = 0; i<v.size(); ++i)
-      v(i) = u(i) - mean;
+
+    v = u;
+    v.add(-mean);
   };
 
   return_op.vmult_add = [](Range &v, const Range &u)
   {
     const auto mean = u.mean_value();
-    for (types::global_dof_index i = 0; i<v.size(); ++i)
-      v(i) += u(i) - mean;
+
+    v += u;
+    v.add(-mean);
   };
 
   return_op.Tvmult = return_op.vmult_add;