alahd.f

Go to the documentation of this file.

      SUBROUTINE alahd( IOUNIT, PATH )
*     
*     -- LAPACK test routine (version 3.1) --
*     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
*     November 2006
*     
*     .. Scalar Arguments ..
      CHARACTER*3        path
      INTEGER            iounit
*     ..
*     
*     Purpose
*     =======
*     
*     ALAHD prints header information for the different test paths.
*     
*     Arguments
*     =========
*     
*     IOUNIT  (input) INTEGER
*     The unit number to which the header information should be
*     printed.
*     
*     PATH    (input) CHARACTER*3
*     The name of the path for which the header information is to
*     be printed.  Current paths are
*     _GE:  General matrices
*     _GB:  General band
*     _GT:  General Tridiagonal
*     _PO:  Symmetric or Hermitian positive definite
*     _PS:  Symmetric or Hermitian positive semi-definite
*     _PP:  Symmetric or Hermitian positive definite packed
*     _PB:  Symmetric or Hermitian positive definite band
*     _PT:  Symmetric or Hermitian positive definite tridiagonal
*     _SY:  Symmetric indefinite
*     _SP:  Symmetric indefinite packed
*     _HE:  (complex) Hermitian indefinite
*     _HP:  (complex) Hermitian indefinite packed
*     _TR:  Triangular
*     _TP:  Triangular packed
*     _TB:  Triangular band
*     _QR:  QR (general matrices)
*     _LQ:  LQ (general matrices)
*     _QL:  QL (general matrices)
*     _RQ:  RQ (general matrices)
*     _QP:  QR with column pivoting
*     _TZ:  Trapezoidal
*     _LS:  Least Squares driver routines
*     _LU:  LU variants
*     _CH:  Cholesky variants
*     _QS:  QR variants
*     The first character must be one of S, D, C, or Z (C or Z only
*     if complex).
*     
*     =====================================================================
*     
*     .. Local Scalars ..
      LOGICAL            corz, sord
      CHARACTER          c1, c3
      CHARACTER*2        p2
      CHARACTER*4        eigcnm
      CHARACTER*32       subnam
      CHARACTER*9        sym
*     ..
*     .. External Functions ..
      LOGICAL            lsame, lsamen
      EXTERNAL           lsame, lsamen
*     ..
*     .. Intrinsic Functions ..
      INTRINSIC          len_trim
*     ..
*     .. Executable Statements ..
*     
      IF( iounit.LE.0 )
     $     return
      c1 = path( 1: 1 )
      c3 = path( 3: 3 )
      p2 = path( 2: 3 )
      sord = lsame( c1, 'S' ) .OR. lsame( c1, 'D' )
      corz = lsame( c1, 'C' ) .OR. lsame( c1, 'Z' )
      IF( .NOT.( sord .OR. corz ) )
     $     return
*     
      IF( lsamen( 2, p2, 'GE' ) ) THEN
*     
*     GE: General dense
*     
         WRITE( iounit, fmt = 9999 )path
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9979 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9962 )1
         WRITE( iounit, fmt = 9961 )2
         WRITE( iounit, fmt = 9960 )3
         WRITE( iounit, fmt = 9959 )4
         WRITE( iounit, fmt = 9958 )5
         WRITE( iounit, fmt = 9957 )6
         WRITE( iounit, fmt = 9956 )7
         WRITE( iounit, fmt = 9955 )8
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'GB' ) ) THEN
*     
*     GB: General band
*     
         WRITE( iounit, fmt = 9998 )path
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9978 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9962 )1
         WRITE( iounit, fmt = 9960 )2
         WRITE( iounit, fmt = 9959 )3
         WRITE( iounit, fmt = 9958 )4
         WRITE( iounit, fmt = 9957 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9955 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'GT' ) ) THEN
*     
*     GT: General tridiagonal
*     
         WRITE( iounit, fmt = 9997 )path
         WRITE( iounit, fmt = 9977 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9962 )1
         WRITE( iounit, fmt = 9960 )2
         WRITE( iounit, fmt = 9959 )3
         WRITE( iounit, fmt = 9958 )4
         WRITE( iounit, fmt = 9957 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9955 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'PO' ) .OR. lsamen( 2, p2, 'PP' ) ) THEN
