From b6ef5221f6e6f4f5dfb4f1d607112b884faba67b Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Wed, 27 Jul 2016 17:46:00 +0200
Subject: [PATCH] Fix initialisation
---
dune/solvers/iterationsteps/choleskystep.hh | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/dune/solvers/iterationsteps/choleskystep.hh b/dune/solvers/iterationsteps/choleskystep.hh
index 2a2157b..b556c5d 100644
--- a/dune/solvers/iterationsteps/choleskystep.hh
+++ b/dune/solvers/iterationsteps/choleskystep.hh
@@ -301,7 +301,7 @@ public:
1. Gy = b
2. G'x = y
*/
- Vector y(G.N());
+ Vector y = *this->x_;
BitVector const *ignore = this->hasIgnore() ? &this->ignore() : nullptr;
Dune::MatrixVector::lowerTriangularSolve(G, *this->rhs_, y, ignore);
Dune::MatrixVector::upperTriangularSolve(G, y, *this->x_, ignore, true);
--
GitLab