diff options
Diffstat (limited to 'TESTING/MATGEN/dlahilb.f')
-rw-r--r-- | TESTING/MATGEN/dlahilb.f | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/TESTING/MATGEN/dlahilb.f b/TESTING/MATGEN/dlahilb.f index 7b2badab..7e30b3bc 100644 --- a/TESTING/MATGEN/dlahilb.f +++ b/TESTING/MATGEN/dlahilb.f @@ -1,4 +1,4 @@ -C> \brief \b DLAHILB +*> \brief \b DLAHILB * * =========== DOCUMENTATION =========== * @@ -140,7 +140,7 @@ C> \brief \b DLAHILB INTEGER TM, TI, R INTEGER M INTEGER I, J - +* .. * .. Parameters .. * NMAX_EXACT the largest dimension where the generated data is * exact. @@ -177,7 +177,7 @@ C> \brief \b DLAHILB IF (N .GT. NMAX_EXACT) THEN INFO = 1 END IF - +* * Compute M = the LCM of the integers [1, 2*N-1]. The largest * reasonable N is small enough that integers suffice (up to N = 11). M = 1 @@ -192,14 +192,14 @@ C> \brief \b DLAHILB END DO M = (M / TI) * I END DO - +* * Generate the scaled Hilbert matrix in A DO J = 1, N DO I = 1, N A(I, J) = DBLE(M) / (I + J - 1) END DO END DO - +* * Generate matrix B as simply the first NRHS columns of M * the * identity. CALL DLASET('Full', N, NRHS, 0.0D+0, DBLE(M), B, LDB) @@ -212,12 +212,12 @@ C> \brief \b DLAHILB 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) END DO END DO - +* END |