dlahd2.f

Go to the documentation of this file.

*> \brief \b DLAHD2
*
*  =========== DOCUMENTATION ===========
*
* Online html documentation available at
*            http://www.netlib.org/lapack/explore-html/
*
*  Definition:
*  ===========
*
*       SUBROUTINE DLAHD2( IOUNIT, PATH )
*
*       .. Scalar Arguments ..
*       CHARACTER*3        PATH
*       INTEGER            IOUNIT
*       ..
*
*
*> \par Purpose:
*  =============
*>
*> \verbatim
*>
*> DLAHD2 prints header information for the different test paths.
*> \endverbatim
*
*  Arguments:
*  ==========
*
*> \param[in] IOUNIT
*> \verbatim
*>          IOUNIT is INTEGER.
*>          On entry, IOUNIT specifies the unit number to which the
*>          header information should be printed.
*> \endverbatim
*>
*> \param[in] PATH
*> \verbatim
*>          PATH is CHARACTER*3.
*>          On entry, PATH contains the name of the path for which the
*>          header information is to be printed.  Current paths are
*>
*>             DHS, ZHS:  Non-symmetric eigenproblem.
*>             DST, ZST:  Symmetric eigenproblem.
*>             DSG, ZSG:  Symmetric Generalized eigenproblem.
*>             DBD, ZBD:  Singular Value Decomposition (SVD)
*>             DBB, ZBB:  General Banded reduction to bidiagonal form
*>
*>          These paths also are supplied in double precision (replace
*>          leading S by D and leading C by Z in path names).
*> \endverbatim
*
*  Authors:
*  ========
*
*> \author Univ. of Tennessee
*> \author Univ. of California Berkeley
*> \author Univ. of Colorado Denver
*> \author NAG Ltd.
*
*> \ingroup double_eig
*
*  =====================================================================
      SUBROUTINE dlahd2( IOUNIT, PATH )
*
*  -- LAPACK test routine --
*  -- LAPACK is a software package provided by Univ. of Tennessee,    --
*  -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..--
*
*     .. Scalar Arguments ..
      CHARACTER*3        PATH
      INTEGER            IOUNIT
*     ..
*
*  =====================================================================
*
*     .. Local Scalars ..
      LOGICAL            CORZ, SORD
      CHARACTER*2        C2
      INTEGER            J
*     ..
*     .. External Functions ..
      LOGICAL            LSAME, LSAMEN
      EXTERNAL           lsame, lsamen
*     ..
*     .. Executable Statements ..
*
      IF( iounit.LE.0 )
     $   RETURN
      sord = lsame( path, 'S' ) .OR. lsame( path, 'D' )
      corz = lsame( path, 'C' ) .OR. lsame( path, 'Z' )
      IF( .NOT.sord .AND. .NOT.corz ) THEN
         WRITE( iounit, fmt = 9999 )path
      END IF
      c2 = path( 2: 3 )
*
      IF( lsamen( 2, c2, 'HS' ) ) THEN
         IF( sord ) THEN
*
*           Real Non-symmetric Eigenvalue Problem:
*
            WRITE( iounit, fmt = 9998 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9988 )
            WRITE( iounit, fmt = 9987 )
            WRITE( iounit, fmt = 9986 )'pairs ', 'pairs ', 'prs.',
     $         'prs.'
            WRITE( iounit, fmt = 9985 )
*
*           Tests performed
*
            WRITE( iounit, fmt = 9984 )'orthogonal', '''=transpose',
     $         ( '''', j = 1, 6 )
*
         ELSE
*
*           Complex Non-symmetric Eigenvalue Problem:
*
            WRITE( iounit, fmt = 9997 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9988 )
            WRITE( iounit, fmt = 9987 )
            WRITE( iounit, fmt = 9986 )'e.vals', 'e.vals', 'e.vs',
     $         'e.vs'
            WRITE( iounit, fmt = 9985 )
*
*           Tests performed
*
            WRITE( iounit, fmt = 9984 )'unitary', '*=conj.transp.',
     $         ( '*', j = 1, 6 )
         END IF
