125 SUBROUTINE slahilb(N, NRHS, A, LDA, X, LDX, B, LDB, WORK, INFO)
133 INTEGER N, NRHS, LDA, LDX, LDB, INFO
135 REAL A(lda, n), X(ldx, nrhs), B(ldb, nrhs), WORK(n)
149 INTEGER NMAX_EXACT, NMAX_APPROX
150 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)
XERBLA
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...