f06qvf (dutsrh) : NAG Library, Mark 27

Integer, Intent (In)	::	n, k1, k2, lda
Real (Kind=nag_wp), Intent (In)	::	c(*)
Real (Kind=nag_wp), Intent (Inout)	::	s(), a(lda,)
Character (1), Intent (In)	::	side

C Header Interface

#include <nag.h>

void	f06qvf_ (const char side, const Integer n, const Integer k1, const Integer k2, const double c[], double s[], double a[], const Integer *lda, const Charlen length_side)

C++ Header Interface

#include <nag.h>

extern "C" {

void	f06qvf_ (const char *side, const Integer &n, const Integer &k1, const Integer &k2, const double c[], double s[], double a[], const Integer &lda, const Charlen length_side)

}

3 Description

f06qvf transforms an

n \times n

real upper triangular matrix

U

to an upper Hessenberg matrix

H

, by applying a given sequence of plane rotations from either the left or the right, in planes

k_{1}

k_{2}

;

H

has nonzero subdiagonal elements

h_{k + 1, k}

, for

k = k_{1}, \dots, k_{2} - 1

only.

side ='L'

, the rotations are applied from the left:

H = P U,

where

P = P_{k_{1}} P_{k_{1} + 1} \dots P_{k_{2} - 1}

side ='R'

, the rotations are applied from the right:

H = U P^{T},

where

P = P_{k_{2} - 1} \dots P_{k_{1} + 1} P_{k_{1}}

The

2 \times 2

plane rotation part of

P_{k}

has the form

(\begin{array}{r} c_{k} & s_{k} \\ - s_{k} & c_{k} \end{array}) .

5 Arguments

side

– Character(1) Input

On entry: specifies whether

U

is operated on from the left or the right.

$side ='L'$: $U$ is pre-multiplied from the left.
$side ='R'$: $U$ is post-multiplied from the right.

Constraint:

side ='L'

'R'

n

– Integer Input

On entry:

n

, the order of the matrices

U

and

H

Constraint:

n \geq 0

k1

– Integer Input

k2

– Integer Input

On entry: the values

k_{1}

and

k_{2}

k1 < 1

k2 \leq k1

k2 > n

, an immediate return is effected.

c (*)

– Real (Kind=nag_wp) array Input

Note: the dimension of the array c must be at least

k2 - 1

On entry:

c (k)

must hold

c_{k}

, the cosine of the rotation

P_{k}

, for

k = k_{1}, \dots, k_{2} - 1

s (*)

– Real (Kind=nag_wp) array Input/Output

Note: the dimension of the array s must be at least

k2 - 1

On entry:

s (k)

must hold

s_{k}

, the sine of the rotation

P_{k}

, for

k = k_{1}, \dots, k_{2} - 1

On exit:

s (k)

holds

h_{k + 1, k}

, the subdiagonal element of

H

, for

k = k_{1}, \dots, k_{2} - 1

a (lda, *)

– Real (Kind=nag_wp) array Input/Output

Note: the second dimension of the array a must be at least

n

On entry: the

n \times n

upper triangular matrix

U

On exit: the upper triangular part of the upper Hessenberg matrix

H

lda

– Integer Input

On entry: the first dimension of the array a as declared in the (sub)program from which f06qvf is called.

Constraint:

lda \geq \max (1, n)

NAG FL Interface
f06qvf (dutsrh)

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Arguments

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

NAG FL Interfacef06qvf (dutsrh)

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Arguments

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

NAG FL Interface
f06qvf (dutsrh)