| ! { 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 |