LAPACK 3.3.1
Linear Algebra PACKage

dlattb.f

Go to the documentation of this file.
00001       SUBROUTINE DLATTB( IMAT, UPLO, TRANS, DIAG, ISEED, N, KD, AB,
00002      $                   LDAB, B, WORK, INFO )
00003 *
00004 *  -- LAPACK test routine (version 3.1) --
00005 *     Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
00006 *     November 2006
00007 *
00008 *     .. Scalar Arguments ..
00009       CHARACTER          DIAG, TRANS, UPLO
00010       INTEGER            IMAT, INFO, KD, LDAB, N
00011 *     ..
00012 *     .. Array Arguments ..
00013       INTEGER            ISEED( 4 )
00014       DOUBLE PRECISION   AB( LDAB, * ), B( * ), WORK( * )
00015 *     ..
00016 *
00017 *  Purpose
00018 *  =======
00019 *
00020 *  DLATTB generates a triangular test matrix in 2-dimensional storage.
00021 *  IMAT and UPLO uniquely specify the properties of the test matrix,
00022 *  which is returned in the array A.
00023 *
00024 *  Arguments
00025 *  =========
00026 *
00027 *  IMAT    (input) INTEGER
00028 *          An integer key describing which matrix to generate for this
00029 *          path.
00030 *
00031 *  UPLO    (input) CHARACTER*1
00032 *          Specifies whether the matrix A will be upper or lower
00033 *          triangular.
00034 *          = 'U':  Upper triangular
00035 *          = 'L':  Lower triangular
00036 *
00037 *  TRANS   (input) CHARACTER*1
00038 *          Specifies whether the matrix or its transpose will be used.
00039 *          = 'N':  No transpose
00040 *          = 'T':  Transpose
00041 *          = 'C':  Conjugate transpose (= transpose)
00042 *
00043 *  DIAG    (output) CHARACTER*1
00044 *          Specifies whether or not the matrix A is unit triangular.
00045 *          = 'N':  Non-unit triangular
00046 *          = 'U':  Unit triangular
00047 *
00048 *  ISEED   (input/output) INTEGER array, dimension (4)
00049 *          The seed vector for the random number generator (used in
00050 *          DLATMS).  Modified on exit.
00051 *
00052 *  N       (input) INTEGER
00053 *          The order of the matrix to be generated.
00054 *
00055 *  KD      (input) INTEGER
00056 *          The number of superdiagonals or subdiagonals of the banded
00057 *          triangular matrix A.  KD >= 0.
00058 *
00059 *  AB      (output) DOUBLE PRECISION array, dimension (LDAB,N)
00060 *          The upper or lower triangular banded matrix A, stored in the
00061 *          first KD+1 rows of AB.  Let j be a column of A, 1<=j<=n.
00062 *          If UPLO = 'U', AB(kd+1+i-j,j) = A(i,j) for max(1,j-kd)<=i<=j.
00063 *          If UPLO = 'L', AB(1+i-j,j)    = A(i,j) for j<=i<=min(n,j+kd).
00064 *
00065 *  LDAB    (input) INTEGER
00066 *          The leading dimension of the array AB.  LDAB >= KD+1.
00067 *
00068 *  B       (workspace) DOUBLE PRECISION array, dimension (N)
00069 *
00070 *  WORK    (workspace) DOUBLE PRECISION array, dimension (2*N)
00071 *
00072 *  INFO    (output) INTEGER
00073 *          = 0:  successful exit
00074 *          < 0: if INFO = -k, the k-th argument had an illegal value
00075 *
00076 *  =====================================================================
00077 *
00078 *     .. Parameters ..
00079       DOUBLE PRECISION   ONE, TWO, ZERO
00080       PARAMETER          ( ONE = 1.0D+0, TWO = 2.0D+0, ZERO = 0.0D+0 )
00081 *     ..
00082 *     .. Local Scalars ..
00083       LOGICAL            UPPER
00084       CHARACTER          DIST, PACKIT, TYPE
00085       CHARACTER*3        PATH
00086       INTEGER            I, IOFF, IY, J, JCOUNT, KL, KU, LENJ, MODE
00087       DOUBLE PRECISION   ANORM, BIGNUM, BNORM, BSCAL, CNDNUM, PLUS1,
00088      $                   PLUS2, REXP, SFAC, SMLNUM, STAR1, TEXP, TLEFT,
00089      $                   TNORM, TSCAL, ULP, UNFL
00090 *     ..
00091 *     .. External Functions ..
00092       LOGICAL            LSAME
00093       INTEGER            IDAMAX
00094       DOUBLE PRECISION   DLAMCH, DLARND
00095       EXTERNAL           LSAME, IDAMAX, DLAMCH, DLARND
00096 *     ..
00097 *     .. External Subroutines ..
00098       EXTERNAL           DCOPY, DLABAD, DLARNV, DLATB4, DLATMS, DSCAL,
00099      $                   DSWAP
00100 *     ..
00101 *     .. Intrinsic Functions ..
00102       INTRINSIC          ABS, DBLE, MAX, MIN, SIGN, SQRT
00103 *     ..
00104 *     .. Executable Statements ..
00105 *
00106       PATH( 1: 1 ) = 'Double precision'
00107       PATH( 2: 3 ) = 'TB'
00108       UNFL = DLAMCH( 'Safe minimum' )
00109       ULP = DLAMCH( 'Epsilon' )*DLAMCH( 'Base' )
00110       SMLNUM = UNFL
00111       BIGNUM = ( ONE-ULP ) / SMLNUM
00112       CALL DLABAD( SMLNUM, BIGNUM )
00113       IF( ( IMAT.GE.6 .AND. IMAT.LE.9 ) .OR. IMAT.EQ.17 ) THEN
00114          DIAG = 'U'
00115       ELSE
00116          DIAG = 'N'
00117       END IF
00118       INFO = 0
00119 *
00120 *     Quick return if N.LE.0.
00121 *
00122       IF( N.LE.0 )
00123      $   RETURN
00124 *
00125 *     Call DLATB4 to set parameters for SLATMS.
00126 *
00127       UPPER = LSAME( UPLO, 'U' )
00128       IF( UPPER ) THEN
00129          CALL DLATB4( PATH, IMAT, N, N, TYPE, KL, KU, ANORM, MODE,
00130      $                CNDNUM, DIST )
00131          KU = KD
00132          IOFF = 1 + MAX( 0, KD-N+1 )
00133          KL = 0
00134          PACKIT = 'Q'
00135       ELSE
00136          CALL DLATB4( PATH, -IMAT, N, N, TYPE, KL, KU, ANORM, MODE,
00137      $                CNDNUM, DIST )
00138          KL = KD
00139          IOFF = 1
00140          KU = 0
00141          PACKIT = 'B'
00142       END IF
00143 *
00144 *     IMAT <= 5:  Non-unit triangular matrix
00145 *
00146       IF( IMAT.LE.5 ) THEN
00147          CALL DLATMS( N, N, DIST, ISEED, TYPE, B, MODE, CNDNUM, ANORM,
00148      $                KL, KU, PACKIT, AB( IOFF, 1 ), LDAB, WORK, INFO )
00149 *
00150 *     IMAT > 5:  Unit triangular matrix
00151 *     The diagonal is deliberately set to something other than 1.
00152 *
00153 *     IMAT = 6:  Matrix is the identity
00154 *
00155       ELSE IF( IMAT.EQ.6 ) THEN
00156          IF( UPPER ) THEN
00157             DO 20 J = 1, N
00158                DO 10 I = MAX( 1, KD+2-J ), KD
00159                   AB( I, J ) = ZERO
00160    10          CONTINUE
00161                AB( KD+1, J ) = J
00162    20       CONTINUE
00163          ELSE
00164             DO 40 J = 1, N
00165                AB( 1, J ) = J
00166                DO 30 I = 2, MIN( KD+1, N-J+1 )
00167                   AB( I, J ) = ZERO
00168    30          CONTINUE
00169    40       CONTINUE
00170          END IF
00171 *
00172 *     IMAT > 6:  Non-trivial unit triangular matrix
00173 *
00174 *     A unit triangular matrix T with condition CNDNUM is formed.
00175 *     In this version, T only has bandwidth 2, the rest of it is zero.
00176 *
00177       ELSE IF( IMAT.LE.9 ) THEN
00178          TNORM = SQRT( CNDNUM )
00179 *
00180 *        Initialize AB to zero.
00181 *
00182          IF( UPPER ) THEN
00183             DO 60 J = 1, N
00184                DO 50 I = MAX( 1, KD+2-J ), KD
00185                   AB( I, J ) = ZERO
00186    50          CONTINUE
00187                AB( KD+1, J ) = DBLE( J )
00188    60       CONTINUE
00189          ELSE
00190             DO 80 J = 1, N
00191                DO 70 I = 2, MIN( KD+1, N-J+1 )
00192                   AB( I, J ) = ZERO
00193    70          CONTINUE
00194                AB( 1, J ) = DBLE( J )
00195    80       CONTINUE
00196          END IF
00197 *
00198 *        Special case:  T is tridiagonal.  Set every other offdiagonal
00199 *        so that the matrix has norm TNORM+1.
00200 *
00201          IF( KD.EQ.1 ) THEN
00202             IF( UPPER ) THEN
00203                AB( 1, 2 ) = SIGN( TNORM, DLARND( 2, ISEED ) )
00204                LENJ = ( N-3 ) / 2
00205                CALL DLARNV( 2, ISEED, LENJ, WORK )
00206                DO 90 J = 1, LENJ
00207                   AB( 1, 2*( J+1 ) ) = TNORM*WORK( J )
00208    90          CONTINUE
00209             ELSE
00210                AB( 2, 1 ) = SIGN( TNORM, DLARND( 2, ISEED ) )
00211                LENJ = ( N-3 ) / 2
00212                CALL DLARNV( 2, ISEED, LENJ, WORK )
00213                DO 100 J = 1, LENJ
00214                   AB( 2, 2*J+1 ) = TNORM*WORK( J )
00215   100          CONTINUE
00216             END IF
00217          ELSE IF( KD.GT.1 ) THEN
00218 *
00219 *           Form a unit triangular matrix T with condition CNDNUM.  T is
00220 *           given by
00221 *                   | 1   +   *                      |
00222 *                   |     1   +                      |
00223 *               T = |         1   +   *              |
00224 *                   |             1   +              |
00225 *                   |                 1   +   *      |
00226 *                   |                     1   +      |
00227 *                   |                          . . . |
00228 *        Each element marked with a '*' is formed by taking the product
00229 *        of the adjacent elements marked with '+'.  The '*'s can be
00230 *        chosen freely, and the '+'s are chosen so that the inverse of
00231 *        T will have elements of the same magnitude as T.
00232 *
00233 *        The two offdiagonals of T are stored in WORK.
00234 *
00235             STAR1 = SIGN( TNORM, DLARND( 2, ISEED ) )
00236             SFAC = SQRT( TNORM )
00237             PLUS1 = SIGN( SFAC, DLARND( 2, ISEED ) )
00238             DO 110 J = 1, N, 2
00239                PLUS2 = STAR1 / PLUS1
00240                WORK( J ) = PLUS1
00241                WORK( N+J ) = STAR1
00242                IF( J+1.LE.N ) THEN
00243                   WORK( J+1 ) = PLUS2
00244                   WORK( N+J+1 ) = ZERO
00245                   PLUS1 = STAR1 / PLUS2
00246 *
00247 *                 Generate a new *-value with norm between sqrt(TNORM)
00248 *                 and TNORM.
00249 *
00250                   REXP = DLARND( 2, ISEED )
00251                   IF( REXP.LT.ZERO ) THEN
00252                      STAR1 = -SFAC**( ONE-REXP )
00253                   ELSE
00254                      STAR1 = SFAC**( ONE+REXP )
00255                   END IF
00256                END IF
00257   110       CONTINUE
00258 *
00259 *           Copy the tridiagonal T to AB.
00260 *
00261             IF( UPPER ) THEN
00262                CALL DCOPY( N-1, WORK, 1, AB( KD, 2 ), LDAB )
00263                CALL DCOPY( N-2, WORK( N+1 ), 1, AB( KD-1, 3 ), LDAB )
00264             ELSE
00265                CALL DCOPY( N-1, WORK, 1, AB( 2, 1 ), LDAB )
00266                CALL DCOPY( N-2, WORK( N+1 ), 1, AB( 3, 1 ), LDAB )
00267             END IF
00268          END IF
00269 *
00270 *     IMAT > 9:  Pathological test cases.  These triangular matrices
00271 *     are badly scaled or badly conditioned, so when used in solving a
00272 *     triangular system they may cause overflow in the solution vector.
00273 *
00274       ELSE IF( IMAT.EQ.10 ) THEN
00275 *
00276 *        Type 10:  Generate a triangular matrix with elements between
00277 *        -1 and 1. Give the diagonal norm 2 to make it well-conditioned.
00278 *        Make the right hand side large so that it requires scaling.
00279 *
00280          IF( UPPER ) THEN
00281             DO 120 J = 1, N
00282                LENJ = MIN( J, KD+1 )
00283                CALL DLARNV( 2, ISEED, LENJ, AB( KD+2-LENJ, J ) )
00284                AB( KD+1, J ) = SIGN( TWO, AB( KD+1, J ) )
00285   120       CONTINUE
00286          ELSE
00287             DO 130 J = 1, N
00288                LENJ = MIN( N-J+1, KD+1 )
00289                IF( LENJ.GT.0 )
00290      $            CALL DLARNV( 2, ISEED, LENJ, AB( 1, J ) )
00291                AB( 1, J ) = SIGN( TWO, AB( 1, J ) )
00292   130       CONTINUE
00293          END IF
00294 *
00295 *        Set the right hand side so that the largest value is BIGNUM.
00296 *
00297          CALL DLARNV( 2, ISEED, N, B )
00298          IY = IDAMAX( N, B, 1 )
00299          BNORM = ABS( B( IY ) )
00300          BSCAL = BIGNUM / MAX( ONE, BNORM )
00301          CALL DSCAL( N, BSCAL, B, 1 )
00302 *
00303       ELSE IF( IMAT.EQ.11 ) THEN
00304 *
00305 *        Type 11:  Make the first diagonal element in the solve small to
00306 *        cause immediate overflow when dividing by T(j,j).
00307 *        In type 11, the offdiagonal elements are small (CNORM(j) < 1).
00308 *
00309          CALL DLARNV( 2, ISEED, N, B )
00310          TSCAL = ONE / DBLE( KD+1 )
00311          IF( UPPER ) THEN
00312             DO 140 J = 1, N
00313                LENJ = MIN( J, KD+1 )
00314                CALL DLARNV( 2, ISEED, LENJ, AB( KD+2-LENJ, J ) )
00315                CALL DSCAL( LENJ-1, TSCAL, AB( KD+2-LENJ, J ), 1 )
00316                AB( KD+1, J ) = SIGN( ONE, AB( KD+1, J ) )
00317   140       CONTINUE
00318             AB( KD+1, N ) = SMLNUM*AB( KD+1, N )
00319          ELSE
00320             DO 150 J = 1, N
00321                LENJ = MIN( N-J+1, KD+1 )
00322                CALL DLARNV( 2, ISEED, LENJ, AB( 1, J ) )
00323                IF( LENJ.GT.1 )
00324      $            CALL DSCAL( LENJ-1, TSCAL, AB( 2, J ), 1 )
00325                AB( 1, J ) = SIGN( ONE, AB( 1, J ) )
00326   150       CONTINUE
00327             AB( 1, 1 ) = SMLNUM*AB( 1, 1 )
00328          END IF
00329 *
00330       ELSE IF( IMAT.EQ.12 ) THEN
00331 *
00332 *        Type 12:  Make the first diagonal element in the solve small to
00333 *        cause immediate overflow when dividing by T(j,j).
00334 *        In type 12, the offdiagonal elements are O(1) (CNORM(j) > 1).
00335 *
00336          CALL DLARNV( 2, ISEED, N, B )
00337          IF( UPPER ) THEN
00338             DO 160 J = 1, N
00339                LENJ = MIN( J, KD+1 )
00340                CALL DLARNV( 2, ISEED, LENJ, AB( KD+2-LENJ, J ) )
00341                AB( KD+1, J ) = SIGN( ONE, AB( KD+1, J ) )
00342   160       CONTINUE
00343             AB( KD+1, N ) = SMLNUM*AB( KD+1, N )
00344          ELSE
00345             DO 170 J = 1, N
00346                LENJ = MIN( N-J+1, KD+1 )
00347                CALL DLARNV( 2, ISEED, LENJ, AB( 1, J ) )
00348                AB( 1, J ) = SIGN( ONE, AB( 1, J ) )
00349   170       CONTINUE
00350             AB( 1, 1 ) = SMLNUM*AB( 1, 1 )
00351          END IF
00352 *
00353       ELSE IF( IMAT.EQ.13 ) THEN
00354 *
00355 *        Type 13:  T is diagonal with small numbers on the diagonal to
00356 *        make the growth factor underflow, but a small right hand side
00357 *        chosen so that the solution does not overflow.
00358 *
00359          IF( UPPER ) THEN
00360             JCOUNT = 1
00361             DO 190 J = N, 1, -1
00362                DO 180 I = MAX( 1, KD+1-( J-1 ) ), KD
00363                   AB( I, J ) = ZERO
00364   180          CONTINUE
00365                IF( JCOUNT.LE.2 ) THEN
00366                   AB( KD+1, J ) = SMLNUM
00367                ELSE
00368                   AB( KD+1, J ) = ONE
00369                END IF
00370                JCOUNT = JCOUNT + 1
00371                IF( JCOUNT.GT.4 )
00372      $            JCOUNT = 1
00373   190       CONTINUE
00374          ELSE
00375             JCOUNT = 1
00376             DO 210 J = 1, N
00377                DO 200 I = 2, MIN( N-J+1, KD+1 )
00378                   AB( I, J ) = ZERO
00379   200          CONTINUE
00380                IF( JCOUNT.LE.2 ) THEN
00381                   AB( 1, J ) = SMLNUM
00382                ELSE
00383                   AB( 1, J ) = ONE
00384                END IF
00385                JCOUNT = JCOUNT + 1
00386                IF( JCOUNT.GT.4 )
00387      $            JCOUNT = 1
00388   210       CONTINUE
00389          END IF
00390 *
00391 *        Set the right hand side alternately zero and small.
00392 *
00393          IF( UPPER ) THEN
00394             B( 1 ) = ZERO
00395             DO 220 I = N, 2, -2
00396                B( I ) = ZERO
00397                B( I-1 ) = SMLNUM
00398   220       CONTINUE
00399          ELSE
00400             B( N ) = ZERO
00401             DO 230 I = 1, N - 1, 2
00402                B( I ) = ZERO
00403                B( I+1 ) = SMLNUM
00404   230       CONTINUE
00405          END IF
00406 *
00407       ELSE IF( IMAT.EQ.14 ) THEN
00408 *
00409 *        Type 14:  Make the diagonal elements small to cause gradual
00410 *        overflow when dividing by T(j,j).  To control the amount of
00411 *        scaling needed, the matrix is bidiagonal.
00412 *
00413          TEXP = ONE / DBLE( KD+1 )
00414          TSCAL = SMLNUM**TEXP
00415          CALL DLARNV( 2, ISEED, N, B )
00416          IF( UPPER ) THEN
00417             DO 250 J = 1, N
00418                DO 240 I = MAX( 1, KD+2-J ), KD
00419                   AB( I, J ) = ZERO
00420   240          CONTINUE
00421                IF( J.GT.1 .AND. KD.GT.0 )
00422      $            AB( KD, J ) = -ONE
00423                AB( KD+1, J ) = TSCAL
00424   250       CONTINUE
00425             B( N ) = ONE
00426          ELSE
00427             DO 270 J = 1, N
00428                DO 260 I = 3, MIN( N-J+1, KD+1 )
00429                   AB( I, J ) = ZERO
00430   260          CONTINUE
00431                IF( J.LT.N .AND. KD.GT.0 )
00432      $            AB( 2, J ) = -ONE
00433                AB( 1, J ) = TSCAL
00434   270       CONTINUE
00435             B( 1 ) = ONE
00436          END IF
00437 *
00438       ELSE IF( IMAT.EQ.15 ) THEN
00439 *
00440 *        Type 15:  One zero diagonal element.
00441 *
00442          IY = N / 2 + 1
00443          IF( UPPER ) THEN
00444             DO 280 J = 1, N
00445                LENJ = MIN( J, KD+1 )
00446                CALL DLARNV( 2, ISEED, LENJ, AB( KD+2-LENJ, J ) )
00447                IF( J.NE.IY ) THEN
00448                   AB( KD+1, J ) = SIGN( TWO, AB( KD+1, J ) )
00449                ELSE
00450                   AB( KD+1, J ) = ZERO
00451                END IF
00452   280       CONTINUE
00453          ELSE
00454             DO 290 J = 1, N
00455                LENJ = MIN( N-J+1, KD+1 )
00456                CALL DLARNV( 2, ISEED, LENJ, AB( 1, J ) )
00457                IF( J.NE.IY ) THEN
00458                   AB( 1, J ) = SIGN( TWO, AB( 1, J ) )
00459                ELSE
00460                   AB( 1, J ) = ZERO
00461                END IF
00462   290       CONTINUE
00463          END IF
00464          CALL DLARNV( 2, ISEED, N, B )
00465          CALL DSCAL( N, TWO, B, 1 )
00466 *
00467       ELSE IF( IMAT.EQ.16 ) THEN
00468 *
00469 *        Type 16:  Make the offdiagonal elements large to cause overflow
00470 *        when adding a column of T.  In the non-transposed case, the
00471 *        matrix is constructed to cause overflow when adding a column in
00472 *        every other step.
00473 *
00474          TSCAL = UNFL / ULP
00475          TSCAL = ( ONE-ULP ) / TSCAL
00476          DO 310 J = 1, N
00477             DO 300 I = 1, KD + 1
00478                AB( I, J ) = ZERO
00479   300       CONTINUE
00480   310    CONTINUE
00481          TEXP = ONE
00482          IF( KD.GT.0 ) THEN
00483             IF( UPPER ) THEN
00484                DO 330 J = N, 1, -KD
00485                   DO 320 I = J, MAX( 1, J-KD+1 ), -2
00486                      AB( 1+( J-I ), I ) = -TSCAL / DBLE( KD+2 )
00487                      AB( KD+1, I ) = ONE
00488                      B( I ) = TEXP*( ONE-ULP )
00489                      IF( I.GT.MAX( 1, J-KD+1 ) ) THEN
00490                         AB( 2+( J-I ), I-1 ) = -( TSCAL / DBLE( KD+2 ) )
00491      $                                          / DBLE( KD+3 )
00492                         AB( KD+1, I-1 ) = ONE
00493                         B( I-1 ) = TEXP*DBLE( ( KD+1 )*( KD+1 )+KD )
00494                      END IF
00495                      TEXP = TEXP*TWO
00496   320             CONTINUE
00497                   B( MAX( 1, J-KD+1 ) ) = ( DBLE( KD+2 ) /
00498      $                                    DBLE( KD+3 ) )*TSCAL
00499   330          CONTINUE
00500             ELSE
00501                DO 350 J = 1, N, KD
00502                   TEXP = ONE
00503                   LENJ = MIN( KD+1, N-J+1 )
00504                   DO 340 I = J, MIN( N, J+KD-1 ), 2
00505                      AB( LENJ-( I-J ), J ) = -TSCAL / DBLE( KD+2 )
00506                      AB( 1, J ) = ONE
00507                      B( J ) = TEXP*( ONE-ULP )
00508                      IF( I.LT.MIN( N, J+KD-1 ) ) THEN
00509                         AB( LENJ-( I-J+1 ), I+1 ) = -( TSCAL /
00510      $                     DBLE( KD+2 ) ) / DBLE( KD+3 )
00511                         AB( 1, I+1 ) = ONE
00512                         B( I+1 ) = TEXP*DBLE( ( KD+1 )*( KD+1 )+KD )
00513                      END IF
00514                      TEXP = TEXP*TWO
00515   340             CONTINUE
00516                   B( MIN( N, J+KD-1 ) ) = ( DBLE( KD+2 ) /
00517      $                                    DBLE( KD+3 ) )*TSCAL
00518   350          CONTINUE
00519             END IF
00520          ELSE
00521             DO 360 J = 1, N
00522                AB( 1, J ) = ONE
00523                B( J ) = DBLE( J )
00524   360       CONTINUE
00525          END IF
00526 *
00527       ELSE IF( IMAT.EQ.17 ) THEN
00528 *
00529 *        Type 17:  Generate a unit triangular matrix with elements
00530 *        between -1 and 1, and make the right hand side large so that it
00531 *        requires scaling.
00532 *
00533          IF( UPPER ) THEN
00534             DO 370 J = 1, N
00535                LENJ = MIN( J-1, KD )
00536                CALL DLARNV( 2, ISEED, LENJ, AB( KD+1-LENJ, J ) )
00537                AB( KD+1, J ) = DBLE( J )
00538   370       CONTINUE
00539          ELSE
00540             DO 380 J = 1, N
00541                LENJ = MIN( N-J, KD )
00542                IF( LENJ.GT.0 )
00543      $            CALL DLARNV( 2, ISEED, LENJ, AB( 2, J ) )
00544                AB( 1, J ) = DBLE( J )
00545   380       CONTINUE
00546          END IF
00547 *
00548 *        Set the right hand side so that the largest value is BIGNUM.
00549 *
00550          CALL DLARNV( 2, ISEED, N, B )
00551          IY = IDAMAX( N, B, 1 )
00552          BNORM = ABS( B( IY ) )
00553          BSCAL = BIGNUM / MAX( ONE, BNORM )
00554          CALL DSCAL( N, BSCAL, B, 1 )
00555 *
00556       ELSE IF( IMAT.EQ.18 ) THEN
00557 *
00558 *        Type 18:  Generate a triangular matrix with elements between
00559 *        BIGNUM/KD and BIGNUM so that at least one of the column
00560 *        norms will exceed BIGNUM.
00561 *
00562          TLEFT = BIGNUM / MAX( ONE, DBLE( KD ) )
00563          TSCAL = BIGNUM*( DBLE( KD ) / DBLE( KD+1 ) )
00564          IF( UPPER ) THEN
00565             DO 400 J = 1, N
00566                LENJ = MIN( J, KD+1 )
00567                CALL DLARNV( 2, ISEED, LENJ, AB( KD+2-LENJ, J ) )
00568                DO 390 I = KD + 2 - LENJ, KD + 1
00569                   AB( I, J ) = SIGN( TLEFT, AB( I, J ) ) +
00570      $                         TSCAL*AB( I, J )
00571   390          CONTINUE
00572   400       CONTINUE
00573          ELSE
00574             DO 420 J = 1, N
00575                LENJ = MIN( N-J+1, KD+1 )
00576                CALL DLARNV( 2, ISEED, LENJ, AB( 1, J ) )
00577                DO 410 I = 1, LENJ
00578                   AB( I, J ) = SIGN( TLEFT, AB( I, J ) ) +
00579      $                         TSCAL*AB( I, J )
00580   410          CONTINUE
00581   420       CONTINUE
00582          END IF
00583          CALL DLARNV( 2, ISEED, N, B )
00584          CALL DSCAL( N, TWO, B, 1 )
00585       END IF
00586 *
00587 *     Flip the matrix if the transpose will be used.
00588 *
00589       IF( .NOT.LSAME( TRANS, 'N' ) ) THEN
00590          IF( UPPER ) THEN
00591             DO 430 J = 1, N / 2
00592                LENJ = MIN( N-2*J+1, KD+1 )
00593                CALL DSWAP( LENJ, AB( KD+1, J ), LDAB-1,
00594      $                     AB( KD+2-LENJ, N-J+1 ), -1 )
00595   430       CONTINUE
00596          ELSE
00597             DO 440 J = 1, N / 2
00598                LENJ = MIN( N-2*J+1, KD+1 )
00599                CALL DSWAP( LENJ, AB( 1, J ), 1, AB( LENJ, N-J+2-LENJ ),
00600      $                     -LDAB+1 )
00601   440       CONTINUE
00602          END IF
00603       END IF
00604 *
00605       RETURN
00606 *
00607 *     End of DLATTB
00608 *
00609       END
 All Files Functions