NAG Library Routine Document

f06ttf (zutsqr)

 Contents

    1  Purpose
    7  Accuracy
    10  Example

1
Purpose

f06ttf performs a QR or RQ factorization of the product of a complex upper triangular matrix and a complex matrix of plane rotations.

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

3
Description

f06ttf performs one of the transformations
RPUQH   or   RQUPH ,  
where U is a given n by 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 k1 to k2.
If side='L',
RPUQH ,  
where P= Pk2-1 Pk1+1 Pk1  and Q= Qk2-1 Qk1+1 Qk1 .
If side='R',
RQUPH ,  
where P=Pk1 Pk1+1 Pk2-1  and Q= Qk1 Qk1+1 Qk2-1 .
In either case Pk and Qk are rotations in the k,k+1 plane.
The 2 by 2 rotation part of Pk or Qk has the form
ck s-k -sk ck  
with ck real.

4
References

None.

5
Arguments

1:     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' or 'R'.
2:     n – IntegerInput
On entry: n, the order of the matrices U and R.
Constraint: n0.
3:     k1 – IntegerInput
4:     k2 – IntegerInput
On entry: the values k1 and k2.
If k1<1 or k2k1 or k2>n, an immediate return is effected.
5:     c* – Real (Kind=nag_wp) arrayInput/Output
Note: the dimension of the array c must be at least k2-1.
On entry: ck must hold the cosine of the rotation Pk, for k=k1,,k2-1.
On exit: ck holds the cosine of the rotation Qk, for k=k1,,k2-1.
6:     s* – Complex (Kind=nag_wp) arrayInput/Output
Note: the dimension of the array s must be at least k2-1.
On entry: sk must hold the sine of the rotation Pk, for k=k1,,k2-1.
On exit: sk holds the sine of the rotation Qk, for k=k1,,k2-1.
7:     alda* – Complex (Kind=nag_wp) arrayInput/Output
Note: the second dimension of the array a must be at least n.
On entry: the n by n upper triangular matrix U.
On exit: the upper triangular matrix R.
8:     lda – IntegerInput
On entry: the first dimension of the array a as declared in the (sub)program from which f06ttf is called.
Constraint: lda max1,n .

6
Error Indicators and Warnings

None.

7
Accuracy

Not applicable.

8
Parallelism and Performance

f06ttf is not threaded in any implementation.

9
Further Comments

None.

10
Example

None.
© The Numerical Algorithms Group Ltd, Oxford, UK. 2017