1c { dg-do compile } 2 subroutine saxpy(n,sa,sx,incx,sy,incy) 3C 4C constant times a vector plus a vector. 5C uses unrolled loop for increments equal to one. 6C jack dongarra, linpack, 3/11/78. 7C modified 12/3/93, array(1) declarations changed to array(*) 8C 9 real sx(*),sy(*),sa 10 integer i,incx,incy,ix,iy,m,mp1,n 11C 12C -ffast-math ICE provoked by this conditional 13 if(sa /= 0.0)then 14C 15C code for both increments equal to 1 16C 17 do i= 1,n 18 sy(i)= sy(i)+sa*sx(i) 19 enddo 20 endif 21 return 22 end 23