NAG Library Routine Document

F08LSF (ZGBBRD)

INTEGER	M, N, NCC, KL, KU, LDAB, LDQ, LDPT, LDC, INFO
REAL (KIND=nag_wp)	D(min(M,N)), E(min(M,N)-1), RWORK(max(M,N))
COMPLEX (KIND=nag_wp)	AB(LDAB,), Q(LDQ,), PT(LDPT,), C(LDC,), WORK(max(M,N))
CHARACTER(1)	VECT

The routine may be called by its LAPACK name zgbbrd.

3 Description

F08LSF (ZGBBRD) reduces a complex

m

n

band matrix to real upper bidiagonal form

B

by a unitary transformation:

A = Q B P^{H}

. The unitary matrices

Q

and

P^{H}

, of order

m

and

n

respectively, are determined as a product of Givens rotation matrices, and may be formed explicitly by the routine if required. A matrix

C

may also be updated to give

\tilde{C} = Q^{H} C

The routine uses a vectorizable form of the reduction.

4 References

None.

5 Parameters

1: VECT – CHARACTER(1)Input

On entry: indicates whether the matrices

Q

and/or

P^{H}

are generated.

$VECT ='N'$: Neither $Q$ nor $P^{H}$ is generated.
$VECT ='Q'$: $Q$ is generated.
$VECT ='P'$: $P^{H}$ is generated.
$VECT ='B'$: Both $Q$ and $P^{H}$ are generated.

Constraint:

VECT ='N'

'Q'

'P'

'B'

2: M – INTEGERInput

On entry:

m

, the number of rows of the matrix

A

Constraint:

M \geq 0

3: N – INTEGERInput

On entry:

n

, the number of columns of the matrix

A

Constraint:

N \geq 0

4: NCC – INTEGERInput

On entry:

n_{C}

, the number of columns of the matrix

C

Constraint:

NCC \geq 0

5: KL – INTEGERInput

On entry: the number of subdiagonals,

k_{l}

, within the band of

A

Constraint:

KL \geq 0

6: KU – INTEGERInput

On entry: the number of superdiagonals,

k_{u}

, within the band of

A

Constraint:

KU \geq 0

7: AB(LDAB, $*$ ) – COMPLEX (KIND=nag_wp) arrayInput/Output

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

\max (1, N)

On entry: the original

m

n

band matrix

A

The matrix is stored in rows

1

k_{l} + k_{u} + 1

, more precisely, the element

A_{i j}

must be stored in

AB (k_{u} + 1 + i - j, j) for ​ \max (1, j - k_{u}) \leq i \leq \min (m, j + k_{l}) .

On exit: AB is overwritten by values generated during the reduction.

8: LDAB – INTEGERInput

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

Constraint:

LDAB \geq KL + KU + 1

9: D( $\min (M, N)$ ) – REAL (KIND=nag_wp) arrayOutput

On exit: the diagonal elements of the bidiagonal matrix

B

10: E( $\min (M, N) - 1$ ) – REAL (KIND=nag_wp) arrayOutput

On exit: the superdiagonal elements of the bidiagonal matrix

B

11: Q(LDQ, $*$ ) – COMPLEX (KIND=nag_wp) arrayOutput

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

\max (1, M)

VECT ='Q'

'B'

, and at least

1

otherwise.

On exit: if

VECT ='Q'

'B'

, contains the

m

m

unitary matrix

Q

VECT ='N'

'P'

, Q is not referenced.

12: LDQ – INTEGERInput

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

Constraints:

if $VECT ='Q'$ or $'B'$ , $LDQ \geq \max (1, M)$ ;
otherwise $LDQ \geq 1$ .

13: PT(LDPT, $*$ ) – COMPLEX (KIND=nag_wp) arrayOutput

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

\max (1, N)

VECT ='P'

'B'

, and at least

1

otherwise.

On exit: the

n

n

unitary matrix

P^{H}

, if

VECT ='P'

'B'

. If

VECT ='N'

'Q'

, PT is not referenced.

14: LDPT – INTEGERInput

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

Constraints:

if $VECT ='P'$ or $'B'$ , $LDPT \geq \max (1, N)$ ;
otherwise $LDPT \geq 1$ .

15: C(LDC, $*$ ) – COMPLEX (KIND=nag_wp) arrayInput/Output

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

\max (1, NCC)

On entry: an

m

n_{C}

matrix

C

On exit: C is overwritten by

Q^{H} C

. If

NCC = 0

, C is not referenced.

16: LDC – INTEGERInput

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

Constraints:

if $NCC > 0$ , $LDC \geq \max (1, M)$ ;
if $NCC = 0$ , $LDC \geq 1$ .

17: WORK( $\max (M, N)$ ) – COMPLEX (KIND=nag_wp) arrayWorkspace

18: RWORK( $\max (M, N)$ ) – REAL (KIND=nag_wp) arrayWorkspace

19: INFO – INTEGEROutput

On exit:

INFO = 0

unless the routine detects an error (see Section 6).

6 Error Indicators and Warnings

$INFO < 0$: If $INFO = - i$ , argument $i$ had an illegal value. An explanatory message is output, and execution of the program is terminated.

7 Accuracy

The computed bidiagonal form

B

satisfies

Q B P^{H} = A + E

, where

{‖E‖}_{2} \leq c (n) ε {‖A‖}_{2},

c (n)

is a modestly increasing function of

n

, and

ε

is the machine precision.

The elements of

B

themselves may be sensitive to small perturbations in

A

or to rounding errors in the computation, but this does not affect the stability of the singular values and vectors.

The computed matrix

Q

differs from an exactly unitary matrix by a matrix

F

such that

{‖F‖}_{2} = O (ε) .

A similar statement holds for the computed matrix

P^{H}

8 Further Comments

The total number of real floating point operations is approximately the sum of:

$20 n^{2} k$ , if $VECT ='N'$ and $NCC = 0$ , and
$10 n^{2} n_{C} (k - 1) / k$ , if $C$ is updated, and
$10 n^{3} (k - 1) / k$ , if either $Q$ or $P^{H}$ is generated (double this if both),

where

k = k_{l} + k_{u}

, assuming

n ≫ k

. For this section we assume that

m = n

The real analogue of this routine is F08LEF (DGBBRD).

9 Example

This example reduces the matrix

A

to upper bidiagonal form, where

A = (\begin{array}{r} 0.96 - 0.81 i & - 0.03 + 0.96 i & 0.00 + 0.00 i & 0.00 + 0.00 i \\ - 0.98 + 1.98 i & - 1.20 + 0.19 i & - 0.66 + 0.42 i & 0.00 + 0.00 i \\ 0.62 - 0.46 i & 1.01 + 0.02 i & 0.63 - 0.17 i & - 1.11 + 0.60 i \\ 0.00 + 0.00 i & 0.19 - 0.54 i & - 0.98 - 0.36 i & 0.22 - 0.20 i \\ 0.00 + 0.00 i & 0.00 + 0.00 i & - 0.17 - 0.46 i & 1.47 + 1.59 i \\ 0.00 + 0.00 i & 0.00 + 0.00 i & 0.00 + 0.00 i & 0.26 + 0.26 i \end{array}) .

NAG Library Routine DocumentF08LSF (ZGBBRD)

+− Contents

1 Purpose

2 Specification

3 Description

4 References

5 Parameters

6 Error Indicators and Warnings

7 Accuracy

8 Further Comments

9 Example

9.1 Program Text

9.2 Program Data

9.3 Program Results

NAG Library Routine Document

F08LSF (ZGBBRD)