| !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 |