diff --git a/dune/matrix-vector/axpy.hh b/dune/matrix-vector/axpy.hh index 12b54d389c63496b267b29a9c1db52fb19c1a9b3..5bcdd028097f97755573ca8a2992fb5ea7017a32 100644 --- a/dune/matrix-vector/axpy.hh +++ b/dune/matrix-vector/axpy.hh @@ -158,12 +158,10 @@ namespace MatrixVector { typedef Dune::FieldMatrix<T, n, m> B; typedef Dune::FieldMatrix<T, m, k> C; static void addProduct(A& a, const B& b, const C& c) { - for (size_t row = 0; row < n; ++row) { - for (size_t col = 0; col < k; ++col) { + for (size_t row = 0; row < n; ++row) + for (size_t col = 0; col < k; ++col) for (size_t i = 0; i < m; ++i) a[row][col] += b[row][i] * c[i][col]; - } - } } }; @@ -175,12 +173,10 @@ namespace MatrixVector { typedef Dune::FieldMatrix<T, n, m> C; typedef Dune::FieldMatrix<T, m, k> D; static void addProduct(A& a, const T& b, const C& c, const D& d) { - for (size_t row = 0; row < n; ++row) { - for (size_t col = 0; col < k; ++col) { + for (size_t row = 0; row < n; ++row) + for (size_t col = 0; col < k; ++col) for (size_t i = 0; i < m; ++i) a[row][col] += b * c[row][i] * d[i][col]; - } - } } };