From 810be78778f47d595e931f4fae0caf01bca049e8 Mon Sep 17 00:00:00 2001
From: Oliver Sander <oliver.sander@tu-dresden.de>
Date: Wed, 23 Jan 2019 22:22:52 +0100
Subject: [PATCH] Fix a bug in 'pre' regarding the handling of scalar-valued
matrices
This bug is really new -- it only came with the introduction of
support for scalar-valued matrices.
---
dune/istl/paamg/amg.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/istl/paamg/amg.hh b/dune/istl/paamg/amg.hh
index 7551b8dda..249c0e093 100644
--- a/dune/istl/paamg/amg.hh
+++ b/dune/istl/paamg/amg.hh
@@ -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()]);
--
GitLab