123 SUBROUTINE slahilb( N, NRHS, A, LDA, X, LDX, B, LDB, WORK, INFO)
130 INTEGER N, NRHS, LDA, LDX, LDB, INFO
132 REAL A(LDA, N), X(LDX, NRHS), B(LDB, NRHS), WORK(N)
146 INTEGER NMAX_EXACT, NMAX_APPROX
147 parameter(nmax_exact = 6, nmax_approx = 11)
161 IF (n .LT. 0 .OR. n .GT. nmax_approx)
THEN
163 ELSE IF (nrhs .LT. 0)
THEN
165 ELSE IF (lda .LT. n)
THEN
167 ELSE IF (ldx .LT. n)
THEN
169 ELSE IF (ldb .LT. n)
THEN
172 IF (info .LT. 0)
THEN
173 CALL xerbla(
'SLAHILB', -info)
176 IF (n .GT. nmax_exact)
THEN
198 a(i, j) = real(m) / (i + j - 1)
204 CALL slaset(
'Full', n, nrhs, 0.0, real(m), b, ldb)
211 work(j) = ( ( (work(j-1)/(j-1)) * (j-1 - n) ) /(j-1) )
217 x(i, j) = (work(i)*work(j)) / (i + j - 1)
subroutine slahilb(n, nrhs, a, lda, x, ldx, b, ldb, work, info)
SLAHILB
subroutine xerbla(srname, info)
subroutine slaset(uplo, m, n, alpha, beta, a, lda)
SLASET initializes the off-diagonal elements and the diagonal elements of a matrix to given values.