1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
|
!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)) call abort()
n = 8
x = 1024.0
y = set_exponent (x, n)
if (exponent (y) .ne. n) call abort()
n = 128
i = 8388607
x = transfer (i, x) ! z'007fffff' Positive denormalized floating-point.
y = set_exponent (x, n)
if (exponent (y) .ne. n) call abort()
n = -148
x = -1024.0
y = set_exponent (x, n)
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
n = 8
x = -1024.0
y = set_exponent (x, n)
if (y .ne. -128.0) call abort()
if (exponent (y) .ne. n) call abort()
n = 128
i = -2139095041
x = transfer (i, x) ! z'807fffff' Negative denormalized floating-point.
y = set_exponent (x, n)
if (exponent (y) .ne. n) call abort()
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)) call abort()
n = 8
x = 1024.0_8
y = set_exponent (x, n)
if (y .ne. 128.0) call abort()
if (exponent (y) .ne. n) call abort()
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) call abort()
n = -1073
x = -1024.0
y = set_exponent (x, n)
if ((y .ne. 0.0) .and. (exponent (y) .ne. n)) call abort()
n = 8
x = -1024.0
y = set_exponent (x, n)
if (y .ne. -128.0) call abort()
if (exponent (y) .ne. n) call abort()
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) call abort()
end
|