1! PR middle-end/57393 2! { dg-do compile } 3! { dg-options "-g -O2 -ffast-math" } 4 5SUBROUTINE pr57393(nn,e,g,t0,t1,t2,t3,t4,t5,t6,t7,& 6 t8,t9,t10,t11,t12,t13,t14,t15,& 7 t16,t17,t18,t19,t20,t21,t22,t23,& 8 t24,t25,t26,t27,t28,t29,t30,& 9 t31,t32,t33,t34,t35,t36,t37,t38,& 10 t39,t40,t41,t42,t43,t44,t45,t46,t47) 11 IMPLICIT REAL*8 (t) 12 INTEGER, PARAMETER :: dp=8 13 REAL(kind=dp) :: e(nn) 14 DO ii=1,nn 15 t48 = 0.1955555555e2_dp * t1 * t2 + & 16 0.6000000000e1_dp * t3 * t4 * t5 17 t49 = 0.1620000000e3_dp * t6 * t7 * t8 + & 18 0.1080000000e3_dp * t6 * t9 * t5 - & 19 0.6000000000e1_dp * t10 * t20 * t21 * t55 - & 20 0.2400000000e2_dp * t10 * t11 * t12 - & 21 0.1200000000e2_dp * t13 * t14 * t15 22 t50 = t49 + t16 23 t51 = (3 * t17 * t18 * t19) + & 24 (t22 * t23 * t19) + (t50 * t19) - & 25 0.3333333336e0_dp * t24 * t25 26 t52 = 0.1555555556e1_dp * t26 * t27 * t12 + & 27 (t51 + t28 + t29 + t30) * & 28 0.3125000000e0_dp * t31 * t32 * t33 * t34 29 t53 = -0.1000000001e1_dp * t35 * t36 * t5 - & 30 (t37 + t38 + t39 + t52) - & 31 0.8333333340e-1_dp * t40 * t41 * t42 32 t54 = -0.1000000001e1_dp * t43 * t44 * t45 - & 33 t47 * (t46 + t53) 34 IF (g >= 3 .OR. g == -3) THEN 35 e(ii) = e(ii) + t54 * t0 36 END IF 37 END DO 38END SUBROUTINE pr57393 39