1 ! { dg
-additional
-options
"-Ofast -floop-parallelize-all" }
3 subroutine subsm
( n
, x
, xp
, xx
)
4 integer n
, m
, x
(n
),xp
(n
), xx
(n
), gg
(n
), dd_p
6 dd_p
= dd_p
+ (x
(i
) - xx
(i
))*gg
(i
)
8 if ( dd_p
.gt
. 0 ) then
9 call dcopy
( n
, xp
, 1, x
, 1 )