home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / library / linpack / srot.for < prev    next >
Encoding:
Text File  |  1984-01-06  |  869 b   |  37 lines
  1.       SUBROUTINE SROT (N,SX,INCX,SY,INCY,C,S)
  2. C
  3. C     APPLIES A PLANE ROTATION.
  4. C     JACK DONGARRA, LINPACK, 3/11/78.
  5. C
  6.       REAL SX(1),SY(1),STEMP,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.         STEMP = C*SX(IX) + S*SY(IY)
  21.         SY(IY) = C*SY(IY) - S*SX(IX)
  22.         SX(IX) = STEMP
  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.         STEMP = C*SX(I) + S*SY(I)
  32.         SY(I) = C*SY(I) - S*SX(I)
  33.         SX(I) = STEMP
  34.    30 CONTINUE
  35.       RETURN
  36.       END
  37.