Subversion Repositories openrisc_2011-10-31

! { 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