home *** CD-ROM | disk | FTP | other *** search
/ The Fred Fish Collection 1.5 / ffcollection-1-5-1992-11.iso / ff_progs / libs / matlab.lzh / MATLAB / MATLAB.LZH / Source / MatLab / WDOTCI.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  625 b   |  31 lines
  1.       DOUBLE PRECISION FUNCTION WDOTCI (N, XR, XI, INCX, YR, YI, INCY)
  2.       IMPLICIT NONE
  3. C
  4.       INTEGER N, INCX, INCY
  5.       DOUBLE PRECISION XR(*), XI(*), YR(*), YI(*)
  6. C
  7.       INTEGER I, IX, IY
  8.       DOUBLE PRECISION S
  9. C
  10.       DOUBLE PRECISION FLOP
  11. C
  12. C
  13.       S = 0.0D0
  14.       IF (N.LE.0) GO TO 20
  15.       IX = 1
  16.       IY = 1
  17.       IF (INCX.LT.0) IX = (-N+1)*INCX+1
  18.       IF (INCY.LT.0) IY = (-N+1)*INCY+1
  19.       DO 10 I = 1, N
  20.          S = S+XR(IX)*YI(IY)-XI(IX)*YR(IY)
  21.          IF (S.NE.0.0D0) S = FLOP (S)
  22.          IX = IX+INCX
  23.          IY = IY+INCY
  24. 10    CONTINUE
  25. C
  26. 20    CONTINUE
  27.       WDOTCI = S
  28. C
  29.       RETURN
  30.       END
  31.