naginterfaces.library.blas.dutsrs

naginterfaces.library.blas.dutsrs(side, k1, k2, c, s, a)

dutsrs transforms a real upper triangular matrix to an upper spiked matrix by applying a given sequence of plane rotations.

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

https://support.nag.com/numeric/nl/nagdoc_30.2/flhtml/f06/f06qwf.html

Parameters

sidestr, length 1

Specifies whether $U$ is operated on from the left or the right.

$side=\text{'L'}$

$U$ is pre-multiplied from the left.

$side=\text{'R'}$

$U$ is post-multiplied from the right.

k1int

The value $k_1$.

k2int

The value $k_2$.

cfloat, array-like, shape $(k2-1)$

$c\left[\text{k}\right]$ must hold $c_{\text{k}}$, the cosine of the rotation $P_{\text{k}}$, for $\text{k}={\text{k}}_1,\dots ,{\text{k}}_2-1$.

sfloat, array-like, shape $(k2-1)$

$s\left[\text{k}\right]$ must hold $s_{\text{k}}$, the sine of the rotation $P_{\text{k}}$, for $\text{k}={\text{k}}_1,\dots ,{\text{k}}_2-1$.

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

The $n\times n$ upper triangular matrix $U$.

Returns

sfloat, ndarray, shape $(k2-1)$

$s\left[\text{k}\right]$ holds a nonzero element of the spike of $H$: $h_{{\text{k}}_2,\text{k}}$ if $side=\text{'L'}$, or $h_{\text{k}+1,{\text{k}}_1}$ if $side=\text{'R'}$, for $\text{k}={\text{k}}_1,\dots ,{\text{k}}_2-1$.

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

The upper triangular part of the upper spiked matrix $H$.

Raises

NagValueError

(errno $1$)

On entry, error in parameter $side$.

Constraint: $side=\text{'L'}$ or $\text{'R'}$.

(errno $2$)

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

Constraint: $n\ge 0$.

Notes

No equivalent traditional C interface for this routine exists in the NAG Library.

dutsrs transforms an $n\times n$ real upper triangular matrix $U$ to an upper spiked matrix $H$, by applying a given sequence of plane rotations from either the left or the right, in planes $k_1$ to $k_2$.

If $side=\text{'L'}$, $H$ has a row spike, with nonzero elements $h_{k_2,k}$, for $k=k_1,k_1+1,\dots ,k_2-1$. The rotations are applied from the left:

$$H=PU\text{,}$$

where $P=P_{k_1}{P}_{k_1+1}\cdots P_{k_2-1}$ and $P_k$ is a rotation in the $(k,k_2)$ plane.

If $side=\text{'R'}$, $H$ has a column spike, with nonzero elements $h_{k+1,k_1}$, for $k=k_1,k_1+1,\dots ,k_2-1$. The rotations are applied from the right:

$$H{P}^T=R\text{,}$$

where $P=P_{k_2-1}\cdots P_{k_1+1}{P}_{k_1}$ and $P_k$ is a rotation in the $(k_1,k+1)$ plane.

The $2\times 2$ plane rotation part of $P_k$ has the form

$$\begin{array}{c}\left(\begin{array}{cc}{c}_k& s_k\\ -s_k& c_k\end{array}\right)\text{.}\end{array}$$