subroutinecaxpycf(n,ca,cx,incx,cy,incy)
c
cconstanttimesavectorplusavector.
cjackdongarra,linpack,3/11/78.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
complexcx(*),cy(*),ca
integeri,incx,incy,ix,iy,n
INTRINSICconjg
c
if(n.le.0)return
if(abs(real(ca))+abs(aimag(ca)).eq.0.0)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
cy(iy)=cy(iy)+ca*conjg(cx(ix))
ix=ix+incx
iy=iy+incy
10continue
return
c
ccodeforbothincrementsequalto1
c
20do30i=1,n
cy(i)=cy(i)+ca*conjg(cx(i))
30continue
return
end