blob: 6866470dd6a4e4d464a278d61cb757aa99d958ea (
plain)
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
|
! Program to test the power (**) operator
program testpow
implicit none
real(kind=4) r, s, two
real(kind=8) :: q
complex(kind=4) :: c, z
real, parameter :: del = 0.0001
integer i, j
i = 2
j = i ** 10
if (abs (j - 1024) .gt. del) call abort
j = i ** (-10)
if (abs (j - 0) .gt. del) call abort
j = i ** 0
if (abs (j - 1) .gt. del) call abort
i = 1
j = i ** 10
if (abs (j - 1) .gt. del) call abort
j = i ** (-10)
if (abs (j - 1) .gt. del) call abort
j = i ** 0
if (abs (j - 1) .gt. del) call abort
i = -1
j = i ** 10
if (abs (j - 1) .gt. del) call abort
j = i ** (-10)
if (abs (j - 1) .gt. del) call abort
j = i ** 0
if (abs (j - 1) .gt. del) call abort
j = i ** 11
if (abs (j - (-1)) .gt. del) call abort
j = i ** (-11)
if (abs (j - (-1)) .gt. del) call abort
c = (2.0, 3.0)
z = c ** 2
if (abs(z - (-5.0, 12.0)) .gt. del) call abort
z = c ** 7
if (abs(z - (6554.0, 4449.0)) .gt. del) call abort
two = 2.0
r = two ** 1
if (abs (r - 2.0) .gt. del) call abort
r = two ** 2
if (abs (r - 4.0) .gt. del) call abort
r = two ** 3
if (abs (r - 8.0) .gt. del) call abort
r = two ** 4
if (abs (r - 16.0) .gt. del) call abort
r = two ** 0
if (abs (r - 1.0) .gt. del) call abort
r = two ** (-1)
if (abs (r - 0.5) .gt. del) call abort
r = two ** (-2)
if (abs (r - 0.25) .gt. del) call abort
r = two ** (-4)
if (abs (r - 0.0625) .gt. del) call abort
s = 3.0
r = two ** s
if (abs (r - 8.0) .gt. del) call abort
s = -3.0
r = two ** s
if (abs (r - 0.125) .gt. del) call abort
i = 3
r = two ** i
if (abs (r - 8.0) .gt. del) call abort
i = -3
r = two ** i
if (abs (r - 0.125) .gt. del) call abort
c = (2.0, 3.0)
c = c ** two
if (abs(c - (-5.0, 12.0)) .gt. del) call abort
end program
|
