blob: c617a7021fcc1bebb8f5ac1378d30cb5614c4224 [file] [log] [blame]
 ! Test that NINT gives right results even in corner cases ! ! PR 31202 ! http://gcc.gnu.org/ml/fortran/2005-04/msg00139.html ! ! { dg-do run } ! { dg-xfail-run-if "PR 33271, math library bug" { powerpc-ibm-aix powerpc-*-linux* powerpc64-*-linux* *-*-mingw* } { "-O0" } { "" } } ! Note that this doesn't fail on powerpc64le-*-linux*. real(kind=8) :: a integer(kind=8) :: i1, i2 real :: b integer :: j1, j2 a = nearest(0.5_8,-1.0_8) i2 = nint(nearest(0.5_8,-1.0_8)) i1 = nint(a) if (i1 /= 0 .or. i2 /= 0) call abort a = 0.5_8 i2 = nint(0.5_8) i1 = nint(a) if (i1 /= 1 .or. i2 /= 1) call abort a = nearest(0.5_8,1.0_8) i2 = nint(nearest(0.5_8,1.0_8)) i1 = nint(a) if (i1 /= 1 .or. i2 /= 1) call abort b = nearest(0.5,-1.0) j2 = nint(nearest(0.5,-1.0)) j1 = nint(b) if (j1 /= 0 .or. j2 /= 0) call abort b = 0.5 j2 = nint(0.5) j1 = nint(b) if (j1 /= 1 .or. j2 /= 1) call abort b = nearest(0.5,1.0) j2 = nint(nearest(0.5,1.0)) j1 = nint(b) if (j1 /= 1 .or. j2 /= 1) call abort a = 4503599627370497.0_8 i1 = nint(a,kind=8) i2 = nint(4503599627370497.0_8,kind=8) if (i1 /= i2 .or. i1 /= 4503599627370497_8) call abort a = -4503599627370497.0_8 i1 = nint(a,kind=8) i2 = nint(-4503599627370497.0_8,kind=8) if (i1 /= i2 .or. i1 /= -4503599627370497_8) call abort end