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 / caxpy.f next >
Encoding:
Text File  |  1994-08-02  |  1.2 KB  |  51 lines
  1.       SUBROUTINE CAXPY( N, CA, CX, INCX, CY, INCY )
  2. *
  3. *     constant times a vector plus a vector.
  4. *     jack dongarra, linpack, 3/11/78.
  5. *
  6. *     .. Scalar Arguments ..
  7.       INTEGER           INCX, INCY, N
  8.       COMPLEX           CA
  9. *     ..
  10. *     .. Array Arguments ..
  11.       COMPLEX           CX( 1 ), CY( 1 )
  12. *     ..
  13. *     .. Local Scalars ..
  14.       INTEGER           I, IX, IY
  15. *     ..
  16. *     .. Intrinsic Functions ..
  17.       INTRINSIC         ABS, AIMAG, REAL
  18. *     ..
  19. *     .. Executable Statements ..
  20. *
  21.       IF( N.LE.0 )
  22.      $   RETURN
  23.       IF( ABS( REAL( CA ) )+ABS( AIMAG( CA ) ).EQ.0.0 )
  24.      $   RETURN
  25.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  26.      $   GO TO 20
  27. *
  28. *        code for unequal increments or equal increments
  29. *          not equal to 1
  30. *
  31.       IX = 1
  32.       IY = 1
  33.       IF( INCX.LT.0 )
  34.      $   IX = ( -N+1 )*INCX + 1
  35.       IF( INCY.LT.0 )
  36.      $   IY = ( -N+1 )*INCY + 1
  37.       DO 10 I = 1, N
  38.          CY( IY ) = CY( IY ) + CA*CX( IX )
  39.          IX = IX + INCX
  40.          IY = IY + INCY
  41.    10 CONTINUE
  42.       RETURN
  43. *
  44. *        code for both increments equal to 1
  45. *
  46.    20 DO 30 I = 1, N
  47.          CY( I ) = CY( I ) + CA*CX( I )
  48.    30 CONTINUE
  49.       RETURN
  50.       END
  51.