LAPACK  3.9.1
LAPACK: Linear Algebra PACKage

◆ zdrscl()

subroutine zdrscl ( integer  N,
double precision  SA,
complex*16, dimension( * )  SX,
integer  INCX 
)

ZDRSCL multiplies a vector by the reciprocal of a real scalar.

Download ZDRSCL + dependencies [TGZ] [ZIP] [TXT]

Purpose:
 ZDRSCL multiplies an n-element complex vector x by the real scalar
 1/a.  This is done without overflow or underflow as long as
 the final result x/a does not overflow or underflow.
Parameters
[in]N
          N is INTEGER
          The number of components of the vector x.
[in]SA
          SA is DOUBLE PRECISION
          The scalar a which is used to divide each component of x.
          SA must be >= 0, or the subroutine will divide by zero.
[in,out]SX
          SX is COMPLEX*16 array, dimension
                         (1+(N-1)*abs(INCX))
          The n-element vector x.
[in]INCX
          INCX is INTEGER
          The increment between successive values of the vector SX.
          > 0:  SX(1) = X(1) and SX(1+(i-1)*INCX) = x(i),     1< i<= n
Author
Univ. of Tennessee
Univ. of California Berkeley
Univ. of Colorado Denver
NAG Ltd.

Definition at line 83 of file zdrscl.f.

84 *
85 * -- LAPACK auxiliary routine --
86 * -- LAPACK is a software package provided by Univ. of Tennessee, --
87 * -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
88 *
89 * .. Scalar Arguments ..
90  INTEGER INCX, N
91  DOUBLE PRECISION SA
92 * ..
93 * .. Array Arguments ..
94  COMPLEX*16 SX( * )
95 * ..
96 *
97 * =====================================================================
98 *
99 * .. Parameters ..
100  DOUBLE PRECISION ZERO, ONE
101  parameter( zero = 0.0d+0, one = 1.0d+0 )
102 * ..
103 * .. Local Scalars ..
104  LOGICAL DONE
105  DOUBLE PRECISION BIGNUM, CDEN, CDEN1, CNUM, CNUM1, MUL, SMLNUM
106 * ..
107 * .. External Functions ..
108  DOUBLE PRECISION DLAMCH
109  EXTERNAL dlamch
110 * ..
111 * .. External Subroutines ..
112  EXTERNAL dlabad, zdscal
113 * ..
114 * .. Intrinsic Functions ..
115  INTRINSIC abs
116 * ..
117 * .. Executable Statements ..
118 *
119 * Quick return if possible
120 *
121  IF( n.LE.0 )
122  $ RETURN
123 *
124 * Get machine parameters
125 *
126  smlnum = dlamch( 'S' )
127  bignum = one / smlnum
128  CALL dlabad( smlnum, bignum )
129 *
130 * Initialize the denominator to SA and the numerator to 1.
131 *
132  cden = sa
133  cnum = one
134 *
135  10 CONTINUE
136  cden1 = cden*smlnum
137  cnum1 = cnum / bignum
138  IF( abs( cden1 ).GT.abs( cnum ) .AND. cnum.NE.zero ) THEN
139 *
140 * Pre-multiply X by SMLNUM if CDEN is large compared to CNUM.
141 *
142  mul = smlnum
143  done = .false.
144  cden = cden1
145  ELSE IF( abs( cnum1 ).GT.abs( cden ) ) THEN
146 *
147 * Pre-multiply X by BIGNUM if CDEN is small compared to CNUM.
148 *
149  mul = bignum
150  done = .false.
151  cnum = cnum1
152  ELSE
153 *
154 * Multiply X by CNUM / CDEN and return.
155 *
156  mul = cnum / cden
157  done = .true.
158  END IF
159 *
160 * Scale the vector X by MUL
161 *
162  CALL zdscal( n, mul, sx, incx )
163 *
164  IF( .NOT.done )
165  $ GO TO 10
166 *
167  RETURN
168 *
169 * End of ZDRSCL
170 *
double precision function dlamch(CMACH)
DLAMCH
Definition: dlamch.f:69
subroutine dlabad(SMALL, LARGE)
DLABAD
Definition: dlabad.f:74
subroutine zdscal(N, DA, ZX, INCX)
ZDSCAL
Definition: zdscal.f:78
Here is the call graph for this function:
Here is the caller graph for this function: