gcc/testsuite/gfortran.dg/nint_2.f90

   1 ! Test that NINT gives right results even in corner cases
   2 !
   3 ! PR 31202
   4 ! http://gcc.gnu.org/ml/fortran/2005-04/msg00139.html
   5 !
   6 ! { dg-do run { xfail powerpc-ibm-aix* powerpc*-*-linux* } }
   7   real(kind=8) :: a
   8   integer(kind=8) :: i1, i2
   9   real :: b
  10   integer :: j1, j2
  11
  12   a = nearest(0.5_8,-1.0_8)
  13   i2 = nint(nearest(0.5_8,-1.0_8))
  14   i1 = nint(a)
  15   if (i1 /= 0 .or. i2 /= 0) call abort
  16
  17   a = 0.5_8
  18   i2 = nint(0.5_8)
  19   i1 = nint(a)
  20   if (i1 /= 1 .or. i2 /= 1) call abort
  21
  22   a = nearest(0.5_8,1.0_8)
  23   i2 = nint(nearest(0.5_8,1.0_8))
  24   i1 = nint(a)
  25   if (i1 /= 1 .or. i2 /= 1) call abort
  26
  27   b = nearest(0.5,-1.0)
  28   j2 = nint(nearest(0.5,-1.0))
  29   j1 = nint(b)
  30   if (j1 /= 0 .or. j2 /= 0) call abort
  31
  32   b = 0.5
  33   j2 = nint(0.5)
  34   j1 = nint(b)
  35   if (j1 /= 1 .or. j2 /= 1) call abort
  36
  37   b = nearest(0.5,1.0)
  38   j2 = nint(nearest(0.5,1.0))
  39   j1 = nint(b)
  40   if (j1 /= 1 .or. j2 /= 1) call abort
  41
  42   a = 4503599627370497.0_8
  43   i1 = nint(a,kind=8)
  44   i2 = nint(4503599627370497.0_8,kind=8)
  45   if (i1 /= i2 .or. i1 /= 4503599627370497_8) call abort
  46
  47   a = -4503599627370497.0_8
  48   i1 = nint(a,kind=8)
  49   i2 = nint(-4503599627370497.0_8,kind=8)
  50   if (i1 /= i2 .or. i1 /= -4503599627370497_8) call abort
  51   end