1! { dg-do compile } 2! { dg-options "-O3" } 3 subroutine hesfcn(n, x, h, ldh) 4 integer n,ldh 5 double precision x(n), h(ldh) 6 7 integer i,j,k,kj 8 double precision th,u1,u2,v2 9 10 kj = 0 11 do 770 j = 1, n 12 kj = kj - j 13 do 760 k = 1, j 14 kj = kj + 1 15 v2 = 2 * x(k) - 1 16 u1 = 0 17 u2 = 2 18 do 750 i = 1, n 19 h(kj) = h(kj) + u2 20 th = 4 * v2 + u2 - u1 21 u1 = u2 22 u2 = th 23 th = v2 - 1 24 750 continue 25 760 continue 26 770 continue 27 28 end 29