122 SUBROUTINE slatm7( MODE, COND, IRSIGN, IDIST, ISEED, D, N,
132 INTEGER IDIST, INFO, IRSIGN, MODE, N, RANK
143 parameter ( one = 1.0e0 )
145 parameter ( zero = 0.0e0 )
147 parameter ( half = 0.5e0 )
161 INTRINSIC abs, exp, log, real
176 IF( mode.LT.-6 .OR. mode.GT.6 )
THEN
178 ELSE IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
179 $ ( irsign.NE.0 .AND. irsign.NE.1 ) )
THEN
181 ELSE IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
184 ELSE IF( ( mode.EQ.6 .OR. mode.EQ.-6 ) .AND.
185 $ ( idist.LT.1 .OR. idist.GT.3 ) )
THEN
187 ELSE IF( n.LT.0 )
THEN
192 CALL xerbla(
'SLATM7', -info )
199 GO TO ( 100, 130, 160, 190, 210, 230 )abs( mode )
207 DO 120 i = rank + 1, n
216 DO 140 i = 1, rank - 1
219 DO 150 i = rank + 1, n
222 d( rank ) = one / cond
229 IF( n.GT.1 .AND. rank.GT.1 )
THEN
230 alpha = cond**( -one /
REAL( RANK-1 ) )
232 d( i ) = alpha**( i-1 )
234 DO 180 i = rank + 1, n
246 alpha = ( one-temp ) /
REAL( n-1 )
248 d( i ) =
REAL( n-i )*ALPHA + TEMP
256 alpha = log( one / cond )
258 d( i ) = exp( alpha*slaran( iseed ) )
265 CALL slarnv( idist, iseed, n, d )
272 IF( ( mode.NE.-6 .AND. mode.NE.0 .AND. mode.NE.6 ) .AND.
275 temp = slaran( iseed )
subroutine xerbla(SRNAME, INFO)
XERBLA
subroutine slarnv(IDIST, ISEED, N, X)
SLARNV returns a vector of random numbers from a uniform or normal distribution.
subroutine slatm7(MODE, COND, IRSIGN, IDIST, ISEED, D, N, RANK, INFO)
SLATM7