home *** CD-ROM | disk | FTP | other *** search
/ Geek Gadgets 1 / ADE-1.bin / ade-dist / octave-1.1.1p1-src.tgz / tar.out / fsf / octave / libcruft / blas / drot.f < prev    next >
Text File  |  1996-09-28  |  845b  |  37 lines
  1.       SUBROUTINE  DROT (N,DX,INCX,DY,INCY,C,S)
  2. C
  3. C     APPLIES A PLANE ROTATION.
  4. C     JACK DONGARRA, LINPACK, 3/11/78.
  5. C
  6.       DOUBLE PRECISION DX(1),DY(1),DTEMP,C,S
  7.       INTEGER I,INCX,INCY,IX,IY,N
  8. C
  9.       IF(N.LE.0)RETURN
  10.       IF(INCX.EQ.1.AND.INCY.EQ.1)GO TO 20
  11. C
  12. C       CODE FOR UNEQUAL INCREMENTS OR EQUAL INCREMENTS NOT EQUAL
  13. C         TO 1
  14. C
  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.         DTEMP = C*DX(IX) + S*DY(IY)
  21.         DY(IY) = C*DY(IY) - S*DX(IX)
  22.         DX(IX) = DTEMP
  23.         IX = IX + INCX
  24.         IY = IY + INCY
  25.    10 CONTINUE
  26.       RETURN
  27. C
  28. C       CODE FOR BOTH INCREMENTS EQUAL TO 1
  29. C
  30.    20 DO 30 I = 1,N
  31.         DTEMP = C*DX(I) + S*DY(I)
  32.         DY(I) = C*DY(I) - S*DX(I)
  33.         DX(I) = DTEMP
  34.    30 CONTINUE
  35.       RETURN
  36.       END
  37.