NAG Library Routine Document

F08QVF (ZTRSYL)

+− 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

1 Purpose

F08QVF (ZTRSYL) solves the complex triangular Sylvester matrix equation.

2 Specification

SUBROUTINE F08QVF (

TRANA, TRANB, ISGN, M, N, A, LDA, B, LDB, C, LDC, SCAL, INFO)

INTEGER	ISGN, M, N, LDA, LDB, LDC, INFO
REAL (KIND=nag_wp)	SCAL
COMPLEX (KIND=nag_wp)	A(LDA,), B(LDB,), C(LDC,*)
CHARACTER(1)	TRANA, TRANB

The routine may be called by its LAPACK name ztrsyl.

3 Description

F08QVF (ZTRSYL) solves the complex Sylvester matrix equation

op (A) X \pm X op (B) = α C,

where

op (A) = A

A^{H}

, and the matrices

A

and

B

are upper triangular;

α

is a scale factor (

\leq 1

) determined by the routine to avoid overflow in

X

;

A

m

m

and

B

n

n

while the right-hand side matrix

C

and the solution matrix

X

are both

m

n

. The matrix

X

is obtained by a straightforward process of back-substitution (see Golub and Van Loan (1996)).

Note that the equation has a unique solution if and only if

α_{i} \pm β_{j} \neq 0

, where

\{α_{i}\}

and

\{β_{j}\}

are the eigenvalues of

A

and

B

respectively and the sign (

+

-

) is the same as that used in the equation to be solved.

4 References

Golub G H and Van Loan C F (1996) Matrix Computations (3rd Edition) Johns Hopkins University Press, Baltimore

Higham N J (1992) Perturbation theory and backward error for

A X - X B = C

Numerical Analysis Report University of Manchester

5 Parameters

1: TRANA – CHARACTER(1)Input

On entry: specifies the option

op (A)

$TRANA ='N'$: $op (A) = A$ .
$TRANA ='C'$: $op (A) = A^{H}$ .

Constraint:

TRANA ='N'

'C'

2: TRANB – CHARACTER(1)Input

On entry: specifies the option

op (B)

$TRANB ='N'$: $op (B) = B$ .
$TRANB ='C'$: $op (B) = B^{H}$ .

Constraint:

TRANB ='N'

'C'

3: ISGN – INTEGERInput

On entry: indicates the form of the Sylvester equation.

$ISGN = +1$: The equation is of the form $op (A) X + X op (B) = α C$ .
$ISGN = -1$: The equation is of the form $op (A) X - X op (B) = α C$ .

Constraint:

ISGN = +1

-1

4: M – INTEGERInput

On entry:

m

, the order of the matrix

A

, and the number of rows in the matrices

X

and

C

Constraint:

M \geq 0

5: N – INTEGERInput

On entry:

n

, the order of the matrix

B

, and the number of columns in the matrices

X

and

C

Constraint:

N \geq 0

6: A(LDA, $*$ ) – COMPLEX (KIND=nag_wp) arrayInput

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

\max (1, M)

On entry: the

m

m

upper triangular matrix

A

7: LDA – INTEGERInput

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

Constraint:

LDA \geq \max (1, M)

8: B(LDB, $*$ ) – COMPLEX (KIND=nag_wp) arrayInput

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

\max (1, N)

On entry: the

n

n

upper triangular matrix

B

9: LDB – INTEGERInput

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

Constraint:

LDB \geq \max (1, N)

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

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

\max (1, N)

On entry: the

m

n

right-hand side matrix

C

On exit: C is overwritten by the solution matrix

X

11: LDC – INTEGERInput

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

Constraint:

LDC \geq \max (1, M)

12: SCAL – REAL (KIND=nag_wp)Output

On exit: the value of the scale factor

α

13: 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.

$INFO = 1$: $A$ and $B$ have common or close eigenvalues, perturbed values of which were used to solve the equation.

7 Accuracy

Consider the equation

A X - X B = C

