1 ! { dg-options "-O2 -g" }
4 INTEGER, PARAMETER :: dp
=8
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
10 INTEGER, INTENT(in
) :: npt
14 sum
=sum
+bmat(jp
,k
)*d(k
)