*     
*     PO: Positive definite full
*     PP: Positive definite packed
*     
         IF( sord ) THEN
            sym = 'Symmetric'
         ELSE
            sym = 'Hermitian'
         END IF
         IF( lsame( c3, 'O' ) ) THEN
            WRITE( iounit, fmt = 9996 )path, sym
         ELSE
            WRITE( iounit, fmt = 9995 )path, sym
         END IF
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9975 )path
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9954 )1
         WRITE( iounit, fmt = 9961 )2
         WRITE( iounit, fmt = 9960 )3
         WRITE( iounit, fmt = 9959 )4
         WRITE( iounit, fmt = 9958 )5
         WRITE( iounit, fmt = 9957 )6
         WRITE( iounit, fmt = 9956 )7
         WRITE( iounit, fmt = 9955 )8
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'PS' ) ) THEN
*     
*     PS: Positive semi-definite full
*     
         IF( sord ) THEN
            sym = 'Symmetric'
         ELSE
            sym = 'Hermitian'
         END IF
         IF( lsame( c1, 'S' ) .OR. lsame( c1, 'C' ) ) THEN
            eigcnm = '1E04'
         ELSE
            eigcnm = '1D12'
         END IF
         WRITE( iounit, fmt = 9995 )path, sym
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 8973 )eigcnm, eigcnm, eigcnm
         WRITE( iounit, fmt = '( '' Difference:'' )' )
         WRITE( iounit, fmt = 8972 )c1
         WRITE( iounit, fmt = '( '' Test ratio:'' )' )
         WRITE( iounit, fmt = 8950 )
         WRITE( iounit, fmt = '( '' Messages:'' )' )
      ELSE IF( lsamen( 2, p2, 'PB' ) ) THEN
*     
*     PB: Positive definite band
*     
         IF( sord ) THEN
            WRITE( iounit, fmt = 9994 )path, 'Symmetric'
         ELSE
            WRITE( iounit, fmt = 9994 )path, 'Hermitian'
         END IF
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9973 )path
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9954 )1
         WRITE( iounit, fmt = 9960 )2
         WRITE( iounit, fmt = 9959 )3
         WRITE( iounit, fmt = 9958 )4
         WRITE( iounit, fmt = 9957 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9955 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'PT' ) ) THEN
*     
*     PT: Positive definite tridiagonal
*     
         IF( sord ) THEN
            WRITE( iounit, fmt = 9993 )path, 'Symmetric'
         ELSE
            WRITE( iounit, fmt = 9993 )path, 'Hermitian'
         END IF
         WRITE( iounit, fmt = 9976 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9952 )1
         WRITE( iounit, fmt = 9960 )2
         WRITE( iounit, fmt = 9959 )3
         WRITE( iounit, fmt = 9958 )4
         WRITE( iounit, fmt = 9957 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9955 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'SY' ) .OR. lsamen( 2, p2, 'SP' ) ) THEN
