1 ! { dg-options "-O2 -fgraphite-identity -fno-loop-block -fno-loop-interchange -fno-loop-strip-mine" }
4 INTEGER, PARAMETER :: dp
=KIND(0.0D0
)
6 FUNCTION basym (a
, b
, lambda
, eps
) RESULT(fn_val
)
7 REAL(dp
) :: a0(21), b0(21), bsum
, c(21), d(21), dsum
, &
8 j0
, j1
, r
, r0
, r1
, s
, sum
, t
, t0
, t1
, &
9 u
, w
, w0
, z
, z0
, z2
, zn
, znm1
18 bsum
= bsum
+ (j
*r
- mmj
)*a0(j
)*b0(mmj
)
20 b0(m
) = r
*a0(m
) + bsum
/m
22 c(i
) = b0(i
)/(i
+ 1.0e0_dp
)
30 END MODULE beta_gamma_psi