1 subroutine saxpy
(n
,sa
,sx
,incx
,sy
,incy
)
3 C constant times a vector plus a vector.
4 C uses unrolled loop for increments equal to one.
5 C jack dongarra, linpack, 3/11/78.
6 C modified 12/3/93, array(1) declarations changed to array(*)
9 integer i
,incx
,incy
,ix
,iy
,m
,mp1
,n
11 C -ffast-math ICE provoked by this conditional
14 C code for both increments equal to 1