Skip to content
Snippets Groups Projects
Commit c86f1cd4 authored by Max Kahnt's avatar Max Kahnt
Browse files

Fix: variabe name was confusing.

parent 0f38f2d3
No related branches found
No related tags found
No related merge requests found
Pipeline #
......@@ -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];
}
}
}
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment