code and comments fixed: subroutine name

Author: jip

TESTING/EIG/cchkhb2stg.f

@@ -1,4 +1,4 @@
-*> \brief \b CCHKHBSTG
+*> \brief \b CCHKHB2STG
 *
 *  =========== DOCUMENTATION ===========
 *
@@ -8,9 +8,9 @@
 *  Definition:
 *  ===========
 *
-*       SUBROUTINE CCHKHBSTG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
+*       SUBROUTINE CCHKHB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
 *                          ISEED, THRESH, NOUNIT, A, LDA, SD, SE, D1,
-*                          D2, D3, U, LDU, WORK, LWORK, RWORK RESULT, 
+*                          D2, D3, U, LDU, WORK, LWORK, RWORK RESULT,
 *                          INFO )
 *
 *       .. Scalar Arguments ..
@@ -31,18 +31,18 @@
 *>
 *> \verbatim
 *>
-*> CCHKHBSTG tests the reduction of a Hermitian band matrix to tridiagonal
+*> CCHKHB2STG tests the reduction of a Hermitian band matrix to tridiagonal
 *> from, used with the Hermitian eigenvalue problem.
 *>
 *> CHBTRD factors a Hermitian band matrix A as  U S U* , where * means
 *> conjugate transpose, S is symmetric tridiagonal, and U is unitary.
 *> CHBTRD can use either just the lower or just the upper triangle
-*> of A; CCHKHBSTG checks both cases.
+*> of A; CCHKHB2STG checks both cases.
 *>
 *> CHETRD_HB2ST factors a Hermitian band matrix A as  U S U* , 
 *> where * means conjugate transpose, S is symmetric tridiagonal, and U is
 *> unitary. CHETRD_HB2ST can use either just the lower or just
-*> the upper triangle of A; CCHKHBSTG checks both cases.
+*> the upper triangle of A; CCHKHB2STG checks both cases.
 *>
 *> DSTEQR factors S as  Z D1 Z'.  
 *> D1 is the matrix of eigenvalues computed when Z is not computed
@@ -52,7 +52,7 @@
 *> D3 is the matrix of eigenvalues computed when Z is not computed
 *> and from the S resulting of DSYTRD_SB2ST "L".
 *>
-*> When CCHKHBSTG is called, a number of matrix "sizes" ("n's"), a number
+*> When CCHKHB2STG is called, a number of matrix "sizes" ("n's"), a number
 *> of bandwidths ("k's"), and a number of matrix "types" are
 *> specified.  For each size ("n"), each bandwidth ("k") less than or
 *> equal to "n", and each type of matrix, one matrix will be generated
@@ -126,7 +126,7 @@
 *> \verbatim
 *>          NSIZES is INTEGER
 *>          The number of sizes of matrices to use.  If it is zero,
-*>          CCHKHBSTG does nothing.  It must be at least zero.
+*>          CCHKHB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] NN
@@ -141,7 +141,7 @@
 *> \verbatim
 *>          NWDTHS is INTEGER
 *>          The number of bandwidths to use.  If it is zero,
-*>          CCHKHBSTG does nothing.  It must be at least zero.
+*>          CCHKHB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] KK
@@ -154,7 +154,7 @@
 *> \param[in] NTYPES
 *> \verbatim
 *>          NTYPES is INTEGER
-*>          The number of elements in DOTYPE.   If it is zero, CCHKHBSTG
+*>          The number of elements in DOTYPE.   If it is zero, CCHKHB2STG
 *>          does nothing.  It must be at least zero.  If it is MAXTYP+1
 *>          and NSIZES is 1, then an additional type, MAXTYP+1 is
 *>          defined, which is to use whatever matrix is in A.  This
@@ -184,7 +184,7 @@
 *>          congruential sequence limited to small integers, and so
 *>          should produce machine independent random numbers. The
 *>          values of ISEED are changed on exit, and can be used in the
-*>          next call to CCHKHBSTG to continue the same random number
+*>          next call to CCHKHB2STG to continue the same random number
 *>          sequence.
 *> \endverbatim
 *>
@@ -432,7 +432,7 @@ SUBROUTINE CCHKHB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       END IF
 *
       IF( INFO.NE.0 ) THEN
-         CALL XERBLA( 'CCHKHBSTG', -INFO )
+         CALL XERBLA( 'CCHKHB2STG', -INFO )
          RETURN
       END IF
 *
@@ -837,7 +837,7 @@ SUBROUTINE CCHKHB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       CALL SLASUM( 'CHB', NOUNIT, NERRS, NTESTT )
       RETURN
 *
