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