!Program to test SET_EXPONENT intrinsic function.

program test_set_exponent
  call test_real4()
  call test_real8()
end

subroutine test_real4()
  real*4 x,y
  integer*4 i,n
  equivalence(x, i)

  n = -148
  x = 1024.0
  y = set_exponent (x, n)
  if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) STOP 1

  n = 8
  x = 1024.0
  y = set_exponent (x, n)
  if (exponent (y) .ne. n) STOP 2

  n = 128
  i = 8388607
  x = transfer (i, x) ! z'007fffff' Positive denormalized floating-point.
  y = set_exponent (x, n)
  if (exponent (y) .ne. n) STOP 3

  n = -148
  x = -1024.0
  y = set_exponent (x, n)
  if  ((y .ne. 0.0) .and. (exponent (y) .ne. n)) STOP 4

  n = 8
  x = -1024.0
  y = set_exponent (x, n)
  if (y .ne. -128.0) STOP 5
  if (exponent (y) .ne. n) STOP 6

  n = 128
  i = -2139095041
  x = transfer (i, x) ! z'807fffff' Negative denormalized floating-point.
  y = set_exponent (x, n)
  if (exponent (y) .ne. n) STOP 7

end

subroutine test_real8()
  implicit none
  real*8 x, y
  integer*8 i, n
  equivalence(x, i)

  n = -1073
  x = 1024.0_8
  y = set_exponent (x, n)
  if  ((y .ne. 0.0_8) .and. (exponent (y) .ne. n)) STOP 8

  n = 8
  x = 1024.0_8
  y = set_exponent (x, n)
  if (y .ne. 128.0) STOP 9
  if (exponent (y) .ne. n) STOP 10

  n = 1024
  i = 4503599627370495_8
  x = transfer (i, x) !z'000fffffffffffff' Positive denormalized floating-point.
  y = set_exponent (x, n)
  if (exponent (y) .ne. n) STOP 11

  n = -1073
  x = -1024.0
  y = set_exponent (x, n)
  if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) STOP 12

  n = 8
  x = -1024.0
  y = set_exponent (x, n)
  if (y .ne. -128.0) STOP 13
  if (exponent (y) .ne. n) STOP 14

  n = 1024
  i = -9218868437227405313_8
  x = transfer (i, x)!z'800fffffffffffff' Negative denormalized floating-point.
  y = set_exponent (x, n)
  if (exponent (y) .ne. n) STOP 15
end