134 SUBROUTINE dlatm1( MODE, COND, IRSIGN, IDIST, ISEED, D, N, INFO )
141 INTEGER IDIST, INFO, IRSIGN, MODE, N
142 DOUBLE PRECISION COND
146 DOUBLE PRECISION D( * )
153 parameter( one = 1.0d0 )
154 DOUBLE PRECISION HALF
155 parameter( half = 0.5d0 )
159 DOUBLE PRECISION ALPHA, TEMP
162 DOUBLE PRECISION DLARAN
169 INTRINSIC abs, dble, exp, log
184 IF( mode.LT.-6 .OR. mode.GT.6 )
THEN
186 ELSE IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
187 $ ( irsign.NE.0 .AND. irsign.NE.1 ) )
THEN
189 ELSE IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
192 ELSE IF( ( mode.EQ.6 .OR. mode.EQ.-6 ) .AND.
193 $ ( idist.LT.1 .OR. idist.GT.3 ) )
THEN
195 ELSE IF( n.LT.0 )
THEN
200 CALL xerbla(
'DLATM1', -info )
207 GO TO ( 10, 30, 50, 70, 90, 110 )abs( mode )
232 alpha = cond**( -one / dble( n-1 ) )
234 d( i ) = alpha**( i-1 )
245 alpha = ( one-temp ) / dble( n-1 )
247 d( i ) = dble( n-i )*alpha + temp
255 alpha = log( one / cond )
257 d( i ) = exp( alpha*dlaran( iseed ) )
264 CALL dlarnv( idist, iseed, n, d )
271 IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
274 temp = dlaran( iseed )
subroutine xerbla(srname, info)
subroutine dlatm1(mode, cond, irsign, idist, iseed, d, n, info)
DLATM1
subroutine dlarnv(idist, iseed, n, x)
DLARNV returns a vector of random numbers from a uniform or normal distribution.