diff --git a/dune/solvers/common/staticmatrixtools.hh b/dune/solvers/common/staticmatrixtools.hh index 576ea7cfc3b394c45828336c47b23730b5873e36..70bd60caa9eee6b690ebbe1baa25545cc364672f 100644 --- a/dune/solvers/common/staticmatrixtools.hh +++ b/dune/solvers/common/staticmatrixtools.hh @@ -150,15 +150,15 @@ class StaticMatrix const typename Dune::FieldMatrix<K3, n, n>& B, const typename Dune::FieldMatrix<K2, n, m>& T2) { - typename Dune::FieldMatrix<K1, m, n> T2transposedB; - T2transposedB = 0; + typename Dune::FieldMatrix<K1, m, n> T1transposedB; + T1transposedB = 0; for (size_t i=0; i<A.N(); ++i) { for (size_t k=0; k<B.N(); ++k) { if (T1[k][i]!=0) for (size_t l=0; l<B.M(); ++l) - T2transposedB[i][l] += T1[k][i] * B[k][l]; + T1transposedB[i][l] += T1[k][i] * B[k][l]; } } for (size_t k=0; k<T2.N(); ++k) @@ -167,7 +167,7 @@ class StaticMatrix { if (T2[k][l]!=0) for (size_t i=0; i<A.N(); ++i) - A[i][l] += T2transposedB[i][k] * T2[k][l]; + A[i][l] += T1transposedB[i][k] * T2[k][l]; } } }