This manual relates to an old release of the Library.
The documentation for the current release is also available on this site.

NAG FL Interface
f06ttf (zutsqr)

Settings help

FL Name Style:


FL Specification Language:


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 <nag.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)
The routine may be called by the names f06ttf or nagf_blas_zutsqr.

3 Description

f06ttf performs one of the transformations
RPUQH   or   RQUPH ,  
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 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×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 Integer Input
On entry: n, the order of the matrices U and R.
Constraint: n0.
3: k1 Integer Input
4: k2 Integer Input
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) array Input/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 Pk, for k=k1,,k2-1.
On exit: c(k) holds the cosine of the rotation Qk, for k=k1,,k2-1.
6: s(*) Complex (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 the sine of the rotation Pk, for k=k1,,k2-1.
On exit: s(k) holds the sine of the rotation Qk, for k=k1,,k2-1.
7: a(lda,*) Complex (Kind=nag_wp) array Input/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.
8: lda Integer Input
On entry: the first dimension of the array a as declared in the (sub)program from which f06ttf is called.
Constraint: lda max(1,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.