changed comments in LAPACK testing routines (s,d,c,z)chksy.f and (s,d,c,z)chksy_rook...

diff --git a/TESTING/LIN/dchksy.f b/TESTING/LIN/dchksy.f
index 3295f13579b36a860908dcbc054feeb30a57d096..894219f9ba4552babd706cdee1f06358855271c7 100644 (file)
--- a/TESTING/LIN/dchksy.f
+++ b/TESTING/LIN/dchksy.f
@@ -134,14 +134,12 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension
-*>                      (NMAX*max(3,NSMAX))
+*>          WORK is DOUBLE PRECISION array, dimension (NMAX*max(3,NSMAX))
 *> \endverbatim
 *>
 *> \param[out] RWORK
 *> \verbatim
-*>          RWORK is DOUBLE PRECISION array, dimension
-*>                      (max(NMAX,2*NSMAX))
+*>          RWORK is DOUBLE PRECISION array, dimension (max(NMAX,2*NSMAX))
 *> \endverbatim
 *>
 *> \param[out] IWORK
@@ -297,6 +295,7 @@
 *
 *              Begin generate the test matrix A.
 *
+*
 *              Set up parameters with DLATB4 for the matrix generator
 *              based on the type of matrix to be generated.
 *
@@ -315,6 +314,9 @@
                IF( INFO.NE.0 ) THEN
                   CALL ALAERH( PATH, 'DLATMS', INFO, 0, UPLO, N, N, -1,
      $                         -1, -1, IMAT, NFAIL, NERRS, NOUT )
+*
+*                    Skip all tests for this generated matrix
+*
                   GO TO 160
                END IF
 *
@@ -357,11 +359,11 @@
    50                   CONTINUE
                      END IF
                   ELSE
-                     IOFF = 0
                      IF( IUPLO.EQ.1 ) THEN
 *
 *                       Set the first IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 70 J = 1, N
                            I2 = MIN( J, IZERO )
                            DO 60 I = 1, I2
@@ -373,6 +375,7 @@
 *
 *                       Set the last IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 90 J = 1, N
                            I1 = MAX( J, IZERO )
                            DO 80 I = I1, N
@@ -506,6 +509,8 @@
                      RCONDC = ZERO
                      GO TO 140
                   END IF
+*
+*                 Do for each value of NRHS in NSVAL.
 *
                   DO 130 IRHS = 1, NNS
                      NRHS = NSVAL( IRHS )
@@ -586,7 +591,7 @@
      $                            RWORK( NRHS+1 ), WORK, IWORK( N+1 ),
      $                            INFO )
 *
