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
|
