naginterfaces.library.lapackeig.zunghr

naginterfaces.library.lapackeig.zunghr(ilo, ihi, a, tau)

zunghr generates the complex unitary matrix $Q$ which was determined by zgehrd() when reducing a complex general matrix $A$ to Hessenberg form.

For full information please refer to the NAG Library document for f08nt

https://support.nag.com/numeric/nl/nagdoc_30.1/flhtml/f08/f08ntf.html

Parameters

iloint

These must be the same arguments $ilo$ and $ihi$, respectively, as supplied to zgehrd().

ihiint

These must be the same arguments $ilo$ and $ihi$, respectively, as supplied to zgehrd().

acomplex, array-like, shape $(n,n)$

Details of the vectors which define the elementary reflectors, as returned by zgehrd().

taucomplex, array-like, shape $(max(1,n-1))$

Further details of the elementary reflectors, as returned by zgehrd().

Returns

acomplex, ndarray, shape $(n,n)$

The $n\times n$ unitary matrix $Q$.

Raises

NagValueError

(errno $-1$)

On entry, error in parameter $\text{n}$.

Constraint: $n\ge 0$.

(errno $-2$)

On entry, error in parameter $ilo$.

(errno $-3$)

On entry, error in parameter $ihi$.

Notes

zunghr is intended to be used following a call to zgehrd(), which reduces a complex general matrix $A$ to upper Hessenberg form $H$ by a unitary similarity transformation: $A=QH{Q}^H$. zgehrd() represents the matrix $Q$ as a product of $i_{hi}-i_{lo}$ elementary reflectors. Here $i_{lo}$ and $i_{hi}$ are values determined by zgebal() when balancing the matrix; if the matrix has not been balanced, $i_{lo}=1$ and $i_{hi}=n$.

This function may be used to generate $Q$ explicitly as a square matrix. $Q$ has the structure:

$$\begin{array}{c}Q=\left(\begin{array}{ccc}I& 0& 0\\ 0& Q_{22}& 0\\ 0& 0& I\end{array}\right)\end{array}$$

where $Q_{22}$ occupies rows and columns $i_{lo}$ to $i_{hi}$.

References

Golub, G H and Van Loan, C F, 1996, Matrix Computations, (3rd Edition), Johns Hopkins University Press, Baltimore