gcc/testsuite/gfortran.dg/pr92094.f90

   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