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