| ! { dg-options "-O1 -fgraphite" } | |
| MODULE powell | |
| INTEGER, PARAMETER :: dp=8 | |
| CONTAINS | |
| SUBROUTINE newuob (n,npt,x,rhobeg,rhoend,maxfun,xbase,& | |
| xopt,xnew,xpt,fval,gq,hq,pq,bmat,zmat,ndim,d,vlag,w,opt) | |
| REAL(dp), DIMENSION(npt, *), & | |
| INTENT(inout) :: xpt | |
| REAL(dp), DIMENSION(*), INTENT(inout) :: fval, gq, hq, pq | |
| 120 IF (dsq <= 1.0e-3_dp*xoptsq) THEN | |
| DO k=1,npt | |
| DO i=1,n | |
| gq(i)=gq(i)+temp*xpt(k,i) | |
| END DO | |
| END DO | |
| END IF | |
| END SUBROUTINE newuob | |
| END MODULE powell |