subroutinezaxpyf(n,za,zx,incx,zy,incy)
c
cconstanttimesavectorplusavector.
cjackdongarra,3/11/78.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
doublecomplexzx(*),zy(*),za
integeri,incx,incy,ix,iy,n
double precisiondcabs1
if(n.le.0)return
if(dcabs1(za).eq.0.0d0)return
if(incx.eq.1.and.incy.eq.1)go to20
c
ccodeforunequalincrementsorequalincrements
cnotequalto1
c
ix=1
iy=1
if(incx.lt.0)ix=(-n+1)*incx+1
if(incy.lt.0)iy=(-n+1)*incy+1
do10i=1,n
zy(iy)=zy(iy)+za*zx(ix)
ix=ix+incx
iy=iy+incy
10continue
return
c
ccodeforbothincrementsequalto1
c
20do30i=1,n
zy(i)=zy(i)+za*zx(i)
30continue
return
end