home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / octave-1.1.1p1-src.tgz / tar.out / fsf / octave / libcruft / blas / dcopy.f < prev    next >
Text File  |  1996-09-28  |  1KB  |  50 lines
  1.       SUBROUTINE  DCOPY(N,DX,INCX,DY,INCY)
  2. C
  3. C     COPIES A VECTOR, X, TO A VECTOR, Y.
  4. C     USES UNROLLED LOOPS FOR INCREMENTS EQUAL TO ONE.
  5. C     JACK DONGARRA, LINPACK, 3/11/78.
  6. C
  7.       DOUBLE PRECISION DX(1),DY(1)
  8.       INTEGER I,INCX,INCY,IX,IY,M,MP1,N
  9. C
  10.       IF(N.LE.0)RETURN
  11.       IF(INCX.EQ.1.AND.INCY.EQ.1)GO TO 20
  12. C
  13. C        CODE FOR UNEQUAL INCREMENTS OR EQUAL INCREMENTS
  14. C          NOT EQUAL TO 1
  15. C
  16.       IX = 1
  17.       IY = 1
  18.       IF(INCX.LT.0)IX = (-N+1)*INCX + 1
  19.       IF(INCY.LT.0)IY = (-N+1)*INCY + 1
  20.       DO 10 I = 1,N
  21.         DY(IY) = DX(IX)
  22.         IX = IX + INCX
  23.         IY = IY + INCY
  24.    10 CONTINUE
  25.       RETURN
  26. C
  27. C        CODE FOR BOTH INCREMENTS EQUAL TO 1
  28. C
  29. C
  30. C        CLEAN-UP LOOP
  31. C
  32.    20 M = MOD(N,7)
  33.       IF( M .EQ. 0 ) GO TO 40
  34.       DO 30 I = 1,M
  35.         DY(I) = DX(I)
  36.    30 CONTINUE
  37.       IF( N .LT. 7 ) RETURN
  38.    40 MP1 = M + 1
  39.       DO 50 I = MP1,N,7
  40.         DY(I) = DX(I)
  41.         DY(I + 1) = DX(I + 1)
  42.         DY(I + 2) = DX(I + 2)
  43.         DY(I + 3) = DX(I + 3)
  44.         DY(I + 4) = DX(I + 4)
  45.         DY(I + 5) = DX(I + 5)
  46.         DY(I + 6) = DX(I + 6)
  47.    50 CONTINUE
  48.       RETURN
  49.       END
  50.