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 / zdotc.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.3 KB  |  52 lines
  1.       DOUBLE COMPLEX   FUNCTION ZDOTC( N, ZX, INCX, ZY, INCY )
  2. *
  3. *     forms the dot product of a vector.
  4. *     jack dongarra, 3/11/78.
  5. *
  6. *     .. Scalar Arguments ..
  7.       INTEGER                          INCX, INCY, N
  8. *     ..
  9. *     .. Array Arguments ..
  10.       COMPLEX*16                       ZX( 1 ), ZY( 1 )
  11. *     ..
  12. *     .. Local Scalars ..
  13.       INTEGER                          I, IX, IY
  14.       COMPLEX*16                       ZTEMP
  15. *     ..
  16. *     .. Intrinsic Functions ..
  17.       INTRINSIC                        DCONJG
  18. *     ..
  19. *     .. Executable Statements ..
  20.       ZTEMP = ( 0.0D0, 0.0D0 )
  21.       ZDOTC = ( 0.0D0, 0.0D0 )
  22.       IF( N.LE.0 )
  23.      $   RETURN
  24.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  25.      $   GO TO 20
  26. *
  27. *        code for unequal increments or equal increments
  28. *          not equal to 1
  29. *
  30.       IX = 1
  31.       IY = 1
  32.       IF( INCX.LT.0 )
  33.      $   IX = ( -N+1 )*INCX + 1
  34.       IF( INCY.LT.0 )
  35.      $   IY = ( -N+1 )*INCY + 1
  36.       DO 10 I = 1, N
  37.          ZTEMP = ZTEMP + DCONJG( ZX( IX ) )*ZY( IY )
  38.          IX = IX + INCX
  39.          IY = IY + INCY
  40.    10 CONTINUE
  41.       ZDOTC = ZTEMP
  42.       RETURN
  43. *
  44. *        code for both increments equal to 1
  45. *
  46.    20 DO 30 I = 1, N
  47.          ZTEMP = ZTEMP + DCONJG( ZX( I ) )*ZY( I )
  48.    30 CONTINUE
  49.       ZDOTC = ZTEMP
  50.       RETURN
  51.       END
  52.