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 / WDOTUR.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  591 b   |  30 lines
  1.       DOUBLE PRECISION FUNCTION WDOTUR (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 = FLOP (S+XR(IX)*YR(IY)-XI(IX)*YI(IY))
  21.         IX = IX+INCX
  22.         IY = IY+INCY
  23. 10    CONTINUE
  24. C
  25. 20    CONTINUE
  26.       WDOTUR = S
  27. C
  28.       RETURN
  29.       END
  30.