gcc/testsuite/gfortran.dg/fold_nearest.f90

   1 ! { dg-do run }
   2 ! Tests for the constant folding of the NEAREST intrinsic
   3 ! We compare against the results of the runtime implementation,
   4 ! thereby making sure that they remain consistent
   5 REAL, PARAMETER :: x(10) = (/ 1., 0.49999997, 0.5, 8388609.0, -1., &
   6                                       -0.49999997, -0.5, -8388609.0, &
   7                                       0., 0. /), &
   8                  dir(10) = (/ -1.,       +1., -1.,       -1., +1., &
   9                                              -1.,  +1.,        +1., &
  10                                      +1.,-1./)
  11 REAL :: a(10)
  12
  13 a = x
  14 if (nearest (x(1), dir(1)) /= nearest (a(1), dir(1))) STOP 1
  15 if (nearest (x(2), dir(2)) /= nearest (a(2), dir(2))) STOP 2
  16 if (nearest (x(3), dir(3)) /= nearest (a(3), dir(3))) STOP 3
  17 if (nearest (x(4), dir(4)) /= nearest (a(4), dir(4))) STOP 4
  18 if (nearest (x(5), dir(5)) /= nearest (a(5), dir(5))) STOP 5
  19 if (nearest (x(6), dir(6)) /= nearest (a(6), dir(6))) STOP 6
  20 if (nearest (x(7), dir(7)) /= nearest (a(7), dir(7))) STOP 7
  21 if (nearest (x(8), dir(8)) /= nearest (a(8), dir(8))) STOP 8
  22 ! These last two tests are commented out because mpfr provides no support
  23 ! for denormals, and therefore we get TINY instead of the correct result.
  24 !if (nearest (x(9), dir(9)) /= nearest (a(9), dir(9))) STOP 9
  25 !if (nearest (x(10), dir(10)) /= nearest (a(10), dir(10))) STOP 10
  26
  27 end