revise infnan

equipez · Oct 6, 2023 · 1c0205f · 1c0205f
1 parent 17323f5
commit 1c0205f
Showing 1 changed file with 4 additions and 4 deletions.
diff --git a/infnan.f90 b/infnan.f90
@@ -46,21 +46,21 @@ elemental pure function is_nan_sp(x) result(y)
 implicit none
 real(SP), intent(in) :: x
 logical :: y
-!y = (.not. (x <= huge_value(x) .and. x >= -huge_value(x))) .and. (.not. abs(x) > huge_value(x))
+!y = ((.not. (x <= huge_value(x) .and. x >= -huge_value(x)))) .and. (.not. abs(x) > huge_value(x))
 !y = (.not. is_finite(x) .and. .not. (abs(x) > huge_value(x))) .or. y
 y = ((.not. is_finite(x)) .and. (.not. is_inf(x)))
-y = (.not. is_inf(x) .and. .not. (x <= huge_value(x) .and. x >= -huge_value(x))) .or. y
+y = ((.not. is_inf(x)) .and. (.not. (x <= huge_value(x) .and. x >= -huge_value(x)))) .or. y
 end function is_nan_sp
 
 elemental pure function is_nan_dp(x) result(y)
 use consts_mod, only : DP
 implicit none
 real(DP), intent(in) :: x
 logical :: y
-!y = (.not. (x <= huge_value(x) .and. x >= -huge_value(x))) .and. (.not. abs(x) > huge_value(x))
+!y = ((.not. (x <= huge_value(x) .and. x >= -huge_value(x)))) .and. (.not. abs(x) > huge_value(x))
 !y = (.not. is_finite(x) .and. .not. (abs(x) > huge_value(x))) .or. y
 y = ((.not. is_finite(x)) .and. (.not. is_inf(x)))
-y = (.not. is_inf(x) .and. .not. (x <= huge_value(x) .and. x >= -huge_value(x))) .or. y
+y = ((.not. is_inf(x)) .and. (.not. (x <= huge_value(x) .and. x >= -huge_value(x)))) .or. y
 end function is_nan_dp