1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
|
! 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
|
