Merge pull request #511 from MESAHub/change_safe_logs

change safe_logs test from IEEE_NORMAL to IEEE_FINITE

Author: evbauer

math/public/math_lib_crmath.f90

@@ -153,9 +153,9 @@ 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
+       log_x = -99._dp
 
     else
 
@@ -172,9 +172,9 @@ 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
+       log10_x = -99._dp
 
     else
 

math/public/math_lib_intrinsic.f90

@@ -171,9 +171,9 @@ 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
+       log_x = -99._dp
 
     else
 
@@ -190,9 +190,9 @@ 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
+       log10_x = -99._dp
 
     else