File: pr87217.f

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 (32 lines) | stat: -rw-r--r-- 791 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
29
30
31
32
 
! { dg-do compile }
! { dg-options "-O3" }
      implicit real*8 (a-h,o-z)
      common clop6(3),dps(6),aml6(6,6)
      dimension y1(3)
      dimension dclo(3)
      dimension dx(3),dy(3)
      save
      do 80 ii=1,itco
        y1(3)=dps(1)
        do 40 l=1,3
          dy(l)=clop6(l)-y1(l)
   40   continue
        dczp=abs(dy(3))
        if(dcx.le.c1m10.and.dcz.le.c1m10.and.dcxp.le.c1m10.and.dczp
     +  .le.c1m10.and.dcy.le.c1m10.and.dcyp.le.c1m10) goto 90
   80 continue
      write(6) itco
      ii=itco
   90 continue
      if(ii.ne.itco) then
        do 65 k=1,3
          do 55 j=1,3
            jj=2*j
            kk=2*k
            dclo(k)=aml6(kk-1,jj-1)*dx(j)+dclo(k)
            dclo(k)=aml6(kk-1,jj)*dy(j)+dclo(k)
   55     continue
   65   continue
      endif
      end