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
|
pure function linspace_n_1_cdp_cdp(x, y, n) result(res)
integer, intent(in) :: n
complex :: res(n)
real, intent(in) :: x(n)
real, intent(in) :: y(n)
res = cmplx(x, y)
end function linspace_n_1_cdp_cdp
program complex_17
implicit none
integer :: i
real :: x(10), y(10)
complex :: z(10)
interface
pure function linspace_n_1_cdp_cdp(x, y, n) result(res)
integer, intent(in) :: n
complex :: res(n)
real, intent(in) :: x(n)
real, intent(in) :: y(n)
end function linspace_n_1_cdp_cdp
end interface
x = 1.0
y = 2.0
z = linspace_n_1_cdp_cdp(x, y, 10)
print *, z
print *, abs(z)
do i = 1, 10
if (abs(abs(z(i)) - 2.236068) > 1e-8) error stop
end do
end program
|
