REAL*8functiondzaminf(n,zx,incx)
c
cfindstheindexofelementhavingmin.absolutevalue.
cjackdongarra,1/15/85.
cmodified3/93toreturnifincx.le.0.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
COMPLEX*16zx(*)
integeri,incx,ix,n
double precisiondcabs1
c
dzaminf=0.
if(n.lt.1.or.incx.le.0)return
dzaminf=dcabs1(zx(1))
if(n.eq.1)return
if(incx.eq.1)go to20
c
ccodeforincrementnotequalto1
c
ix=1
dzaminf=dcabs1(zx(1))
ix=ix+incx
do10i=2,n
if(dcabs1(zx(ix)).ge.dzaminf)go to5
dzaminf=dcabs1(zx(ix))
5ix=ix+incx
10continue
return
c
ccodeforincrementequalto1
c
20dzaminf=dcabs1(zx(1))
do30i=2,n
if(dcabs1(zx(i)).ge.dzaminf)go to30
dzaminf=dcabs1(zx(i))
30continue
return
end