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
|
module testmod
implicit none
contains
subroutine foo(n)
integer, intent(in) :: n
real :: r(0:n,-n:n), a(0:n,-n:n), dj
integer :: k, j
! initialize with some dummy values
do j = -n, n
a(:, j) = j
r(:,j) = j + 1
end do
! here be dragons
do k = 0, n
dj = r(k, k - 2) * a(k, k - 2)
r(k,k) = a(k, k - 1) * dj
enddo
if (r(0,0) .ne. -2.) call abort
end subroutine
end module
program test
use testmod
implicit none
call foo(5)
end program
|
