1 ! Program to test the power (**) operator
6 complex(kind
=4) :: c
, z
7 real, parameter :: del
= 0.0001
12 if (abs (j
- 1024) .gt
. del
) call abort
14 if (abs (j
- 0) .gt
. del
) call abort
16 if (abs (j
- 1) .gt
. del
) call abort
19 if (abs (j
- 1) .gt
. del
) call abort
21 if (abs (j
- 1) .gt
. del
) call abort
23 if (abs (j
- 1) .gt
. del
) call abort
26 if (abs (j
- 1) .gt
. del
) call abort
28 if (abs (j
- 1) .gt
. del
) call abort
30 if (abs (j
- 1) .gt
. del
) call abort
32 if (abs (j
- (-1)) .gt
. del
) call abort
34 if (abs (j
- (-1)) .gt
. del
) call abort
38 if (abs(z
- (-5.0, 12.0)) .gt
. del
) call abort
40 if (abs(z
- (6554.0, 4449.0)) .gt
. del
) call abort
45 if (abs (r
- 2.0) .gt
. del
) call abort
47 if (abs (r
- 4.0) .gt
. del
) call abort
49 if (abs (r
- 8.0) .gt
. del
) call abort
51 if (abs (r
- 16.0) .gt
. del
) call abort
53 if (abs (r
- 1.0) .gt
. del
) call abort
55 if (abs (r
- 0.5) .gt
. del
) call abort
57 if (abs (r
- 0.25) .gt
. del
) call abort
59 if (abs (r
- 0.0625) .gt
. del
) call abort
62 if (abs (r
- 8.0) .gt
. del
) call abort
65 if (abs (r
- 0.125) .gt
. del
) call abort
68 if (abs (r
- 8.0) .gt
. del
) call abort
71 if (abs (r
- 0.125) .gt
. del
) call abort
74 if (abs(c
- (-5.0, 12.0)) .gt
. del
) call abort