From c86f1cd45974b08f23d595fc69770a61ede48dea Mon Sep 17 00:00:00 2001 From: Max Kahnt <max.kahnt@fu-berlin.de> Date: Thu, 25 Feb 2016 11:57:47 +0100 Subject: [PATCH] Fix: variabe name was confusing. --- dune/solvers/common/staticmatrixtools.hh | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/dune/solvers/common/staticmatrixtools.hh b/dune/solvers/common/staticmatrixtools.hh index 576ea7cf..70bd60ca 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]; } } } -- GitLab