double precisionfunctiondzasumf(n,zx,incx)
c
ctakesthesumoftheabsolutevalues.
cjackdongarra,3/11/78.
cmodified3/93toreturnifincx.le.0.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
doublecomplexzx(*)
double precisionstemp,dcabs1
integeri,incx,ix,n
c
dzasumf=0.0d0
stemp=0.0d0
if(n.le.0.or.incx.le.0)return
if(incx.eq.1)go to20
c
ccodeforincrementnotequalto1
c
ix=1
do10i=1,n
stemp=stemp+dcabs1(zx(ix))
ix=ix+incx
10continue
dzasumf=stemp
return
c
ccodeforincrementequalto1
c
20do30i=1,n
stemp=stemp+dcabs1(zx(i))
30continue
dzasumf=stemp
return
end