home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / linpack / scopy.for < prev    next >
Encoding:
Text File  |  1984-01-05  |  1.1 KB  |  50 lines
  1.       SUBROUTINE SCOPY(N,SX,INCX,SY,INCY)
  2. C
  3. C     COPIES A VECTOR, X, TO A VECTOR, Y.
  4. C     USES UNROLLED LOOPS FOR INCREMENTS EQUAL TO 1.
  5. C     JACK DONGARRA, LINPACK, 3/11/78.
  6. C
  7.       REAL SX(1),SY(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.         SY(IY) = SX(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.         SY(I) = SX(I)
  36.    30 CONTINUE
  37.       IF( N .LT. 7 ) RETURN
  38.    40 MP1 = M + 1
  39.       DO 50 I = MP1,N,7
  40.         SY(I) = SX(I)
  41.         SY(I + 1) = SX(I + 1)
  42.         SY(I + 2) = SX(I + 2)
  43.         SY(I + 3) = SX(I + 3)
  44.         SY(I + 4) = SX(I + 4)
  45.         SY(I + 5) = SX(I + 5)
  46.         SY(I + 6) = SX(I + 6)
  47.    50 CONTINUE
  48.       RETURN
  49.       END
  50.