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)
162 IF (n .LT. 0 .OR. n .GT. nmax_approx)
THEN
164 ELSE IF (nrhs .LT. 0)
THEN
166 ELSE IF (lda .LT. n)
THEN
168 ELSE IF (ldx .LT. n)
THEN
170 ELSE IF (ldb .LT. n)
THEN
173 IF (info .LT. 0)
THEN
174 CALL xerbla(
'SLAHILB', -info)
177 IF (n .GT. nmax_exact)
THEN
199 a(i, j) =
REAL(M) / (i + j - 1)
205 CALL slaset(
'Full', n, nrhs, 0.0,
REAL(M), B, LDB)
212 work(j) = ( ( (work(j-1)/(j-1)) * (j-1 - n) ) /(j-1) )
218 x(i, j) = (work(i)*work(j)) / (i + j - 1)
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...
subroutine slahilb(N, NRHS, A, LDA, X, LDX, B, LDB, WORK, INFO)
SLAHILB