1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
! { dg-compile }
! { dg-options "-fgraphite -O -ffast-math" }
MODULE powell
INTEGER, PARAMETER :: dp=8
CONTAINS
SUBROUTINE trsapp (n,npt,xopt,xpt,gq,hq,pq,delta,step,d,g,hd,hs,crvmin)
REAL(dp), DIMENSION(*), INTENT(INOUT) :: step, d, g, hd, hs
LOGICAL :: jump1, jump2
REAL(dp) :: alpha, angle, angtest, bstep, cf, cth, dd, delsq, dg, dhd, &
reduc, sg, sgk, shs, ss, sth, temp, tempa, tempb
DO i=1,n
dd=dd+d(i)**2
END DO
mainloop : DO
IF ( .NOT. jump2 ) THEN
IF ( .NOT. jump1 ) THEN
bstep=temp/(ds+SQRT(ds*ds+dd*temp))
IF (alpha < bstep) THEN
IF (ss < delsq) CYCLE mainloop
END IF
IF (gg <= 1.0e-4_dp*ggbeg) EXIT mainloop
END IF
END IF
END DO mainloop
END SUBROUTINE trsapp
END MODULE powell
! { dg-final { cleanup-modules "powell" } }
|
