3 ! { dg-options "-O3 -ffast-math" }
5 SUBROUTINE urep_egr(erep
,derep
,surr
)
6 INTEGER, PARAMETER :: dp
=8
7 REAL(dp
), INTENT(inout
) :: erep
, derep(3)
8 REAL(dp
), INTENT(in
) :: surr(2)
10 INTEGER :: isp
,spdim
,jsp
,nsp
11 IF (n_urpoly
> 0) THEN
12 IF (r
< spxr(1,1)) THEN
13 ispg
: DO isp
= 1,spdim
! condition ca)
14 IF (isp
/= spdim
) THEN
15 nsp
= 5 ! condition cb
19 erep
= erep
+ surr(jsp
-3)*rz
**(jsp
)
26 END SUBROUTINE urep_egr