zlaqz3()

subroutine zlaqz3	(	logical, intent(in)	ILSCHUR,
		logical, intent(in)	ILQ,
		logical, intent(in)	ILZ,
		integer, intent(in)	N,
		integer, intent(in)	ILO,
		integer, intent(in)	IHI,
		integer, intent(in)	NSHIFTS,
		integer, intent(in)	NBLOCK_DESIRED,
		complex*16, dimension( * ), intent(inout)	ALPHA,
		complex*16, dimension( * ), intent(inout)	BETA,
		complex*16, dimension( lda, * ), intent(inout)	A,
		integer, intent(in)	LDA,
		complex*16, dimension( ldb, * ), intent(inout)	B,
		integer, intent(in)	LDB,
		complex*16, dimension( ldq, * ), intent(inout)	Q,
		integer, intent(in)	LDQ,
		complex*16, dimension( ldz, * ), intent(inout)	Z,
		integer, intent(in)	LDZ,
		complex*16, dimension( ldqc, * ), intent(inout)	QC,
		integer, intent(in)	LDQC,
		complex*16, dimension( ldzc, * ), intent(inout)	ZC,
		integer, intent(in)	LDZC,
		complex*16, dimension( * ), intent(inout)	WORK,
		integer, intent(in)	LWORK,
		integer, intent(out)	INFO
	)

ZLAQZ3

Download ZLAQZ3 + dependencies [TGZ] [ZIP] [TXT]

Purpose:

 ZLAQZ3 Executes a single multishift QZ sweep

Parameters

[in]	ILSCHUR	ILSCHUR is LOGICAL Determines whether or not to update the full Schur form
[in]	ILQ	ILQ is LOGICAL Determines whether or not to update the matrix Q
[in]	ILZ	ILZ is LOGICAL Determines whether or not to update the matrix Z
[in]	N	N is INTEGER The order of the matrices A, B, Q, and Z. N >= 0.
[in]	ILO	ILO is INTEGER
[in]	IHI	IHI is INTEGER
[in]	NSHIFTS	NSHIFTS is INTEGER The desired number of shifts to use
[in]	NBLOCK_DESIRED	NBLOCK_DESIRED is INTEGER The desired size of the computational windows
[in]	ALPHA	ALPHA is COMPLEX*16 array. SR contains the alpha parts of the shifts to use.
[in]	BETA	BETA is COMPLEX*16 array. SS contains the scale of the shifts to use.
[in,out]	A	A is COMPLEX*16 array, dimension (LDA, N)
[in]	LDA	LDA is INTEGER The leading dimension of the array A. LDA >= max( 1, N ).
[in,out]	B	B is COMPLEX*16 array, dimension (LDB, N)
[in]	LDB	LDB is INTEGER The leading dimension of the array B. LDB >= max( 1, N ).
[in,out]	Q	Q is COMPLEX*16 array, dimension (LDQ, N)
[in]	LDQ	LDQ is INTEGER
[in,out]	Z	Z is COMPLEX*16 array, dimension (LDZ, N)
[in]	LDZ	LDZ is INTEGER
[in,out]	QC	QC is COMPLEX*16 array, dimension (LDQC, NBLOCK_DESIRED)
[in]	LDQC	LDQC is INTEGER
[in,out]	ZC	ZC is COMPLEX*16 array, dimension (LDZC, NBLOCK_DESIRED)
[in]	LDZC	LDZ is INTEGER
[out]	WORK	WORK is COMPLEX*16 array, dimension (MAX(1,LWORK)) On exit, if INFO >= 0, WORK(1) returns the optimal LWORK.
[in]	LWORK	LWORK is INTEGER The dimension of the array WORK. LWORK >= max(1,N). If LWORK = -1, then a workspace query is assumed; the routine only calculates the optimal size of the WORK array, returns this value as the first entry of the WORK array, and no error message related to LWORK is issued by XERBLA.
[out]	INFO	INFO is INTEGER = 0: successful exit < 0: if INFO = -i, the i-th argument had an illegal value

Author

Thijs Steel, KU Leuven

Date

May 2020