- 9999 FORMAT( ' CCHKHBSTG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
+ 9999 FORMAT( ' CCHKHB2STG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
      $      I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
  9998 FORMAT( / 1X, A3,
      $     ' -- Complex Hermitian Banded Tridiagonal Reduction Routines'

TESTING/EIG/dchksb2stg.f

@@ -1,4 +1,4 @@
-*> \brief \b DCHKSBSTG
+*> \brief \b DCHKSB2STG
 *
 *  =========== DOCUMENTATION ===========
 *
@@ -30,18 +30,18 @@
 *>
 *> \verbatim
 *>
-*> DCHKSBSTG tests the reduction of a symmetric band matrix to tridiagonal
+*> DCHKSB2STG tests the reduction of a symmetric band matrix to tridiagonal
 *> form, used with the symmetric eigenvalue problem.
 *>
 *> DSBTRD factors a symmetric band matrix A as  U S U' , where ' means
 *> transpose, S is symmetric tridiagonal, and U is orthogonal.
 *> DSBTRD can use either just the lower or just the upper triangle
-*> of A; DCHKSBSTG checks both cases.
+*> of A; DCHKSB2STG checks both cases.
 *>
 *> DSYTRD_SB2ST factors a symmetric band matrix A as  U S U' , 
 *> where ' means transpose, S is symmetric tridiagonal, and U is
 *> orthogonal. DSYTRD_SB2ST can use either just the lower or just
-*> the upper triangle of A; DCHKSBSTG checks both cases.
+*> the upper triangle of A; DCHKSB2STG checks both cases.
 *>
 *> DSTEQR factors S as  Z D1 Z'.  
 *> D1 is the matrix of eigenvalues computed when Z is not computed
@@ -51,7 +51,7 @@
 *> D3 is the matrix of eigenvalues computed when Z is not computed
 *> and from the S resulting of DSYTRD_SB2ST "L".
 *>
-*> When DCHKSBSTG is called, a number of matrix "sizes" ("n's"), a number
+*> When DCHKSB2STG is called, a number of matrix "sizes" ("n's"), a number
 *> of bandwidths ("k's"), and a number of matrix "types" are
 *> specified.  For each size ("n"), each bandwidth ("k") less than or
 *> equal to "n", and each type of matrix, one matrix will be generated
@@ -125,7 +125,7 @@
 *> \verbatim
 *>          NSIZES is INTEGER
 *>          The number of sizes of matrices to use.  If it is zero,
-*>          DCHKSBSTG does nothing.  It must be at least zero.
+*>          DCHKSB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] NN
@@ -140,7 +140,7 @@
 *> \verbatim
 *>          NWDTHS is INTEGER
 *>          The number of bandwidths to use.  If it is zero,
-*>          DCHKSBSTG does nothing.  It must be at least zero.
+*>          DCHKSB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] KK
@@ -153,7 +153,7 @@
 *> \param[in] NTYPES
 *> \verbatim
 *>          NTYPES is INTEGER
-*>          The number of elements in DOTYPE.   If it is zero, DCHKSBSTG
+*>          The number of elements in DOTYPE.   If it is zero, DCHKSB2STG
 *>          does nothing.  It must be at least zero.  If it is MAXTYP+1
 *>          and NSIZES is 1, then an additional type, MAXTYP+1 is
 *>          defined, which is to use whatever matrix is in A.  This
@@ -183,7 +183,7 @@
 *>          congruential sequence limited to small integers, and so
 *>          should produce machine independent random numbers. The
 *>          values of ISEED are changed on exit, and can be used in the
-*>          next call to DCHKSBSTG to continue the same random number
+*>          next call to DCHKSB2STG to continue the same random number
 *>          sequence.
 *> \endverbatim
 *>
@@ -422,7 +422,7 @@ SUBROUTINE DCHKSB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       END IF
 *
       IF( INFO.NE.0 ) THEN
-         CALL XERBLA( 'DCHKSBSTG', -INFO )
+         CALL XERBLA( 'DCHKSB2STG', -INFO )
          RETURN
       END IF
 *
@@ -827,12 +827,12 @@ SUBROUTINE DCHKSB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       CALL DLASUM( 'DSB', NOUNIT, NERRS, NTESTT )
       RETURN
 *
- 9999 FORMAT( ' DCHKSBSTG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
+ 9999 FORMAT( ' DCHKSB2STG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
      $      I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
 *
  9998 FORMAT( / 1X, A3,
      $      ' -- Real Symmetric Banded Tridiagonal Reduction Routines' )
- 9997 FORMAT( ' Matrix types (see DCHKSBSTG for details): ' )
+ 9997 FORMAT( ' Matrix types (see DCHKSB2STG for details): ' )
 *
  9996 FORMAT( / ' Special Matrices:',
      $      / '  1=Zero matrix.                        ',

TESTING/EIG/schksb2stg.f

@@ -1,4 +1,4 @@
-*> \brief \b SCHKSBSTG
+*> \brief \b SCHKSB2STG
 *
 *  =========== DOCUMENTATION ===========
 *
@@ -30,18 +30,18 @@
 *>
 *> \verbatim
 *>
-*> SCHKSBSTG tests the reduction of a symmetric band matrix to tridiagonal
+*> SCHKSB2STG tests the reduction of a symmetric band matrix to tridiagonal
 *> form, used with the symmetric eigenvalue problem.
 *>
 *> SSBTRD factors a symmetric band matrix A as  U S U' , where ' means
 *> transpose, S is symmetric tridiagonal, and U is orthogonal.
 *> SSBTRD can use either just the lower or just the upper triangle
-*> of A; SCHKSBSTG checks both cases.
+*> of A; SCHKSB2STG checks both cases.
 *>
 *> SSYTRD_SB2ST factors a symmetric band matrix A as  U S U' , 
 *> where ' means transpose, S is symmetric tridiagonal, and U is
 *> orthogonal. SSYTRD_SB2ST can use either just the lower or just
-*> the upper triangle of A; SCHKSBSTG checks both cases.
+*> the upper triangle of A; SCHKSB2STG checks both cases.
 *>
 *> SSTEQR factors S as  Z D1 Z'.  
 *> D1 is the matrix of eigenvalues computed when Z is not computed
@@ -51,7 +51,7 @@
 *> D3 is the matrix of eigenvalues computed when Z is not computed
 *> and from the S resulting of SSYTRD_SB2ST "L".
 *>
-*> When SCHKSBSTG is called, a number of matrix "sizes" ("n's"), a number
+*> When SCHKSB2STG is called, a number of matrix "sizes" ("n's"), a number
 *> of bandwidths ("k's"), and a number of matrix "types" are
 *> specified.  For each size ("n"), each bandwidth ("k") less than or
 *> equal to "n", and each type of matrix, one matrix will be generated
@@ -125,7 +125,7 @@
 *> \verbatim
 *>          NSIZES is INTEGER
 *>          The number of sizes of matrices to use.  If it is zero,
-*>          SCHKSBSTG does nothing.  It must be at least zero.
+*>          SCHKSB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] NN
@@ -140,7 +140,7 @@
 *> \verbatim
 *>          NWDTHS is INTEGER
 *>          The number of bandwidths to use.  If it is zero,
-*>          SCHKSBSTG does nothing.  It must be at least zero.
+*>          SCHKSB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] KK
@@ -153,7 +153,7 @@
 *> \param[in] NTYPES
 *> \verbatim
 *>          NTYPES is INTEGER
-*>          The number of elements in DOTYPE.   If it is zero, SCHKSBSTG
+*>          The number of elements in DOTYPE.   If it is zero, SCHKSB2STG
 *>          does nothing.  It must be at least zero.  If it is MAXTYP+1
 *>          and NSIZES is 1, then an additional type, MAXTYP+1 is
 *>          defined, which is to use whatever matrix is in A.  This
@@ -183,7 +183,7 @@
 *>          congruential sequence limited to small integers, and so
 *>          should produce machine independent random numbers. The
 *>          values of ISEED are changed on exit, and can be used in the
-*>          next call to SCHKSBSTG to continue the same random number
+*>          next call to SCHKSB2STG to continue the same random number
 *>          sequence.
 *> \endverbatim
 *>
@@ -422,7 +422,7 @@ SUBROUTINE SCHKSB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       END IF
 *
       IF( INFO.NE.0 ) THEN
-         CALL XERBLA( 'SCHKSBSTG', -INFO )
+         CALL XERBLA( 'SCHKSB2STG', -INFO )
          RETURN
       END IF
 *
@@ -827,12 +827,12 @@ SUBROUTINE SCHKSB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       CALL SLASUM( 'SSB', NOUNIT, NERRS, NTESTT )
       RETURN
 *
- 9999 FORMAT( ' SCHKSBSTG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
+ 9999 FORMAT( ' SCHKSB2STG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
      $      I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
 *
  9998 FORMAT( / 1X, A3,
      $      ' -- Real Symmetric Banded Tridiagonal Reduction Routines' )
- 9997 FORMAT( ' Matrix types (see SCHKSBSTG for details): ' )
+ 9997 FORMAT( ' Matrix types (see SCHKSB2STG for details): ' )
 *
  9996 FORMAT( / ' Special Matrices:',
      $      / '  1=Zero matrix.                        ',

TESTING/EIG/zchkhb2stg.f

@@ -1,4 +1,4 @@
-*> \brief \b ZCHKHBSTG
+*> \brief \b ZCHKHB2STG
 *
 *  =========== DOCUMENTATION ===========
 *
@@ -8,7 +8,7 @@
 *  Definition:
 *  ===========
 *
-*       SUBROUTINE ZCHKHBSTG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
+*       SUBROUTINE ZCHKHB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
 *                          ISEED, THRESH, NOUNIT, A, LDA, SD, SE, D1,
 *                          D2, D3, U, LDU, WORK, LWORK, RWORK RESULT, 
 *                          INFO )
@@ -31,18 +31,18 @@
 *>
 *> \verbatim
 *>
-*> ZCHKHBSTG tests the reduction of a Hermitian band matrix to tridiagonal
+*> ZCHKHB2STG tests the reduction of a Hermitian band matrix to tridiagonal
 *> from, used with the Hermitian eigenvalue problem.
 *>
 *> ZHBTRD factors a Hermitian band matrix A as  U S U* , where * means
 *> conjugate transpose, S is symmetric tridiagonal, and U is unitary.
 *> ZHBTRD can use either just the lower or just the upper triangle
-*> of A; ZCHKHBSTG checks both cases.
+*> of A; ZCHKHB2STG checks both cases.
 *>
 *> ZHETRD_HB2ST factors a Hermitian band matrix A as  U S U* , 
 *> where * means conjugate transpose, S is symmetric tridiagonal, and U is
 *> unitary. ZHETRD_HB2ST can use either just the lower or just
-*> the upper triangle of A; ZCHKHBSTG checks both cases.
+*> the upper triangle of A; ZCHKHB2STG checks both cases.
 *>
 *> DSTEQR factors S as  Z D1 Z'.  
 *> D1 is the matrix of eigenvalues computed when Z is not computed
@@ -52,7 +52,7 @@
 *> D3 is the matrix of eigenvalues computed when Z is not computed
 *> and from the S resulting of DSYTRD_SB2ST "L".
 *>
-*> When ZCHKHBSTG is called, a number of matrix "sizes" ("n's"), a number
+*> When ZCHKHB2STG is called, a number of matrix "sizes" ("n's"), a number
 *> of bandwidths ("k's"), and a number of matrix "types" are
 *> specified.  For each size ("n"), each bandwidth ("k") less than or
 *> equal to "n", and each type of matrix, one matrix will be generated
@@ -126,7 +126,7 @@
 *> \verbatim
 *>          NSIZES is INTEGER
 *>          The number of sizes of matrices to use.  If it is zero,
-*>          ZCHKHBSTG does nothing.  It must be at least zero.
+*>          ZCHKHB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] NN
@@ -141,7 +141,7 @@
 *> \verbatim
 *>          NWDTHS is INTEGER
 *>          The number of bandwidths to use.  If it is zero,
-*>          ZCHKHBSTG does nothing.  It must be at least zero.
+*>          ZCHKHB2STG does nothing.  It must be at least zero.
 *> \endverbatim
 *>
 *> \param[in] KK
@@ -154,7 +154,7 @@
 *> \param[in] NTYPES
 *> \verbatim
 *>          NTYPES is INTEGER
-*>          The number of elements in DOTYPE.   If it is zero, ZCHKHBSTG
+*>          The number of elements in DOTYPE.   If it is zero, ZCHKHB2STG
 *>          does nothing.  It must be at least zero.  If it is MAXTYP+1
 *>          and NSIZES is 1, then an additional type, MAXTYP+1 is
 *>          defined, which is to use whatever matrix is in A.  This
@@ -184,7 +184,7 @@
 *>          congruential sequence limited to small integers, and so
 *>          should produce machine independent random numbers. The
 *>          values of ISEED are changed on exit, and can be used in the
-*>          next call to ZCHKHBSTG to continue the same random number
+*>          next call to ZCHKHB2STG to continue the same random number
 *>          sequence.
 *> \endverbatim
 *>
@@ -432,7 +432,7 @@ SUBROUTINE ZCHKHB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       END IF
 *
       IF( INFO.NE.0 ) THEN
-         CALL XERBLA( 'ZCHKHBSTG', -INFO )
+         CALL XERBLA( 'ZCHKHB2STG', -INFO )
          RETURN
       END IF
 *
@@ -837,7 +837,7 @@ SUBROUTINE ZCHKHB2STG( NSIZES, NN, NWDTHS, KK, NTYPES, DOTYPE,
       CALL DLASUM( 'ZHB', NOUNIT, NERRS, NTESTT )
       RETURN
 *
- 9999 FORMAT( ' ZCHKHBSTG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
+ 9999 FORMAT( ' ZCHKHB2STG: ', A, ' returned INFO=', I6, '.', / 9X, 'N=',
      $      I6, ', JTYPE=', I6, ', ISEED=(', 3( I5, ',' ), I5, ')' )
  9998 FORMAT( / 1X, A3,
      $     ' -- Complex Hermitian Banded Tridiagonal Reduction Routines'