1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
! { 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
|