*
      ELSE IF( lsamen( 2, c2, 'ST' ) ) THEN
*
         IF( sord ) THEN
*
*           Real Symmetric Eigenvalue Problem:
*
            WRITE( iounit, fmt = 9996 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9983 )
            WRITE( iounit, fmt = 9982 )
            WRITE( iounit, fmt = 9981 )'Symmetric'
*
*           Tests performed
*
            WRITE( iounit, fmt = 9968 )
*
         ELSE
*
*           Complex Hermitian Eigenvalue Problem:
*
            WRITE( iounit, fmt = 9995 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9983 )
            WRITE( iounit, fmt = 9982 )
            WRITE( iounit, fmt = 9981 )'Hermitian'
*
*           Tests performed
*
            WRITE( iounit, fmt = 9967 )
         END IF
*
      ELSE IF( lsamen( 2, c2, 'SG' ) ) THEN
*
         IF( sord ) THEN
*
*           Real Symmetric Generalized Eigenvalue Problem:
*
            WRITE( iounit, fmt = 9992 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9980 )
            WRITE( iounit, fmt = 9979 )
            WRITE( iounit, fmt = 9978 )'Symmetric'
*
*           Tests performed
*
            WRITE( iounit, fmt = 9977 )
            WRITE( iounit, fmt = 9976 )
*
         ELSE
*
*           Complex Hermitian Generalized Eigenvalue Problem:
*
            WRITE( iounit, fmt = 9991 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9980 )
            WRITE( iounit, fmt = 9979 )
            WRITE( iounit, fmt = 9978 )'Hermitian'
*
*           Tests performed
*
            WRITE( iounit, fmt = 9975 )
            WRITE( iounit, fmt = 9974 )
*
         END IF
*
      ELSE IF( lsamen( 2, c2, 'BD' ) ) THEN
*
         IF( sord ) THEN
*
*           Real Singular Value Decomposition:
*
            WRITE( iounit, fmt = 9994 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9973 )
*
*           Tests performed
*
            WRITE( iounit, fmt = 9972 )'orthogonal'
            WRITE( iounit, fmt = 9971 )
         ELSE
*
*           Complex Singular Value Decomposition:
*
            WRITE( iounit, fmt = 9993 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9973 )
*
*           Tests performed
*
            WRITE( iounit, fmt = 9972 )'unitary   '
            WRITE( iounit, fmt = 9971 )
         END IF
*
      ELSE IF( lsamen( 2, c2, 'BB' ) ) THEN
*
         IF( sord ) THEN
*
*           Real General Band reduction to bidiagonal form:
*
            WRITE( iounit, fmt = 9990 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9970 )
*
*           Tests performed
*
            WRITE( iounit, fmt = 9969 )'orthogonal'
         ELSE
*
*           Complex Band reduction to bidiagonal form:
*
            WRITE( iounit, fmt = 9989 )path
*
*           Matrix types
*
            WRITE( iounit, fmt = 9970 )
*
*           Tests performed
*
            WRITE( iounit, fmt = 9969 )'unitary   '
         END IF
*
      ELSE
*
         WRITE( iounit, fmt = 9999 )path
         RETURN
      END IF
*
      RETURN
*
 9999 FORMAT( 1x, a3, ':  no header available' )
 9998 FORMAT( / 1x, a3, ' -- Real Non-symmetric eigenvalue problem' )
 9997 FORMAT( / 1x, a3, ' -- Complex Non-symmetric eigenvalue problem' )
 9996 FORMAT( / 1x, a3, ' -- Real Symmetric eigenvalue problem' )
 9995 FORMAT( / 1x, a3, ' -- Complex Hermitian eigenvalue problem' )
 9994 FORMAT( / 1x, a3, ' -- Real Singular Value Decomposition' )
 9993 FORMAT( / 1x, a3, ' -- Complex Singular Value Decomposition' )
 9992 FORMAT( / 1x, a3, ' -- real symmetric generalized eigenvalue ',
     $      'problem' )
 9991 FORMAT( / 1X, A3, ' -- Complex Hermitian Generalized eigenvalue ',
     $      'problem' )
 9990 FORMAT( / 1X, A3, ' -- Real Band reduc. to bidiagonal form' )
 9989 FORMAT( / 1X, A3, ' -- Complex Band reduc. to bidiagonal form' )
