!Program to test SET_EXPONENT intrinsic function.
|
!Program to test SET_EXPONENT intrinsic function.
|
|
|
program test_set_exponent
|
program test_set_exponent
|
call test_real4()
|
call test_real4()
|
call test_real8()
|
call test_real8()
|
end
|
end
|
|
|
subroutine test_real4()
|
subroutine test_real4()
|
real*4 x,y
|
real*4 x,y
|
integer*4 i,n
|
integer*4 i,n
|
equivalence(x, i)
|
equivalence(x, i)
|
|
|
n = -148
|
n = -148
|
x = 1024.0
|
x = 1024.0
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
|
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
|
|
|
n = 8
|
n = 8
|
x = 1024.0
|
x = 1024.0
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
n = 128
|
n = 128
|
i = 8388607
|
i = 8388607
|
x = transfer (i, x) ! z'007fffff' Positive denormalized floating-point.
|
x = transfer (i, x) ! z'007fffff' Positive denormalized floating-point.
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
n = -148
|
n = -148
|
x = -1024.0
|
x = -1024.0
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
|
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
|
|
|
n = 8
|
n = 8
|
x = -1024.0
|
x = -1024.0
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (y .ne. -128.0) call abort()
|
if (y .ne. -128.0) call abort()
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
n = 128
|
n = 128
|
i = -2139095041
|
i = -2139095041
|
x = transfer (i, x) ! z'807fffff' Negative denormalized floating-point.
|
x = transfer (i, x) ! z'807fffff' Negative denormalized floating-point.
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
end
|
end
|
|
|
subroutine test_real8()
|
subroutine test_real8()
|
implicit none
|
implicit none
|
real*8 x, y
|
real*8 x, y
|
integer*8 i, n
|
integer*8 i, n
|
equivalence(x, i)
|
equivalence(x, i)
|
|
|
n = -1073
|
n = -1073
|
x = 1024.0_8
|
x = 1024.0_8
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if ((y .ne. 0.0_8) .and. (exponent (y) .ne. n)) call abort()
|
if ((y .ne. 0.0_8) .and. (exponent (y) .ne. n)) call abort()
|
|
|
n = 8
|
n = 8
|
x = 1024.0_8
|
x = 1024.0_8
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (y .ne. 128.0) call abort()
|
if (y .ne. 128.0) call abort()
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
n = 1024
|
n = 1024
|
i = 4503599627370495_8
|
i = 4503599627370495_8
|
x = transfer (i, x) !z'000fffffffffffff' Positive denormalized floating-point.
|
x = transfer (i, x) !z'000fffffffffffff' Positive denormalized floating-point.
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
n = -1073
|
n = -1073
|
x = -1024.0
|
x = -1024.0
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
|
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
|
|
|
n = 8
|
n = 8
|
x = -1024.0
|
x = -1024.0
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (y .ne. -128.0) call abort()
|
if (y .ne. -128.0) call abort()
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
|
|
n = 1024
|
n = 1024
|
i = -9218868437227405313_8
|
i = -9218868437227405313_8
|
x = transfer (i, x)!z'800fffffffffffff' Negative denormalized floating-point.
|
x = transfer (i, x)!z'800fffffffffffff' Negative denormalized floating-point.
|
y = set_exponent (x, n)
|
y = set_exponent (x, n)
|
if (exponent (y) .ne. n) call abort()
|
if (exponent (y) .ne. n) call abort()
|
end
|
end
|
|
|