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