100 DOUBLE PRECISION ZERO, ONE
101 parameter( zero = 0.0d+0, one = 1.0d+0 )
105 DOUBLE PRECISION BIGNUM, CDEN, CDEN1, CNUM, CNUM1, MUL, SMLNUM
108 DOUBLE PRECISION DLAMCH
126 smlnum = dlamch(
'S' )
127 bignum = one / smlnum
136 cnum1 = cnum / bignum
137 IF( abs( cden1 ).GT.abs( cnum ) .AND. cnum.NE.zero )
THEN
144 ELSE IF( abs( cnum1 ).GT.abs( cden ) )
THEN
161 CALL zdscal( n, mul, sx, incx )
subroutine zdrscl(n, sa, sx, incx)
ZDRSCL multiplies a vector by the reciprocal of a real scalar.
subroutine zdscal(n, da, zx, incx)
ZDSCAL