summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gfortran.dg/graphite/pr42185.f90
blob: 6fa6e303677133952ecb9ebe746d4c313f99f496 (plain)
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-do 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" } }