1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
! { dg-do compile }
! { dg-options "-O3" }
function pr52701 (x, z, e, f, g, l)
integer a, b, c, d, e, f, g, i, j, l, pr52701
double precision x(e), z(e*e)
do i = l, f
do j = l, i
d = 0
do a = 1, g
c = a - g
do b = 1, g
d = d + 1
c = c + g
z(d) = z(d) / (x(i) + x(j) - x(f + a) - x(f + b))
end do
end do
end do
end do
pr52701 = c
end
|