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
|
program flip_sign
implicit none
real :: rxsp = 5.5, epsrsp = 1e-6
real(8) :: rxdp = 5.5, epsrdp = 1e-6
integer :: ixsp = 5, epsisp = 16
integer(8) :: ixdp = 5, epsidp = 16
rxsp = rxsp * sign(1._4, epsrsp)
print *, rxsp
if (abs(rxsp - 5.5) > epsrsp) error stop
rxsp = rxsp * sign(1._4, -epsrsp)
print *, rxsp
if (abs(rxsp + 5.5) > epsrsp) error stop
rxdp = rxdp * sign(1._8, epsrdp)
print *, rxdp
if (abs(rxdp - 5.5) > epsrdp) error stop
rxdp = rxdp * sign(1._8, -epsrdp)
print *, rxdp
if (abs(rxdp + 5.5) > epsrdp) error stop
ixsp = ixsp * sign(1_4, epsisp)
print *, ixsp
if (ixsp /= 5) error stop
ixdp = ixdp * sign(1_8, epsidp)
print *, ixdp
if (ixdp /= 5) error stop
end program
|
