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.