Skip to content
Snippets Groups Projects

Fix a bug in 'pre' regarding the handling of scalar-valued matrices

Merged Oliver Sander requested to merge bugfix-for-amg-with-scalar-matrix into master
1 file
+ 1
1
Compare changes
  • Side-by-side
  • Inline
+ 1
1
@@ -587,7 +587,7 @@ namespace Dune
{
Hybrid::ifElse(IsNumber<Block>(),
[&](auto id) {
x[row.index()] = id(diagonal) / b[row.index()];
x[row.index()] = b[row.index()] / id(diagonal);
},
[&](auto id) {
id(diagonal).solve(x[row.index()], b[row.index()]);
Loading