3 ! { dg-options "-g -O2 -ffast-math" }
5 SUBROUTINE 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
)
12 INTEGER, PARAMETER :: dp
=8
13 REAL(kind
=dp
) :: e(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
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
- &
34 IF (g
>= 3 .OR
. g
== -3) THEN
35 e(ii
) = e(ii
) + t54
* t0
38 END SUBROUTINE pr57393