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 / WAXPY.FOR < prev    next >
Encoding:
Text File  |  1991-04-13  |  660 b   |  28 lines
  1.       SUBROUTINE WAXPY (N, SR, SI, XR, XI, INCX, YR, YI, INCY)
  2.       IMPLICIT NONE
  3. C
  4.       INTEGER N, INCX, INCY
  5.       DOUBLE PRECISION SR, SI, XR(*), XI(*), YR(*), YI(*)
  6. C
  7.       INTEGER I, IX, IY
  8. C
  9.       DOUBLE PRECISION FLOP
  10. C
  11. C
  12.       IF (N.LE.0) RETURN
  13.       IF (SR.EQ.0.0D0 .AND. SI.EQ.0.0D0) RETURN
  14.       IX = 1
  15.       IY = 1
  16.       IF (INCX.LT.0) IX = (-N+1)*INCX+1
  17.       IF (INCY.LT.0) IY = (-N+1)*INCY+1
  18.       DO 10 I = 1, N
  19.         YR(IY) = FLOP (YR(IY)+SR*XR(IX)-SI*XI(IX))
  20.         YI(IY) = YI(IY)+SR*XI(IX)+SI*XR(IX)
  21.         IF (YI(IY).NE.0.0D0) YI(IY) = FLOP (YI(IY))
  22.         IX = IX+INCX
  23.         IY = IY+INCY
  24. 10    CONTINUE
  25. C
  26.       RETURN
  27.       END
  28.