home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / daxpy.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.5 KB  |  67 lines
  1.       SUBROUTINE DAXPY( N, DA, DX, INCX, DY, INCY )
  2. *
  3. *     constant times a vector plus a vector.
  4. *     uses unrolled loops for increments equal to one.
  5. *     jack dongarra, linpack, 3/11/78.
  6. *
  7. *     .. Scalar Arguments ..
  8.       INTEGER           INCX, INCY, N
  9.       DOUBLE PRECISION  DA
  10. *     ..
  11. *     .. Array Arguments ..
  12.       DOUBLE PRECISION  DX( 1 ), DY( 1 )
  13. *     ..
  14. *     .. Local Scalars ..
  15.       INTEGER           I, IX, IY, M, MP1
  16. *     ..
  17. *     .. Intrinsic Functions ..
  18.       INTRINSIC         MOD
  19. *     ..
  20. *     .. Executable Statements ..
  21. *
  22.       IF( N.LE.0 )
  23.      $   RETURN
  24.       IF( DA.EQ.0.0D0 )
  25.      $   RETURN
  26.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  27.      $   GO TO 20
  28. *
  29. *        code for unequal increments or equal increments
  30. *          not equal to 1
  31. *
  32.       IX = 1
  33.       IY = 1
  34.       IF( INCX.LT.0 )
  35.      $   IX = ( -N+1 )*INCX + 1
  36.       IF( INCY.LT.0 )
  37.      $   IY = ( -N+1 )*INCY + 1
  38.       DO 10 I = 1, N
  39.          DY( IY ) = DY( IY ) + DA*DX( IX )
  40.          IX = IX + INCX
  41.          IY = IY + INCY
  42.    10 CONTINUE
  43.       RETURN
  44. *
  45. *        code for both increments equal to 1
  46. *
  47. *
  48. *        clean-up loop
  49. *
  50.    20 M = MOD( N, 4 )
  51.       IF( M.EQ.0 )
  52.      $   GO TO 40
  53.       DO 30 I = 1, M
  54.          DY( I ) = DY( I ) + DA*DX( I )
  55.    30 CONTINUE
  56.       IF( N.LT.4 )
  57.      $   RETURN
  58.    40 MP1 = M + 1
  59.       DO 50 I = MP1, N, 4
  60.          DY( I ) = DY( I ) + DA*DX( I )
  61.          DY( I+1 ) = DY( I+1 ) + DA*DX( I+1 )
  62.          DY( I+2 ) = DY( I+2 ) + DA*DX( I+2 )
  63.          DY( I+3 ) = DY( I+3 ) + DA*DX( I+3 )
  64.    50 CONTINUE
  65.       RETURN
  66.       END
  67.