ZGLMTS tests ZGGGLM - a subroutine for solving the generalized
 linear model problem.

N

          N is INTEGER
          The number of rows of the matrices A and B.  N >= 0.

M

          M is INTEGER
          The number of columns of the matrix A.  M >= 0.

P

          P is INTEGER
          The number of columns of the matrix B.  P >= 0.

A

          A is COMPLEX*16 array, dimension (LDA,M)
          The N-by-M matrix A.

AF

          AF is COMPLEX*16 array, dimension (LDA,M)

LDA

          LDA is INTEGER
          The leading dimension of the arrays A, AF. LDA >= max(M,N).

B

          B is COMPLEX*16 array, dimension (LDB,P)
          The N-by-P matrix A.

BF

          BF is COMPLEX*16 array, dimension (LDB,P)

LDB

          LDB is INTEGER
          The leading dimension of the arrays B, BF. LDB >= max(P,N).

D

          D is COMPLEX*16 array, dimension( N )
          On input, the left hand side of the GLM.

DF

          DF is COMPLEX*16 array, dimension( N )

X

          X is COMPLEX*16 array, dimension( M )
          solution vector X in the GLM problem.

U

          U is COMPLEX*16 array, dimension( P )
          solution vector U in the GLM problem.

WORK

          WORK is COMPLEX*16 array, dimension (LWORK)

LWORK

          LWORK is INTEGER
          The dimension of the array WORK.

RWORK

          RWORK is DOUBLE PRECISION array, dimension (M)

RESULT

          RESULT is DOUBLE PRECISION
          The test ratio:
                           norm( d - A*x - B*u )
            RESULT = -----------------------------------------
                     (norm(A)+norm(B))*(norm(x)+norm(u))*EPS

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.