nag_zungtr (f08ftc) generates the complex unitary matrix
, which was determined by
nag_zhetrd (f08fsc) when reducing a Hermitian matrix to tridiagonal form.
nag_zungtr (f08ftc) is intended to be used after a call to
nag_zhetrd (f08fsc), which reduces a complex Hermitian matrix
to real symmetric tridiagonal form
by a unitary similarity transformation:
.
nag_zhetrd (f08fsc) represents the unitary matrix
as a product of
elementary reflectors.
The computed matrix
differs from an exactly unitary matrix by a matrix
such that
where
is the
machine precision.
The real analogue of this function is
nag_dorgtr (f08ffc).
This example computes all the eigenvalues and eigenvectors of the matrix
, where
Here
is Hermitian and must first be reduced to tridiagonal form by
nag_zhetrd (f08fsc). The program then calls nag_zungtr (f08ftc) to form
, and passes this matrix to
nag_zsteqr (f08jsc) which computes the eigenvalues and eigenvectors of
.