*     
*     SY: Symmetric indefinite full
*     SP: Symmetric indefinite packed
*     
         IF( lsame( c3, 'Y' ) ) THEN
            WRITE( iounit, fmt = 9992 )path, 'Symmetric'
         ELSE
            WRITE( iounit, fmt = 9991 )path, 'Symmetric'
         END IF
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         IF( sord ) THEN
            WRITE( iounit, fmt = 9972 )
         ELSE
            WRITE( iounit, fmt = 9971 )
         END IF
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9953 )1
         WRITE( iounit, fmt = 9961 )2
         WRITE( iounit, fmt = 9960 )3
         WRITE( iounit, fmt = 9959 )4
         WRITE( iounit, fmt = 9958 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9957 )7
         WRITE( iounit, fmt = 9955 )8
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'HE' ) .OR. lsamen( 2, p2, 'HP' ) ) THEN
*     
*     HE: Hermitian indefinite full
*     HP: Hermitian indefinite packed
*     
         IF( lsame( c3, 'E' ) ) THEN
            WRITE( iounit, fmt = 9992 )path, 'Hermitian'
         ELSE
            WRITE( iounit, fmt = 9991 )path, 'Hermitian'
         END IF
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9972 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9953 )1
         WRITE( iounit, fmt = 9961 )2
         WRITE( iounit, fmt = 9960 )3
         WRITE( iounit, fmt = 9959 )4
         WRITE( iounit, fmt = 9958 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9957 )7
         WRITE( iounit, fmt = 9955 )8
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'TR' ) .OR. lsamen( 2, p2, 'TP' ) ) THEN
*     
*     TR: Triangular full
*     TP: Triangular packed
*     
         IF( lsame( c3, 'R' ) ) THEN
            WRITE( iounit, fmt = 9990 )path
            subnam = path( 1: 1 ) // 'LATRS'
         ELSE
            WRITE( iounit, fmt = 9989 )path
            subnam = path( 1: 1 ) // 'LATPS'
         END IF
         WRITE( iounit, fmt = 9966 )path
         WRITE( iounit, fmt = 9965 )subnam(1:len_trim( subnam ))
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9961 )1
         WRITE( iounit, fmt = 9960 )2
         WRITE( iounit, fmt = 9959 )3
         WRITE( iounit, fmt = 9958 )4
         WRITE( iounit, fmt = 9957 )5
         WRITE( iounit, fmt = 9956 )6
         WRITE( iounit, fmt = 9955 )7
         WRITE( iounit, fmt = 9951 )subnam(1:len_trim( subnam )), 8
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'TB' ) ) THEN
*     
*     TB: Triangular band
*     
         WRITE( iounit, fmt = 9988 )path
         subnam = path( 1: 1 ) // 'LATBS'
         WRITE( iounit, fmt = 9964 )path
         WRITE( iounit, fmt = 9963 )subnam(1:len_trim( subnam ))
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9960 )1
         WRITE( iounit, fmt = 9959 )2
         WRITE( iounit, fmt = 9958 )3
         WRITE( iounit, fmt = 9957 )4
         WRITE( iounit, fmt = 9956 )5
         WRITE( iounit, fmt = 9955 )6
         WRITE( iounit, fmt = 9951 )subnam(1:len_trim( subnam )), 7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'QR' ) ) THEN
*     
*     QR decomposition of rectangular matrices
*     
         WRITE( iounit, fmt = 9987 )path, 'QR'
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9970 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9950 )1
         WRITE( iounit, fmt = 9946 )2
         WRITE( iounit, fmt = 9944 )3, 'M'
         WRITE( iounit, fmt = 9943 )4, 'M'
         WRITE( iounit, fmt = 9942 )5, 'M'
         WRITE( iounit, fmt = 9941 )6, 'M'
         WRITE( iounit, fmt = 9960 )7
         WRITE( iounit, fmt = 6660 )8
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'LQ' ) ) THEN
*     
*     LQ decomposition of rectangular matrices
*     
         WRITE( iounit, fmt = 9987 )path, 'LQ'
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9970 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9949 )1
         WRITE( iounit, fmt = 9945 )2
         WRITE( iounit, fmt = 9944 )3, 'N'
         WRITE( iounit, fmt = 9943 )4, 'N'
         WRITE( iounit, fmt = 9942 )5, 'N'
         WRITE( iounit, fmt = 9941 )6, 'N'
         WRITE( iounit, fmt = 9960 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'QL' ) ) THEN
*     
*     QL decomposition of rectangular matrices
*     
         WRITE( iounit, fmt = 9987 )path, 'QL'
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9970 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9948 )1
         WRITE( iounit, fmt = 9946 )2
         WRITE( iounit, fmt = 9944 )3, 'M'
         WRITE( iounit, fmt = 9943 )4, 'M'
         WRITE( iounit, fmt = 9942 )5, 'M'
         WRITE( iounit, fmt = 9941 )6, 'M'
         WRITE( iounit, fmt = 9960 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'RQ' ) ) THEN
