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
|
subroutine sub(a)
double precision :: a
dimension :: a(16)
print *, a
if (abs(a(1) - 1) > 1d-15) error stop
if (abs(a(2) - 2) > 1d-15) error stop
if (.not. all(abs(a(3:) - 1) < 1d-15)) error stop
end subroutine
program main
double precision :: a(16,3)
integer :: j
j = 3
a = 1
a(2,3) = 2
call sub(a(1,j))
a(2,2) = 2
j = 2
call sub2(a(1,j))
contains
subroutine sub2(a)
real(8) :: a(16)
print *, a
if (abs(a(1) - 1) > 1d-15) error stop
if (abs(a(2) - 2) > 1d-15) error stop
if (.not. all(abs(a(3:) - 1) < 1d-15)) error stop
end subroutine
end program
|
