realfunctionsasumf(n,sx,incx)
c
ctakesthesumoftheabsolutevalues.
cusesunrolledloopsforincrementequaltoone.
cjackdongarra,linpack,3/11/78.
cmodified3/93toreturnifincx.le.0.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
realsx(*),stemp
integeri,incx,m,mp1,n,nincx
c
sasumf=0.0e0
stemp=0.0e0
if(n.le.0.or.incx.le.0)return
if(incx.eq.1)go to20
c
ccodeforincrementnotequalto1
c
nincx=n*incx
do10i=1,nincx,incx
stemp=stemp+abs(sx(i))
10continue
sasumf=stemp
return
c
ccodeforincrementequalto1
c
c
cclean-uploop
c
20m=mod(n,6)
if(m.eq.0)go to40
do30i=1,m
stemp=stemp+abs(sx(i))
30continue
if(n.lt.6)go to60
40mp1=m+1
do50i=mp1,n,6
stemp=stemp+abs(sx(i))+abs(sx(i+1))+abs(sx(i+2))
*+abs(sx(i+3))+abs(sx(i+4))+abs(sx(i+5))
50continue
60sasumf=stemp
return
end