[ci optional] change safe_logs to test for FINITE IEEE values rather than NORMAL

Author: Matthias Fabry

math/public/math_lib_crmath.f90

@@ -153,7 +153,7 @@ elemental function safe_log_ (x) result (log_x)
     real(dp), intent(in) :: x
     real(dp)             :: log_x
 
-    if (.NOT. IEEE_IS_NORMAL(x)) then
+    if (.NOT. IEEE_IS_FINITE(x)) then
 
        log_x = 1.E-99_dp
 
@@ -172,7 +172,7 @@ elemental function safe_log10_ (x) result (log10_x)
     real(dp), intent(in) :: x
     real(dp)             :: log10_x
 
-    if (.NOT. IEEE_IS_NORMAL(x)) then
+    if (.NOT. IEEE_IS_FINITE(x)) then
 
        log10_x = 1.E-99_dp
 

math/public/math_lib_intrinsic.f90

@@ -171,7 +171,7 @@ elemental function safe_log_ (x) result (log_x)
     real(dp), intent(in) :: x
     real(dp)             :: log_x
 
-    if (.NOT. IEEE_IS_NORMAL(x)) then
+    if (.NOT. IEEE_IS_FINITE(x)) then
 
        log_x = 1.E-99_dp
 
@@ -190,7 +190,7 @@ elemental function safe_log10_ (x) result (log10_x)
     real(dp), intent(in) :: x
     real(dp)             :: log10_x
 
-    if (.NOT. IEEE_IS_NORMAL(x)) then
+    if (.NOT. IEEE_IS_FINITE(x)) then
 
        log10_x = 1.E-99_dp