flang/test/Evaluate/folding28.f90

   1 ! RUN: %python %S/test_folding.py %s %flang_fc1
   2 ! Tests folding of SQRT()
   3 module m
   4   implicit none
   5   ! +Inf
   6   real(8), parameter :: inf8 = z'7ff0000000000000'
   7   logical, parameter :: test_inf8 = sqrt(inf8) == inf8
   8   ! max finite
   9   real(8), parameter :: h8 = huge(1.0_8), h8z = z'7fefffffffffffff'
  10   logical, parameter :: test_h8 = h8 == h8z
  11   real(8), parameter :: sqrt_h8 = sqrt(h8), sqrt_h8z = z'5fefffffffffffff'
  12   logical, parameter :: test_sqrt_h8 = sqrt_h8 == sqrt_h8z
  13   real(8), parameter :: sqr_sqrt_h8 = sqrt_h8 * sqrt_h8, sqr_sqrt_h8z = z'7feffffffffffffe'
  14   logical, parameter :: test_sqr_sqrt_h8 = sqr_sqrt_h8 == sqr_sqrt_h8z
  15   ! -0 (sqrt is -0)
  16   real(8), parameter :: n08 = z'8000000000000000'
  17   real(8), parameter :: sqrt_n08 = sqrt(n08)
  18   !WARN: warning: division by zero
  19   real(8), parameter :: inf_n08 = 1.0_8 / sqrt_n08, inf_n08z = z'fff0000000000000'
  20   logical, parameter :: test_n08 = inf_n08 == inf_n08z
  21   ! min normal
  22   real(8), parameter :: t8 = tiny(1.0_8), t8z = z'0010000000000000'
  23   logical, parameter :: test_t8 = t8 == t8z
  24   real(8), parameter :: sqrt_t8 = sqrt(t8), sqrt_t8z = z'2000000000000000'
  25   logical, parameter :: test_sqrt_t8 = sqrt_t8 == sqrt_t8z
  26   real(8), parameter :: sqr_sqrt_t8 = sqrt_t8 * sqrt_t8
  27   logical, parameter :: test_sqr_sqrt_t8 = sqr_sqrt_t8 == t8
  28   ! max subnormal
  29   real(8), parameter :: maxs8 = z'000fffffffffffff'
  30   real(8), parameter :: sqrt_maxs8 = sqrt(maxs8), sqrt_maxs8z = z'1fffffffffffffff'
  31   logical, parameter :: test_sqrt_maxs8 = sqrt_maxs8 == sqrt_maxs8z
  32   ! min subnormal
  33   real(8), parameter :: mins8 = z'1'
  34   real(8), parameter :: sqrt_mins8 = sqrt(mins8), sqrt_mins8z = z'1e60000000000000'
  35   logical, parameter :: test_sqrt_mins8 = sqrt_mins8 == sqrt_mins8z
  36   real(8), parameter :: sqr_sqrt_mins8 = sqrt_mins8 * sqrt_mins8
  37   logical, parameter :: test_sqr_sqrt_mins8 = sqr_sqrt_mins8 == mins8
  38   ! regression tests: cases near 1.
  39   real(4), parameter :: sqrt_under1 = sqrt(.96875)
  40   logical, parameter :: test_sqrt_under1 = sqrt_under1 == .984250962734222412109375
  41   ! oddball case: the value before 1. is also its own sqrt, but not its own square
  42   real(4), parameter :: before_1 = z'3f7fffff' ! .999999940395355224609375
  43   real(4), parameter :: sqrt_before_1 = sqrt(before_1)
  44   logical, parameter :: test_before_1 = sqrt_before_1 == before_1
  45   real(4), parameter :: sq_sqrt_before_1 = sqrt_before_1 * sqrt_before_1
  46   logical, parameter :: test_sq_before_1 = sq_sqrt_before_1 < before_1
  47 end module