Skip to content
Snippets Groups Projects
Commit 6feff92b authored by Oliver Sander's avatar Oliver Sander Committed by sander@FU-BERLIN.DE
Browse files

bugfix for the vector valued case: '1' has to be the identity matrix, not the 'all-1' matrix

[[Imported from SVN: r3322]]
parent f64473fc
Branches
Tags
No related merge requests found
#include <dune/istl/btdmatrix.hh> #include <dune/istl/btdmatrix.hh>
#include <dune/istl/scaledidmatrix.hh>
template<class OperatorType, class DiscFuncType, class BitVectorType> template<class OperatorType, class DiscFuncType, class BitVectorType>
inline inline
...@@ -62,7 +62,7 @@ void LineGSStep<OperatorType, DiscFuncType, BitVectorType >::iterate() ...@@ -62,7 +62,7 @@ void LineGSStep<OperatorType, DiscFuncType, BitVectorType >::iterate()
tridiagonalMatrix[j][j-1] = 0; tridiagonalMatrix[j][j-1] = 0;
// diagonal // diagonal
tridiagonalMatrix[j][j] = 1; tridiagonalMatrix[j][j] = Dune::ScaledIdentityMatrix<double,BlockSize>(1);
// left off-diagonal: // left off-diagonal:
if (j<current_block_size-1) if (j<current_block_size-1)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment