! { 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 |