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 54