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 / crotg.f < prev    next >
Encoding:
Text File  |  1994-08-02  |  764 b   |  31 lines
  1.       SUBROUTINE CROTG( CA, CB, C, S )
  2. *     .. Scalar Arguments ..
  3.       REAL              C
  4.       COMPLEX           CA, CB, S
  5. *     ..
  6. *     .. Local Scalars ..
  7.       REAL              NORM, SCALE
  8.       COMPLEX           ALPHA
  9. *     ..
  10. *     .. Intrinsic Functions ..
  11.       INTRINSIC         CABS, CONJG, SQRT
  12. *     ..
  13. *     .. Executable Statements ..
  14.       IF( CABS( CA ).NE.0. )
  15.      $   GO TO 10
  16.       C = 0.
  17.       S = ( 1., 0. )
  18.       CA = CB
  19.       GO TO 20
  20.    10 CONTINUE
  21.       SCALE = CABS( CA ) + CABS( CB )
  22.       NORM = SCALE*SQRT( ( CABS( CA / SCALE ) )**2+
  23.      $       ( CABS( CB / SCALE ) )**2 )
  24.       ALPHA = CA / CABS( CA )
  25.       C = CABS( CA ) / NORM
  26.       S = ALPHA*CONJG( CB ) / NORM
  27.       CA = ALPHA*NORM
  28.    20 CONTINUE
  29.       RETURN
  30.       END
  31.