home *** CD-ROM | disk | FTP | other *** search
/ SGI Developer Toolbox 6.1 / SGI Developer Toolbox 6.1 - Disc 4.iso / lib / mathlib / libblas / src_original / drot.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  1.1 KB  |  51 lines
  1.       SUBROUTINE DROT( N, DX, INCX, DY, INCY, C, S )
  2. *
  3. *     applies a plane rotation.
  4. *     jack dongarra, linpack, 3/11/78.
  5. *
  6. *     .. Scalar Arguments ..
  7.       INTEGER          INCX, INCY, N
  8.       DOUBLE PRECISION C, S
  9. *     ..
  10. *     .. Array Arguments ..
  11.       DOUBLE PRECISION DX( 1 ), DY( 1 )
  12. *     ..
  13. *     .. Local Scalars ..
  14.       INTEGER          I, IX, IY
  15.       DOUBLE PRECISION DTEMP
  16. *     ..
  17. *     .. Executable Statements ..
  18. *
  19.       IF( N.LE.0 )
  20.      $   RETURN
  21.       IF( INCX.EQ.1 .AND. INCY.EQ.1 )
  22.      $   GO TO 20
  23. *
  24. *        code for unequal increments or equal increments
  25. *          not equal to 1
  26. *
  27.       IX = 1
  28.       IY = 1
  29.       IF( INCX.LT.0 )
  30.      $   IX = ( -N+1 )*INCX + 1
  31.       IF( INCY.LT.0 )
  32.      $   IY = ( -N+1 )*INCY + 1
  33.       DO 10 I = 1, N
  34.          DTEMP = C*DX( IX ) + S*DY( IY )
  35.          DY( IY ) = C*DY( IY ) - S*DX( IX )
  36.          DX( IX ) = DTEMP
  37.          IX = IX + INCX
  38.          IY = IY + INCY
  39.    10 CONTINUE
  40.       RETURN
  41. *
  42. *        code for both increments equal to 1
  43. *
  44.    20 DO 30 I = 1, N
  45.          DTEMP = C*DX( I ) + S*DY( I )
  46.          DY( I ) = C*DY( I ) - S*DX( I )
  47.          DX( I ) = DTEMP
  48.    30 CONTINUE
  49.       RETURN
  50.       END
  51.