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
|
! { dg-do compile }
! { dg-options "-O3" }
subroutine hesfcn(n, x, h, ldh)
integer n,ldh
double precision x(n), h(ldh)
integer i,j,k,kj
double precision th,u1,u2,v2
kj = 0
do 770 j = 1, n
kj = kj - j
do 760 k = 1, j
kj = kj + 1
v2 = 2 * x(k) - 1
u1 = 0
u2 = 2
do 750 i = 1, n
h(kj) = h(kj) + u2
th = 4 * v2 + u2 - u1
u1 = u2
u2 = th
th = v2 - 1
750 continue
760 continue
770 continue
end
|