136 SUBROUTINE zlatm1( MODE, COND, IRSIGN, IDIST, ISEED, D, N,
144 INTEGER IDIST, INFO, IRSIGN, MODE, N
145 DOUBLE PRECISION COND
156 parameter( one = 1.0d0 )
160 DOUBLE PRECISION ALPHA, TEMP
164 DOUBLE PRECISION DLARAN
166 EXTERNAL dlaran, zlarnd
172 INTRINSIC abs, dble, exp, log
187 IF( mode.LT.-6 .OR. mode.GT.6 )
THEN
189 ELSE IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
190 $ ( irsign.NE.0 .AND. irsign.NE.1 ) )
THEN
192 ELSE IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
195 ELSE IF( ( mode.EQ.6 .OR. mode.EQ.-6 ) .AND.
196 $ ( idist.LT.1 .OR. idist.GT.4 ) )
THEN
198 ELSE IF( n.LT.0 )
THEN
203 CALL xerbla(
'ZLATM1', -info )
210 GO TO ( 10, 30, 50, 70, 90, 110 )abs( mode )
235 alpha = cond**( -one / dble( n-1 ) )
237 d( i ) = alpha**( i-1 )
248 alpha = ( one-temp ) / dble( n-1 )
250 d( i ) = dble( n-i )*alpha + temp
258 alpha = log( one / cond )
260 d( i ) = exp( alpha*dlaran( iseed ) )
267 CALL zlarnv( idist, iseed, n, d )
274 IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
277 ctemp = zlarnd( 3, iseed )
278 d( i ) = d( i )*( ctemp / abs( ctemp ) )
subroutine zlarnv(idist, iseed, n, x)
ZLARNV returns a vector of random numbers from a uniform or normal distribution.
subroutine zlatm1(mode, cond, irsign, idist, iseed, d, n, info)
ZLATM1