blob: b623ac1079405d0deb7e40b9159234df9f39af8f [file] [log] [blame]
! { dg-do compile }
! { dg-options "-Ofast -std=legacy" }
! { dg-additional-options "-march=haswell" { target x86_64-*-* i?86-*-* } }
COMMON /JMSG80/ PI4,PIF,P120,R12,P340,R34,FCS(4,3),E34MAX,
7 IJSAME,KLSAME,IKSMJL
DIMENSION P1(3),FQ(0:5),F1(0:4),F2(0:4),WS(8),WP(8)
DIMENSION VEA(12),VES(9),WES(6)
DIMENSION T(0:20),U(0:20)
DIMENSION T3R(6,3,3,3),T9B(0:20,3,3,3)
DIMENSION F5X(0:12,3,3,3),F6X(0: 6,3,3,3,3)
DIMENSION A31(0:3,0:3),C31(2,0:3),A32(0:3,0:3),C32(2,0:3)
DIMENSION A41(0:3,0:3),C41(2,0:3),A42(0:3,0:3),C42(2,0:3)
DIMENSION A33(16),FIJ00(2),A43(16),FI0K0(2)
DIMENSION SEJJK0( 3),A54(16, 3),C54(2, 3)
DIMENSION A56(0:22,3,0:3),C56(2,0:3)
DIMENSION A60(0:3,0:3),C60(2,0:3),A61(0:3,0:3),C61(2,0:3)
DIMENSION A62(16),FI00L(2),A63(16),F0J0L(2)
DIMENSION A64(0:3,0:3),C64(2, 3),A65(0:3,0:3),C65(2, 3)
DIMENSION A69(0:3, 3),C69(2,0:3),A70(0:3, 3),C70(2,0:3)
DIMENSION A71(18, 3),C71(2, 3)
DIMENSION A72(18, 3),C72(2, 3)
DIMENSION A73(18,0:3),C73(2,0:3)
DIMENSION SE0LKL( 3),A75(16,3),C75(2,0:3)
DIMENSION SE0JLL( 3),A76(16,3),C76(2,0:3)
DIMENSION A77(0:25,3,0:3),C77(2,0:3),A78(0:31,3,0:3),C78(2,0:3)
DIMENSION A79(0:31,3,0:3),C79(2,0:3)
DIMENSION A80(0: 2,2),A81(0:24,3),A82(0:31,2),A83(0:22,2)
DIMENSION A84(0:13,2),A85(0:13,2),A86(0: 6)
DIMENSION S4(0:14),Q4(0:4),FIJKL(2)
IF(XVA.LT.CUG) THEN
ENDIF
F1(M)= FQ0*TMP
F2(M)= FQ0*TMP
XX1=-X12*X43
IF(JI.EQ.1) THEN
DO 255 J=1,3
255 CONTINUE
DO 268 K=1,3
SEJJK00= E0+E(2,2,K,0)+E(3,3,K,0)
A54( 5,K)= A540
268 CONTINUE
297 F5X(3+M,I,I,I)=-R3(M,I,I,I)
DO 299 J=1,3
F5X(3+M,I,I,J)=-R3(M,J,I,I)
299 CONTINUE
DO 300 L=0,M56
DO 300 M=1,3
300 A56(N,M,L)= ZER
A60(2,L)= A600+P34(I,3)*E(I,0,0,L)
A61(0,L)= A610+D1I *E(L,0,0,I)
A61(1,L)= A610+P12(I,3)*E(L,0,0,I)
SEL00L= E(1,0,0,1)+E(2,0,0,2)+E(3,0,0,3)
IF(I.NE.J) THEN
K=6-I-J
F6X(0,J,I,I,I)= ZER
F6X(0,I,J,I,I)= ZER
F6X(0,I,I,J,I)= ZER
F6X(0,I,I,I,J)= ZER
F6X(M,I,I,K,J)= R2(M,K,J)
ENDIF
391 A82( M,N)= ZER
392 A83( M,N)= ZER
A84(M,N)= ZER
A85(M,N)= ZER
397 A86( M)= ZER
DO 399 K=1,3
DO 399 J=1,3
DO 398 M=1,6
T9B(M+ 2,I,J,K)= T3R0
T9B(M+ 8,I,J,K)= T1R(M,I,J,K)
T9B(M+14,I,J,K)= T3R0
398 CONTINUE
399 CONTINUE
417 A77( M,3,K)= A770+F5X0*GEIJKL
445 A81( M,3) = A81( M,3)+T( M)*TMP
IF(K.EQ.L)A81( 5,3)=A81( 5,3)+TMP
IF(I.EQ.J) THEN
DO 447 M=6,11
447 A81( M,3) = A81( M,3)+T( M)*GEIJKL
ENDIF
ENDIF
IF(LK.EQ.1) THEN
IF(JTYPE.NE.4) THEN
DO 510 J=0,3
A31(3,J)= A310+ A310*Y02
A32(3,J)= A320+ A320*Y02
510 CONTINUE
A33( 6)=-AEIJ00*Y1Y+T01
A33( 7)= A330-0*Y01+T01
A33( 8)= A330- A330*Y01
A33(15)= A330+0*Y02
A33(16)= A330+ A330*Y02
ENDIF
A84(12,N)= A84( 7,N)+ A84( 8,N)*Y02
A84(13,N)= A84( 9,N)
A85(10,2)= A85(10,2)- A85(10,1)+ A850
A85(11,2)= A85(11,2)- A85(11,1)+ A850
A85(12,2)= A85(12,2)- A85(12,1)+ A850
A85(13,2)= A85(13,2)- A85(13,1)
Q4(0)= S4( 0)+(S4( 1)+(S4( 2)+(S4( 3)+S4( 4)*Y02)*Y02)*Y02)*Y02
Q4(1)= S4( 5)+(S4( 6)+(S4( 7)+ S4( 8)*Y02)*Y02)*Y02
Q4(2)= S4( 9)+(S4(10)+ S4(11)*Y02)*Y02
Q4(3)= S4(12)+ S4(13)*Y02
Q4(4)= S4(14)
ENDIF
IF(JTYPE.NE.4) THEN
ENDIF
C42(1,M)= T0*F10-T0*F10
C42(2,M)= T0*F20-T0*F20
T(1)= T01+(A43( 4)- A43( 5)*Y04)*Y04
F0J0L(2)= T(0)*F20-T0*F20+T(2)*F20
DO 660 N=1,3
T(0)= A64(3,N)- A64(0,N)- A64(1,N)
T(1)= A640- A640*Y04
C64(1,N)= T0*F10-T0*F10
C64(2,N)= T0*F20-T0*F2(1)
C65(1,N)= T0*F10-T0*F10
C65(2,N)= T0*F20-T0*F2(1)
C70(2,N)= T0*F20-T0*F20
T(2)=(A71(17,N)-(A71(18,N)- A71(16,N)*Y04)*Y04)*XX1
C71(1,N)= T0*F10-T0*F10+T0*F10
C71(2,N)= T0*F20-T0*F20+T0*F20
T(1)=(A720+ A720- A720-T0)*XX1
C72(1,N)= T0*F10-T0*F10+T0*F10
C75(1,N)= T(0)*F10-T0*F1(1)+T(2)*F1(2)
C75(2,N)= T(0)*F20-T0*F2(1)+T(2)*F2(2)
T01 = A76( 6,N)*XX1
T(1)=(T01- A760-(A760- A76( 7,N)- A76( 8,N)
2 -(A760+ A76( 3,N))*Y04)*Y04)*XX1
660 CONTINUE
C73(2,M)= T0*F20+T0*F20+T(2)*F20
2 +(A77(23,1,M)+ A77(24,1,M)*Y04)*Y03
T(2)=(A77(21,2,M)+(A77(22,2,M)+ A77(23,2,M)*Y04)*Y04)*XX1
2 -(A77(24,2,M)-(A77(25,2,M)+ A77(20,2,M)*Y04)*Y04)*Y03
T(3)=(A77(21,3,M)+(A77(22,3,M)+(A77(23,3,M)
2 + A77(24,3,M)*Y04)*Y04)*Y04)*XX1
C77(1,M)= T0*F10-T0*F10-T0*F10+T0*F10
C77(2,M)= T(0)*F20-T(1)*F20-T(2)*F20+T(3)*F20
T(0)=(A780+ A78(24,1,M))*Y3Y+ A780*XX1
T(1)=(A78(23,1,M)+(A78(21,1,M)+A78(22,1,M)*Y04)*Y04)*XX1
2 +(A78(25,1,M)+ A78(26,1,M)*Y04)*Y3Y- A78(27,1,M)*XX2
T(2)=(A78(21,2,M)+(A78(22,2,M)+ A78(28,2,M)*Y04)*Y04)*XX2
2 +(A78(23,2,M)-(A78(24,2,M)+ A78(25,2,M)*Y04)*Y04)*XX1
3 -(A78(29,2,M)-(A78(30,2,M)+ A78(31,2,M)*Y04)*Y04)*Y41
T(3)=(A78(21,3,M)+(A78(22,3,M)+(A78(23,3,M)
2 + A78(24,3,M)*Y04)*Y04)*Y04)*XX2
C78(1,M)= T0*F10-T0*F1(1)-T(2)*F1(2)+T(3)*F1(3)
C78(2,M)= T(0)*F2(0)-T(1)*F2(1)-T(2)*F2(2)+T(3)*F2(3)
T(0)=-A79(24,1,M)*Y3Y+ A79(30,1,M)*H43
T(1)=(A79(21,1,M)-(A79(23,1,M)+ A79(22,1,M)*Y04)*Y04)*XX1
2 +(A79(25,1,M)+ A79(26,1,M)*Y04)*Y3Y- A79(29,1,M)*XX2
T(2)=(A79(21,2,M)+(A79(22,2,M)- A79(28,2,M)*Y04)*Y04)*XX2
2 +(A79(23,2,M)-(A79(24,2,M)+ A79(25,2,M)*Y04)*Y04)*XX1
3 -(A79(29,2,M)-(A79(30,2,M)+ A79(31,2,M)*Y04)*Y04)*Y41
T(3)=(A79(21,3,M)+(A79(22,3,M)+(A79(23,3,M)
2 + A79(24,3,M)*Y04)*Y04)*Y04)*XX2
C79(1,M)= T(0)*F1(0)-T(1)*F1(1)-T(2)*F1(2)+T(3)*F1(3)
C79(2,M)= T(0)*F2(0)-T(1)*F2(1)-T(2)*F2(2)+T(3)*F2(3)
T(0)= A80( 2,1)*Y3Y+ A80( 2,2)*H43
T(1)=(A81(16,1)+(A81(14,1)+A81(15,1)*Y04)*Y04)*XX1
2 +(A81(18,1)+ A81(19,1)*Y04)*Y3Y- A81(20,1)*XX2
T(2)=(A82(21,1)+(A82(22,1)+ A82(28,1)*Y04)*Y04)*XX2
2 +(A82(23,1)-(A82(24,1)+ A82(25,1)*Y04)*Y04)*XX1
3 +(A83(15,2)+(A83(16,2)+ A82(31,1)*Y04)*Y04)*Y41
4 -(A83(17,2)-(A83(18,2)- A83(19,2)*Y04)*Y04)*Y3Y
T(3)=(A84(10,1)+(A84(11,1)+(A84(12,1)+A84(13,1)*Y04)*Y04)*Y04)*XX2
2 +(A85(10,1)+(A85(10,2)+(A85(11,2)+(A85(12,2)
3 + A85(13,2)*Y04)*Y04)*Y04)*Y04)*XX1
T(4)=(Q4(0)+(Q4(1)+(Q4(2)+(Q4(3)+Q4(4)*Y04)*Y04)*Y04)*Y04)*XX2
FIJKL(1)= T(0)*F1(0)-T(1)*F1(1)-T(2)*F1(2)+T(3)*F1(3)+T(4)*F1(4)
FIJKL(2)= T(0)*F2(0)-T(1)*F2(1)-T(2)*F2(2)+T(3)*F2(3)+T(4)*F2(4)
DO 800 ICP=1,2
VE0= VE0+C61(ICP,0)*WP(1)+FI00L(ICP)*WP(2)
2 +F0J0L(ICP)*WP(3)+C77(ICP,0)*WP(4)
3 +C73(ICP,0)*WP(5)+C78(ICP,0)*WP(6)
4 +C79(ICP,0)*WP(7)+FIJKL(ICP)*WP(8)
IF(ICP.EQ.1) THEN
DO 720 M=1,3
VES( M)= VES( M)+VEA( M)
VES(3+M)= VES(3+M)+VEA(3+M)
T01 = VEA(6+M)
T02 = VEA(9+M)
WES( M)=-T01+(T01+T02)*Y03
720 CONTINUE
FE1= X24*VE0
DO 730 M=1,3
T01 = VEA( M)+VEA(3+M)
T02 = VEA(6+M)+VEA(9+M)
WES(3+M)=-T01+(T01+T02)*Z02
730 CONTINUE
ENDIF
800 CONTINUE
WES( M)= WES( M)+P34(M,3)*FE0
WES(3+M)= WES(3+M)-P1(M)*FE1
VES(6+M)= VES(6+M)+WES(3+M)
FCS(3,M)= FCS(3,M)-WES(3+M)*Y03+WES( M)
DO 925 M=1,3
T01= VES( M)
T02= VES(3+M)+VES(6+M)
T01=-T01+(T01+T02)*Y01+P12(M,3)*TMP
FCS(2,M)= FCS(2,M)-T01+VES(6+M)
925 CONTINUE
END