naginterfaces.library.lapackeig.dgebrd

naginterfaces.library.lapackeig.dgebrd(a)

dgebrd reduces a real $m\times n$ matrix to bidiagonal form.

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

https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/f08/f08kef.html

Parameters

afloat, array-like, shape $(m,n)$

The $m\times n$ matrix $A$.

Returns

afloat, ndarray, shape $(m,n)$

If $m\ge n$, the diagonal and first superdiagonal are overwritten by the upper bidiagonal matrix $B$, elements below the diagonal are overwritten by details of the orthogonal matrix $Q$ and elements above the first superdiagonal are overwritten by details of the orthogonal matrix $P$.

If $m<n$, the diagonal and first subdiagonal are overwritten by the lower bidiagonal matrix $B$, elements below the first subdiagonal are overwritten by details of the orthogonal matrix $Q$ and elements above the diagonal are overwritten by details of the orthogonal matrix $P$.

dfloat, ndarray, shape $(min(m,n))$

The diagonal elements of the bidiagonal matrix $B$.

efloat, ndarray, shape $(min(m,n)-1)$

The off-diagonal elements of the bidiagonal matrix $B$.

tauqfloat, ndarray, shape $(min(m,n))$

Further details of the orthogonal matrix $Q$.

taupfloat, ndarray, shape $(min(m,n))$

Further details of the orthogonal matrix $P$.

Raises

NagValueError

(errno $-1$)

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

Constraint: $m\ge 0$.

(errno $-2$)

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

Constraint: $n\ge 0$.

Notes

dgebrd reduces a real $m\times n$ matrix $A$ to bidiagonal form $B$ by an orthogonal transformation: $A=QB{P}^T$, where $Q$ and $P^T$ are orthogonal matrices of order $m$ and $n$ respectively.

If $m\ge n$, the reduction is given by:

$$\begin{array}{c}A=Q\left(\begin{array}{c}{B}_1\\ 0\end{array}\right)P^T=Q_1{B}_1{P}^T\text{,}\end{array}$$

where $B_1$ is an $n\times n$ upper bidiagonal matrix and $Q_1$ consists of the first $n$ columns of $Q$.

If $m<n$, the reduction is given by

$$A=Q\left(\begin{array}{cc}{B}_1& 0\end{array}\right)P^T=Q{B}_1{P}_1^T\text{,}$$

where $B_1$ is an $m\times m$ lower bidiagonal matrix and $P_1^T$ consists of the first $m$ rows of $P^T$.

The orthogonal matrices $Q$ and $P$ are not formed explicitly but are represented as products of elementary reflectors (see the F08 Introduction for details). Functions are provided to work with $Q$ and $P$ in this representation (see Further Comments).

References

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