F06 Chapter Contents
F06 Chapter Introduction
NAG Library Manual

# NAG Library Routine DocumentF06TSF

Note:  before using this routine, please read the Users' Note for your implementation to check the interpretation of bold italicised terms and other implementation-dependent details.

## 1  Purpose

F06TSF performs a $QR$ or $RQ$ factorization (as a sequence of plane rotations) of a complex upper spiked matrix.

## 2  Specification

 SUBROUTINE F06TSF ( SIDE, N, K1, K2, C, S, A, LDA)
 INTEGER N, K1, K2, LDA REAL (KIND=nag_wp) C(K2-1) COMPLEX (KIND=nag_wp) S(*), A(LDA,*) CHARACTER(1) SIDE

## 3  Description

F06TSF transforms an $n$ by $n$ complex upper spiked matrix $H$ to upper triangular form $R$ by applying a complex unitary matrix $P$ from the left or the right. $H$ is assumed to have real diagonal elements except where the spike joins the diagonal; $R$ has real diagonal elements. $P$ is formed as a sequence of plane rotations in planes ${k}_{1}$ to ${k}_{2}$.
If ${\mathbf{SIDE}}=\text{'L'}$, $H$ is assumed to have a row spike, with nonzero elements ${h}_{{\mathit{k}}_{2},\mathit{k}}$, for $\mathit{k}={\mathit{k}}_{1},\dots ,{\mathit{k}}_{2}-1$. The rotations are applied from the left:
 $PH=R ,$
where $P=D{P}_{{k}_{2}-1}\cdots {P}_{{k}_{1}+1}{P}_{{k}_{1}}$, ${P}_{k}$ is a rotation in the $\left(k,{k}_{2}\right)$ plane and $D=\mathrm{diag}\left(1,\dots ,1,{d}_{{k}_{2}},1,\dots ,1\right)$ with $\left|{d}_{{k}_{2}}\right|=1$.
If ${\mathbf{SIDE}}=\text{'R'}$, $H$ is assumed to have a column spike, with nonzero elements ${h}_{\mathit{k}+1,{\mathit{k}}_{1}}$, for $\mathit{k}={\mathit{k}}_{1},\dots ,{\mathit{k}}_{2}-1$. The rotations are applied from the right:
 $HPH=R ,$
where $P=D{P}_{{k}_{1}}{P}_{{k}_{1}+1}\cdots {P}_{{k}_{2}-1}$, ${P}_{k}$ is a rotation in the $\left({k}_{1},k+1\right)$ plane and $D=\mathrm{diag}\left(1,\dots ,1,{d}_{{k}_{1}},1,\dots ,1\right)$ with $\left|{d}_{{k}_{1}}\right|=1$.
The $2$ by $2$ plane rotation part of ${P}_{k}$ has the form
 $ck s-k -sk ck$
with ${c}_{k}$ real.
None.

## 5  Parameters

1:     SIDE – CHARACTER(1)Input
On entry: specifies whether $H$ is operated on from the left or the right.
${\mathbf{SIDE}}=\text{'L'}$
$H$ is pre-multiplied from the left.
${\mathbf{SIDE}}=\text{'R'}$
$H$ is post-multiplied from the right.
Constraint: ${\mathbf{SIDE}}=\text{'L'}$ or $\text{'R'}$.
2:     N – INTEGERInput
On entry: $n$, the order of the matrix $H$.
Constraint: ${\mathbf{N}}\ge 0$.
3:     K1 – INTEGERInput
4:     K2 – INTEGERInput
On entry: the values ${k}_{1}$ and ${k}_{2}$.
If ${\mathbf{K1}}<1$ or ${\mathbf{K2}}\le {\mathbf{K1}}$ or ${\mathbf{K2}}>{\mathbf{N}}$, an immediate return is effected.
5:     C(${\mathbf{K2}}-1$) – REAL (KIND=nag_wp) arrayOutput
On exit: ${\mathbf{C}}\left(\mathit{k}\right)$ holds ${c}_{\mathit{k}}$, the cosine of the rotation ${P}_{\mathit{k}}$, for $\mathit{k}={\mathit{k}}_{1},\dots ,{\mathit{k}}_{2}-1$.
6:     S($*$) – COMPLEX (KIND=nag_wp) arrayInput/Output
Note: the dimension of the array S must be at least ${\mathbf{K2}}-{\mathbf{K1}}$.
On entry: the nonzero elements of the spike of $H$: ${\mathbf{S}}\left(\mathit{k}\right)$ must hold ${h}_{{\mathit{k}}_{2},\mathit{k}}$ if ${\mathbf{SIDE}}=\text{'L'}$, and ${h}_{\mathit{k}+1,{\mathit{k}}_{1}}$ if ${\mathbf{SIDE}}=\text{'R'}$, for $\mathit{k}={\mathit{k}}_{1},\dots ,{\mathit{k}}_{2}-1$.
On exit: ${\mathbf{S}}\left(\mathit{k}\right)$ holds ${s}_{\mathit{k}}$, the sine of the rotation ${P}_{\mathit{k}}$, for $\mathit{k}={\mathit{k}}_{1},\dots ,{\mathit{k}}_{2}-1$; ${\mathbf{S}}\left({\mathit{k}}_{2}\right)$ holds ${d}_{{k}_{2}}$, the ${k}_{2}$th diagonal element of $D$, if ${\mathbf{SIDE}}=\text{'L'}$, or ${d}_{{k}_{1}}$, the ${k}_{1}$th diagonal element of $D$, if ${\mathbf{SIDE}}=\text{'R'}$.
7:     A(LDA,$*$) – COMPLEX (KIND=nag_wp) arrayInput/Output
Note: the second dimension of the array A must be at least ${\mathbf{N}}$.
On entry: the upper triangular part of the $n$ by $n$ upper spiked matrix $H$. The imaginary parts of the diagonal elements must be zero, except for the $\left({k}_{2},{k}_{2}\right)$ element if ${\mathbf{SIDE}}=\text{'L'}$, or the $\left({k}_{1},{k}_{1}\right)$ element if ${\mathbf{SIDE}}=\text{'R'}$.
On exit: the upper triangular matrix $R$. The imaginary parts of the diagonal elements are set to zero.
8:     LDA – INTEGERInput
On entry: the first dimension of the array A as declared in the (sub)program from which F06TSF is called.
Constraint: ${\mathbf{LDA}}\ge \mathrm{max}\phantom{\rule{0.125em}{0ex}}\left(1,{\mathbf{N}}\right)$.

None.

Not applicable.