From 055d2006cde87fbf4fdde7840b27606f6d12459e Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Fri, 9 Feb 2018 12:35:55 +0100
Subject: [PATCH] [NL] Explicit quiet NaN in DeltaX crit. output.
---
NumLib/ODESolver/ConvergenceCriterionDeltaX.cpp | 4 +++-
1 file changed, 3 insertions(+), 1 deletion(-)
diff --git a/NumLib/ODESolver/ConvergenceCriterionDeltaX.cpp b/NumLib/ODESolver/ConvergenceCriterionDeltaX.cpp
index 68008f193bb..4cff22004c2 100644
--- a/NumLib/ODESolver/ConvergenceCriterionDeltaX.cpp
+++ b/NumLib/ODESolver/ConvergenceCriterionDeltaX.cpp
@@ -36,7 +36,9 @@ void ConvergenceCriterionDeltaX::checkDeltaX(const GlobalVector& minus_delta_x,
auto norm_x = MathLib::LinAlg::norm(x, _norm_type);
INFO("Convergence criterion: |dx|=%.4e, |x|=%.4e, |dx|/|x|=%.4e", error_dx,
- norm_x, error_dx / norm_x);
+ norm_x,
+ (norm_x == 0. ? std::numeric_limits<double>::quiet_NaN()
+ : (error_dx / norm_x)));
bool satisfied_abs = false;
bool satisfied_rel = false;
--
GitLab