DGET31 tests DLALN2, a routine for solving
    (ca A - w D)X = sB
 where A is an NA by NA matrix (NA=1 or 2 only), w is a real (NW=1) or
 complex (NW=2) constant, ca is a real constant, D is an NA by NA real
 diagonal matrix, and B is an NA by NW matrix (when NW=2 the second
 column of B contains the imaginary part of the solution).  The code
 returns X and s, where s is a scale factor, less than or equal to 1,
 which is chosen to avoid overflow in X.
 If any singular values of ca A-w D are less than another input
 parameter SMIN, they are perturbed up to SMIN.
 The test condition is that the scaled residual
     norm( (ca A-w D)*X - s*B ) /
           ( max( ulp*norm(ca A-w D), SMIN )*norm(X) )
 should be on the order of 1.  Here, ulp is the machine precision.
 Also, it is verified that SCALE is less than or equal to 1, and that
 XNORM = infinity-norm(X).

RMAX

          RMAX is DOUBLE PRECISION
          Value of the largest test ratio.

LMAX

          LMAX is INTEGER
          Example number where largest test ratio achieved.

NINFO

          NINFO is INTEGER array, dimension (2)
          NINFO(1) = number of examples with INFO less than 0
          NINFO(2) = number of examples with INFO greater than 0

KNT

          KNT is INTEGER
          Total number of examples tested.

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.