DQRT17 computes the ratio
    norm(R**T * op(A)) / ( norm(A) * alpha * max(M,N,NRHS) * EPS ),
 where R = B - op(A)*X, op(A) is A or A**T, depending on TRANS, EPS
 is the machine epsilon, and
    alpha = norm(B) if IRESID = 1 (zero-residual problem)
    alpha = norm(R) if IRESID = 2 (otherwise).
 The norm used is the 1-norm.

TRANS

          TRANS is CHARACTER*1
          Specifies whether or not the transpose of A is used.
          = 'N':  No transpose, op(A) = A.
          = 'T':  Transpose, op(A) = A**T.

IRESID

          IRESID is INTEGER
          IRESID = 1 indicates zero-residual problem.
          IRESID = 2 indicates non-zero residual.

M

          M is INTEGER
          The number of rows of the matrix A.
          If TRANS = 'N', the number of rows of the matrix B.
          If TRANS = 'T', the number of rows of the matrix X.

N

          N is INTEGER
          The number of columns of the matrix  A.
          If TRANS = 'N', the number of rows of the matrix X.
          If TRANS = 'T', the number of rows of the matrix B.

NRHS

          NRHS is INTEGER
          The number of columns of the matrices X and B.

A

          A is DOUBLE PRECISION array, dimension (LDA,N)
          The m-by-n matrix A.

LDA

          LDA is INTEGER
          The leading dimension of the array A. LDA >= M.

X

          X is DOUBLE PRECISION array, dimension (LDX,NRHS)
          If TRANS = 'N', the n-by-nrhs matrix X.
          If TRANS = 'T', the m-by-nrhs matrix X.

LDX

          LDX is INTEGER
          The leading dimension of the array X.
          If TRANS = 'N', LDX >= N.
          If TRANS = 'T', LDX >= M.

B

          B is DOUBLE PRECISION array, dimension (LDB,NRHS)
          If TRANS = 'N', the m-by-nrhs matrix B.
          If TRANS = 'T', the n-by-nrhs matrix B.

LDB

          LDB is INTEGER
          The leading dimension of the array B.
          If TRANS = 'N', LDB >= M.
          If TRANS = 'T', LDB >= N.

C

          C is DOUBLE PRECISION array, dimension (LDB,NRHS)

WORK

          WORK is DOUBLE PRECISION array, dimension (LWORK)

LWORK

          LWORK is INTEGER
          The length of the array WORK.  LWORK >= NRHS*(M+N).

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.