diff options
| author | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2011-12-25 04:50:57 +0000 |
|---|---|---|
| committer | igor175 <igor175@8a072113-8704-0410-8d35-dd094bca7971> | 2011-12-25 04:50:57 +0000 |
| commit | acad7414deabcaa751530d3d35f91c613935b28d (patch) | |
| tree | 5e7b8f651dd0f7e92b116332eb666afe01ecdf5c /TESTING/LIN/schksy.f | |
| parent | ecfc5e7828f1d654da14886ebe5b5ef97a14faac (diff) | |
| download | lapack-acad7414deabcaa751530d3d35f91c613935b28d.tar.gz lapack-acad7414deabcaa751530d3d35f91c613935b28d.tar.bz2 lapack-acad7414deabcaa751530d3d35f91c613935b28d.zip | |
added a few comments in routines schksy.f dchksy.f cchksy.f zchksy.f
Diffstat (limited to 'TESTING/LIN/schksy.f')
| -rw-r--r-- | TESTING/LIN/schksy.f | 20 |
1 files changed, 17 insertions, 3 deletions
diff --git a/TESTING/LIN/schksy.f b/TESTING/LIN/schksy.f index c03b3a9c..a5f322aa 100644 --- a/TESTING/LIN/schksy.f +++ b/TESTING/LIN/schksy.f @@ -513,6 +513,9 @@ *+ TEST 3 (Using DSYTRS) * Solve and compute residual for A * X = B. * +* Choose a set of NRHS random solution vectors +* stored in XACT and set up the right hand side B +* SRNAMT = 'SLARHS' CALL SLARHS( PATH, XTYPE, UPLO, ' ', N, N, KL, KU, $ NRHS, A, LDA, XACT, LDA, B, LDA, @@ -523,7 +526,7 @@ CALL SSYTRS( UPLO, N, NRHS, AFAC, LDA, IWORK, X, $ LDA, INFO ) * -* Check error code from SSYTRS. +* Check error code from SSYTRS and handle error. * IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'SSYTRS', INFO, 0, UPLO, N, @@ -531,12 +534,18 @@ $ NERRS, NOUT ) * CALL SLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA ) +* +* Compute the residual for the solution +* CALL SPOT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK, $ LDA, RWORK, RESULT( 3 ) ) * *+ TEST 4 (Using DSYTRS2) * Solve and compute residual for A * X = B. * +* Choose a set of NRHS random solution vectors +* stored in XACT and set up the right hand side B +* SRNAMT = 'SLARHS' CALL SLARHS( PATH, XTYPE, UPLO, ' ', N, N, KL, KU, $ NRHS, A, LDA, XACT, LDA, B, LDA, @@ -547,7 +556,7 @@ CALL SSYTRS2( UPLO, N, NRHS, AFAC, LDA, IWORK, X, $ LDA, WORK, INFO ) * -* Check error code from SSYTRS2. +* Check error code from SSYTRS2 and handle error. * IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'SSYTRS2', INFO, 0, UPLO, N, @@ -555,6 +564,9 @@ $ NERRS, NOUT ) * CALL SLACPY( 'Full', N, NRHS, B, LDA, WORK, LDA ) +* +* Compute the residual for the solution +* CALL SPOT02( UPLO, N, NRHS, A, LDA, X, LDA, WORK, $ LDA, RWORK, RESULT( 4 ) ) * @@ -610,12 +622,14 @@ CALL SSYCON( UPLO, N, AFAC, LDA, IWORK, ANORM, RCOND, $ WORK, IWORK( N+1 ), INFO ) * -* Check error code from SSYCON. +* Check error code from SSYCON and handle error. * IF( INFO.NE.0 ) $ CALL ALAERH( PATH, 'SSYCON', INFO, 0, UPLO, N, N, $ -1, -1, -1, IMAT, NFAIL, NERRS, NOUT ) * +* Compute the test ratio to compare to values of RCOND +* RESULT( 9 ) = SGET06( RCOND, RCONDC ) * * Print information about the tests that did not pass |