From aa0e7c9f674dfcf6748b4a687467de955b0b9132 Mon Sep 17 00:00:00 2001
From: Elias Pipping <elias.pipping@fu-berlin.de>
Date: Mon, 17 Nov 2014 13:31:11 +0100
Subject: [PATCH] [Output] Special handling for lambda = 0
---
src/fixedpointiterator.cc | 3 ++-
1 file changed, 2 insertions(+), 1 deletion(-)
diff --git a/src/fixedpointiterator.cc b/src/fixedpointiterator.cc
index 8a73ba48..4bc8dee7 100644
--- a/src/fixedpointiterator.cc
+++ b/src/fixedpointiterator.cc
@@ -72,7 +72,8 @@ FixedPointIterator<Factory, StateUpdater, VelocityUpdater, ErrorNorm>::run(
ScalarVector newAlpha;
stateUpdater->extractAlpha(newAlpha);
- if (errorNorm_.diff(alpha, newAlpha) < fixedPointTolerance_) {
+ if (lambda_ < 1e-12 or
+ errorNorm_.diff(alpha, newAlpha) < fixedPointTolerance_) {
fixedPointIteration++;
break;
}
--
GitLab