*     
*     RQ decomposition of rectangular matrices
*     
         WRITE( iounit, fmt = 9987 )path, 'RQ'
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9970 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9947 )1
         WRITE( iounit, fmt = 9945 )2
         WRITE( iounit, fmt = 9944 )3, 'N'
         WRITE( iounit, fmt = 9943 )4, 'N'
         WRITE( iounit, fmt = 9942 )5, 'N'
         WRITE( iounit, fmt = 9941 )6, 'N'
         WRITE( iounit, fmt = 9960 )7
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'QP' ) ) THEN
*     
*     QR decomposition with column pivoting
*     
         WRITE( iounit, fmt = 9986 )path
         WRITE( iounit, fmt = 9969 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9940 )1
         WRITE( iounit, fmt = 9939 )2
         WRITE( iounit, fmt = 9938 )3
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'TZ' ) ) THEN
*     
*     TZ:  Trapezoidal
*     
         WRITE( iounit, fmt = 9985 )path
         WRITE( iounit, fmt = 9968 )
         WRITE( iounit, fmt = 9929 )c1, c1
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
         WRITE( iounit, fmt = 9940 )1
         WRITE( iounit, fmt = 9937 )2
         WRITE( iounit, fmt = 9938 )3
         WRITE( iounit, fmt = 9940 )4
         WRITE( iounit, fmt = 9937 )5
         WRITE( iounit, fmt = 9938 )6
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'LS' ) ) THEN
*     
*     LS:  Least Squares driver routines for
*     LS, LSD, LSS, LSX and LSY.
*     
         WRITE( iounit, fmt = 9984 )path
         WRITE( iounit, fmt = 9967 )
         WRITE( iounit, fmt = 9921 )c1, c1, c1, c1, c1
         WRITE( iounit, fmt = 9935 )1
         WRITE( iounit, fmt = 9931 )2
         WRITE( iounit, fmt = 9933 )3
         WRITE( iounit, fmt = 9935 )4
         WRITE( iounit, fmt = 9934 )5
         WRITE( iounit, fmt = 9932 )6
         WRITE( iounit, fmt = 9920 )
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'LU' ) ) THEN
*     
*     LU factorization variants
*     
         WRITE( iounit, fmt = 9983 )path
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9979 )
         WRITE( iounit, fmt = '( '' Test ratio:'' )' )
         WRITE( iounit, fmt = 9962 )1
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'CH' ) ) THEN
*     
*     Cholesky factorization variants
*     
         WRITE( iounit, fmt = 9982 )path
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9974 )
         WRITE( iounit, fmt = '( '' Test ratio:'' )' )
         WRITE( iounit, fmt = 9954 )1
         WRITE( iounit, fmt = '( '' Messages:'' )' )
*     
      ELSE IF( lsamen( 2, p2, 'QS' ) ) THEN
*     
*     QR factorization variants
*     
         WRITE( iounit, fmt = 9981 )path
         WRITE( iounit, fmt = '( '' Matrix types:'' )' )
         WRITE( iounit, fmt = 9970 )
         WRITE( iounit, fmt = '( '' Test ratios:'' )' )
*     
      ELSE
*     
*     Print error message if no header is available.
*     
         WRITE( iounit, fmt = 9980 )path
      END IF
*     
*     First line of header
*     
 9999 format( / 1x, a3, ':  General dense matrices' )
 9998 format( / 1x, a3, ':  General band matrices' )
 9997 format( / 1x, a3, ':  General tridiagonal' )
 9996 format( / 1x, a3, ':  ', a9, ' positive definite matrices' )
 9995 format( / 1x, a3, ':  ', a9, ' positive definite packed matrices'
     $     )
 9994 format( / 1x, a3, ':  ', a9, ' positive definite band matrices' )
 9993 format( / 1x, a3, ':  ', a9, ' positive definite tridiagonal' )
 9992 format( / 1x, a3, ':  ', a9, ' indefinite matrices' )
 9991 format( / 1x, a3, ':  ', a9, ' indefinite packed matrices' )
 9990 format( / 1x, a3, ':  Triangular matrices' )
 9989 format( / 1x, a3, ':  Triangular packed matrices' )
 9988 format( / 1x, a3, ':  Triangular band matrices' )
 9987 format( / 1x, a3, ':  ', a2, ' factorization of general matrices'
     $     )
 9986 format( / 1x, a3, ':  QR factorization with column pivoting' )
 9985 format( / 1x, a3, ':  RQ factorization of trapezoidal matrix' )
 9984 format( / 1x, a3, ':  Least squares driver routines' )
 9983 format( / 1x, a3, ':  LU factorization variants' )
 9982 format( / 1x, a3, ':  Cholesky factorization variants' )
 9981 format( / 1x, a3, ':  QR factorization variants' )
 9980 format( / 1x, a3, ':  No header available' )
