1! { dg-do compile }
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
7    DO k = 0, n-2
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
11       END DO
12    END DO
13  END SUBROUTINE influence_factor
14