From 2820a31d6c10cd63ccaaab2fffc811f12d1d94e8 Mon Sep 17 00:00:00 2001
From: Max Kahnt <max.kahnt@fu-berlin.de>
Date: Tue, 26 Jul 2016 13:43:36 +0200
Subject: [PATCH] Avoid changing the result in ignored nodes.
---
dune/matrix-vector/triangularsolve.hh | 2 --
1 file changed, 2 deletions(-)
diff --git a/dune/matrix-vector/triangularsolve.hh b/dune/matrix-vector/triangularsolve.hh
index 7bf57e4..a9ab4c8 100644
--- a/dune/matrix-vector/triangularsolve.hh
+++ b/dune/matrix-vector/triangularsolve.hh
@@ -12,7 +12,6 @@ namespace MatrixVector {
static void lowerTriangularSolve(Matrix const& L, Vector b, Vector& x,
BitVector const* ignore,
bool transpose = false) {
- x = 0;
if (transpose) {
for (auto it = L.begin(); it != L.end(); ++it) {
const size_t i = it.index();
@@ -55,7 +54,6 @@ namespace MatrixVector {
static void upperTriangularSolve(Matrix const& U, Vector b, Vector& x,
BitVector const* ignore,
bool transpose = false) {
- x = 0;
if (transpose) {
for (auto it = U.beforeEnd(); it != U.beforeBegin(); --it) {
const size_t i = it.index();
--
GitLab