*     
*     GE matrix types
*     
 9979 format( 4x, '1. Diagonal', 24x, '7. Last n/2 columns zero', / 4x,
     $     '2. Upper triangular', 16x,
     $     '8. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '3. Lower triangular', 16x, '9. Random, CNDNUM = 0.1/EPS',
     $     / 4x, '4. Random, CNDNUM = 2', 13x,
     $     '10. Scaled near underflow', / 4x, '5. First column zero',
     $     14x, '11. Scaled near overflow', / 4x,
     $     '6. Last column zero' )
*     
*     GB matrix types
*     
 9978 format( 4x, '1. Random, CNDNUM = 2', 14x,
     $     '5. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '2. First column zero', 15x, '6. Random, CNDNUM = .01/EPS',
     $     / 4x, '3. Last column zero', 16x,
     $     '7. Scaled near underflow', / 4x,
     $     '4. Last n/2 columns zero', 11x, '8. Scaled near overflow' )
*     
*     GT matrix types
*     
 9977 format( ' Matrix types (1-6 have specified condition numbers):',
     $     / 4x, '1. Diagonal', 24x, '7. Random, unspecified CNDNUM',
     $     / 4x, '2. Random, CNDNUM = 2', 14x, '8. First column zero',
     $     / 4x, '3. Random, CNDNUM = sqrt(0.1/EPS)', 2x,
     $     '9. Last column zero', / 4x, '4. Random, CNDNUM = 0.1/EPS',
     $     7x, '10. Last n/2 columns zero', / 4x,
     $     '5. Scaled near underflow', 10x,
     $     '11. Scaled near underflow', / 4x,
     $     '6. Scaled near overflow', 11x, '12. Scaled near overflow' )
*     
*     PT matrix types
*     
 9976 format( ' Matrix types (1-6 have specified condition numbers):',
     $     / 4x, '1. Diagonal', 24x, '7. Random, unspecified CNDNUM',
     $     / 4x, '2. Random, CNDNUM = 2', 14x,
     $     '8. First row and column zero', / 4x,
     $     '3. Random, CNDNUM = sqrt(0.1/EPS)', 2x,
     $     '9. Last row and column zero', / 4x,
     $     '4. Random, CNDNUM = 0.1/EPS', 7x,
     $     '10. Middle row and column zero', / 4x,
     $     '5. Scaled near underflow', 10x,
     $     '11. Scaled near underflow', / 4x,
     $     '6. Scaled near overflow', 11x, '12. Scaled near overflow' )
*     
*     PO, PP matrix types
*     
 9975 format( 4x, '1. Diagonal', 24x,
     $     '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '2. Random, CNDNUM = 2', 14x, '7. Random, CNDNUM = 0.1/EPS',
     $     / 3x, '*3. First row and column zero', 7x,
     $     '8. Scaled near underflow', / 3x,
     $     '*4. Last row and column zero', 8x,
     $     '9. Scaled near overflow', / 3x,
     $     '*5. Middle row and column zero', / 3x,
     $     '(* - tests error exits from ', a3,
     $     'TRF, no test ratios are computed)' )
*     
*     CH matrix types
*     
 9974 format( 4x, '1. Diagonal', 24x,
     $     '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '2. Random, CNDNUM = 2', 14x, '7. Random, CNDNUM = 0.1/EPS',
     $     / 3x, '*3. First row and column zero', 7x,
     $     '8. Scaled near underflow', / 3x,
     $     '*4. Last row and column zero', 8x,
     $     '9. Scaled near overflow', / 3x,
     $     '*5. Middle row and column zero', / 3x,
     $     '(* - tests error exits, no test ratios are computed)' )
