F06TWF : NAG Library, Mark 25

NAG Library Routine Document

F06TWF

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

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

2 Specification

3 Description

F06TWF transforms an

n

n

complex upper triangular matrix

U

with real diagonal elements, 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}

k_{2}

H

has real diagonal elements except where the spike joins the diagonal.

SIDE ='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 = P U,

where

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

and

P_{k}

is a rotation in the

(k, k_{2})

plane.

SIDE ='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^{H} = R,

where

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

and

P_{k}

is a rotation in the

(k_{1}, k + 1)

plane.

The

2

2

plane rotation part of

P_{k}

has the form

(\begin{array}{r} c_{k} & {\bar{s}}_{k} \\ - s_{k} & c_{k} \end{array})

with

c_{k}

real.

5 Parameters

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

– INTEGERInput

On entry:

n

, the order of the matrices

U

and

H

Constraint:

N \geq 0

K1

– INTEGERInput

K2

– INTEGERInput

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) arrayInput

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

K2 - K1

On entry:

C (k)

must hold

c_{k}

, the cosine of the rotation

P_{k}

, for

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

S (*)

– COMPLEX (KIND=nag_wp) arrayInput/Output

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

K2 - K1

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 a nonzero element of the spike of

H

h_{k_{2}, k}

SIDE ='L'

, or

h_{k + 1, k_{1}}

SIDE ='R'

, for

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

A (LDA *)

– COMPLEX (KIND=nag_wp) arrayInput/Output

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

N

On entry: the

n

n

upper triangular matrix

U

. The imaginary parts of the diagonal elements must be zero.

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

H

. The imaginary parts of the diagonal elements are set to zero except for the

(k_{2}, k_{2})

element if

SIDE ='L'

, or the

(k_{1}, k_{1})

element if

SIDE ='R'

LDA

– INTEGERInput

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

Constraint:

LDA \geq \max (1, N)

NAG Library Routine Document

F06TWF

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

NAG Library Routine DocumentF06TWF

▸▿ Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Parallelism and Performance

9 Further Comments

10 Example

NAG Library Routine Document

F06TWF