- naginterfaces.library.matop.complex_trapez_rq(m, a)¶
complex_trapez_rqreduces the complex () upper trapezoidal matrix to upper triangular form by means of unitary transformations.
For full information please refer to the NAG Library document for f01rg
, the number of rows of the matrix .
When then an immediate return is effected.
- acomplex, array-like, shape
The leading upper trapezoidal part of the array must contain the matrix to be factorized.
- acomplex, ndarray, shape
The upper triangular part of will contain the upper triangular matrix , and the upper trapezoidal part of will contain details of the factorization as described in Notes.
- thetacomplex, ndarray, shape
contains the scalar for the th transformation. If then ; if
then , otherwise contains as described in Notes and is always in the range .
- (errno )
On entry, and .
- (errno )
On entry, .
No equivalent traditional C interface for this routine exists in the NAG Library.
The upper trapezoidal matrix given by
where is an upper triangular matrix, is factorized as
where is an unitary matrix and is an upper triangular matrix.
is given as a sequence of Householder transformation matrices
the th transformation matrix, , being used to introduce zeros into the th row of . has the form
is a scalar for which , is a real scalar and is an element vector. , and are chosen to annihilate the elements of the th row of and to make the diagonal elements of real.
The scalar and the vector are returned in the th element of the array and in the th row of , such that , given by
is in and the elements of are in . The elements of are returned in the upper triangular part of .
For further information on this factorization and its use see Section 6.5 of Golub and Van Loan (1996).
Golub, G H and Van Loan, C F, 1996, Matrix Computations, (3rd Edition), Johns Hopkins University Press, Baltimore
Wilkinson, J H, 1965, The Algebraic Eigenvalue Problem, Oxford University Press, Oxford