complexfunctioncdotcf(n,cx,incx,cy,incy)
c
cformsthedotproductoftwovectors,conjugatingthefirst
cvector.
cjackdongarra,linpack,3/11/78.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
complexcx(*),cy(*),ctemp
integeri,incx,incy,ix,iy,n
c
ctemp=(0.0,0.0)
cdotcf=(0.0,0.0)
if(n.le.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
ctemp=ctemp+conjg(cx(ix))*cy(iy)
ix=ix+incx
iy=iy+incy
10continue
cdotcf=ctemp
return
c
ccodeforbothincrementsequalto1
c
20do30i=1,n
ctemp=ctemp+conjg(cx(i))*cy(i)
30continue
cdotcf=ctemp
return
end