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
|
program complex_07
implicit none
integer, parameter :: real_part = 42
real, parameter :: img_part = 3.14
complex :: x = cmplx(real_part)
complex :: y = cmplx(img_part)
complex :: z = cmplx(real_part, img_part)
double complex :: c1, x2, y2, z2
integer :: i = 42
real :: x1 = 3.14
complex :: z1
z1 = cmplx(i, x1)
c1 = (12, 24)
print *, c1
x2 = dcmplx(real_part)
y2 = dcmplx(img_part)
z2 = dcmplx(real_part, img_part)
if (abs(x - (real_part, 0)) > 1e-5) error stop
if (abs(y - (img_part, 0)) > 1e-5) error stop
if (abs(z - (real_part, img_part)) > 1e-5) error stop
if (abs(x2 - (real_part, 0)) > 1e-5) error stop
if (abs(y2 - (img_part, 0)) > 1e-5) error stop
if (abs(z2 - (real_part, img_part)) > 1e-5) error stop
print *, dcmplx(i)
if (abs(dcmplx(i) - (42.000000,0.000000)) > 1e-5) error stop
print *, dcmplx(x1)
if (abs(dcmplx(x1) - (3.140000,0.000000)) > 1e-5) error stop
print *, dcmplx(z1)
if (abs(dcmplx(z1) - (42.000000,3.140000)) > 1e-5) error stop
print *, dcmplx(x1, i)
if (abs(dcmplx(x1, i) - (3.140000,42.000000)) > 1e-5) error stop
print *, x, y, z
end program
|