*     
*     PS matrix types
*     
 8973 format( 4x, '1. Diagonal', / 4x, '2. Random, CNDNUM = 2', 14x,
     $     / 3x, '*3. Nonzero eigenvalues of: D(1:RANK-1)=1 and ',
     $     'D(RANK) = 1.0/', a4, / 3x,
     $     '*4. Nonzero eigenvalues of: D(1)=1 and ',
     $     ' D(2:RANK) = 1.0/', a4, / 3x,
     $     '*5. Nonzero eigenvalues of: D(I) = ', a4,
     $     '**(-(I-1)/(RANK-1)) ', ' I=1:RANK', / 4x,
     $     '6. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '7. Random, CNDNUM = 0.1/EPS', / 4x,
     $     '8. Scaled near underflow', / 4x, '9. Scaled near overflow',
     $     / 3x, '(* - Semi-definite tests )' )
 8972 format( 3x, 'RANK minus computed rank, returned by ', a, 'PSTRF' )
*     
*     PB matrix types
*     
 9973 format( 4x, '1. Random, CNDNUM = 2', 14x,
     $     '5. Random, CNDNUM = sqrt(0.1/EPS)', / 3x,
     $     '*2. First row and column zero', 7x,
     $     '6. Random, CNDNUM = 0.1/EPS', / 3x,
     $     '*3. Last row and column zero', 8x,
     $     '7. Scaled near underflow', / 3x,
     $     '*4. Middle row and column zero', 6x,
     $     '8. Scaled near overflow', / 3x,
     $     '(* - tests error exits from ', a3,
     $     'TRF, no test ratios are computed)' )
*     
*     SSY, SSP, CHE, CHP matrix types
*     
 9972 format( 4x, '1. Diagonal', 24x,
     $     '6. Last n/2 rows and columns zero', / 4x,
     $     '2. Random, CNDNUM = 2', 14x,
     $     '7. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '3. First row and column zero', 7x,
     $     '8. Random, CNDNUM = 0.1/EPS', / 4x,
     $     '4. Last row and column zero', 8x,
     $     '9. Scaled near underflow', / 4x,
     $     '5. Middle row and column zero', 5x,
     $     '10. Scaled near overflow' )
*     
*     CSY, CSP matrix types
*     
 9971 format( 4x, '1. Diagonal', 24x,
     $     '7. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '2. Random, CNDNUM = 2', 14x, '8. Random, CNDNUM = 0.1/EPS',
     $     / 4x, '3. First row and column zero', 7x,
     $     '9. Scaled near underflow', / 4x,
     $     '4. Last row and column zero', 7x,
     $     '10. Scaled near overflow', / 4x,
     $     '5. Middle row and column zero', 5x,
     $     '11. Block diagonal matrix', / 4x,
     $     '6. Last n/2 rows and columns zero' )
*     
*     QR matrix types
*     
 9970 format( 4x, '1. Diagonal', 24x,
     $     '5. Random, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '2. Upper triangular', 16x, '6. Random, CNDNUM = 0.1/EPS',
     $     / 4x, '3. Lower triangular', 16x,
     $     '7. Scaled near underflow', / 4x, '4. Random, CNDNUM = 2',
     $     14x, '8. Scaled near overflow' )
*     
*     QP matrix types
*     
 9969 format( ' Matrix types (2-6 have condition 1/EPS):', / 4x,
     $     '1. Zero matrix', 21x, '4. First n/2 columns fixed', / 4x,
     $     '2. One small eigenvalue', 12x, '5. Last n/2 columns fixed',
     $     / 4x, '3. Geometric distribution', 10x,
     $     '6. Every second column fixed' )
*     
*     TZ matrix types
*     
 9968 format( ' Matrix types (2-3 have condition 1/EPS):', / 4x,
     $     '1. Zero matrix', / 4x, '2. One small eigenvalue', / 4x,
     $     '3. Geometric distribution' )
*     
*     LS matrix types
*     
 9967 format( ' Matrix types (1-3: full rank, 4-6: rank deficient):',
     $     / 4x, '1 and 4. Normal scaling', / 4x,
     $     '2 and 5. Scaled near overflow', / 4x,
     $     '3 and 6. Scaled near underflow' )
