| ! PR middle-end/57393 |
| ! { dg-do compile } |
| ! { dg-options "-g -O2 -ffast-math" } |
| |
| SUBROUTINE pr57393(nn,e,g,t0,t1,t2,t3,t4,t5,t6,t7,& |
| t8,t9,t10,t11,t12,t13,t14,t15,& |
| t16,t17,t18,t19,t20,t21,t22,t23,& |
| t24,t25,t26,t27,t28,t29,t30,& |
| t31,t32,t33,t34,t35,t36,t37,t38,& |
| t39,t40,t41,t42,t43,t44,t45,t46,t47) |
| IMPLICIT REAL*8 (t) |
| INTEGER, PARAMETER :: dp=8 |
| REAL(kind=dp) :: e(nn) |
| DO ii=1,nn |
| t48 = 0.1955555555e2_dp * t1 * t2 + & |
| 0.6000000000e1_dp * t3 * t4 * t5 |
| t49 = 0.1620000000e3_dp * t6 * t7 * t8 + & |
| 0.1080000000e3_dp * t6 * t9 * t5 - & |
| 0.6000000000e1_dp * t10 * t20 * t21 * t55 - & |
| 0.2400000000e2_dp * t10 * t11 * t12 - & |
| 0.1200000000e2_dp * t13 * t14 * t15 |
| t50 = t49 + t16 |
| t51 = (3 * t17 * t18 * t19) + & |
| (t22 * t23 * t19) + (t50 * t19) - & |
| 0.3333333336e0_dp * t24 * t25 |
| t52 = 0.1555555556e1_dp * t26 * t27 * t12 + & |
| (t51 + t28 + t29 + t30) * & |
| 0.3125000000e0_dp * t31 * t32 * t33 * t34 |
| t53 = -0.1000000001e1_dp * t35 * t36 * t5 - & |
| (t37 + t38 + t39 + t52) - & |
| 0.8333333340e-1_dp * t40 * t41 * t42 |
| t54 = -0.1000000001e1_dp * t43 * t44 * t45 - & |
| t47 * (t46 + t53) |
| IF (g >= 3 .OR. g == -3) THEN |
| e(ii) = e(ii) + t54 * t0 |
| END IF |
| END DO |
| END SUBROUTINE pr57393 |