zget54()

subroutine zget54	(	integer	N,
		complex*16, dimension( lda, * )	A,
		integer	LDA,
		complex*16, dimension( ldb, * )	B,
		integer	LDB,
		complex*16, dimension( lds, * )	S,
		integer	LDS,
		complex*16, dimension( ldt, * )	T,
		integer	LDT,
		complex*16, dimension( ldu, * )	U,
		integer	LDU,
		complex*16, dimension( ldv, * )	V,
		integer	LDV,
		complex*16, dimension( * )	WORK,
		double precision	RESULT
	)

ZGET54

Purpose:

 ZGET54 checks a generalized decomposition of the form

          A = U*S*V'  and B = U*T* V'

 where ' means conjugate transpose and U and V are unitary.

 Specifically,

   RESULT = ||( A - U*S*V', B - U*T*V' )|| / (||( A, B )||*n*ulp )

Parameters

[in]	N	N is INTEGER The size of the matrix. If it is zero, DGET54 does nothing. It must be at least zero.
[in]	A	A is COMPLEX*16 array, dimension (LDA, N) The original (unfactored) matrix A.
[in]	LDA	LDA is INTEGER The leading dimension of A. It must be at least 1 and at least N.
[in]	B	B is COMPLEX*16 array, dimension (LDB, N) The original (unfactored) matrix B.
[in]	LDB	LDB is INTEGER The leading dimension of B. It must be at least 1 and at least N.
[in]	S	S is COMPLEX*16 array, dimension (LDS, N) The factored matrix S.
[in]	LDS	LDS is INTEGER The leading dimension of S. It must be at least 1 and at least N.
[in]	T	T is COMPLEX*16 array, dimension (LDT, N) The factored matrix T.
[in]	LDT	LDT is INTEGER The leading dimension of T. It must be at least 1 and at least N.
[in]	U	U is COMPLEX*16 array, dimension (LDU, N) The orthogonal matrix on the left-hand side in the decomposition.
[in]	LDU	LDU is INTEGER The leading dimension of U. LDU must be at least N and at least 1.
[in]	V	V is COMPLEX*16 array, dimension (LDV, N) The orthogonal matrix on the left-hand side in the decomposition.
[in]	LDV	LDV is INTEGER The leading dimension of V. LDV must be at least N and at least 1.
[out]	WORK	WORK is COMPLEX*16 array, dimension (3*N**2)
[out]	RESULT	RESULT is DOUBLE PRECISION The value RESULT, It is currently limited to 1/ulp, to avoid overflow. Errors are flagged by RESULT=10/ulp.

Author

Univ. of Tennessee

Univ. of California Berkeley

Univ. of Colorado Denver

NAG Ltd.