naginterfaces.library.blas.dsrotg

naginterfaces.library.blas.dsrotg(pivot, direct, alpha, x)

dsrotg generates a sequence of real plane rotations.

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

https://www.nag.com/numeric/nl/nagdoc_29.3/flhtml/f06/f06fqf.html

Parameters

pivotstr, length 1

Specifies the plane rotated by $P_k$.

$pivot=\text{'V'}$ (variable pivot)

$P_k$ rotates the $(k,k+1)$ plane.

$pivot=\text{'F'}$ (fixed pivot)

$P_k$ rotates the $(1,k+1)$ plane if $direct=\text{'F'}$, or the $(k,n+1)$ plane if $direct=\text{'B'}$.

directstr, length 1

Specifies the sequence direction.

$direct=\text{'F'}$ (forward sequence)

$P=P_n\cdots P_2{P}_1$.

$direct=\text{'B'}$ (backward sequence)

$P=P_1{P}_2\cdots P_n$.

alphafloat

The scalar $\alpha$.

xfloat, array-like, shape $\left(n\right)$

The $n$-element vector $x$.

Returns

alphafloat

The scalar $\beta$.

xfloat, ndarray, shape $\left(n\right)$

The referenced elements are overwritten by details of the sequence of plane rotations.

cfloat, ndarray, shape $\left(n\right)$

The values $c_k$, the cosines of the rotations.

sfloat, ndarray, shape $\left(n\right)$

The values $s_k$, the sines of the rotations.

Raises

NagValueError

(errno $1$)

On entry, error in parameter $pivot$.

Constraint: $pivot=\text{'V'}$ or $\text{'F'}$.

(errno $2$)

On entry, error in parameter $direct$.

Constraint: $direct=\text{'F'}$ or $\text{'B'}$.

(errno $6$)

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

Constraint: $\text{incx}>0$.

Notes

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

dsrotg generates the parameters of a real orthogonal matrix $P$, of order $n+1$, chosen so as to set to zero the elements of a supplied $n$-element real vector $x$.

If $pivot=\text{'F'}$ and $direct=\text{'F'}$, or if $pivot=\text{'V'}$ and $direct=\text{'B'}$,

$$\begin{array}{c}P\left(\begin{array}{c}\alpha \\ x\end{array}\right)=\left(\begin{array}{c}\beta \\ 0\end{array}\right)\text{;}\end{array}$$

If $pivot=\text{'F'}$ and $direct=\text{'B'}$, or if $pivot=\text{'V'}$ and $direct=\text{'F'}$,

$$\begin{array}{c}P\left(\begin{array}{c}x\\ \alpha \end{array}\right)=\left(\begin{array}{c}0\\ \beta \end{array}\right)\text{.}\end{array}$$

Here $\alpha$ and $\beta$ are real scalars.

$P$ is represented as a sequence of $n$ plane rotations $P_k$, as specified by $pivot$ and $direct$; $P_k$ is chosen to annihilate $x_k$, and its $2\times 2$ plane rotation part has the form

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

The tangent of the rotation, $t_k$, is overwritten on $x_k$.