diff options
Diffstat (limited to 'TESTING/MATGEN/slahilb.f')
-rw-r--r-- | TESTING/MATGEN/slahilb.f | 26 |
1 files changed, 13 insertions, 13 deletions
diff --git a/TESTING/MATGEN/slahilb.f b/TESTING/MATGEN/slahilb.f index fb9c837c..ebcdeff4 100644 --- a/TESTING/MATGEN/slahilb.f +++ b/TESTING/MATGEN/slahilb.f @@ -2,20 +2,20 @@ * * =========== DOCUMENTATION =========== * -* Online html documentation available at -* http://www.netlib.org/lapack/explore-html/ +* Online html documentation available at +* http://www.netlib.org/lapack/explore-html/ * * Definition: * =========== * * SUBROUTINE SLAHILB( N, NRHS, A, LDA, X, LDX, B, LDB, WORK, INFO) -* +* * .. Scalar Arguments .. * INTEGER N, NRHS, LDA, LDX, LDB, INFO * .. Array Arguments .. * REAL A(LDA, N), X(LDX, NRHS), B(LDB, NRHS), WORK(N) * .. -* +* * *> \par Purpose: * ============= @@ -26,8 +26,8 @@ *> NRHS right-hand sides in B and solutions in X such that A*X=B. *> *> The Hilbert matrix is scaled by M = LCM(1, 2, ..., 2*N-1) so that all -*> entries are integers. The right-hand sides are the first NRHS -*> columns of M * the identity matrix, and the solutions are the +*> entries are integers. The right-hand sides are the first NRHS +*> columns of M * the identity matrix, and the solutions are the *> first NRHS columns of the inverse Hilbert matrix. *> *> The condition number of the Hilbert matrix grows exponentially with @@ -36,7 +36,7 @@ *> generated exactly without extra precision. Precision is exhausted *> when the largest entry in the inverse Hilbert matrix is greater than *> 2 to the power of the number of bits in the fraction of the data type -*> used plus one, which is 24 for single precision. +*> used plus one, which is 24 for single precision. *> *> In single, the generated solution is exact for N <= 6 and has *> small componentwise error for 7 <= N <= 11. @@ -50,7 +50,7 @@ *> N is INTEGER *> The dimension of the matrix A. *> \endverbatim -*> +*> *> \param[in] NRHS *> \verbatim *> NRHS is INTEGER @@ -112,10 +112,10 @@ * Authors: * ======== * -*> \author Univ. of Tennessee -*> \author Univ. of California Berkeley -*> \author Univ. of Colorado Denver -*> \author NAG Ltd. +*> \author Univ. of Tennessee +*> \author Univ. of California Berkeley +*> \author Univ. of Colorado Denver +*> \author NAG Ltd. * *> \date November 2015 * @@ -212,7 +212,7 @@ WORK(J) = ( ( (WORK(J-1)/(J-1)) * (J-1 - N) ) /(J-1) ) $ * (N +J -1) END DO - + DO J = 1, NRHS DO I = 1, N X(I, J) = (WORK(I)*WORK(J)) / (I + J - 1) |