f06ttf : NAG Library, Mark 26

f06ttf performs a

Q R

R Q

factorization of the product of a complex upper triangular matrix and a complex matrix of plane rotations.

Specification

Fortran Interface

Subroutine f06ttf (

side, n, k1, k2, c, s, a, lda)

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

C Header Interface

#include nagmk26.h

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

Description

f06ttf performs one of the transformations

R \leftarrow P U Q^{H} or R \leftarrow Q U P^{H},

where

U

is a given

n

n

complex upper triangular matrix,

P

is a given complex unitary matrix, and

Q

is a complex unitary matrix chosen to make

R

upper triangular. Both

P

and

Q

are represented as sequences of plane rotations in planes

k_{1}

k_{2}

side ='L'

R \leftarrow P U Q^{H},

where

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

and

Q = Q_{k_{2} - 1} \dots Q_{k_{1} + 1} Q_{k_{1}}

side ='R'

R \leftarrow Q U P^{H},

where

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

and

Q = Q_{k_{1}} Q_{k_{1} + 1} \dots Q_{k_{2} - 1}

In either case

P_{k}

and

Q_{k}

are rotations in the

(k, k + 1)

plane.

The

2

2

rotation part of

P_{k}

Q_{k}

has the form

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

with

c_{k}

real.

Arguments

side

– Character(1)Input

On entry: specifies whether

P

is applied from the left or the right in the transformation.

$side ='L'$: $P$ is applied from the left.
$side ='R'$: $P$ is applied from the right.

Constraint:

side ='L'

'R'

n

– IntegerInput

On entry:

n

, the order of the matrices

U

and

R

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/Output

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

k2 - 1

On entry:

c (k)

must hold the cosine of the rotation

P_{k}

, for

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

On exit:

c (k)

holds the cosine of the rotation

Q_{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 - 1

On entry:

s (k)

must hold the sine of the rotation

P_{k}

, for

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

On exit:

s (k)

holds the sine of the rotation

Q_{k}

, 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

On exit: the upper triangular matrix

R

lda

– IntegerInput

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

Constraint:

lda \geq \max (1, n)

NAG Library Routine Document

f06ttf (zutsqr)

▸▿ 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 Library Routine Document

f06ttf (zutsqr)

▸▿ 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

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