!>
!> 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 COMPLEX*16, 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 COMPLEX*16, dimension ( LDA, N )
!> 

D

!>          D is COMPLEX*16, dimension ( LDA, M )
!> 

E

!>          E is COMPLEX*16, dimension ( LDA, N )
!>
!>          The matrices used in forming the output matrix Z.
!> 

Z

!>          Z is COMPLEX*16, 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.