*
 9988 FORMAT( ' Matrix types (see xCHKHS for details): ' )
*
 9987 FORMAT( / ' Special Matrices:', / '  1=zero matrix.             ',
     $      '           ', '  5=diagonal: geometr. spaced entries.',
     $      / '  2=identity matrix.                    ', '  6=diagona',
     $      'l: clustered entries.', / '  3=transposed jordan block.  ',
     $      '          ', '  7=diagonal: large, evenly spaced.', / '  ',
     $      '4=diagonal: evenly spaced entries.    ', '  8=diagonal: s',
     $      'mall, evenly spaced.' )
 9986 FORMAT( ' dense, non-symmetric matrices:', / '  9=well-cond., ev',
     $      'enly spaced eigenvals.', ' 14=ill-cond., geomet. spaced e',
     $      'igenals.', / ' 10=well-cond., geom. spaced eigenvals. ',
     $      ' 15=ill-conditioned, clustered e.vals.', / ' 11=well-cond',
     $      'itioned, clustered e.vals. ', ' 16=ill-cond., random comp',
     $      'lex ', A6, / ' 12=well-cond., random complex ', A6, '   ',
     $      ' 17=ill-cond., large rand. complx ', A4, / ' 13=ill-condi',
     $      'tioned, evenly spaced.     ', ' 18=ill-cond., small rand.',
     $      ' complx ', A4 )
 9985 FORMAT( ' 19=matrix with random o(1) entries.    ', ' 21=matrix ',
     $      'with small random entries.', / ' 20=matrix with large ran',
     $      'dom entries.   ' )
 9984 FORMAT( / ' tests performed:   ', '(h is hessenberg, t is schur,',
     $      ' u and z are ', A, ',', / 20X, A, ', w is a diagonal matr',
     $      'ix of eigenvalues,', / 20X, 'l and r are the left and rig',
     $      'ht eigenvector matrices)', / '  1 = | a - u h u', A1, ' |',
     $      ' / ( |a| n ulp )         ', '  2 = | i - u u', A1, ' | / ',
     $      '( n ulp )', / '  3 = | h - z t z', A1, ' | / ( |h| n ulp ',
     $      ')         ', '  4 = | i - z z', A1, ' | / ( n ulp )',
     $      / '  5 = | a - uz t(uz)', A1, ' | / ( |a| n ulp )     ',
     $      '  6 = | i - uz(uz)', A1, ' | / ( n ulp )', / '  7 = | t(',
     $      'e.vects.) - t(no e.vects.) | / ( |t| ulp )', / '  8 = | w',
     $      '(e.vects.) - w(no e.vects.) | / ( |w| ulp )', / '  9 = | ',
     $      'tr - rw | / ( |t| |r| ulp )     ', ' 10 = | lt - wl | / (',
     $      ' |t| |l| ulp )', / ' 11= |hx - xw| / (|h| |x| ulp)  (inv.',
     $      'it)', ' 12= |yh - wy| / (|h| |y| ulp)  (inv.it)' )
*
*     Symmetric/Hermitian eigenproblem
*
 9983 FORMAT( ' matrix types(see xdrvst for details): ' )
*
 9982 FORMAT( / ' special matrices:', / '  1=zero matrix.             ',
     $      '           ', '  5=diagonal: clustered entries.', / '  2=',
     $      'identity matrix.                    ', '  6=diagonal: lar',
     $      'ge, evenly spaced.', / '  3=diagonal: evenly spaced entri',
     $      'es.    ', '  7=diagonal: small, evenly spaced.', / '  4=d',
     $      'iagonal: geometr. spaced entries.' )
 9981 FORMAT( ' dense ', A, ' matrices:', / '  8=evenly spaced eigen',
     $      'vals.            ', ' 12=small, evenly spaced eigenvals.',
     $      / '  9=geometrically spaced eigenvals.     ', ' 13=matrix ',
     $      'with random o(1) entries.', / ' 10=clustered eigenvalues.',
     $      '              ', ' 14=matrix with large random entries.',
     $      / ' 11=large, evenly spaced eigenvals.     ', ' 15=matrix ',
     $      'with small random entries.' )
