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 / cdotc.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.3 KB  |  54 lines
  1.       COMPLEX          FUNCTION CDOTC( N, CX, INCX, CY, INCY )
  2. *
  3. *     forms the dot product of two vectors, conjugating the first
  4. *     vector.
  5. *     jack dongarra, linpack,  3/11/78.
  6. *
  7. *     .. Scalar Arguments ..
  8.       INTEGER                          INCX, INCY, N
  9. *     ..
  10. *     .. Array Arguments ..
  11.       COMPLEX                          CX( 1 ), CY( 1 )
  12. *     ..
  13. *     .. Local Scalars ..
  14.       INTEGER                          I, IX, IY
  15.       COMPLEX                          CTEMP
  16. *     ..
  17. *     .. Intrinsic Functions ..
  18.       INTRINSIC                        CONJG
  19. *     ..
  20. *     .. Executable Statements ..
  21. *
  22.       CTEMP = ( 0.0, 0.0 )
  23.       CDOTC = ( 0.0, 0.0 )
  24.       IF( N.LE.0 )
  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.          CTEMP = CTEMP + CONJG( CX( IX ) )*CY( IY )
  40.          IX = IX + INCX
  41.          IY = IY + INCY
  42.    10 CONTINUE
  43.       CDOTC = CTEMP
  44.       RETURN
  45. *
  46. *        code for both increments equal to 1
  47. *
  48.    20 DO 30 I = 1, N
  49.          CTEMP = CTEMP + CONJG( CX( I ) )*CY( I )
  50.    30 CONTINUE
  51.       CDOTC = CTEMP
  52.       RETURN
  53.       END
  54.