dlakf2()

subroutine dlakf2	(	integer	M,
		integer	N,
		double precision, dimension( lda, * )	A,
		integer	LDA,
		double precision, dimension( lda, * )	B,
		double precision, dimension( lda, * )	D,
		double precision, dimension( lda, * )	E,
		double precision, dimension( ldz, * )	Z,
		integer	LDZ
	)

DLAKF2

Purpose:

 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.

Parameters

[in]	M	M is INTEGER Size of matrix, must be >= 1.
[in]	N	N is INTEGER Size of matrix, must be >= 1.
[in]	A	A is DOUBLE PRECISION, dimension ( LDA, M ) The matrix A in the output matrix Z.
[in]	LDA	LDA is INTEGER The leading dimension of A, B, D, and E. ( LDA >= M+N )
[in]	B	B is DOUBLE PRECISION, dimension ( LDA, N )
[in]	D	D is DOUBLE PRECISION, dimension ( LDA, M )
[in]	E	E is DOUBLE PRECISION, dimension ( LDA, N ) The matrices used in forming the output matrix Z.
[out]	Z	Z is DOUBLE PRECISION, dimension ( LDZ, 2*M*N ) The resultant Kronecker M*N*2 by M*N*2 matrix (see above.)
[in]	LDZ	LDZ is INTEGER The leading dimension of Z. ( LDZ >= 2*M*N )

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.