subroutinescopyf(n,sx,incx,sy,incy)
c
ccopiesavector,x,toavector,y.
cusesunrolledloopsforincrementsequalto1.
cjackdongarra,linpack,3/11/78.
cmodified12/3/93,array(1)declarationschangedtoarray(*)
c
realsx(*),sy(*)
integeri,incx,incy,ix,iy,m,mp1,n
c
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
sy(iy)=sx(ix)
ix=ix+incx
iy=iy+incy
10continue
return
c
ccodeforbothincrementsequalto1
c
c
cclean-uploop
c
20m=mod(n,7)
if(m.eq.0)go to40
do30i=1,m
sy(i)=sx(i)
30continue
if(n.lt.7)return
40mp1=m+1
do50i=mp1,n,7
sy(i)=sx(i)
sy(i+1)=sx(i+1)
sy(i+2)=sx(i+2)
sy(i+3)=sx(i+3)
sy(i+4)=sx(i+4)
sy(i+5)=sx(i+5)
sy(i+6)=sx(i+6)
50continue
return
end