Form the 2*M*N by 2*M*N matrix
        Z = [ kron(In, A)  -kron(B', Im) ]
            [ kron(In, D)  -kron(E', Im) ],
 where In is the identity matrix of size n and X' is the transpose
 of X. kron(X, Y) is the Kronecker product between the matrices X
 and Y.

M

          M is INTEGER
          Size of matrix, must be >= 1.

N

          N is INTEGER
          Size of matrix, must be >= 1.

A

          A is DOUBLE PRECISION, dimension ( LDA, M )
          The matrix A in the output matrix Z.

LDA

          LDA is INTEGER
          The leading dimension of A, B, D, and E. ( LDA >= M+N )

B

          B is DOUBLE PRECISION, dimension ( LDA, N )

D

          D is DOUBLE PRECISION, dimension ( LDA, M )

E

          E is DOUBLE PRECISION, dimension ( LDA, N )
          The matrices used in forming the output matrix Z.

Z

          Z is DOUBLE PRECISION, dimension ( LDZ, 2*M*N )
          The resultant Kronecker M*N*2 by M*N*2 matrix (see above.)

LDZ

          LDZ is INTEGER
          The leading dimension of Z. ( LDZ >= 2*M*N )

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.