File: pr92094.f90

package info (click to toggle)
gcc-arm-none-eabi 15%3A12.2.rel1-1
  • links: PTS, VCS
  • area: main
  • in suites: bookworm
  • size: 959,712 kB
  • sloc: cpp: 3,275,382; ansic: 2,061,766; ada: 840,956; f90: 208,513; makefile: 76,132; asm: 73,433; xml: 50,448; exp: 34,146; sh: 32,436; objc: 15,637; fortran: 14,012; python: 11,991; pascal: 6,787; awk: 4,779; perl: 3,054; yacc: 338; ml: 285; lex: 201; haskell: 122
file content (28 lines) | stat: -rw-r--r-- 608 bytes parent folder | download | duplicates (2) 
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