la_xisnan.F90

Go to the documentation of this file.

module la_xisnan
   interface la_isnan
 
   module procedure sisnan
   module procedure disnan
 
   end interface
 
contains
   
   logical function sisnan( x )
   use la_constants, only: wp=>sp
#ifdef USE_IEEE_INTRINSIC
   use, intrinsic :: ieee_arithmetic
#elif USE_ISNAN
   intrinsic :: isnan
#endif
   real(wp) :: x
#ifdef USE_IEEE_INTRINSIC
   sisnan = ieee_is_nan(x)
#elif USE_ISNAN
   sisnan = isnan(x)
#else
   sisnan = slaisnan(x,x)
 
   contains
   logical function slaisnan( x, y )
   use la_constants, only: wp=>sp
   real(wp) :: x, y
   slaisnan = ( x.ne.y )
   end function slaisnan
#endif
   end function sisnan
 
   logical function disnan( x )
   use la_constants, only: wp=>dp
#ifdef USE_IEEE_INTRINSIC
   use, intrinsic :: ieee_arithmetic
#elif USE_ISNAN
   intrinsic :: isnan
#endif
   real(wp) :: x
#ifdef USE_IEEE_INTRINSIC
   disnan = ieee_is_nan(x)
#elif USE_ISNAN
   disnan = isnan(x)
#else
   disnan = dlaisnan(x,x)
 
   contains
   logical function dlaisnan( x, y )
   use la_constants, only: wp=>dp
   real(wp) :: x, y
   dlaisnan = ( x.ne.y )
   end function dlaisnan
#endif
   end function disnan
 
end module la_xisnan