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