diff --git a/dune/matrix-vector/triangularsolve.hh b/dune/matrix-vector/triangularsolve.hh index e21426ba0d5d4698832297c0ac3868ca757a74d1..1237abb776efdb9e2a86e435bd68331356dc0fc1 100644 --- a/dune/matrix-vector/triangularsolve.hh +++ b/dune/matrix-vector/triangularsolve.hh @@ -4,13 +4,13 @@ namespace Dune { namespace MatrixVector { template <class Matrix, class Vector, class BitVector> - static void lowerTriangularSolve(Matrix const& A, Vector const& b, Vector& x, + static void lowerTriangularSolve(Matrix const& L, Vector const& b, Vector& x, BitVector const* ignore, bool transpose = false) { x = 0; Vector r = b; if (transpose) { - for (auto it = A.begin(); it != A.end(); ++it) { + for (auto it = L.begin(); it != L.end(); ++it) { const size_t i = it.index(); if (ignore != nullptr and (*ignore)[i].all()) continue; @@ -23,7 +23,7 @@ namespace MatrixVector { } } } else { - for (auto it = A.begin(); it != A.end(); ++it) { + for (auto it = L.begin(); it != L.end(); ++it) { size_t i = it.index(); if (ignore != nullptr and (*ignore)[i].all()) continue;