From: Wolfgang Bangerth <bangerth@math.tamu.edu>
Date: Thu, 28 May 2015 20:37:37 +0000 (-0500)
Subject: Use a standard method to find whether a number is normal.
X-Git-Tag: v8.3.0-rc1~142^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=ba81b62200bf17280130713a2064c15033afe123;p=dealii.git

Use a standard method to find whether a number is normal.

In the current context, this slightly changes the behavior in that the new
condition also catches infinities, rather than just NaNs as before, but the
downstream code will still yield the same answers.
---

diff --git a/include/deal.II/base/logstream.h b/include/deal.II/base/logstream.h
index 6d6720ffad..577f9dc691 100644
--- a/include/deal.II/base/logstream.h
+++ b/include/deal.II/base/logstream.h
@@ -1,6 +1,6 @@
 // ---------------------------------------------------------------------
 //
-// Copyright (C) 1998 - 2014 by the deal.II authors
+// Copyright (C) 1998 - 2015 by the deal.II authors
 //
 // This file is part of the deal.II library.
 //
@@ -594,11 +594,10 @@ LogStream::operator<< (const double t)
 {
   std::ostringstream &stream = get_stream();
 
+  // drop small numbers or skew them away from zero.
   // we have to make sure that we don't catch NaN's and +-Inf's with the
   // test, because for these denormals all comparisons are always false.
-  // thus, for a NaN, both t<=0 and t>=0 are false at the same time, which
-  // can't be said for any other number
-  if (! (t<=0) && !(t>=0))
+  if (! numbers::is_finite(t))
     stream << t;
   else if (std::fabs(t) < double_threshold)
     stream << '0';