LAPACK 3.3.0

dlahd2.f

Go to the documentation of this file.
00001       SUBROUTINE DLAHD2( IOUNIT, PATH )
00002 *
00003 *  -- LAPACK auxiliary test routine (version 2.0) --
00004 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
00005 *     November 2006
00006 *
00007 *     .. Scalar Arguments ..
00008       CHARACTER*3        PATH
00009       INTEGER            IOUNIT
00010 *     ..
00011 *
00012 *  Purpose
00013 *  =======
00014 *
00015 *  DLAHD2 prints header information for the different test paths.
00016 *
00017 *  Arguments
00018 *  =========
00019 *
00020 *  IOUNIT  (input) INTEGER.
00021 *          On entry, IOUNIT specifies the unit number to which the
00022 *          header information should be printed.
00023 *
00024 *  PATH    (input) CHARACTER*3.
00025 *          On entry, PATH contains the name of the path for which the
00026 *          header information is to be printed.  Current paths are
00027 *
00028 *             DHS, ZHS:  Non-symmetric eigenproblem.
00029 *             DST, ZST:  Symmetric eigenproblem.
00030 *             DSG, ZSG:  Symmetric Generalized eigenproblem.
00031 *             DBD, ZBD:  Singular Value Decomposition (SVD)
00032 *             DBB, ZBB:  General Banded reduction to bidiagonal form
00033 *
00034 *          These paths also are supplied in double precision (replace
00035 *          leading S by D and leading C by Z in path names).
00036 *
00037 *  =====================================================================
00038 *
00039 *     .. Local Scalars ..
00040       LOGICAL            CORZ, SORD
00041       CHARACTER*2        C2
00042       INTEGER            J
00043 *     ..
00044 *     .. External Functions ..
00045       LOGICAL            LSAME, LSAMEN
00046       EXTERNAL           LSAME, LSAMEN
00047 *     ..
00048 *     .. Executable Statements ..
00049 *
00050       IF( IOUNIT.LE.0 )
00051      $   RETURN
00052       SORD = LSAME( PATH, 'S' ) .OR. LSAME( PATH, 'D' )
00053       CORZ = LSAME( PATH, 'C' ) .OR. LSAME( PATH, 'Z' )
00054       IF( .NOT.SORD .AND. .NOT.CORZ ) THEN
00055          WRITE( IOUNIT, FMT = 9999 )PATH
00056       END IF
00057       C2 = PATH( 2: 3 )
00058 *
00059       IF( LSAMEN( 2, C2, 'HS' ) ) THEN
00060          IF( SORD ) THEN
00061 *
00062 *           Real Non-symmetric Eigenvalue Problem:
00063 *
00064             WRITE( IOUNIT, FMT = 9998 )PATH
00065 *
00066 *           Matrix types
00067 *
00068             WRITE( IOUNIT, FMT = 9988 )
00069             WRITE( IOUNIT, FMT = 9987 )
00070             WRITE( IOUNIT, FMT = 9986 )'pairs ', 'pairs ', 'prs.',
00071      $         'prs.'
00072             WRITE( IOUNIT, FMT = 9985 )
00073 *
00074 *           Tests performed
00075 *
00076             WRITE( IOUNIT, FMT = 9984 )'orthogonal', '''=transpose',
00077      $         ( '''', J = 1, 6 )
00078 *
00079          ELSE
00080 *
00081 *           Complex Non-symmetric Eigenvalue Problem:
00082 *
00083             WRITE( IOUNIT, FMT = 9997 )PATH
00084 *
00085 *           Matrix types
00086 *
00087             WRITE( IOUNIT, FMT = 9988 )
00088             WRITE( IOUNIT, FMT = 9987 )
00089             WRITE( IOUNIT, FMT = 9986 )'e.vals', 'e.vals', 'e.vs',
00090      $         'e.vs'
00091             WRITE( IOUNIT, FMT = 9985 )
00092 *
00093 *           Tests performed
00094 *
00095             WRITE( IOUNIT, FMT = 9984 )'unitary', '*=conj.transp.',
00096      $         ( '*', J = 1, 6 )
00097          END IF
00098 *
00099       ELSE IF( LSAMEN( 2, C2, 'ST' ) ) THEN
00100 *
00101          IF( SORD ) THEN
00102 *
00103 *           Real Symmetric Eigenvalue Problem:
00104 *
00105             WRITE( IOUNIT, FMT = 9996 )PATH
00106 *
00107 *           Matrix types
00108 *
00109             WRITE( IOUNIT, FMT = 9983 )
00110             WRITE( IOUNIT, FMT = 9982 )
00111             WRITE( IOUNIT, FMT = 9981 )'Symmetric'
00112 *
00113 *           Tests performed
00114 *
00115             WRITE( IOUNIT, FMT = 9968 )
00116 *
00117          ELSE
00118 *
00119 *           Complex Hermitian Eigenvalue Problem:
00120 *
00121             WRITE( IOUNIT, FMT = 9995 )PATH
00122 *
00123 *           Matrix types
00124 *
00125             WRITE( IOUNIT, FMT = 9983 )
00126             WRITE( IOUNIT, FMT = 9982 )
00127             WRITE( IOUNIT, FMT = 9981 )'Hermitian'
00128 *
00129 *           Tests performed
00130 *
00131             WRITE( IOUNIT, FMT = 9967 )
00132          END IF
00133 *
00134       ELSE IF( LSAMEN( 2, C2, 'SG' ) ) THEN
00135 *
00136          IF( SORD ) THEN
00137 *
00138 *           Real Symmetric Generalized Eigenvalue Problem:
00139 *
00140             WRITE( IOUNIT, FMT = 9992 )PATH
00141 *
00142 *           Matrix types
00143 *
00144             WRITE( IOUNIT, FMT = 9980 )
00145             WRITE( IOUNIT, FMT = 9979 )
00146             WRITE( IOUNIT, FMT = 9978 )'Symmetric'
00147 *
00148 *           Tests performed
00149 *
00150             WRITE( IOUNIT, FMT = 9977 )
00151             WRITE( IOUNIT, FMT = 9976 )
00152 *
00153          ELSE
00154 *
00155 *           Complex Hermitian Generalized Eigenvalue Problem:
00156 *
00157             WRITE( IOUNIT, FMT = 9991 )PATH
00158 *
00159 *           Matrix types
00160 *
00161             WRITE( IOUNIT, FMT = 9980 )
00162             WRITE( IOUNIT, FMT = 9979 )
00163             WRITE( IOUNIT, FMT = 9978 )'Hermitian'
00164 *
00165 *           Tests performed
00166 *
00167             WRITE( IOUNIT, FMT = 9975 )
00168             WRITE( IOUNIT, FMT = 9974 )
00169 *
00170          END IF
00171 *
00172       ELSE IF( LSAMEN( 2, C2, 'BD' ) ) THEN
00173 *
00174          IF( SORD ) THEN
00175 *
00176 *           Real Singular Value Decomposition:
00177 *
00178             WRITE( IOUNIT, FMT = 9994 )PATH
00179 *
00180 *           Matrix types
00181 *
00182             WRITE( IOUNIT, FMT = 9973 )
00183 *
00184 *           Tests performed
00185 *
00186             WRITE( IOUNIT, FMT = 9972 )'orthogonal'
00187             WRITE( IOUNIT, FMT = 9971 )
00188          ELSE
00189 *
00190 *           Complex Singular Value Decomposition:
00191 *
00192             WRITE( IOUNIT, FMT = 9993 )PATH
00193 *
00194 *           Matrix types
00195 *
00196             WRITE( IOUNIT, FMT = 9973 )
00197 *
00198 *           Tests performed
00199 *
00200             WRITE( IOUNIT, FMT = 9972 )'unitary   '
00201             WRITE( IOUNIT, FMT = 9971 )
00202          END IF
00203 *
00204       ELSE IF( LSAMEN( 2, C2, 'BB' ) ) THEN
00205 *
00206          IF( SORD ) THEN
00207 *
00208 *           Real General Band reduction to bidiagonal form:
00209 *
00210             WRITE( IOUNIT, FMT = 9990 )PATH
00211 *
00212 *           Matrix types
00213 *
00214             WRITE( IOUNIT, FMT = 9970 )
00215 *
00216 *           Tests performed
00217 *
00218             WRITE( IOUNIT, FMT = 9969 )'orthogonal'
00219          ELSE
00220 *
00221 *           Complex Band reduction to bidiagonal form:
00222 *
00223             WRITE( IOUNIT, FMT = 9989 )PATH
00224 *
00225 *           Matrix types
00226 *
00227             WRITE( IOUNIT, FMT = 9970 )
00228 *
00229 *           Tests performed
00230 *
00231             WRITE( IOUNIT, FMT = 9969 )'unitary   '
00232          END IF
00233 *
00234       ELSE
00235 *
00236          WRITE( IOUNIT, FMT = 9999 )PATH
00237          RETURN
00238       END IF
00239 *
00240       RETURN
00241 *
00242  9999 FORMAT( 1X, A3, ':  no header available' )
00243  9998 FORMAT( / 1X, A3, ' -- Real Non-symmetric eigenvalue problem' )
00244  9997 FORMAT( / 1X, A3, ' -- Complex Non-symmetric eigenvalue problem' )
00245  9996 FORMAT( / 1X, A3, ' -- Real Symmetric eigenvalue problem' )
00246  9995 FORMAT( / 1X, A3, ' -- Complex Hermitian eigenvalue problem' )
00247  9994 FORMAT( / 1X, A3, ' -- Real Singular Value Decomposition' )
00248  9993 FORMAT( / 1X, A3, ' -- Complex Singular Value Decomposition' )
00249  9992 FORMAT( / 1X, A3, ' -- Real Symmetric Generalized eigenvalue ',
00250      $      'problem' )
00251  9991 FORMAT( / 1X, A3, ' -- Complex Hermitian Generalized eigenvalue ',
00252      $      'problem' )
00253  9990 FORMAT( / 1X, A3, ' -- Real Band reduc. to bidiagonal form' )
00254  9989 FORMAT( / 1X, A3, ' -- Complex Band reduc. to bidiagonal form' )
00255 *
00256  9988 FORMAT( ' Matrix types (see xCHKHS for details): ' )
00257 *
00258  9987 FORMAT( / ' Special Matrices:', / '  1=Zero matrix.             ',
00259      $      '           ', '  5=Diagonal: geometr. spaced entries.',
00260      $      / '  2=Identity matrix.                    ', '  6=Diagona',
00261      $      'l: clustered entries.', / '  3=Transposed Jordan block.  ',
00262      $      '          ', '  7=Diagonal: large, evenly spaced.', / '  ',
00263      $      '4=Diagonal: evenly spaced entries.    ', '  8=Diagonal: s',
00264      $      'mall, evenly spaced.' )
00265  9986 FORMAT( ' Dense, Non-Symmetric Matrices:', / '  9=Well-cond., ev',
00266      $      'enly spaced eigenvals.', ' 14=Ill-cond., geomet. spaced e',
00267      $      'igenals.', / ' 10=Well-cond., geom. spaced eigenvals. ',
00268      $      ' 15=Ill-conditioned, clustered e.vals.', / ' 11=Well-cond',
00269      $      'itioned, clustered e.vals. ', ' 16=Ill-cond., random comp',
00270      $      'lex ', A6, / ' 12=Well-cond., random complex ', A6, '   ',
00271      $      ' 17=Ill-cond., large rand. complx ', A4, / ' 13=Ill-condi',
00272      $      'tioned, evenly spaced.     ', ' 18=Ill-cond., small rand.',
00273      $      ' complx ', A4 )
00274  9985 FORMAT( ' 19=Matrix with random O(1) entries.    ', ' 21=Matrix ',
00275      $      'with small random entries.', / ' 20=Matrix with large ran',
00276      $      'dom entries.   ' )
00277  9984 FORMAT( / ' Tests performed:   ', '(H is Hessenberg, T is Schur,',
00278      $      ' U and Z are ', A, ',', / 20X, A, ', W is a diagonal matr',
00279      $      'ix of eigenvalues,', / 20X, 'L and R are the left and rig',
00280      $      'ht eigenvector matrices)', / '  1 = | A - U H U', A1, ' |',
00281      $      ' / ( |A| n ulp )         ', '  2 = | I - U U', A1, ' | / ',
00282      $      '( n ulp )', / '  3 = | H - Z T Z', A1, ' | / ( |H| n ulp ',
00283      $      ')         ', '  4 = | I - Z Z', A1, ' | / ( n ulp )',
00284      $      / '  5 = | A - UZ T (UZ)', A1, ' | / ( |A| n ulp )     ',
00285      $      '  6 = | I - UZ (UZ)', A1, ' | / ( n ulp )', / '  7 = | T(',
00286      $      'e.vects.) - T(no e.vects.) | / ( |T| ulp )', / '  8 = | W',
00287      $      '(e.vects.) - W(no e.vects.) | / ( |W| ulp )', / '  9 = | ',
00288      $      'TR - RW | / ( |T| |R| ulp )     ', ' 10 = | LT - WL | / (',
00289      $      ' |T| |L| ulp )', / ' 11= |HX - XW| / (|H| |X| ulp)  (inv.',
00290      $      'it)', ' 12= |YH - WY| / (|H| |Y| ulp)  (inv.it)' )
00291 *
00292 *     Symmetric/Hermitian eigenproblem
00293 *
00294  9983 FORMAT( ' Matrix types (see xDRVST for details): ' )
00295 *
00296  9982 FORMAT( / ' Special Matrices:', / '  1=Zero matrix.             ',
00297      $      '           ', '  5=Diagonal: clustered entries.', / '  2=',
00298      $      'Identity matrix.                    ', '  6=Diagonal: lar',
00299      $      'ge, evenly spaced.', / '  3=Diagonal: evenly spaced entri',
00300      $      'es.    ', '  7=Diagonal: small, evenly spaced.', / '  4=D',
00301      $      'iagonal: geometr. spaced entries.' )
00302  9981 FORMAT( ' Dense ', A, ' Matrices:', / '  8=Evenly spaced eigen',
00303      $      'vals.            ', ' 12=Small, evenly spaced eigenvals.',
00304      $      / '  9=Geometrically spaced eigenvals.     ', ' 13=Matrix ',
00305      $      'with random O(1) entries.', / ' 10=Clustered eigenvalues.',
00306      $      '              ', ' 14=Matrix with large random entries.',
00307      $      / ' 11=Large, evenly spaced eigenvals.     ', ' 15=Matrix ',
00308      $      'with small random entries.' )
00309 *
00310 *     Symmetric/Hermitian Generalized eigenproblem
00311 *
00312  9980 FORMAT( ' Matrix types (see xDRVSG for details): ' )
00313 *
00314  9979 FORMAT( / ' Special Matrices:', / '  1=Zero matrix.             ',
00315      $      '           ', '  5=Diagonal: clustered entries.', / '  2=',
00316      $      'Identity matrix.                    ', '  6=Diagonal: lar',
00317      $      'ge, evenly spaced.', / '  3=Diagonal: evenly spaced entri',
00318      $      'es.    ', '  7=Diagonal: small, evenly spaced.', / '  4=D',
00319      $      'iagonal: geometr. spaced entries.' )
00320  9978 FORMAT( ' Dense or Banded ', A, ' Matrices: ',
00321      $      / '  8=Evenly spaced eigenvals.         ',
00322      $      ' 15=Matrix with small random entries.',
00323      $      / '  9=Geometrically spaced eigenvals.  ',
00324      $      ' 16=Evenly spaced eigenvals, KA=1, KB=1.',
00325      $      / ' 10=Clustered eigenvalues.           ',
00326      $      ' 17=Evenly spaced eigenvals, KA=2, KB=1.',
00327      $      / ' 11=Large, evenly spaced eigenvals.  ',
00328      $      ' 18=Evenly spaced eigenvals, KA=2, KB=2.',
00329      $      / ' 12=Small, evenly spaced eigenvals.  ',
00330      $      ' 19=Evenly spaced eigenvals, KA=3, KB=1.',
00331      $      / ' 13=Matrix with random O(1) entries. ',
00332      $      ' 20=Evenly spaced eigenvals, KA=3, KB=2.',
00333      $      / ' 14=Matrix with large random entries.',
00334      $      ' 21=Evenly spaced eigenvals, KA=3, KB=3.' )
00335  9977 FORMAT( / ' Tests performed:   ',
00336      $      / '( For each pair (A,B), where A is of the given type ',
00337      $      / ' and B is a random well-conditioned matrix. D is ',
00338      $      / ' diagonal, and Z is orthogonal. )',
00339      $      / ' 1 = DSYGV, with ITYPE=1 and UPLO=''U'':',
00340      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00341      $      / ' 2 = DSPGV, with ITYPE=1 and UPLO=''U'':',
00342      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00343      $      / ' 3 = DSBGV, with ITYPE=1 and UPLO=''U'':',
00344      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00345      $      / ' 4 = DSYGV, with ITYPE=1 and UPLO=''L'':',
00346      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00347      $      / ' 5 = DSPGV, with ITYPE=1 and UPLO=''L'':',
00348      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00349      $      / ' 6 = DSBGV, with ITYPE=1 and UPLO=''L'':',
00350      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ' )
00351  9976 FORMAT( ' 7 = DSYGV, with ITYPE=2 and UPLO=''U'':',
00352      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00353      $      / ' 8 = DSPGV, with ITYPE=2 and UPLO=''U'':',
00354      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00355      $      / ' 9 = DSPGV, with ITYPE=2 and UPLO=''L'':',
00356      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00357      $      / '10 = DSPGV, with ITYPE=2 and UPLO=''L'':',
00358      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00359      $      / '11 = DSYGV, with ITYPE=3 and UPLO=''U'':',
00360      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ',
00361      $      / '12 = DSPGV, with ITYPE=3 and UPLO=''U'':',
00362      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ',
00363      $      / '13 = DSYGV, with ITYPE=3 and UPLO=''L'':',
00364      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ',
00365      $      / '14 = DSPGV, with ITYPE=3 and UPLO=''L'':',
00366      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ' )
00367  9975 FORMAT( / ' Tests performed:   ',
00368      $      / '( For each pair (A,B), where A is of the given type ',
00369      $      / ' and B is a random well-conditioned matrix. D is ',
00370      $      / ' diagonal, and Z is unitary. )',
00371      $      / ' 1 = ZHEGV, with ITYPE=1 and UPLO=''U'':',
00372      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00373      $      / ' 2 = ZHPGV, with ITYPE=1 and UPLO=''U'':',
00374      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00375      $      / ' 3 = ZHBGV, with ITYPE=1 and UPLO=''U'':',
00376      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00377      $      / ' 4 = ZHEGV, with ITYPE=1 and UPLO=''L'':',
00378      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00379      $      / ' 5 = ZHPGV, with ITYPE=1 and UPLO=''L'':',
00380      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ',
00381      $      / ' 6 = ZHBGV, with ITYPE=1 and UPLO=''L'':',
00382      $      '  | A Z - B Z D | / ( |A| |Z| n ulp )     ' )
00383  9974 FORMAT( ' 7 = ZHEGV, with ITYPE=2 and UPLO=''U'':',
00384      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00385      $      / ' 8 = ZHPGV, with ITYPE=2 and UPLO=''U'':',
00386      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00387      $      / ' 9 = ZHPGV, with ITYPE=2 and UPLO=''L'':',
00388      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00389      $      / '10 = ZHPGV, with ITYPE=2 and UPLO=''L'':',
00390      $      '  | A B Z - Z D | / ( |A| |Z| n ulp )     ',
00391      $      / '11 = ZHEGV, with ITYPE=3 and UPLO=''U'':',
00392      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ',
00393      $      / '12 = ZHPGV, with ITYPE=3 and UPLO=''U'':',
00394      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ',
00395      $      / '13 = ZHEGV, with ITYPE=3 and UPLO=''L'':',
00396      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ',
00397      $      / '14 = ZHPGV, with ITYPE=3 and UPLO=''L'':',
00398      $      '  | B A Z - Z D | / ( |A| |Z| n ulp )     ' )
00399 *
00400 *     Singular Value Decomposition
00401 *
00402  9973 FORMAT( ' Matrix types (see xCHKBD for details):',
00403      $      / ' Diagonal matrices:', / '   1: Zero', 28X,
00404      $      ' 5: Clustered entries', / '   2: Identity', 24X,
00405      $      ' 6: Large, evenly spaced entries',
00406      $      / '   3: Evenly spaced entries', 11X,
00407      $      ' 7: Small, evenly spaced entries',
00408      $      / '   4: Geometrically spaced entries',
00409      $      / ' General matrices:', / '   8: Evenly spaced sing. vals.',
00410      $      7X, '12: Small, evenly spaced sing vals',
00411      $      / '   9: Geometrically spaced sing vals  ',
00412      $      '13: Random, O(1) entries', / '  10: Clustered sing. vals.',
00413      $      11X, '14: Random, scaled near overflow',
00414      $      / '  11: Large, evenly spaced sing vals  ',
00415      $      '15: Random, scaled near underflow' )
00416 *
00417  9972 FORMAT( / ' Test ratios:  ',
00418      $      '(B: bidiagonal, S: diagonal, Q, P, U, and V: ', A10, / 16X,
00419      $      'X: m x nrhs, Y = Q'' X, and Z = U'' Y)',
00420      $      / '   1: norm( A - Q B P'' ) / ( norm(A) max(m,n) ulp )',
00421      $      / '   2: norm( I - Q'' Q )   / ( m ulp )',
00422      $      / '   3: norm( I - P'' P )   / ( n ulp )',
00423      $      / '   4: norm( B - U S V'' ) / ( norm(B) min(m,n) ulp )', /
00424      $      '   5: norm( Y - U Z )    / ( norm(Z) max(min(m,n),k) ulp )'
00425      $      , / '   6: norm( I - U'' U )   / ( min(m,n) ulp )',
00426      $      / '   7: norm( I - V'' V )   / ( min(m,n) ulp )' )
00427  9971 FORMAT( '   8: Test ordering of S  (0 if nondecreasing, 1/ulp ',
00428      $      ' otherwise)', /
00429      $      '   9: norm( S - S2 )     / ( norm(S) ulp ),',
00430      $      ' where S2 is computed', / 44X,
00431      $      'without computing U and V''',
00432      $      / '  10: Sturm sequence test ',
00433      $      '(0 if sing. vals of B within THRESH of S)',
00434      $      / '  11: norm( A - (QU) S (V'' P'') ) / ',
00435      $      '( norm(A) max(m,n) ulp )', /
00436      $      '  12: norm( X - (QU) Z )         / ( |X| max(M,k) ulp )',
00437      $      / '  13: norm( I - (QU)''(QU) )      / ( M ulp )',
00438      $      / '  14: norm( I - (V'' P'') (P V) )  / ( N ulp )' )
00439 *
00440 *     Band reduction to bidiagonal form
00441 *
00442  9970 FORMAT( ' Matrix types (see xCHKBB for details):',
00443      $      / ' Diagonal matrices:', / '   1: Zero', 28X,
00444      $      ' 5: Clustered entries', / '   2: Identity', 24X,
00445      $      ' 6: Large, evenly spaced entries',
00446      $      / '   3: Evenly spaced entries', 11X,
00447      $      ' 7: Small, evenly spaced entries',
00448      $      / '   4: Geometrically spaced entries',
00449      $      / ' General matrices:', / '   8: Evenly spaced sing. vals.',
00450      $      7X, '12: Small, evenly spaced sing vals',
00451      $      / '   9: Geometrically spaced sing vals  ',
00452      $      '13: Random, O(1) entries', / '  10: Clustered sing. vals.',
00453      $      11X, '14: Random, scaled near overflow',
00454      $      / '  11: Large, evenly spaced sing vals  ',
00455      $      '15: Random, scaled near underflow' )
00456 *
00457  9969 FORMAT( / ' Test ratios:  ', '(B: upper bidiagonal, Q and P: ',
00458      $      A10, / 16X, 'C: m x nrhs, PT = P'', Y = Q'' C)',
00459      $      / ' 1: norm( A - Q B PT ) / ( norm(A) max(m,n) ulp )',
00460      $      / ' 2: norm( I - Q'' Q )   / ( m ulp )',
00461      $      / ' 3: norm( I - PT PT'' )   / ( n ulp )',
00462      $      / ' 4: norm( Y - Q'' C )   / ( norm(Y) max(m,nrhs) ulp )' )
00463  9968 FORMAT( / ' Tests performed:  See sdrvst.f' )
00464  9967 FORMAT( / ' Tests performed:  See cdrvst.f' )
00465 *
00466 *     End of DLAHD2
00467 *
00468       END
 All Files Functions