*     
*     TR, TP matrix types
*     
 9966 format( ' Matrix types for ', a3, ' routines:', / 4x,
     $     '1. Diagonal', 24x, '6. Scaled near overflow', / 4x,
     $     '2. Random, CNDNUM = 2', 14x, '7. Identity', / 4x,
     $     '3. Random, CNDNUM = sqrt(0.1/EPS)  ',
     $     '8. Unit triangular, CNDNUM = 2', / 4x,
     $     '4. Random, CNDNUM = 0.1/EPS', 8x,
     $     '9. Unit, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '5. Scaled near underflow', 10x,
     $     '10. Unit, CNDNUM = 0.1/EPS' )
 9965 format( ' Special types for testing ', a, ':', / 3x,
     $     '11. Matrix elements are O(1), large right hand side', / 3x,
     $     '12. First diagonal causes overflow,',
     $     ' offdiagonal column norms < 1', / 3x,
     $     '13. First diagonal causes overflow,',
     $     ' offdiagonal column norms > 1', / 3x,
     $     '14. Growth factor underflows, solution does not overflow',
     $     / 3x, '15. Small diagonal causes gradual overflow', / 3x,
     $     '16. One zero diagonal element', / 3x,
     $     '17. Large offdiagonals cause overflow when adding a column'
     $     , / 3x, '18. Unit triangular with large right hand side' )
*     
*     TB matrix types
*     
 9964 format( ' Matrix types for ', a3, ' routines:', / 4x,
     $     '1. Random, CNDNUM = 2', 14x, '6. Identity', / 4x,
     $     '2. Random, CNDNUM = sqrt(0.1/EPS)  ',
     $     '7. Unit triangular, CNDNUM = 2', / 4x,
     $     '3. Random, CNDNUM = 0.1/EPS', 8x,
     $     '8. Unit, CNDNUM = sqrt(0.1/EPS)', / 4x,
     $     '4. Scaled near underflow', 11x,
     $     '9. Unit, CNDNUM = 0.1/EPS', / 4x,
     $     '5. Scaled near overflow' )
 9963 format( ' Special types for testing ', a, ':', / 3x,
     $     '10. Matrix elements are O(1), large right hand side', / 3x,
     $     '11. First diagonal causes overflow,',
     $     ' offdiagonal column norms < 1', / 3x,
     $     '12. First diagonal causes overflow,',
     $     ' offdiagonal column norms > 1', / 3x,
     $     '13. Growth factor underflows, solution does not overflow',
     $     / 3x, '14. Small diagonal causes gradual overflow', / 3x,
     $     '15. One zero diagonal element', / 3x,
     $     '16. Large offdiagonals cause overflow when adding a column'
     $     , / 3x, '17. Unit triangular with large right hand side' )
