diff --git a/dumux/boxmodels/1p/1plocalresidual.hh b/dumux/boxmodels/1p/1plocalresidual.hh
index 0454731a4e7fa480a3d49f2b5bfa88e835f0068e..eeb39d9d831fb97be764dc6eb6114d883049e9b7 100644
--- a/dumux/boxmodels/1p/1plocalresidual.hh
+++ b/dumux/boxmodels/1p/1plocalresidual.hh
@@ -118,10 +118,7 @@ public:
         fluxVars.intrinsicPermeability().mv(fluxVars.potentialGrad(),
                                             tmpVec);
 
-        Scalar normalFlux = 0;
-        for (int i = 0; i < Vector::size; ++i)
-            normalFlux += tmpVec[i]*fluxVars.face().normal[i];
-        normalFlux *= -1;
+        Scalar normalFlux = -(tmpVec*fluxVars.face().normal);
 
         const VolumeVariables &up = this->curVolVars_(fluxVars.upstreamIdx(normalFlux));
         const VolumeVariables &dn = this->curVolVars_(fluxVars.downstreamIdx(normalFlux));
diff --git a/dumux/boxmodels/1p2c/1p2cfluxvariables.hh b/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
index 79096bfcec46a160fa173db68356fe91a7a69fa2..c1cf464a23b20b1970190979b86d20d73a85c4bb 100644
--- a/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
+++ b/dumux/boxmodels/1p2c/1p2cfluxvariables.hh
@@ -412,10 +412,7 @@ protected:
                               const ElementVolumeVariables &elemDat)
     {
         K_.mv(potentialGrad_, Kmvp_);
-        KmvpNormal_ = 0;
-        for (int i = 0; i < Vector::size; ++i)
-            KmvpNormal_ += Kmvp_[i] * face().normal[i];
-        KmvpNormal_ *= -1;
+        KmvpNormal_ = -(Kmvp_*face().normal);
 
         // set the upstream and downstream vertices
         upstreamIdx_ = face().i;
diff --git a/dumux/boxmodels/2p/2plocalresidual.hh b/dumux/boxmodels/2p/2plocalresidual.hh
index 1e840efdb0fa0a8aadee12b2b2cce0466a6d9ee8..1c9a2a1f0e366e4bb43bbde84f45b271e99c4627 100644
--- a/dumux/boxmodels/2p/2plocalresidual.hh
+++ b/dumux/boxmodels/2p/2plocalresidual.hh
@@ -153,10 +153,7 @@ public:
             // (the minus comes from the Darcy law which states that
             // the flux is from high to low pressure potentials.)
             fluxVars.intrinsicPermeability().mv(fluxVars.potentialGrad(phaseIdx), tmpVec);
-            Scalar normalFlux = 0;
-            for (int i = 0; i < Vector::size; ++i)
-                normalFlux += tmpVec[i]*fluxVars.face().normal[i];
-            normalFlux *= -1;
+            Scalar normalFlux = -(tmpVec*fluxVars.face().normal);
 
             // data attached to upstream and the downstream vertices
             // of the current phase
diff --git a/dumux/boxmodels/2p2c/2p2cfluxvariables.hh b/dumux/boxmodels/2p2c/2p2cfluxvariables.hh
index 706197bc911d0cd1495c5fe9ae8e0023626bf1c7..38dfe098076c3e716d8bed8039df0c1b91f07a85 100644
--- a/dumux/boxmodels/2p2c/2p2cfluxvariables.hh
+++ b/dumux/boxmodels/2p2c/2p2cfluxvariables.hh
@@ -231,10 +231,7 @@ protected:
                                                                     fvElemGeom_,
                                                                     face().j));
             K.mv(potentialGrad_[phaseIdx], Kmvp_[phaseIdx]);
-            KmvpNormal_[phaseIdx] = 0;
-            for (int i = 0; i < Vector::size; ++i)
-                KmvpNormal_[phaseIdx] += Kmvp_[phaseIdx][i] * face().normal[i];
-            KmvpNormal_[phaseIdx] *= -1;
+            KmvpNormal_[phaseIdx] = -(Kmvp_[phaseIdx]*face().normal);
         }
 
         // set the upstream and downstream vertices
diff --git a/dumux/boxmodels/richards/richardslocalresidual.hh b/dumux/boxmodels/richards/richardslocalresidual.hh
index 4fdc70f55114ae5bdbf7b417ba03b1bdb146f56e..7d5e74e2cbd733d40118b030be24a512183f8cc4 100644
--- a/dumux/boxmodels/richards/richardslocalresidual.hh
+++ b/dumux/boxmodels/richards/richardslocalresidual.hh
@@ -129,10 +129,7 @@ public:
         Vector tmpVec;
         fluxVars.intrinsicPermeability().mv(fluxVars.potentialGradW(),
                                             tmpVec);
-        Scalar normalFlux = 0;
-        for (int i = 0; i < Vector::size; ++ i)
-             normalFlux += tmpVec[i] + fluxVars.face().normal[i];
-        normalFlux *= -1;
+        Scalar normalFlux = -(tmpVec*fluxVars.face().normal);
 
         // data attached to upstream and the downstream vertices
         // of the current phase