LAPACK: TESTING/LIN/alahd.f Source File

Go to the documentation of this file.
00001       SUBROUTINE ALAHD( IOUNIT, PATH )
00002 *
00003 *  -- LAPACK test routine (version 3.3.0) --
00004 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
00005 *     November 2010
00006 *
00007 *     .. Scalar Arguments ..
00008       CHARACTER*3        PATH
00009       INTEGER            IOUNIT
00010 *     ..
00011 *
00012 *  Purpose
00013 *  =======
00014 *
00015 *  ALAHD prints header information for the different test paths.
00016 *
00017 *  Arguments
00018 *  =========
00019 *
00020 *  IOUNIT  (input) INTEGER
00021 *          The unit number to which the header information should be
00022 *          printed.
00023 *
00024 *  PATH    (input) CHARACTER*3
00025 *          The name of the path for which the header information is to
00026 *          be printed.  Current paths are
00027 *             _GE:  General matrices
00028 *             _GB:  General band
00029 *             _GT:  General Tridiagonal
00030 *             _PO:  Symmetric or Hermitian positive definite
00031 *             _PS:  Symmetric or Hermitian positive semi-definite
00032 *             _PP:  Symmetric or Hermitian positive definite packed
00033 *             _PB:  Symmetric or Hermitian positive definite band
00034 *             _PT:  Symmetric or Hermitian positive definite tridiagonal
00035 *             _SY:  Symmetric indefinite
00036 *             _SP:  Symmetric indefinite packed
00037 *             _HE:  (complex) Hermitian indefinite
00038 *             _HP:  (complex) Hermitian indefinite packed
00039 *             _TR:  Triangular
00040 *             _TP:  Triangular packed
00041 *             _TB:  Triangular band
00042 *             _QR:  QR (general matrices)
00043 *             _LQ:  LQ (general matrices)
00044 *             _QL:  QL (general matrices)
00045 *             _RQ:  RQ (general matrices)
00046 *             _QP:  QR with column pivoting
00047 *             _TZ:  Trapezoidal
00048 *             _LS:  Least Squares driver routines
00049 *             _LU:  LU variants
00050 *             _CH:  Cholesky variants
00051 *             _QS:  QR variants
00052 *          The first character must be one of S, D, C, or Z (C or Z only
00053 *          if complex).
00054 *
00055 *  =====================================================================
00056 *
00057 *     .. Local Scalars ..
00058       LOGICAL            CORZ, SORD
00059       CHARACTER          C1, C3
00060       CHARACTER*2        P2
00061       CHARACTER*4        EIGCNM
00062       CHARACTER*32       SUBNAM
00063       CHARACTER*9        SYM
00064 *     ..
00065 *     .. External Functions ..
00066       LOGICAL            LSAME, LSAMEN
00067       EXTERNAL           LSAME, LSAMEN
00068 *     ..
00069 *     .. Intrinsic Functions ..
00070       INTRINSIC          LEN_TRIM
00071 *     ..
00072 *     .. Executable Statements ..
00073 *
00074       IF( IOUNIT.LE.0 )
00075      $   RETURN
00076       C1 = PATH( 1: 1 )
00077       C3 = PATH( 3: 3 )
00078       P2 = PATH( 2: 3 )
00079       SORD = LSAME( C1, 'S' ) .OR. LSAME( C1, 'D' )
00080       CORZ = LSAME( C1, 'C' ) .OR. LSAME( C1, 'Z' )
00081       IF( .NOT.( SORD .OR. CORZ ) )
00082      $   RETURN
00083 *
00084       IF( LSAMEN( 2, P2, 'GE' ) ) THEN
00085 *
00086 *        GE: General dense
00087 *
00088          WRITE( IOUNIT, FMT = 9999 )PATH
00089          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00090          WRITE( IOUNIT, FMT = 9979 )
00091          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00092          WRITE( IOUNIT, FMT = 9962 )1
00093          WRITE( IOUNIT, FMT = 9961 )2
00094          WRITE( IOUNIT, FMT = 9960 )3
00095          WRITE( IOUNIT, FMT = 9959 )4
00096          WRITE( IOUNIT, FMT = 9958 )5
00097          WRITE( IOUNIT, FMT = 9957 )6
00098          WRITE( IOUNIT, FMT = 9956 )7
00099          WRITE( IOUNIT, FMT = 9955 )8
00100          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00101 *
00102       ELSE IF( LSAMEN( 2, P2, 'GB' ) ) THEN
00103 *
00104 *        GB: General band
00105 *
00106          WRITE( IOUNIT, FMT = 9998 )PATH
00107          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00108          WRITE( IOUNIT, FMT = 9978 )
00109          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00110          WRITE( IOUNIT, FMT = 9962 )1
00111          WRITE( IOUNIT, FMT = 9960 )2
00112          WRITE( IOUNIT, FMT = 9959 )3
00113          WRITE( IOUNIT, FMT = 9958 )4
00114          WRITE( IOUNIT, FMT = 9957 )5
00115          WRITE( IOUNIT, FMT = 9956 )6
00116          WRITE( IOUNIT, FMT = 9955 )7
00117          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00118 *
00119       ELSE IF( LSAMEN( 2, P2, 'GT' ) ) THEN
00120 *
00121 *        GT: General tridiagonal
00122 *
00123          WRITE( IOUNIT, FMT = 9997 )PATH
00124          WRITE( IOUNIT, FMT = 9977 )
00125          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00126          WRITE( IOUNIT, FMT = 9962 )1
00127          WRITE( IOUNIT, FMT = 9960 )2
00128          WRITE( IOUNIT, FMT = 9959 )3
00129          WRITE( IOUNIT, FMT = 9958 )4
00130          WRITE( IOUNIT, FMT = 9957 )5
00131          WRITE( IOUNIT, FMT = 9956 )6
00132          WRITE( IOUNIT, FMT = 9955 )7
00133          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00134 *
00135       ELSE IF( LSAMEN( 2, P2, 'PO' ) .OR. LSAMEN( 2, P2, 'PP' ) ) THEN
00136 *
00137 *        PO: Positive definite full
00138 *        PP: Positive definite packed
00139 *
00140          IF( SORD ) THEN
00141             SYM = 'Symmetric'
00142          ELSE
00143             SYM = 'Hermitian'
00144          END IF
00145          IF( LSAME( C3, 'O' ) ) THEN
00146             WRITE( IOUNIT, FMT = 9996 )PATH, SYM
00147          ELSE
00148             WRITE( IOUNIT, FMT = 9995 )PATH, SYM
00149          END IF
00150          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00151          WRITE( IOUNIT, FMT = 9975 )PATH
00152          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00153          WRITE( IOUNIT, FMT = 9954 )1
00154          WRITE( IOUNIT, FMT = 9961 )2
00155          WRITE( IOUNIT, FMT = 9960 )3
00156          WRITE( IOUNIT, FMT = 9959 )4
00157          WRITE( IOUNIT, FMT = 9958 )5
00158          WRITE( IOUNIT, FMT = 9957 )6
00159          WRITE( IOUNIT, FMT = 9956 )7
00160          WRITE( IOUNIT, FMT = 9955 )8
00161          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00162 *
00163       ELSE IF( LSAMEN( 2, P2, 'PS' ) ) THEN
00164 *
00165 *        PS: Positive semi-definite full
00166 *
00167          IF( SORD ) THEN
00168             SYM = 'Symmetric'
00169          ELSE
00170             SYM = 'Hermitian'
00171          END IF
00172          IF( LSAME( C1, 'S' ) .OR. LSAME( C1, 'C' ) ) THEN
00173             EIGCNM = '1E04'
00174          ELSE
00175             EIGCNM = '1D12'
00176          END IF
00177          WRITE( IOUNIT, FMT = 9995 )PATH, SYM
00178          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00179          WRITE( IOUNIT, FMT = 8973 )EIGCNM, EIGCNM, EIGCNM
00180          WRITE( IOUNIT, FMT = '( '' Difference:'' )' )
00181          WRITE( IOUNIT, FMT = 8972 )C1
00182          WRITE( IOUNIT, FMT = '( '' Test ratio:'' )' )
00183          WRITE( IOUNIT, FMT = 8950 )
00184          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00185       ELSE IF( LSAMEN( 2, P2, 'PB' ) ) THEN
00186 *
00187 *        PB: Positive definite band
00188 *
00189          IF( SORD ) THEN
00190             WRITE( IOUNIT, FMT = 9994 )PATH, 'Symmetric'
00191          ELSE
00192             WRITE( IOUNIT, FMT = 9994 )PATH, 'Hermitian'
00193          END IF
00194          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00195          WRITE( IOUNIT, FMT = 9973 )PATH
00196          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00197          WRITE( IOUNIT, FMT = 9954 )1
00198          WRITE( IOUNIT, FMT = 9960 )2
00199          WRITE( IOUNIT, FMT = 9959 )3
00200          WRITE( IOUNIT, FMT = 9958 )4
00201          WRITE( IOUNIT, FMT = 9957 )5
00202          WRITE( IOUNIT, FMT = 9956 )6
00203          WRITE( IOUNIT, FMT = 9955 )7
00204          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00205 *
00206       ELSE IF( LSAMEN( 2, P2, 'PT' ) ) THEN
00207 *
00208 *        PT: Positive definite tridiagonal
00209 *
00210          IF( SORD ) THEN
00211             WRITE( IOUNIT, FMT = 9993 )PATH, 'Symmetric'
00212          ELSE
00213             WRITE( IOUNIT, FMT = 9993 )PATH, 'Hermitian'
00214          END IF
00215          WRITE( IOUNIT, FMT = 9976 )
00216          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00217          WRITE( IOUNIT, FMT = 9952 )1
00218          WRITE( IOUNIT, FMT = 9960 )2
00219          WRITE( IOUNIT, FMT = 9959 )3
00220          WRITE( IOUNIT, FMT = 9958 )4
00221          WRITE( IOUNIT, FMT = 9957 )5
00222          WRITE( IOUNIT, FMT = 9956 )6
00223          WRITE( IOUNIT, FMT = 9955 )7
00224          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00225 *
00226       ELSE IF( LSAMEN( 2, P2, 'SY' )  ) THEN
00227 *
00228 *        SY: Symmetric indefinite full
00229 *
00230          IF( LSAME( C3, 'Y' ) ) THEN
00231             WRITE( IOUNIT, FMT = 9992 )PATH, 'Symmetric'
00232          ELSE
00233             WRITE( IOUNIT, FMT = 9991 )PATH, 'Symmetric'
00234          END IF
00235          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00236          IF( SORD ) THEN
00237             WRITE( IOUNIT, FMT = 9972 )
00238          ELSE
00239             WRITE( IOUNIT, FMT = 9971 )
00240          END IF
00241          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00242          WRITE( IOUNIT, FMT = 9953 )1
00243          WRITE( IOUNIT, FMT = 9961 )2
00244          WRITE( IOUNIT, FMT = 9960 )3
00245          WRITE( IOUNIT, FMT = 9960 )4
00246          WRITE( IOUNIT, FMT = 9959 )5
00247          WRITE( IOUNIT, FMT = 9958 )6
00248          WRITE( IOUNIT, FMT = 9956 )7
00249          WRITE( IOUNIT, FMT = 9957 )8
00250          WRITE( IOUNIT, FMT = 9955 )9
00251          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00252 *
00253       ELSE IF( LSAMEN( 2, P2, 'SP' ) ) THEN
00254 *
00255 *        SP: Symmetric indefinite packed
00256 *
00257          IF( LSAME( C3, 'Y' ) ) THEN
00258             WRITE( IOUNIT, FMT = 9992 )PATH, 'Symmetric'
00259          ELSE
00260             WRITE( IOUNIT, FMT = 9991 )PATH, 'Symmetric'
00261          END IF
00262          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00263          IF( SORD ) THEN
00264             WRITE( IOUNIT, FMT = 9972 )
00265          ELSE
00266             WRITE( IOUNIT, FMT = 9971 )
00267          END IF
00268          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00269          WRITE( IOUNIT, FMT = 9953 )1
00270          WRITE( IOUNIT, FMT = 9961 )2
00271          WRITE( IOUNIT, FMT = 9960 )3
00272          WRITE( IOUNIT, FMT = 9959 )4
00273          WRITE( IOUNIT, FMT = 9958 )5
00274          WRITE( IOUNIT, FMT = 9956 )6
00275          WRITE( IOUNIT, FMT = 9957 )7
00276          WRITE( IOUNIT, FMT = 9955 )8
00277          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00278 *
00279       ELSE IF( LSAMEN( 2, P2, 'HE' )  ) THEN
00280 *
00281 *        HE: Hermitian indefinite full
00282 *
00283          IF( LSAME( C3, 'E' ) ) THEN
00284             WRITE( IOUNIT, FMT = 9992 )PATH, 'Hermitian'
00285          ELSE
00286             WRITE( IOUNIT, FMT = 9991 )PATH, 'Hermitian'
00287          END IF
00288          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00289          IF( SORD ) THEN
00290             WRITE( IOUNIT, FMT = 9972 )
00291          ELSE
00292             WRITE( IOUNIT, FMT = 9971 )
00293          END IF
00294          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00295          WRITE( IOUNIT, FMT = 9953 )1
00296          WRITE( IOUNIT, FMT = 9961 )2
00297          WRITE( IOUNIT, FMT = 9960 )3
00298          WRITE( IOUNIT, FMT = 9960 )4
00299          WRITE( IOUNIT, FMT = 9959 )5
00300          WRITE( IOUNIT, FMT = 9958 )6
00301          WRITE( IOUNIT, FMT = 9956 )7
00302          WRITE( IOUNIT, FMT = 9957 )8
00303          WRITE( IOUNIT, FMT = 9955 )9
00304          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00305 *
00306       ELSE IF( LSAMEN( 2, P2, 'HP' ) ) THEN
00307 *
00308 *        HP: Hermitian indefinite packed
00309 *
00310          IF( LSAME( C3, 'E' ) ) THEN
00311             WRITE( IOUNIT, FMT = 9992 )PATH, 'Hermitian'
00312          ELSE
00313             WRITE( IOUNIT, FMT = 9991 )PATH, 'Hermitian'
00314          END IF
00315          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00316          WRITE( IOUNIT, FMT = 9972 )
00317          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00318          WRITE( IOUNIT, FMT = 9953 )1
00319          WRITE( IOUNIT, FMT = 9961 )2
00320          WRITE( IOUNIT, FMT = 9960 )3
00321          WRITE( IOUNIT, FMT = 9959 )4
00322          WRITE( IOUNIT, FMT = 9958 )5
00323          WRITE( IOUNIT, FMT = 9956 )6
00324          WRITE( IOUNIT, FMT = 9957 )7
00325          WRITE( IOUNIT, FMT = 9955 )8
00326          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00327 *
00328       ELSE IF( LSAMEN( 2, P2, 'TR' ) .OR. LSAMEN( 2, P2, 'TP' ) ) THEN
00329 *
00330 *        TR: Triangular full
00331 *        TP: Triangular packed
00332 *
00333          IF( LSAME( C3, 'R' ) ) THEN
00334             WRITE( IOUNIT, FMT = 9990 )PATH
00335             SUBNAM = PATH( 1: 1 ) // 'LATRS'
00336          ELSE
00337             WRITE( IOUNIT, FMT = 9989 )PATH
00338             SUBNAM = PATH( 1: 1 ) // 'LATPS'
00339          END IF
00340          WRITE( IOUNIT, FMT = 9966 )PATH
00341          WRITE( IOUNIT, FMT = 9965 )SUBNAM(1:LEN_TRIM( SUBNAM ))
00342          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00343          WRITE( IOUNIT, FMT = 9961 )1
00344          WRITE( IOUNIT, FMT = 9960 )2
00345          WRITE( IOUNIT, FMT = 9959 )3
00346          WRITE( IOUNIT, FMT = 9958 )4
00347          WRITE( IOUNIT, FMT = 9957 )5
00348          WRITE( IOUNIT, FMT = 9956 )6
00349          WRITE( IOUNIT, FMT = 9955 )7
00350          WRITE( IOUNIT, FMT = 9951 )SUBNAM(1:LEN_TRIM( SUBNAM )), 8
00351          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00352 *
00353       ELSE IF( LSAMEN( 2, P2, 'TB' ) ) THEN
00354 *
00355 *        TB: Triangular band
00356 *
00357          WRITE( IOUNIT, FMT = 9988 )PATH
00358          SUBNAM = PATH( 1: 1 ) // 'LATBS'
00359          WRITE( IOUNIT, FMT = 9964 )PATH
00360          WRITE( IOUNIT, FMT = 9963 )SUBNAM(1:LEN_TRIM( SUBNAM ))
00361          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00362          WRITE( IOUNIT, FMT = 9960 )1
00363          WRITE( IOUNIT, FMT = 9959 )2
00364          WRITE( IOUNIT, FMT = 9958 )3
00365          WRITE( IOUNIT, FMT = 9957 )4
00366          WRITE( IOUNIT, FMT = 9956 )5
00367          WRITE( IOUNIT, FMT = 9955 )6
00368          WRITE( IOUNIT, FMT = 9951 )SUBNAM(1:LEN_TRIM( SUBNAM )), 7
00369          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00370 *
00371       ELSE IF( LSAMEN( 2, P2, 'QR' ) ) THEN
00372 *
00373 *        QR decomposition of rectangular matrices
00374 *
00375          WRITE( IOUNIT, FMT = 9987 )PATH, 'QR'
00376          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00377          WRITE( IOUNIT, FMT = 9970 )
00378          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00379          WRITE( IOUNIT, FMT = 9950 )1
00380          WRITE( IOUNIT, FMT = 6950 )8
00381          WRITE( IOUNIT, FMT = 9946 )2
00382          WRITE( IOUNIT, FMT = 9944 )3, 'M'
00383          WRITE( IOUNIT, FMT = 9943 )4, 'M'
00384          WRITE( IOUNIT, FMT = 9942 )5, 'M'
00385          WRITE( IOUNIT, FMT = 9941 )6, 'M'
00386          WRITE( IOUNIT, FMT = 9960 )7
00387          WRITE( IOUNIT, FMT = 6660 )9
00388          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00389 *
00390       ELSE IF( LSAMEN( 2, P2, 'LQ' ) ) THEN
00391 *
00392 *        LQ decomposition of rectangular matrices
00393 *
00394          WRITE( IOUNIT, FMT = 9987 )PATH, 'LQ'
00395          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00396          WRITE( IOUNIT, FMT = 9970 )
00397          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00398          WRITE( IOUNIT, FMT = 9949 )1
00399          WRITE( IOUNIT, FMT = 9945 )2
00400          WRITE( IOUNIT, FMT = 9944 )3, 'N'
00401          WRITE( IOUNIT, FMT = 9943 )4, 'N'
00402          WRITE( IOUNIT, FMT = 9942 )5, 'N'
00403          WRITE( IOUNIT, FMT = 9941 )6, 'N'
00404          WRITE( IOUNIT, FMT = 9960 )7
00405          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00406 *
00407       ELSE IF( LSAMEN( 2, P2, 'QL' ) ) THEN
00408 *
00409 *        QL decomposition of rectangular matrices
00410 *
00411          WRITE( IOUNIT, FMT = 9987 )PATH, 'QL'
00412          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00413          WRITE( IOUNIT, FMT = 9970 )
00414          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00415          WRITE( IOUNIT, FMT = 9948 )1
00416          WRITE( IOUNIT, FMT = 9946 )2
00417          WRITE( IOUNIT, FMT = 9944 )3, 'M'
00418          WRITE( IOUNIT, FMT = 9943 )4, 'M'
00419          WRITE( IOUNIT, FMT = 9942 )5, 'M'
00420          WRITE( IOUNIT, FMT = 9941 )6, 'M'
00421          WRITE( IOUNIT, FMT = 9960 )7
00422          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00423 *
00424       ELSE IF( LSAMEN( 2, P2, 'RQ' ) ) THEN
00425 *
00426 *        RQ decomposition of rectangular matrices
00427 *
00428          WRITE( IOUNIT, FMT = 9987 )PATH, 'RQ'
00429          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00430          WRITE( IOUNIT, FMT = 9970 )
00431          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00432          WRITE( IOUNIT, FMT = 9947 )1
00433          WRITE( IOUNIT, FMT = 9945 )2
00434          WRITE( IOUNIT, FMT = 9944 )3, 'N'
00435          WRITE( IOUNIT, FMT = 9943 )4, 'N'
00436          WRITE( IOUNIT, FMT = 9942 )5, 'N'
00437          WRITE( IOUNIT, FMT = 9941 )6, 'N'
00438          WRITE( IOUNIT, FMT = 9960 )7
00439          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00440 *
00441       ELSE IF( LSAMEN( 2, P2, 'QP' ) ) THEN
00442 *
00443 *        QR decomposition with column pivoting
00444 *
00445          WRITE( IOUNIT, FMT = 9986 )PATH
00446          WRITE( IOUNIT, FMT = 9969 )
00447          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00448          WRITE( IOUNIT, FMT = 9940 )1
00449          WRITE( IOUNIT, FMT = 9939 )2
00450          WRITE( IOUNIT, FMT = 9938 )3
00451          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00452 *
00453       ELSE IF( LSAMEN( 2, P2, 'TZ' ) ) THEN
00454 *
00455 *        TZ:  Trapezoidal
00456 *
00457          WRITE( IOUNIT, FMT = 9985 )PATH
00458          WRITE( IOUNIT, FMT = 9968 )
00459          WRITE( IOUNIT, FMT = 9929 )C1, C1
00460          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00461          WRITE( IOUNIT, FMT = 9940 )1
00462          WRITE( IOUNIT, FMT = 9937 )2
00463          WRITE( IOUNIT, FMT = 9938 )3
00464          WRITE( IOUNIT, FMT = 9940 )4
00465          WRITE( IOUNIT, FMT = 9937 )5
00466          WRITE( IOUNIT, FMT = 9938 )6
00467          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00468 *
00469       ELSE IF( LSAMEN( 2, P2, 'LS' ) ) THEN
00470 *
00471 *        LS:  Least Squares driver routines for
00472 *             LS, LSD, LSS, LSX and LSY.
00473 *
00474          WRITE( IOUNIT, FMT = 9984 )PATH
00475          WRITE( IOUNIT, FMT = 9967 )
00476          WRITE( IOUNIT, FMT = 9921 )C1, C1, C1, C1, C1
00477          WRITE( IOUNIT, FMT = 9935 )1
00478          WRITE( IOUNIT, FMT = 9931 )2
00479          WRITE( IOUNIT, FMT = 9933 )3
00480          WRITE( IOUNIT, FMT = 9935 )4
00481          WRITE( IOUNIT, FMT = 9934 )5
00482          WRITE( IOUNIT, FMT = 9932 )6
00483          WRITE( IOUNIT, FMT = 9920 )
00484          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00485 *
00486       ELSE IF( LSAMEN( 2, P2, 'LU' ) ) THEN
00487 *
00488 *        LU factorization variants
00489 *
00490          WRITE( IOUNIT, FMT = 9983 )PATH
00491          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00492          WRITE( IOUNIT, FMT = 9979 )
00493          WRITE( IOUNIT, FMT = '( '' Test ratio:'' )' )
00494          WRITE( IOUNIT, FMT = 9962 )1
00495          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00496 *
00497       ELSE IF( LSAMEN( 2, P2, 'CH' ) ) THEN
00498 *
00499 *        Cholesky factorization variants
00500 *
00501          WRITE( IOUNIT, FMT = 9982 )PATH
00502          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00503          WRITE( IOUNIT, FMT = 9974 )
00504          WRITE( IOUNIT, FMT = '( '' Test ratio:'' )' )
00505          WRITE( IOUNIT, FMT = 9954 )1
00506          WRITE( IOUNIT, FMT = '( '' Messages:'' )' )
00507 *
00508       ELSE IF( LSAMEN( 2, P2, 'QS' ) ) THEN
00509 *
00510 *        QR factorization variants
00511 *
00512          WRITE( IOUNIT, FMT = 9981 )PATH
00513          WRITE( IOUNIT, FMT = '( '' Matrix types:'' )' )
00514          WRITE( IOUNIT, FMT = 9970 )
00515          WRITE( IOUNIT, FMT = '( '' Test ratios:'' )' )
00516 *
00517       ELSE
00518 *
00519 *        Print error message if no header is available.
00520 *
00521          WRITE( IOUNIT, FMT = 9980 )PATH
00522       END IF
00523 *
00524 *     First line of header
00525 *
00526  9999 FORMAT( / 1X, A3, ':  General dense matrices' )
00527  9998 FORMAT( / 1X, A3, ':  General band matrices' )
00528  9997 FORMAT( / 1X, A3, ':  General tridiagonal' )
00529  9996 FORMAT( / 1X, A3, ':  ', A9, ' positive definite matrices' )
00530  9995 FORMAT( / 1X, A3, ':  ', A9, ' positive definite packed matrices'
00531      $       )
00532  9994 FORMAT( / 1X, A3, ':  ', A9, ' positive definite band matrices' )
00533  9993 FORMAT( / 1X, A3, ':  ', A9, ' positive definite tridiagonal' )
00534  9992 FORMAT( / 1X, A3, ':  ', A9, ' indefinite matrices' )
00535  9991 FORMAT( / 1X, A3, ':  ', A9, ' indefinite packed matrices' )
00536  9990 FORMAT( / 1X, A3, ':  Triangular matrices' )
00537  9989 FORMAT( / 1X, A3, ':  Triangular packed matrices' )
00538  9988 FORMAT( / 1X, A3, ':  Triangular band matrices' )
00539  9987 FORMAT( / 1X, A3, ':  ', A2, ' factorization of general matrices'
00540      $       )
00541  9986 FORMAT( / 1X, A3, ':  QR factorization with column pivoting' )
00542  9985 FORMAT( / 1X, A3, ':  RQ factorization of trapezoidal matrix' )
00543  9984 FORMAT( / 1X, A3, ':  Least squares driver routines' )
00544  9983 FORMAT( / 1X, A3, ':  LU factorization variants' )
00545  9982 FORMAT( / 1X, A3, ':  Cholesky factorization variants' )
00546  9981 FORMAT( / 1X, A3, ':  QR factorization variants' )
00547  9980 FORMAT( / 1X, A3, ':  No header available' )
00548 *
00549 *     GE matrix types
00550 *
00551  9979 FORMAT( 4X, '1. Diagonal', 24X, '7. Last n/2 columns zero', / 4X,
00552      $      '2. Upper triangular', 16X,
00553      $      '8. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00554      $      '3. Lower triangular', 16X, '9. Random, CNDNUM = 0.1/EPS',
00555      $      / 4X, '4. Random, CNDNUM = 2', 13X,
00556      $      '10. Scaled near underflow', / 4X, '5. First column zero',
00557      $      14X, '11. Scaled near overflow', / 4X,
00558      $      '6. Last column zero' )
00559 *
00560 *     GB matrix types
00561 *
00562  9978 FORMAT( 4X, '1. Random, CNDNUM = 2', 14X,
00563      $      '5. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00564      $      '2. First column zero', 15X, '6. Random, CNDNUM = .01/EPS',
00565      $      / 4X, '3. Last column zero', 16X,
00566      $      '7. Scaled near underflow', / 4X,
00567      $      '4. Last n/2 columns zero', 11X, '8. Scaled near overflow' )
00568 *
00569 *     GT matrix types
00570 *
00571  9977 FORMAT( ' Matrix types (1-6 have specified condition numbers):',
00572      $      / 4X, '1. Diagonal', 24X, '7. Random, unspecified CNDNUM',
00573      $      / 4X, '2. Random, CNDNUM = 2', 14X, '8. First column zero',
00574      $      / 4X, '3. Random, CNDNUM = sqrt(0.1/EPS)', 2X,
00575      $      '9. Last column zero', / 4X, '4. Random, CNDNUM = 0.1/EPS',
00576      $      7X, '10. Last n/2 columns zero', / 4X,
00577      $      '5. Scaled near underflow', 10X,
00578      $      '11. Scaled near underflow', / 4X,
00579      $      '6. Scaled near overflow', 11X, '12. Scaled near overflow' )
00580 *
00581 *     PT matrix types
00582 *
00583  9976 FORMAT( ' Matrix types (1-6 have specified condition numbers):',
00584      $      / 4X, '1. Diagonal', 24X, '7. Random, unspecified CNDNUM',
00585      $      / 4X, '2. Random, CNDNUM = 2', 14X,
00586      $      '8. First row and column zero', / 4X,
00587      $      '3. Random, CNDNUM = sqrt(0.1/EPS)', 2X,
00588      $      '9. Last row and column zero', / 4X,
00589      $      '4. Random, CNDNUM = 0.1/EPS', 7X,
00590      $      '10. Middle row and column zero', / 4X,
00591      $      '5. Scaled near underflow', 10X,
00592      $      '11. Scaled near underflow', / 4X,
00593      $      '6. Scaled near overflow', 11X, '12. Scaled near overflow' )
00594 *
00595 *     PO, PP matrix types
00596 *
00597  9975 FORMAT( 4X, '1. Diagonal', 24X,
00598      $      '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00599      $      '2. Random, CNDNUM = 2', 14X, '7. Random, CNDNUM = 0.1/EPS',
00600      $      / 3X, '*3. First row and column zero', 7X,
00601      $      '8. Scaled near underflow', / 3X,
00602      $      '*4. Last row and column zero', 8X,
00603      $      '9. Scaled near overflow', / 3X,
00604      $      '*5. Middle row and column zero', / 3X,
00605      $      '(* - tests error exits from ', A3,
00606      $      'TRF, no test ratios are computed)' )
00607 *
00608 *     CH matrix types
00609 *
00610  9974 FORMAT( 4X, '1. Diagonal', 24X,
00611      $      '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00612      $      '2. Random, CNDNUM = 2', 14X, '7. Random, CNDNUM = 0.1/EPS',
00613      $      / 3X, '*3. First row and column zero', 7X,
00614      $      '8. Scaled near underflow', / 3X,
00615      $      '*4. Last row and column zero', 8X,
00616      $      '9. Scaled near overflow', / 3X,
00617      $      '*5. Middle row and column zero', / 3X,
00618      $      '(* - tests error exits, no test ratios are computed)' )
00619 *
00620 *     PS matrix types
00621 *
00622  8973 FORMAT( 4X, '1. Diagonal', / 4X, '2. Random, CNDNUM = 2', 14X,
00623      $      / 3X, '*3. Nonzero eigenvalues of: D(1:RANK-1)=1 and ',
00624      $      'D(RANK) = 1.0/', A4, / 3X,
00625      $      '*4. Nonzero eigenvalues of: D(1)=1 and ',
00626      $      ' D(2:RANK) = 1.0/', A4, / 3X,
00627      $      '*5. Nonzero eigenvalues of: D(I) = ', A4,
00628      $      '**(-(I-1)/(RANK-1)) ', ' I=1:RANK', / 4X,
00629      $      '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00630      $      '7. Random, CNDNUM = 0.1/EPS', / 4X,
00631      $      '8. Scaled near underflow', / 4X, '9. Scaled near overflow',
00632      $      / 3X, '(* - Semi-definite tests )' )
00633  8972 FORMAT( 3X, 'RANK minus computed rank, returned by ', A, 'PSTRF' )
00634 *
00635 *     PB matrix types
00636 *
00637  9973 FORMAT( 4X, '1. Random, CNDNUM = 2', 14X,
00638      $      '5. Random, CNDNUM = sqrt(0.1/EPS)', / 3X,
00639      $      '*2. First row and column zero', 7X,
00640      $      '6. Random, CNDNUM = 0.1/EPS', / 3X,
00641      $      '*3. Last row and column zero', 8X,
00642      $      '7. Scaled near underflow', / 3X,
00643      $      '*4. Middle row and column zero', 6X,
00644      $      '8. Scaled near overflow', / 3X,
00645      $      '(* - tests error exits from ', A3,
00646      $      'TRF, no test ratios are computed)' )
00647 *
00648 *     SSY, SSP, CHE, CHP matrix types
00649 *
00650  9972 FORMAT( 4X, '1. Diagonal', 24X,
00651      $      '6. Last n/2 rows and columns zero', / 4X,
00652      $      '2. Random, CNDNUM = 2', 14X,
00653      $      '7. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00654      $      '3. First row and column zero', 7X,
00655      $      '8. Random, CNDNUM = 0.1/EPS', / 4X,
00656      $      '4. Last row and column zero', 8X,
00657      $      '9. Scaled near underflow', / 4X,
00658      $      '5. Middle row and column zero', 5X,
00659      $      '10. Scaled near overflow' )
00660 *
00661 *     CSY, CSP matrix types
00662 *
00663  9971 FORMAT( 4X, '1. Diagonal', 24X,
00664      $      '7. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00665      $      '2. Random, CNDNUM = 2', 14X, '8. Random, CNDNUM = 0.1/EPS',
00666      $      / 4X, '3. First row and column zero', 7X,
00667      $      '9. Scaled near underflow', / 4X,
00668      $      '4. Last row and column zero', 7X,
00669      $      '10. Scaled near overflow', / 4X,
00670      $      '5. Middle row and column zero', 5X,
00671      $      '11. Block diagonal matrix', / 4X,
00672      $      '6. Last n/2 rows and columns zero' )
00673 *
00674 *     QR matrix types
00675 *
00676  9970 FORMAT( 4X, '1. Diagonal', 24X,
00677      $      '5. Random, CNDNUM = sqrt(0.1/EPS)', / 4X,
00678      $      '2. Upper triangular', 16X, '6. Random, CNDNUM = 0.1/EPS',
00679      $      / 4X, '3. Lower triangular', 16X,
00680      $      '7. Scaled near underflow', / 4X, '4. Random, CNDNUM = 2',
00681      $      14X, '8. Scaled near overflow' )
00682 *
00683 *     QP matrix types
00684 *
00685  9969 FORMAT( ' Matrix types (2-6 have condition 1/EPS):', / 4X,
00686      $      '1. Zero matrix', 21X, '4. First n/2 columns fixed', / 4X,
00687      $      '2. One small eigenvalue', 12X, '5. Last n/2 columns fixed',
00688      $      / 4X, '3. Geometric distribution', 10X,
00689      $      '6. Every second column fixed' )
00690 *
00691 *     TZ matrix types
00692 *
00693  9968 FORMAT( ' Matrix types (2-3 have condition 1/EPS):', / 4X,
00694      $      '1. Zero matrix', / 4X, '2. One small eigenvalue', / 4X,
00695      $      '3. Geometric distribution' )
00696 *
00697 *     LS matrix types
00698 *
00699  9967 FORMAT( ' Matrix types (1-3: full rank, 4-6: rank deficient):',
00700      $      / 4X, '1 and 4. Normal scaling', / 4X,
00701      $      '2 and 5. Scaled near overflow', / 4X,
00702      $      '3 and 6. Scaled near underflow' )
00703 *
00704 *     TR, TP matrix types
00705 *
00706  9966 FORMAT( ' Matrix types for ', A3, ' routines:', / 4X,
00707      $      '1. Diagonal', 24X, '6. Scaled near overflow', / 4X,
00708      $      '2. Random, CNDNUM = 2', 14X, '7. Identity', / 4X,
00709      $      '3. Random, CNDNUM = sqrt(0.1/EPS)  ',
00710      $      '8. Unit triangular, CNDNUM = 2', / 4X,
00711      $      '4. Random, CNDNUM = 0.1/EPS', 8X,
00712      $      '9. Unit, CNDNUM = sqrt(0.1/EPS)', / 4X,
00713      $      '5. Scaled near underflow', 10X,
00714      $      '10. Unit, CNDNUM = 0.1/EPS' )
00715  9965 FORMAT( ' Special types for testing ', A, ':', / 3X,
00716      $      '11. Matrix elements are O(1), large right hand side', / 3X,
00717      $      '12. First diagonal causes overflow,',
00718      $      ' offdiagonal column norms < 1', / 3X,
00719      $      '13. First diagonal causes overflow,',
00720      $      ' offdiagonal column norms > 1', / 3X,
00721      $      '14. Growth factor underflows, solution does not overflow',
00722      $      / 3X, '15. Small diagonal causes gradual overflow', / 3X,
00723      $      '16. One zero diagonal element', / 3X,
00724      $      '17. Large offdiagonals cause overflow when adding a column'
00725      $      , / 3X, '18. Unit triangular with large right hand side' )
00726 *
00727 *     TB matrix types
00728 *
00729  9964 FORMAT( ' Matrix types for ', A3, ' routines:', / 4X,
00730      $      '1. Random, CNDNUM = 2', 14X, '6. Identity', / 4X,
00731      $      '2. Random, CNDNUM = sqrt(0.1/EPS)  ',
00732      $      '7. Unit triangular, CNDNUM = 2', / 4X,
00733      $      '3. Random, CNDNUM = 0.1/EPS', 8X,
00734      $      '8. Unit, CNDNUM = sqrt(0.1/EPS)', / 4X,
00735      $      '4. Scaled near underflow', 11X,
00736      $      '9. Unit, CNDNUM = 0.1/EPS', / 4X,
00737      $      '5. Scaled near overflow' )
00738  9963 FORMAT( ' Special types for testing ', A, ':', / 3X,
00739      $      '10. Matrix elements are O(1), large right hand side', / 3X,
00740      $      '11. First diagonal causes overflow,',
00741      $      ' offdiagonal column norms < 1', / 3X,
00742      $      '12. First diagonal causes overflow,',
00743      $      ' offdiagonal column norms > 1', / 3X,
00744      $      '13. Growth factor underflows, solution does not overflow',
00745      $      / 3X, '14. Small diagonal causes gradual overflow', / 3X,
00746      $      '15. One zero diagonal element', / 3X,
00747      $      '16. Large offdiagonals cause overflow when adding a column'
00748      $      , / 3X, '17. Unit triangular with large right hand side' )
00749 *
00750 *     Test ratios
00751 *
00752  9962 FORMAT( 3X, I2, ': norm( L * U - A )  / ( N * norm(A) * EPS )' )
00753  9961 FORMAT( 3X, I2, ': norm( I - A*AINV ) / ',
00754      $      '( N * norm(A) * norm(AINV) * EPS )' )
00755  9960 FORMAT( 3X, I2, ': norm( B - A * X )  / ',
00756      $      '( norm(A) * norm(X) * EPS )' )
00757  6660 FORMAT( 3X, I2, ': diagonal is not non-negative')
00758  9959 FORMAT( 3X, I2, ': norm( X - XACT )   / ',
00759      $      '( norm(XACT) * CNDNUM * EPS )' )
00760  9958 FORMAT( 3X, I2, ': norm( X - XACT )   / ',
00761      $      '( norm(XACT) * CNDNUM * EPS ), refined' )
00762  9957 FORMAT( 3X, I2, ': norm( X - XACT )   / ',
00763      $      '( norm(XACT) * (error bound) )' )
00764  9956 FORMAT( 3X, I2, ': (backward error)   / EPS' )
00765  9955 FORMAT( 3X, I2, ': RCOND * CNDNUM - 1.0' )
00766  9954 FORMAT( 3X, I2, ': norm( U'' * U - A ) / ( N * norm(A) * EPS )',
00767      $      ', or', / 7X, 'norm( L * L'' - A ) / ( N * norm(A) * EPS )'
00768      $       )
00769  8950 FORMAT( 3X,
00770      $      'norm( P * U'' * U * P'' - A ) / ( N * norm(A) * EPS )',
00771      $      ', or', / 3X,
00772      $      'norm( P * L * L'' * P'' - A ) / ( N * norm(A) * EPS )' )
00773  9953 FORMAT( 3X, I2, ': norm( U*D*U'' - A ) / ( N * norm(A) * EPS )',
00774      $      ', or', / 7X, 'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
00775      $       )
00776  9952 FORMAT( 3X, I2, ': norm( U''*D*U - A ) / ( N * norm(A) * EPS )',
00777      $      ', or', / 7X, 'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
00778      $       )
00779  9951 FORMAT( ' Test ratio for ', A, ':', / 3X, I2,
00780      $      ': norm( s*b - A*x )  / ( norm(A) * norm(x) * EPS )' )
00781  9950 FORMAT( 3X, I2, ': norm( R - Q'' * A ) / ( M * norm(A) * EPS )' )
00782  6950 FORMAT( 3X, I2, ': norm( R - Q'
00783 ' * A ) / ( M * norm(A) * EPS )     $       [RFPG]' )
00784  9949 FORMAT( 3X, I2, ': norm( L - A * Q'' ) / ( N * norm(A) * EPS )' )
00785  9948 FORMAT( 3X, I2, ': norm( L - Q'' * A ) / ( M * norm(A) * EPS )' )
00786  9947 FORMAT( 3X, I2, ': norm( R - A * Q'' ) / ( N * norm(A) * EPS )' )
00787  9946 FORMAT( 3X, I2, ': norm( I - Q''*Q )   / ( M * EPS )' )
00788  9945 FORMAT( 3X, I2, ': norm( I - Q*Q'' )   / ( N * EPS )' )
00789  9944 FORMAT( 3X, I2, ': norm( Q*C - Q*C )  / ', '( ', A1,
00790      $      ' * norm(C) * EPS )' )
00791  9943 FORMAT( 3X, I2, ': norm( C*Q - C*Q )  / ', '( ', A1,
00792      $      ' * norm(C) * EPS )' )
00793  9942 FORMAT( 3X, I2, ': norm( Q''*C - Q''*C )/ ', '( ', A1,
00794      $      ' * norm(C) * EPS )' )
00795  9941 FORMAT( 3X, I2, ': norm( C*Q'' - C*Q'' )/ ', '( ', A1,
00796      $      ' * norm(C) * EPS )' )
00797  9940 FORMAT( 3X, I2, ': norm(svd(A) - svd(R)) / ',
00798      $      '( M * norm(svd(R)) * EPS )' )
00799  9939 FORMAT( 3X, I2, ': norm( A*P - Q*R )     / ( M * norm(A) * EPS )'
00800      $       )
00801  9938 FORMAT( 3X, I2, ': norm( I - Q''*Q )      / ( M * EPS )' )
00802  9937 FORMAT( 3X, I2, ': norm( A - R*Q )       / ( M * norm(A) * EPS )'
00803      $       )
00804  9936 FORMAT( ' Test ratios (1-2: ', A1, 'GELS, 3-6: ', A1,
00805      $      'GELSS, 7-10: ', A1, 'GELSX):' )
00806  9935 FORMAT( 3X, I2, ': norm( B - A * X )   / ',
00807      $      '( max(M,N) * norm(A) * norm(X) * EPS )' )
00808  9934 FORMAT( 3X, I2, ': norm( (A*X-B)'' *A ) / ',
00809      $      '( max(M,N,NRHS) * norm(A) * norm(B) * EPS )' )
00810  9933 FORMAT( 3X, I2, ': norm(svd(A)-svd(R)) / ',
00811      $      '( min(M,N) * norm(svd(R)) * EPS )' )
00812  9932 FORMAT( 3X, I2, ': Check if X is in the row space of A or A''' )
00813  9931 FORMAT( 3X, I2, ': norm( (A*X-B)'' *A ) / ',
00814      $      '( max(M,N,NRHS) * norm(A) * norm(B) * EPS )', / 7X,
00815      $      'if TRANS=''N'' and M.GE.N or TRANS=''T'' and M.LT.N, ',
00816      $      'otherwise', / 7X,
00817      $      'check if X is in the row space of A or A'' ',
00818      $      '(overdetermined case)' )
00819  9930 FORMAT( 3X, ' 7-10: same as 3-6' )
00820  9929 FORMAT( ' Test ratios (1-3: ', A1, 'TZRQF, 4-6: ', A1,
00821      $      'TZRZF):' )
00822  9920 FORMAT( 3X, ' 7-10: same as 3-6', 3X, ' 11-14: same as 3-6',
00823      $      3X, ' 15-18: same as 3-6' )
00824  9921 FORMAT( ' Test ratios:', / '    (1-2: ', A1, 'GELS, 3-6: ', A1,
00825      $      'GELSX, 7-10: ', A1, 'GELSY, 11-14: ', A1, 'GELSS, 15-18: ',
00826      $      A1, 'GELSD)' )
00827 *
00828       RETURN
00829 *
00830 *     End of ALAHD
00831 *
00832       END