summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gfortran.dg/pr42294.f
blob: 9464379083d7b718e94bca3d0cbedec14283c93f (plain)
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
33
34
35
36
37
38
39
40
41
C PR rtl-optimization/42294
C { dg-do compile { target powerpc*-*-* ia64-*-* x86_64-*-* } }
C { dg-options "-O2 -fselective-scheduling2 -fsel-sched-pipelining -funroll-all-loops" }

      SUBROUTINE ORIEN(IW,NATOT,NTOTORB,NATORB,P,T)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      DIMENSION NATORB(NATOT),P(NTOTORB*(NTOTORB+1)/2)
      DIMENSION T(NTOTORB,NTOTORB)
      DO 9000 IATOM=1,NATOT
         ILAST = NTOTORB
         IF (IATOM.NE.NATOT) ILAST=NATORB(IATOM+1)-1
         DO 8000 IAOI=NATORB(IATOM),ILAST
            DO 7000 IAOJ = IAOI+1,ILAST
               R2 = 0.0D+00
               R3 = 0.0D+00
               DO 6000 INOTA=1,NATOT
                  DO 5000 IK=NATORB(INOTA),NTOTORB
                     IMAI=MAX(IK,IAOI)
                     IMII=MIN(IK,IAOI)
                     IMAJ=MAX(IK,IAOJ)
                     IMIJ=MIN(IK,IAOJ)
                     IKI=(IMAI*(IMAI-1))/2 + IMII
                     IKJ=(IMAJ*(IMAJ-1))/2 + IMIJ
                     PIKI=P(IKI)
                     PIKJ=P(IKJ)
                     R2 = R2 + (PIKI**4)-6*(PIKI*PIKI*PIKJ*PIKJ)+(PIKJ)
 5000             CONTINUE
 6000          CONTINUE
               R2 = (R2/4.0D+00)
               Q = SQRT(R2*R2 + R3*R3)
               IF (Q.LT.1.0D-08) GO TO 7000
               A = COS(THETA)
               B = -SIN(THETA)
               CALL ROT1INT(NTOTORB,IAOI,IAOJ,A,B,P)
 7000       CONTINUE
 8000    CONTINUE
 9000 CONTINUE
      RETURN
      END