DBDT02 tests the change of basis C = U**H * B by computing the


 residual


    RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ),


 where B and C are M by N matrices, U is an M by M orthogonal matrix,


 and EPS is the machine precision.

M

          M is INTEGER


          The number of rows of the matrices B and C and the order of


          the matrix Q.

N

          N is INTEGER


          The number of columns of the matrices B and C.

B

          B is DOUBLE PRECISION array, dimension (LDB,N)


          The m by n matrix B.

LDB

          LDB is INTEGER


          The leading dimension of the array B.  LDB >= max(1,M).

C

          C is DOUBLE PRECISION array, dimension (LDC,N)


          The m by n matrix C, assumed to contain U**H * B.

LDC

          LDC is INTEGER


          The leading dimension of the array C.  LDC >= max(1,M).

U

          U is DOUBLE PRECISION array, dimension (LDU,M)


          The m by m orthogonal matrix U.

LDU

          LDU is INTEGER


          The leading dimension of the array U.  LDU >= max(1,M).

WORK

          WORK is DOUBLE PRECISION array, dimension (M)

RESID

          RESID is DOUBLE PRECISION


          RESID = norm(B - U * C) / ( max(m,n) * norm(B) * EPS ),

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.