diff --git a/dune/gfe/symmetricmatrix.hh b/dune/gfe/symmetricmatrix.hh
index 3ce94c909a357274eedb30499c06c0beb51b2294..d6d32d7cddd4cd1ed791ea7eb0b37f0173694aec 100644
--- a/dune/gfe/symmetricmatrix.hh
+++ b/dune/gfe/symmetricmatrix.hh
@@ -72,7 +72,7 @@ public:
       T result = 0;
       for (size_t i=0; i<N; i++)
         for (size_t j=0; j<=i; j++)
-          result += (1+(i!=j)) * operator()(i,j) * v1[i] * v2[j];
+            result += (1-0.5*(i==j)) * operator()(i,j) * (v1[i] * v2[j] + v1[j] * v2[i]);
 
       return result;
     }