NAG FL Interface
f06tmf (zhesrc)

1 Purpose

f06tmf performs a Unitary similarity transformation (as a sequence of plane rotations) of a complex Hermitian matrix.

2 Specification

Fortran Interface
Subroutine f06tmf ( uplo, pivot, direct, n, k1, k2, c, s, a, lda)
Integer, Intent (In) :: n, k1, k2, lda
Real (Kind=nag_wp), Intent (In) :: c(*)
Complex (Kind=nag_wp), Intent (In) :: s(*)
Complex (Kind=nag_wp), Intent (Inout) :: a(lda,*)
Character (1), Intent (In) :: uplo, pivot, direct
C Header Interface
#include <nag.h>
void  f06tmf_ (const char *uplo, const char *pivot, const char *direct, const Integer *n, const Integer *k1, const Integer *k2, const double c[], const Complex s[], Complex a[], const Integer *lda, const Charlen length_uplo, const Charlen length_pivot, const Charlen length_direct)
The routine may be called by the names f06tmf or nagf_blas_zhesrc.

3 Description

f06tmf performs the transformation
APAPH  
where A is an n by n complex Hermitian matrix, and P is a complex unitary matrix defined as a sequence of plane rotations, Pk, applied in planes k1 to k2.
The 2 by 2 plane rotation part of Pk is assumed to have the form
ck s¯k -sk ck  
with ck real.

4 References

None.

5 Arguments

1: uplo Character(1) Input
On entry: specifies whether the upper or lower triangular part of A is stored.
uplo='U'
The upper triangular part of A is stored.
uplo='L'
The lower triangular part of A is stored.
Constraint: uplo='U' or 'L'.
2: pivot Character(1) Input
On entry: specifies the plane rotated by Pk.
pivot='V' (variable pivot)
Pk rotates the k,k+1 plane.
pivot='T' (top pivot)
Pk rotates the k1,k+1 plane.
pivot='B' (bottom pivot)
Pk rotates the k,k2 plane.
Constraint: pivot='V', 'T' or 'B'.
3: direct Character(1) Input
On entry: specifies the sequence direction.
direct='F' (forward sequence)
P=Pk2-1Pk1+1Pk1.
direct='B' (backward sequence)
P=Pk1Pk1+1Pk2-1.
Constraint: direct='F' or 'B'.
4: n Integer Input
On entry: n, the order of the matrix A.
Constraint: n0.
5: k1 Integer Input
6: k2 Integer Input
On entry: the values k1 and k2.
If k1<1 or k2k1 or k2>n, an immediate return is effected.
7: c* Real (Kind=nag_wp) array Input
Note: the dimension of the array c must be at least k2-1.
On entry: ck must hold ck, the cosine of the rotation Pk, for k=k1,,k2-1.
8: s* Complex (Kind=nag_wp) array Input
Note: the dimension of the array s must be at least k2-1.
On entry: sk must hold sk, the sine of the rotation Pk, for k=k1,,k2-1.
9: alda* Complex (Kind=nag_wp) array Input/Output
Note: the second dimension of the array a must be at least max1,n.
On entry: the n by n Hermitian matrix A.
  • If uplo='U', the upper triangular part of A must be stored and the elements of the array below the diagonal are not referenced.
  • If uplo='L', the lower triangular part of A must be stored and the elements of the array above the diagonal are not referenced.
On exit: the transformed matrix A. The imaginary parts of the diagonal elements are set to zero.
10: lda Integer Input
On entry: the first dimension of the array a as declared in the (sub)program from which f06tmf is called.
Constraint: lda max1,n .

6 Error Indicators and Warnings

None.

7 Accuracy

Not applicable.

8 Parallelism and Performance

f06tmf is not threaded in any implementation.

9 Further Comments

None.

10 Example

None.