2 ! { dg-additional-options "-ffast-math" }
3 SUBROUTINE influence_factor ( gftype
, error
)
4 INTEGER, PARAMETER :: dp
=8
5 INTEGER :: k
,n
,lb(3),ub(3),dim
,pt
6 COMPLEX(KIND
=dp
) :: b_m
, exp_m
, sum_m
8 DO pt
= lb (dim
), ub (dim
)
9 sum_m
= CMPLX ( 0.0_dp
, 0.0_dp
,KIND
=dp
)
10 b_m
= exp_m
** ( n
- 1 ) / sum_m
13 END SUBROUTINE influence_factor