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 }
   7 ! { dg-xfail-run-if "PR 33271, math library bug" { powerpc-ibm-aix* powerpc-*-linux* powerpc64-*-linux* *-*-mingw* } { "-O0" } { "" } }
   8 ! Note that this doesn't fail on powerpc64le-*-linux*.
   9   real(kind=8) :: a
  10   integer(kind=8) :: i1, i2
  11   real :: b
  12   integer :: j1, j2
  13
  14   a = nearest(0.5_8,-1.0_8)
  15   i2 = nint(nearest(0.5_8,-1.0_8))
  16   i1 = nint(a)
  17   if (i1 /= 0 .or. i2 /= 0) call abort
  18
  19   a = 0.5_8
  20   i2 = nint(0.5_8)
  21   i1 = nint(a)
  22   if (i1 /= 1 .or. i2 /= 1) call abort
  23
  24   a = nearest(0.5_8,1.0_8)
  25   i2 = nint(nearest(0.5_8,1.0_8))
  26   i1 = nint(a)
  27   if (i1 /= 1 .or. i2 /= 1) call abort
  28
  29   b = nearest(0.5,-1.0)
  30   j2 = nint(nearest(0.5,-1.0))
  31   j1 = nint(b)
  32   if (j1 /= 0 .or. j2 /= 0) call abort
  33
  34   b = 0.5
  35   j2 = nint(0.5)
  36   j1 = nint(b)
  37   if (j1 /= 1 .or. j2 /= 1) call abort
  38
  39   b = nearest(0.5,1.0)
  40   j2 = nint(nearest(0.5,1.0))
  41   j1 = nint(b)
  42   if (j1 /= 1 .or. j2 /= 1) call abort
  43
  44   a = 4503599627370497.0_8
  45   i1 = nint(a,kind=8)
  46   i2 = nint(4503599627370497.0_8,kind=8)
  47   if (i1 /= i2 .or. i1 /= 4503599627370497_8) call abort
  48
  49   a = -4503599627370497.0_8
  50   i1 = nint(a,kind=8)
  51   i2 = nint(-4503599627370497.0_8,kind=8)
  52   if (i1 /= i2 .or. i1 /= -4503599627370497_8) call abort
  53   end