103 DOUBLE PRECISION ZERO, ONE
104 parameter( zero = 0.0d+0, one = 1.0d+0 )
108 DOUBLE PRECISION BIGNUM, CDEN, CDEN1, CNUM, CNUM1, MUL, SMLNUM
111 DOUBLE PRECISION DLAMCH
129 smlnum = dlamch(
'S' )
130 bignum = one / smlnum
131 CALL
dlabad( smlnum, bignum )
140 cnum1 = cnum / bignum
141 IF( abs( cden1 ).GT.abs( cnum ) .AND. cnum.NE.zero )
THEN
148 ELSE IF( abs( cnum1 ).GT.abs( cden ) )
THEN
165 CALL
zdscal( n, mul, sx, incx )
subroutine zdscal(N, DA, ZX, INCX)
ZDSCAL
subroutine dlabad(SMALL, LARGE)
DLABAD
subroutine zdrscl(N, SA, SX, INCX)
ZDRSCL multiplies a vector by the reciprocal of a real scalar.