*
*     Symmetric/Hermitian Generalized eigenproblem
*
 9980 FORMAT( ' matrix types(see xdrvsg for details): ' )
*
 9979 FORMAT( / ' special matrices:', / '  1=zero matrix.             ',
     $      '           ', '  5=diagonal: clustered entries.', / '  2=',
     $      'identity matrix.                    ', '  6=diagonal: lar',
     $      'ge, evenly spaced.', / '  3=diagonal: evenly spaced entri',
     $      'es.    ', '  7=diagonal: small, evenly spaced.', / '  4=d',
     $      'iagonal: geometr. spaced entries.' )
 9978 FORMAT( ' dense or banded ', A, ' matrices: ',
     $      / '  8=evenly spaced eigenvals.         ',
     $      ' 15=matrix with small random entries.',
     $      / '  9=geometrically spaced eigenvals.  ',
     $      ' 16=evenly spaced eigenvals, ka=1, kb=1.',
     $      / ' 10=Clustered eigenvalues.           ',
     $      ' 17=Evenly spaced eigenvals, KA=2, KB=1.',
     $      / ' 11=Large, evenly spaced eigenvals.  ',
     $      ' 18=Evenly spaced eigenvals, KA=2, KB=2.',
     $      / ' 12=Small, evenly spaced eigenvals.  ',
     $      ' 19=Evenly spaced eigenvals, KA=3, KB=1.',
     $      / ' 13=Matrix with random O(1) entries. ',
     $      ' 20=Evenly spaced eigenvals, KA=3, KB=2.',
     $      / ' 14=Matrix with large random entries.',
     $      ' 21=Evenly spaced eigenvals, KA=3, KB=3.' )
 9977 FORMAT( / ' Tests performed:   ',
     $      / '( For each pair (A,B), where A is of the given type ',
     $      / ' and B is a random well-conditioned matrix. D is ',
     $      / ' diagonal, and Z is orthogonal. )',
     $      / ' 1 = DSYGV, with ITYPE=1 and UPLO=''U'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 2 = DSPGV, with ITYPE=1 and UPLO=''U'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 3 = DSBGV, with ITYPE=1 and UPLO=''U'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 4 = DSYGV, with ITYPE=1 and UPLO=''L'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 5 = DSPGV, with ITYPE=1 and UPLO=''L'':',
     $      '  | a z - b z d | / ( |a| |z| n ulp )     ',
     $      / ' 6 = dsbgv, with itype=1 and uplo=''l'':',
     $      '  | a z - b z d | / ( |a| |z| n ulp )     ' )
 9976 FORMAT( ' 7 = dsygv, with itype=2 and uplo=''u'':',
     $      '  | a b z - z d | / ( |a| |z| n ulp )     ',
     $      / ' 8 = dspgv, with itype=2 and uplo=''u'':',
     $      '  | a b z - z d | / ( |a| |z| n ulp )     ',
     $      / ' 9 = dspgv, with itype=2 and uplo=''l'':',
     $      '  | a b z - z d | / ( |a| |z| n ulp )     ',
     $      / '10 = dspgv, with itype=2 and uplo=''l'':',
     $      '  | a b z - z d | / ( |a| |z| n ulp )     ',
     $      / '11 = dsygv, with itype=3 and uplo=''u'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ',
     $      / '12 = dspgv, with itype=3 and uplo=''u'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ',
     $      / '13 = dsygv, with itype=3 and uplo=''l'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ',
     $      / '14 = dspgv, with itype=3 and uplo=''l'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ' )
 9975 FORMAT( / ' tests performed:   ',
     $      / '( for each pair (a,b), where a is of the given type ',
     $      / ' and b is a random well-conditioned matrix. d is ',
     $      / ' diagonal, and z is unitary. )',
     $      / ' 1 = zhegv, with itype=1 and uplo=''u'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 2 = ZHPGV, with ITYPE=1 and UPLO=''U'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 3 = ZHBGV, with ITYPE=1 and UPLO=''U'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 4 = ZHEGV, with ITYPE=1 and UPLO=''L'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 5 = ZHPGV, with ITYPE=1 and UPLO=''L'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 6 = ZHBGV, with ITYPE=1 and UPLO=''L'':',
     $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ' )
 9974 FORMAT( ' 7 = ZHEGV, with ITYPE=2 and UPLO=''U'':',
     $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 8 = ZHPGV, with ITYPE=2 and UPLO=''U'':',
     $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
     $      / ' 9 = ZHPGV, with ITYPE=2 and UPLO=''L'':',
     $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
     $      / '10 = ZHPGV, with ITYPE=2 and UPLO=''L'':',
     $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
     $      / '11 = zhegv, with itype=3 and uplo=''u'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ',
     $      / '12 = zhpgv, with itype=3 and uplo=''u'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ',
     $      / '13 = zhegv, with itype=3 and uplo=''l'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ',
     $      / '14 = zhpgv, with itype=3 and uplo=''l'':',
     $      '  | b a z - z d | / ( |a| |z| n ulp )     ' )
*
*     Singular Value Decomposition
*
 9973 FORMAT( ' matrix types(see xchkbd for details):',
     $      / ' diagonal matrices:', / '   1: zero', 28X,
     $      ' 5: clustered entries', / '   2: identity', 24X,
     $      ' 6: large, evenly spaced entries',
     $      / '   3: evenly spaced entries', 11X,
     $      ' 7: small, evenly spaced entries',
     $      / '   4: geometrically spaced entries',
     $      / ' general matrices:', / '   8: evenly spaced sing. vals.',
     $      7X, '12: small, evenly spaced sing vals',
     $      / '   9: geometrically spaced sing vals  ',
     $      '13: random, o(1) entries', / '  10: clustered sing. vals.',
     $      11X, '14: random, scaled near overflow',
     $      / '  11: large, evenly spaced sing vals  ',
     $      '15: random, scaled near underflow' )
*
 9972 FORMAT( / ' test ratios:  ',
     $      '(b: bidiagonal, s: diagonal, q, p, u, and v: ', A10, / 16X,
     $      'x: m x nrhs, y = q'' x, and z = u'' y)' )
 9971 FORMAT( '   1: norm( a - q b p'' ) / ( norm(a) max(m,n) ulp )',
     $      / '   2: norm( i - q'' q )   / ( m ulp )',
     $      / '   3: norm( i - p'' p )   / ( n ulp )',
     $      / '   4: norm( b - u s v'' ) / ( norm(b) min(m,n) ulp )',
     $      / '   5: norm( y - u z )    / ',
     $        '( norm(z) max(min(m,n),k) ulp )',
     $      / '   6: norm( i - u'' u )   / ( min(m,n) ulp )',
     $      / '   7: norm( i - v'' v )   / ( min(m,n) ulp )',
     $      / '   8: test ordering of s(0 if nondecreasing, 1/ulp ',
     $        ' otherwise)',
     $      / '   9: norm( s - s1 )     / ( norm(s) ulp ),',
     $        ' where s1 is computed', / 43X,
     $        ' without computing u and v''',
     $      / '  10: sturm sequence test ',
     $        '(0 if sing. vals of b within thresh of s)',
     $      / '  11: norm( a - (qu) s(v'' p'') ) / ',
     $        '( norm(a) max(m,n) ulp )',
     $      / '  12: norm( x - (qu) z )         / ( |x| max(m,k) ulp )',
     $      / '  13: norm( i - (qu)''(qu) )      / ( m ulp )',
     $      / '  14: norm( i - (v'' p'') (p v) )  / ( n ulp )',
     $      / '  15: norm( b - u s v'' ) / ( norm(b) min(m,n) ulp )',
     $      / '  16: norm( i - u'' u )   / ( min(m,n) ulp )',
     $      / '  17: norm( i - v'' v )   / ( min(m,n) ulp )',
     $      / '  18: test ordering of s(0 if nondecreasing, 1/ulp ',
     $        ' otherwise)',
     $      / '  19: norm( s - s1 )     / ( norm(s) ulp ),',
     $        ' where s1 is computed', / 43X,
     $        ' without computing u and v''',
     $      / '  20: norm( b - u s v'' )  / ( norm(b) min(m,n) ulp )',
     $        '  dbdsvx(v,a)',
     $      / '  21: norm( i - u'' u )    / ( min(m,n) ulp )',
     $      / '  22: norm( i - v'' v )    / ( min(m,n) ulp )',
     $      / '  23: test ordering of s(0 if nondecreasing, 1/ulp ',
     $        ' otherwise)',
     $      / '  24: norm( s - s1 )      / ( norm(s) ulp ),',
     $        ' where s1 is computed', / 44X,
     $        ' without computing u and v''',
     $      / '  25: norm( s - u'' b v ) / ( norm(b) n ulp )',
     $        '  dbdsvx(v,i)',
     $      / '  26: norm( i - u'' u )    / ( min(m,n) ulp )',
     $      / '  27: norm( i - v'' v )    / ( min(m,n) ulp )',
     $      / '  28: test ordering of s(0 if nondecreasing, 1/ulp ',
     $        ' otherwise)',
     $      / '  29: norm( s - s1 )      / ( norm(s) ulp ),',
     $        ' where s1 is computed', / 44X,
     $        ' without computing u and v''',
     $      / '  30: norm( s - u'' b v ) / ( norm(b) n ulp )',
     $        '  dbdsvx(v,v)',
     $      / '  31: norm( i - u'' u )    / ( min(m,n) ulp )',
     $      / '  32: norm( i - v'' v )    / ( min(m,n) ulp )',
     $      / '  33: test ordering of s(0 if nondecreasing, 1/ulp ',
     $        ' otherwise)',
     $      / '  34: norm( s - s1 )      / ( norm(s) ulp ),',
     $        ' where s1 is computed', / 44X,
     $        ' without computing u and v''' )
*
*     Band reduction to bidiagonal form
*
 9970 FORMAT( ' matrix types(see xchkbb for details):',
     $      / ' diagonal matrices:', / '   1: zero', 28X,
     $      ' 5: clustered entries', / '   2: identity', 24X,
     $      ' 6: large, evenly spaced entries',
     $      / '   3: evenly spaced entries', 11X,
     $      ' 7: small, evenly spaced entries',
     $      / '   4: geometrically spaced entries',
     $      / ' general matrices:', / '   8: evenly spaced sing. vals.',
     $      7X, '12: small, evenly spaced sing vals',
     $      / '   9: geometrically spaced sing vals  ',
     $      '13: random, o(1) entries', / '  10: clustered sing. vals.',
     $      11X, '14: random, scaled near overflow',
     $      / '  11: large, evenly spaced sing vals  ',
     $      '15: random, scaled near underflow' )
*
 9969 FORMAT( / ' test ratios:  ', '(b: upper bidiagonal, q and p: ',
     $      A10, / 16X, 'c: m x nrhs, pt = p'', y = q'' c)',
     $      / ' 1: norm( a - q b pt ) / ( norm(a) max(m,n) ulp )',
     $      / ' 2: norm( i - q'' q )   / ( m ulp )',
     $      / ' 3: norm( i - pt pt'' )   / ( n ulp )',
     $      / ' 4: norm( y - q'' c )   / ( norm(y) max(m,nrhs) ulp )' )
 9968 FORMAT( / ' tests performed:  see sdrvst.f' )
 9967 FORMAT( / ' tests performed:  see cdrvst.f' )
*
*     End of DLAHD2
*
      END