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
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.
double precision function dlamch(CMACH)
DLAMCH