blob: 6d49a6fa6866582f4807a5c308e40f6091e5876f [file] [log] [blame]
! Test if MIN and MAX intrinsics behave correctly when passed NaNs
! as arguments
!
! { dg-do run }
! { dg-add-options ieee }
!
module aux2
interface isnan
module procedure isnan_r
module procedure isnan_d
end interface isnan
interface isinf
module procedure isinf_r
module procedure isinf_d
end interface isinf
contains
pure function isnan_r(x) result (isnan)
logical :: isnan
real, intent(in) :: x
isnan = (.not.(x == x))
end function isnan_r
pure function isnan_d(x) result (isnan)
logical :: isnan
double precision, intent(in) :: x
isnan = (.not.(x == x))
end function isnan_d
pure function isinf_r(x) result (isinf)
logical :: isinf
real, intent(in) :: x
isinf = (x > huge(x)) .or. (x < -huge(x))
end function isinf_r
pure function isinf_d(x) result (isinf)
logical :: isinf
double precision, intent(in) :: x
isinf = (x > huge(x)) .or. (x < -huge(x))
end function isinf_d
end module aux2
program test
use aux2
implicit none
real :: nan, large, inf
! Create a NaN and check it
nan = 0
nan = nan / nan
if (nan == nan .or. nan > nan .or. nan < nan .or. nan >= nan &
.or. nan <= nan) STOP 1
if (isnan (2.d0) .or. (.not. isnan(nan)) .or. &
(.not. isnan(real(nan,kind=kind(2.d0))))) STOP 2
! Create an INF and check it
large = huge(large)
inf = 2 * large
if (isinf(nan) .or. isinf(large) .or. .not. isinf(inf)) STOP 3
if (isinf(-nan) .or. isinf(-large) .or. .not. isinf(-inf)) STOP 4
! Check that MIN and MAX behave correctly
if (.not. isnan(min(nan,nan))) STOP 13
if (.not. isnan(max(nan,nan))) STOP 14
! Same thing, with more arguments
if (.not. isnan(min(nan,nan,nan))) STOP 27
if (.not. isnan(max(nan,nan,nan))) STOP 28
if (.not. isnan(min(nan,nan,nan,nan))) STOP 29
if (.not. isnan(max(nan,nan,nan,nan))) STOP 30
if (.not. isnan(min(nan,nan,nan,nan,nan))) STOP 31
if (.not. isnan(max(nan,nan,nan,nan,nan))) STOP 32
! Large values, INF and NaNs
if (.not. isinf(max(large, inf))) STOP 33
if (isinf(min(large, inf))) STOP 34
if (.not. isinf(min(-large, -inf))) STOP 41
if (isinf(max(-large, -inf))) STOP 42
end program test