From 6feff92bb6cd07d7a1b74ddcfc0152720f00a0ed Mon Sep 17 00:00:00 2001
From: Oliver Sander <sander@igpm.rwth-aachen.de>
Date: Fri, 7 May 2010 15:49:37 +0000
Subject: [PATCH] bugfix for the vector valued case: '1' has to be the identity
matrix, not the 'all-1' matrix
[[Imported from SVN: r3322]]
---
dune/solvers/iterationsteps/linegsstep.cc | 4 ++--
1 file changed, 2 insertions(+), 2 deletions(-)
diff --git a/dune/solvers/iterationsteps/linegsstep.cc b/dune/solvers/iterationsteps/linegsstep.cc
index 47e1e4dd..a28d981d 100755
--- a/dune/solvers/iterationsteps/linegsstep.cc
+++ b/dune/solvers/iterationsteps/linegsstep.cc
@@ -1,5 +1,5 @@
#include <dune/istl/btdmatrix.hh>
-
+#include <dune/istl/scaledidmatrix.hh>
template<class OperatorType, class DiscFuncType, class BitVectorType>
inline
@@ -62,7 +62,7 @@ void LineGSStep<OperatorType, DiscFuncType, BitVectorType >::iterate()
tridiagonalMatrix[j][j-1] = 0;
// diagonal
- tridiagonalMatrix[j][j] = 1;
+ tridiagonalMatrix[j][j] = Dune::ScaledIdentityMatrix<double,BlockSize>(1);
// left off-diagonal:
if (j<current_block_size-1)
--
GitLab