*     
*     Test ratios
*     
 9962 format( 3x, i2, ': norm( L * U - A )  / ( N * norm(A) * EPS )' )
 9961 format( 3x, i2, ': norm( I - A*AINV ) / ',
     $     '( N * norm(A) * norm(AINV) * EPS )' )
 9960 format( 3x, i2, ': norm( B - A * X )  / ',
     $     '( norm(A) * norm(X) * EPS )' )
 6660 format( 3x, i2, ': diagonal is not non-negative')
 9959 format( 3x, i2, ': norm( X - XACT )   / ',
     $     '( norm(XACT) * CNDNUM * EPS )' )
 9958 format( 3x, i2, ': norm( X - XACT )   / ',
     $     '( norm(XACT) * CNDNUM * EPS ), refined' )
 9957 format( 3x, i2, ': norm( X - XACT )   / ',
     $     '( norm(XACT) * (error bound) )' )
 9956 format( 3x, i2, ': (backward error)   / EPS' )
 9955 format( 3x, i2, ': RCOND * CNDNUM - 1.0' )
 9954 format( 3x, i2, ': norm( U'' * U - A ) / ( N * norm(A) * EPS )',
     $     ', or', / 7x, 'norm( L * L'' - A ) / ( N * norm(A) * EPS )'
     $     )
 8950 format( 3x,
     $     'norm( P * U'' * U * P'' - A ) / ( N * norm(A) * EPS )',
     $     ', or', / 3x,
     $     'norm( P * L * L'' * P'' - A ) / ( N * norm(A) * EPS )' )
 9953 format( 3x, i2, ': norm( U*D*U'' - A ) / ( N * norm(A) * EPS )',
     $     ', or', / 7x, 'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
     $     )
 9952 format( 3x, i2, ': norm( U''*D*U - A ) / ( N * norm(A) * EPS )',
     $     ', or', / 7x, 'norm( L*D*L'' - A ) / ( N * norm(A) * EPS )'
     $     )
 9951 format( ' Test ratio for ', a, ':', / 3x, i2,
     $     ': norm( s*b - A*x )  / ( norm(A) * norm(x) * EPS )' )
 9950 format( 3x, i2, ': norm( R - Q'' * A ) / ( M * norm(A) * EPS )' )
 9949 format( 3x, i2, ': norm( L - A * Q'' ) / ( N * norm(A) * EPS )' )
 9948 format( 3x, i2, ': norm( L - Q'' * A ) / ( M * norm(A) * EPS )' )
 9947 format( 3x, i2, ': norm( R - A * Q'' ) / ( N * norm(A) * EPS )' )
 9946 format( 3x, i2, ': norm( I - Q''*Q )   / ( M * EPS )' )
 9945 format( 3x, i2, ': norm( I - Q*Q'' )   / ( N * EPS )' )
 9944 format( 3x, i2, ': norm( Q*C - Q*C )  / ', '( ', a1,
     $     ' * norm(C) * EPS )' )
 9943 format( 3x, i2, ': norm( C*Q - C*Q )  / ', '( ', a1,
     $     ' * norm(C) * EPS )' )
 9942 format( 3x, i2, ': norm( Q''*C - Q''*C )/ ', '( ', a1,
     $     ' * norm(C) * EPS )' )
 9941 format( 3x, i2, ': norm( C*Q'' - C*Q'' )/ ', '( ', a1,
     $     ' * norm(C) * EPS )' )
 9940 format( 3x, i2, ': norm(svd(A) - svd(R)) / ',
     $     '( M * norm(svd(R)) * EPS )' )
 9939 format( 3x, i2, ': norm( A*P - Q*R )     / ( M * norm(A) * EPS )'
     $     )
 9938 format( 3x, i2, ': norm( I - Q''*Q )      / ( M * EPS )' )
 9937 format( 3x, i2, ': norm( A - R*Q )       / ( M * norm(A) * EPS )'
     $     )
 9936 format( ' Test ratios (1-2: ', a1, 'GELS, 3-6: ', a1,
     $     'GELSS, 7-10: ', a1, 'GELSX):' )
 9935 format( 3x, i2, ': norm( B - A * X )   / ',
     $     '( max(M,N) * norm(A) * norm(X) * EPS )' )
 9934 format( 3x, i2, ': norm( (A*X-B)'' *A ) / ',
     $     '( max(M,N,NRHS) * norm(A) * norm(B) * EPS )' )
 9933 format( 3x, i2, ': norm(svd(A)-svd(R)) / ',
     $     '( min(M,N) * norm(svd(R)) * EPS )' )
 9932 format( 3x, i2, ': Check if X is in the row space of A or A''' )
 9931 format( 3x, i2, ': norm( (A*X-B)'' *A ) / ',
     $     '( max(M,N,NRHS) * norm(A) * norm(B) * EPS )', / 7x,
     $     'if TRANS=''N'.GE.' and MN or TRANS=''T'.LT.' and MN, ',
     $     'otherwise', / 7x,
     $     'check if X is in the row space of A or A'' ',
     $     '(overdetermined case)' )
 9930 format( 3x, ' 7-10: same as 3-6' )
 9929 format( ' Test ratios (1-3: ', a1, 'TZRQF, 4-6: ', a1,
     $     'TZRZF):' )
 9920 format( 3x, ' 7-10: same as 3-6', 3x, ' 11-14: same as 3-6',
     $     3x, ' 15-18: same as 3-6' )
 9921 format( ' Test ratios:', / '    (1-2: ', a1, 'GELS, 3-6: ', a1,
     $     'GELSX, 7-10: ', a1, 'GELSY, 11-14: ', a1, 'GELSS, 15-18: ',
     $     a1, 'GELSD)' )
*     
      return
*     
*     End of ALAHD
*     
      END