. (To apply the remarks to the equation

A X + X B = C

, simply replace

B

- B

Let

\tilde{X}

be the computed solution and

R

the residual matrix:

R = C - (A \tilde{X} - \tilde{X} B) .

Then the residual is always small:

{‖R‖}_{F} = O (ε) ({‖A‖}_{F} + {‖B‖}_{F}) {‖\tilde{X}‖}_{F} .

However,

\tilde{X}

is not necessarily the exact solution of a slightly perturbed equation; in other words, the solution is not backwards stable.

For the forward error, the following bound holds:

{‖\tilde{X} - X‖}_{F} \leq \frac{{‖R‖}_{F}}{sep (A, B)}

but this may be a considerable over estimate. See Golub and Van Loan (1996) for a definition of

sep (A, B)

, and Higham (1992) for further details.

These remarks also apply to the solution of a general Sylvester equation, as described in Section 8.

8 Further Comments

The total number of real floating point operations is approximately

4 m n (m + n)

To solve the general complex Sylvester equation

A X \pm X B = C

where

A

and

B

are general matrices,

A

and

B

must first be reduced to Schur form (by calling F08PNF (ZGEES), for example):

A = Q_{1} \tilde{A} Q_{1}^{H} and B = Q_{2} \tilde{B} Q_{2}^{H}

where

\tilde{A}

and

\tilde{B}

are upper triangular and

Q_{1}

and

Q_{2}

are unitary. The original equation may then be transformed to:

\tilde{A} \tilde{X} \pm \tilde{X} \tilde{B} = \tilde{C}

where

\tilde{X} = Q_{1}^{H} X Q_{2}

and

\tilde{C} = Q_{1}^{H} C Q_{2}

\tilde{C}

may be computed by matrix multiplication; F08QVF (ZTRSYL) may be used to solve the transformed equation; and the solution to the original equation can be obtained as

X = Q_{1} \tilde{X} Q_{2}^{H}

The real analogue of this routine is F08QHF (DTRSYL).

9 Example

This example solves the Sylvester equation

A X + X B = C

, where

A = (\begin{array}{r} - 6.00 - 7.00 i & 0.36 - 0.36 i & - 0.19 + 0.48 i & 0.88 - 0.25 i \\ 0.00 + 0.00 i & - 5.00 + 2.00 i & - 0.03 - 0.72 i & - 0.23 + 0.13 i \\ 0.00 + 0.00 i & 0.00 + 0.00 i & 8.00 - 1.00 i & 0.94 + 0.53 i \\ 0.00 + 0.00 i & 0.00 + 0.00 i & 0.00 + 0.00 i & 3.00 - 4.00 i \end{array}),

B = (\begin{array}{r} 0.50 - 0.20 i & - 0.29 - 0.16 i & - 0.37 + 0.84 i & - 0.55 + 0.73 i \\ 0.00 + 0.00 i & - 0.40 + 0.90 i & 0.06 + 0.22 i & - 0.43 + 0.17 i \\ 0.00 + 0.00 i & 0.00 + 0.00 i & - 0.90 - 0.10 i & - 0.89 - 0.42 i \\ 0.00 + 0.00 i & 0.00 + 0.00 i & 0.00 + 0.00 i & 0.30 - 0.70 i \end{array})

and

C = (\begin{array}{r} 0.63 + 0.35 i & 0.45 - 0.56 i & 0.08 - 0.14 i & - 0.17 - 0.23 i \\ - 0.17 + 0.09 i & - 0.07 - 0.31 i & 0.27 - 0.54 i & 0.35 + 1.21 i \\ - 0.93 - 0.44 i & - 0.33 - 0.35 i & 0.41 - 0.03 i & 0.57 + 0.84 i \\ 0.54 + 0.25 i & - 0.62 - 0.05 i & - 0.52 - 0.13 i & 0.11 - 0.08 i \end{array}) .

NAG Library Routine DocumentF08QVF (ZTRSYL)

+− 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

F08QVF (ZTRSYL)