-*                 Check error code from DSYRFS.
+*                    Check error code from DSYRFS and handle error.
 *
                      IF( INFO.NE.0 )
      $                  CALL ALAERH( PATH, 'DSYRFS', INFO, 0, UPLO, N,
@@ -612,6 +617,9 @@
                         END IF
   120                CONTINUE
                      NRUN = NRUN + 6
+*
+*                 End do for each value of NRHS in NSVAL.
+*
   130             CONTINUE
 *
 *+    TEST 9
@@ -629,7 +637,7 @@
      $               CALL ALAERH( PATH, 'DSYCON', INFO, 0, UPLO, N, N,
      $                            -1, -1, -1, IMAT, NFAIL, NERRS, NOUT )
 *
-*                 Compute the test ratio to compare to values of RCOND
+*                 Compute the test ratio to compare values of RCOND
 *
                   RESULT( 9 ) = DGET06( RCOND, RCONDC )
 *

diff --git a/TESTING/LIN/dchksy_rook.f b/TESTING/LIN/dchksy_rook.f
index 5be791ea2117b8ded063bbc1dde8eba80e73acde..c7cbf8fa188292319d449082c45f8c2b9efe5d63 100644 (file)
--- a/TESTING/LIN/dchksy_rook.f
+++ b/TESTING/LIN/dchksy_rook.f
@@ -135,14 +135,12 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension
-*>                      (NMAX*max(3,NSMAX))
+*>          WORK is DOUBLE PRECISION array, dimension (NMAX*max(3,NSMAX))
 *> \endverbatim
 *>
 *> \param[out] RWORK
 *> \verbatim
-*>          RWORK is DOUBLE PRECISION array, dimension
-*>                      (max(NMAX,2*NSMAX))
+*>          RWORK is DOUBLE PRECISION array, dimension (max(NMAX,2*NSMAX))
 *> \endverbatim
 *>
 *> \param[out] IWORK
@@ -375,11 +373,11 @@
    50                   CONTINUE
                      END IF
                   ELSE
-                     IOFF = 0
                      IF( IUPLO.EQ.1 ) THEN
 *
 *                       Set the first IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 70 J = 1, N
                            I2 = MIN( J, IZERO )
                            DO 60 I = 1, I2
@@ -391,6 +389,7 @@
 *
 *                       Set the last IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 90 J = 1, N
                            I1 = MAX( J, IZERO )
                            DO 80 I = I1, N
@@ -406,6 +405,7 @@
 *
 *              End generate the test matrix A.
 *
+*
 *              Do for each value of NB in NBVAL
 *
                DO 240 INB = 1, NNB
@@ -449,7 +449,7 @@
                      END IF
                   END IF
 *
-*                 Check error code from DSYTRF and handle error.
+*                 Check error code from DSYTRF_ROOK and handle error.
 *
                   IF( INFO.NE.K)
      $               CALL ALAERH( PATH, 'DSYTRF_ROOK', INFO, K,
@@ -634,7 +634,7 @@
 *
                         DTEMP = LAM_MAX / LAM_MIN
 *
-*                    DTEMP should be bounded by CONST
+*                       DTEMP should be bounded by CONST
 *
                         DTEMP = ABS( DTEMP ) - CONST + THRESH
                         IF( DTEMP.GT.RESULT( 4 ) )
@@ -675,7 +675,7 @@
 *
                         DTEMP = LAM_MAX / LAM_MIN
 *
-*                    DTEMP should be bounded by CONST
+*                       DTEMP should be bounded by CONST
 *
                         DTEMP = ABS( DTEMP ) - CONST + THRESH
                         IF( DTEMP.GT.RESULT( 4 ) )
@@ -702,7 +702,7 @@
                         NFAIL = NFAIL + 1
                      END IF
   200             CONTINUE
-                  NRUN = NRUN + NT
+                  NRUN = NRUN + 2
 *
 *                 Skip the other tests if this is not the first block
 *                 size.
@@ -716,13 +716,12 @@
                      RCONDC = ZERO
                      GO TO 230
                   END IF
+*
+*                 Do for each value of NRHS in NSVAL.
 *
                   DO 220 IRHS = 1, NNS
                      NRHS = NSVAL( IRHS )
 *
-*                    Begin loop over NRHS values
-*
-*
 *+    TEST 5 ( Using TRS_ROOK)
 *                 Solve and compute residual for  A * X = B.
 *
@@ -773,7 +772,7 @@
   210                CONTINUE
                      NRUN = NRUN + 2
 *
-*                    End loop over NRHS values
+*                 End do for each value of NRHS in NSVAL.
 *
   220             CONTINUE
 *

diff --git a/TESTING/LIN/schksy.f b/TESTING/LIN/schksy.f
index 5c698e5a4c6cc689b27826b0173d0c2cc52e410f..792586cdc1c1384730e01043d9e7a88db2b0034d 100644 (file)
--- a/TESTING/LIN/schksy.f
+++ b/TESTING/LIN/schksy.f
@@ -134,14 +134,12 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension
-*>                      (NMAX*max(3,NSMAX))
+*>          WORK is REAL array, dimension (NMAX*max(3,NSMAX))
 *> \endverbatim
 *>
 *> \param[out] RWORK
 *> \verbatim
-*>          RWORK is REAL array, dimension
-*>                      (max(NMAX,2*NSMAX))
+*>          RWORK is REAL array, dimension (max(NMAX,2*NSMAX))
 *> \endverbatim
 *>
 *> \param[out] IWORK
@@ -315,6 +313,9 @@
                IF( INFO.NE.0 ) THEN
                   CALL ALAERH( PATH, 'SLATMS', INFO, 0, UPLO, N, N, -1,
      $                         -1, -1, IMAT, NFAIL, NERRS, NOUT )
+*
+*                 Skip all tests for this generated matrix
+*
                   GO TO 160
                END IF
 *
@@ -357,11 +358,11 @@
    50                   CONTINUE
                      END IF
                   ELSE
-                     IOFF = 0
                      IF( IUPLO.EQ.1 ) THEN
 *
 *                       Set the first IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 70 J = 1, N
                            I2 = MIN( J, IZERO )
                            DO 60 I = 1, I2
@@ -373,6 +374,7 @@
 *
 *                       Set the last IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 90 J = 1, N
                            I1 = MAX( J, IZERO )
                            DO 80 I = I1, N
@@ -388,6 +390,7 @@
 *
 *              End generate the test matrix A.
 *
+*
 *              Do for each value of NB in NBVAL
 *
                DO 150 INB = 1, NNB
@@ -506,6 +509,8 @@
                      RCONDC = ZERO
                      GO TO 140
                   END IF
+*
+*                 Do for each value of NRHS in NSVAL.
 *
                   DO 130 IRHS = 1, NNS
                      NRHS = NSVAL( IRHS )
@@ -585,7 +590,7 @@
      $                            RWORK( NRHS+1 ), WORK, IWORK( N+1 ),
      $                            INFO )
 *
-*                 Check error code from SSYRFS.
+*                    Check error code from SSYRFS and handle error.
 *
                      IF( INFO.NE.0 )
      $                  CALL ALAERH( PATH, 'SSYRFS', INFO, 0, UPLO, N,
@@ -611,6 +616,9 @@
                         END IF
   120                CONTINUE
                      NRUN = NRUN + 6
+*
+*                 End do for each value of NRHS in NSVAL.
+*
   130             CONTINUE
 *
 *+    TEST 9

diff --git a/TESTING/LIN/schksy_rook.f b/TESTING/LIN/schksy_rook.f
index 27c457f0f953b6b641af5b5674c3dbd40f499695..8db53f73f0e03e70fa6c954c2c4cfb2a28b87f2a 100644 (file)
--- a/TESTING/LIN/schksy_rook.f
+++ b/TESTING/LIN/schksy_rook.f
@@ -135,14 +135,12 @@
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension
-*>                      (NMAX*max(3,NSMAX))
+*>          WORK is REAL array, dimension (NMAX*max(3,NSMAX))
 *> \endverbatim
 *>
 *> \param[out] RWORK
 *> \verbatim
-*>          RWORK is REAL array, dimension
-*>                      (max(NMAX,2*NSMAX))
+*>          RWORK is REAL array, dimension (max(NMAX,2*NSMAX))
 *> \endverbatim
 *>
 *> \param[out] IWORK
@@ -375,11 +373,11 @@
    50                   CONTINUE
                      END IF
                   ELSE
-                     IOFF = 0
                      IF( IUPLO.EQ.1 ) THEN
 *
 *                       Set the first IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 70 J = 1, N
                            I2 = MIN( J, IZERO )
                            DO 60 I = 1, I2
@@ -391,6 +389,7 @@
 *
 *                       Set the last IZERO rows and columns to zero.
 *
+                        IOFF = 0
                         DO 90 J = 1, N
                            I1 = MAX( J, IZERO )
                            DO 80 I = I1, N
@@ -406,6 +405,7 @@
 *
 *              End generate the test matrix A.
 *
+*
 *              Do for each value of NB in NBVAL
 *
                DO 240 INB = 1, NNB
@@ -449,7 +449,7 @@
                      END IF
                   END IF
 *
-*                 Check error code from SSYTRF and handle error.
+*                 Check error code from SSYTRF_ROOK and handle error.
 *
                   IF( INFO.NE.K)
      $               CALL ALAERH( PATH, 'SSYTRF_ROOK', INFO, K,
@@ -702,7 +702,7 @@
                         NFAIL = NFAIL + 1
                      END IF
   200             CONTINUE
-                  NRUN = NRUN + NT
+                  NRUN = NRUN + 2
 *
 *                 Skip the other tests if this is not the first block
 *                 size.
@@ -716,13 +716,12 @@
                      RCONDC = ZERO
                      GO TO 230
                   END IF
+*
+*                 Do for each value of NRHS in NSVAL.
 *
                   DO 220 IRHS = 1, NNS
                      NRHS = NSVAL( IRHS )
 *
-*                    Begin loop over NRHS values
-*
-*
 *+    TEST 5 ( Using TRS_ROOK)
 *                 Solve and compute residual for  A * X = B.
 *
@@ -773,7 +772,7 @@
   210                CONTINUE
                      NRUN = NRUN + 2
 *
-*                    End loop over NRHS values
+*                 End do for each value of NRHS in NSVAL.
 *
   220             CONTINUE
 *
@@ -793,7 +792,7 @@
      $                             UPLO, N, N, -1, -1, -1, IMAT,
      $                             NFAIL, NERRS, NOUT )
 *
-*                 Compute the test ratio to compare to values of RCOND
+*                 Compute the test ratio to compare values of RCOND
 *
                   RESULT( 7 ) = SGET06( RCOND, RCONDC )
 *

diff --git a/TESTING/LIN/zchkhe_rook.f b/TESTING/LIN/zchkhe_rook.f
index fece0e5f62bd24481e39eef1b3a2d84e16d6dc41..17c73d58ba09177ab3e4950d09322c69bb7503c5 100644 (file)
--- a/TESTING/LIN/zchkhe_rook.f
+++ b/TESTING/LIN/zchkhe_rook.f
@@ -199,7 +199,7 @@
       PARAMETER          ( ONEHALF = 0.5D+0 )
       DOUBLE PRECISION   EIGHT, SEVTEN
       PARAMETER          ( EIGHT = 8.0D+0, SEVTEN = 17.0D+0 )
-      COMPLEX            CZERO
+      COMPLEX*16         CZERO
       PARAMETER          ( CZERO = ( 0.0D+0, 0.0D+0 ) )
       INTEGER            NTYPES
       PARAMETER          ( NTYPES = 10 )
@@ -215,7 +215,7 @@
      $                   LWORK, MODE, N, NB, NERRS, NFAIL, NIMAT, NRHS,
      $                   NRUN, NT
       DOUBLE PRECISION   ALPHA, ANORM, CNDNUM, CONST, LAM_MAX, LAM_MIN,
-     $                   RCOND, RCONDC, STEMP
+     $                   RCOND, RCONDC, DTEMP
 *     ..
 *     .. Local Arrays ..
       CHARACTER          UPLOS( 2 )
@@ -523,7 +523,7 @@
 *                 Compute largest element in U or L
 *
                   RESULT( 3 ) = ZERO
-                  STEMP = ZERO
+                  DTEMP = ZERO
 *
                   CONST = ( ( ALPHA**2-ONE ) / ( ALPHA**2-ONEHALF ) ) /
      $                    ( ONE-ALPHA )
@@ -542,24 +542,24 @@
 *                       Get max absolute value from elements
 *                       in column k in in U
 *
-                        STEMP = ZLANGE( 'M', K-1, 1,
+                        DTEMP = ZLANGE( 'M', K-1, 1,
      $                          AFAC( ( K-1 )*LDA+1 ), LDA, RWORK )
                      ELSE
 *
 *                       Get max absolute value from elements
 *                       in columns k and k-1 in U
 *
-                        STEMP = ZLANGE( 'M', K-2, 2,
+                        DTEMP = ZLANGE( 'M', K-2, 2,
      $                          AFAC( ( K-2 )*LDA+1 ), LDA, RWORK )
                         K = K - 1
 *
                      END IF
 *
-*                    STEMP should be bounded by CONST
+*                    DTEMP should be bounded by CONST
 *
-                     STEMP = STEMP - CONST + THRESH
-                     IF( STEMP.GT.RESULT( 3 ) )
-     $                  RESULT( 3 ) = STEMP
+                     DTEMP = DTEMP - CONST + THRESH
+                     IF( DTEMP.GT.RESULT( 3 ) )
+     $                  RESULT( 3 ) = DTEMP
 *
                      K = K - 1
 *
@@ -580,24 +580,24 @@
 *                       Get max absolute value from elements
 *                       in column k in L
 *
-                        STEMP = ZLANGE( 'M', N-K, 1,
+                        DTEMP = ZLANGE( 'M', N-K, 1,
      $                          AFAC( ( K-1 )*LDA+K+1 ), LDA, RWORK )
                      ELSE
 *
 *                       Get max absolute value from elements
 *                       in columns k and k+1 in L
 *
-                        STEMP = ZLANGE( 'M', N-K-1, 2,
+                        DTEMP = ZLANGE( 'M', N-K-1, 2,
      $                          AFAC( ( K-1 )*LDA+K+2 ), LDA, RWORK )
                         K = K + 1
 *
                      END IF
 *
-*                    STEMP should be bounded by CONST
+*                    DTEMP should be bounded by CONST
 *
-                     STEMP = STEMP - CONST + THRESH
-                     IF( STEMP.GT.RESULT( 3 ) )
-     $                  RESULT( 3 ) = STEMP
+                     DTEMP = DTEMP - CONST + THRESH
+                     IF( DTEMP.GT.RESULT( 3 ) )
+     $                  RESULT( 3 ) = DTEMP
 *
                      K = K + 1
 *
@@ -610,7 +610,7 @@
 *                 Compute largest 2-Norm of 2-by-2 diag blocks
 *
                   RESULT( 4 ) = ZERO
-                  STEMP = ZERO
+                  DTEMP = ZERO
 *
                   CONST = ( ( ALPHA**2-ONE ) / ( ALPHA**2-ONEHALF ) )*
      $                    ( ( ONE + ALPHA ) / ( ONE - ALPHA ) )
@@ -636,7 +636,7 @@
 *
                         CALL ZHEEVX( 'N', 'N', 'N', 'N', 2, BLOCK,
      $                               2, WORK, CDUMMY, 1, CDUMMY, 1,
-     $                               ITEMP, ITEMP2, RWORK, STEMP,
+     $                               ITEMP, ITEMP2, RWORK, DTEMP,
      $                               RWORK( 3 ), RWORK( 5 ), WORK( 3 ),
      $                               4, RWORK( 7 ), INFO )
 *
@@ -645,13 +645,13 @@
                         LAM_MIN = MIN( ABS( WORK( 1 ) ),
      $                            ABS( WORK( 2 ) ) )
 *
-                        STEMP = LAM_MAX / LAM_MIN
+                        DTEMP = LAM_MAX / LAM_MIN
 *
-*                       STEMP should be bounded by CONST
+*                       DTEMP should be bounded by CONST
 *
-                        STEMP = ABS( STEMP ) - CONST + THRESH
-                        IF( STEMP.GT.RESULT( 4 ) )
-     $                     RESULT( 4 ) = STEMP
+                        DTEMP = ABS( DTEMP ) - CONST + THRESH
+                        IF( DTEMP.GT.RESULT( 4 ) )
+     $                     RESULT( 4 ) = DTEMP
                         K = K - 1
 *
                      END IF
@@ -682,7 +682,7 @@
 *
                         CALL ZHEEVX( 'N', 'N', 'N', 'N', 2, BLOCK,
      $                               2, WORK, CDUMMY, 1, CDUMMY, 1,
-     $                               ITEMP, ITEMP2, RWORK, STEMP,
+     $                               ITEMP, ITEMP2, RWORK, DTEMP,
      $                               RWORK( 3 ), RWORK( 5 ), WORK( 3 ),
      $                               4, RWORK( 7 ), INFO )
 *
@@ -691,13 +691,13 @@
                         LAM_MIN = MIN( ABS( WORK( 1 ) ),
      $                            ABS( WORK( 2 ) ) )
 *
-                        STEMP = LAM_MAX / LAM_MIN
+                        DTEMP = LAM_MAX / LAM_MIN
 *
-*                       STEMP should be bounded by CONST
+*                       DTEMP should be bounded by CONST
 *
-                        STEMP = ABS( STEMP ) - CONST + THRESH
-                        IF( STEMP.GT.RESULT( 4 ) )
-     $                     RESULT( 4 ) = STEMP
+                        DTEMP = ABS( DTEMP ) - CONST + THRESH
+                        IF( DTEMP.GT.RESULT( 4 ) )
+     $                     RESULT( 4 ) = DTEMP
                         K = K + 1
 *
                      END IF