1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
! { dg-options "-O2 -g" }
MODULE powell
INTEGER, PARAMETER :: dp=8
CONTAINS
SUBROUTINE newuob (n, bmat, ndim, d, vlag, w, npt)
REAL(dp), DIMENSION(ndim, *), INTENT(inout) :: bmat
REAL(dp), DIMENSION(*), INTENT(inout) :: d, vlag, w
REAL(dp) :: sum
INTEGER, INTENT(in) :: npt
DO j=1,n
jp=npt+j
DO k=1,n
sum=sum+bmat(jp,k)*d(k)
END DO
vlag(jp)=sum
END DO
END SUBROUTINE newuob
END MODULE powell
|
