Open AccessArticle

A Globally Convergent Iterative Method for Matrix Sign Function and Its Application for Determining the Eigenvalues of a Matrix Pencil

Department of Mathematics, Thapar Institute of Engineering and Technology, Patiala 147004, India

School of Commerce, SVKM’s Narsee Monjee Institute of Management Studies, Chandigarh 160014, India

Department of Physics and Chemistry, Technical University of Cluj-Napoca, 400114 Cluj-Napoca, Romania

Author to whom correspondence should be addressed.

Symmetry 2024, 16(4), 481; https://doi.org/10.3390/sym16040481

Submission received: 11 March 2024 / Revised: 8 April 2024 / Accepted: 10 April 2024 / Published: 16 April 2024

(This article belongs to the Special Issue The Advances of Nonlinear Equations: Mathematical Models, Symmetry and Applications)

Download

Browse Figures

Versions Notes

Abstract

In this research article, we propose a new matrix iterative method with a convergence order of five for computing the sign of a complex matrix by examining the different patterns and symmetry of existing methods. Analysis of the convergence of the method was explored on a global scale, and attraction basins were demonstrated. In addition to this, the asymptotic stability of the scheme was explored.Then, an algorithm for determing thegeneralized eigenvalues for the case of regular matrix pencils was investigated using the matrix sign computation. We performed a series of numerical experiments using numerous matrices to confirm the usefulness and superiority of the proposed method.

Keywords:

iterative method; matrix sign function; basins of attraction; generalized eigenvalue problem; matrix pencil

MSC:

65F60; 15A24

1. Introduction

Over the past several years, the efficient computing of matrix functions has emerged as a prominent area of research. Symmetry plays a significant role in the study of matrix functions, particularly in understanding their properties, behavior, and computational methods. Symmetry principles are widely used in physics and engineering to analyze and solve problems involving matrix functions. For instance, in quantum mechanics, symmetric matrices often arise in the context of Hamiltonians, and efficient computation of matrix functions is crucial for simulating quantum systems and predicting their behavior. In order to compute matrix functions, various iterative schemes have been suggested, which include methods for the evaluation of matrix sign function. The scalar sign function defined for

z \in C

such that

R e (z) \neq 0

which is given as,

s i g n (z) = \{\begin{matrix} 1, R e (z) > 0, \\ - 1, R e (z) < 0 . \end{matrix}

(1)

is extended to form the matrix sign function. Roberts [1] introduced the Cauchy integral representation of the matrix sign function in 1980, given by

s i g n (A) = \frac{2}{π} A \int_{0}^{\infty} {(t^{2} I + A^{2})}^{- 1} d t,

(2)

which can be applied in model reduction and in solving algebraic Riccati and Lyapunov equations. Recall the definition of matrix sign function

s i g n (A)

via the Jordan canonical form [2], in which if

A = P J P^{- 1}

is the Jordan decomposition of A such that

J = d i a g (J_{1}, J_{2})

, then

s i g n (A) = P [\begin{matrix} - I_{p} & 0 \\ 0 & I_{n - p} \end{matrix}] P^{- 1},

(3)

where

A \in C^{n \times n}

has no eigenvalue on the imaginary axis, P is a nonsingular matrix, and p and

n - p

are sizes of Jordan blocks

J_{1}

J_{2}

corresponding to eigenvalues lying in left half plane

(C^{-})

and right half plane

(C^{+})

, respectively.

In the pursuit of evaluating

s i g n (A)

as a solution of

W^{2} - I = 0,

(4)

Roberts applied Newton’s method (

N M

) resulting into the following iterative formula:

W_{l + 1} = \frac{1}{2} (W_{l} + W_{l}^{- 1}), l = 0, 1, 2, \dots,

(5)

with the initial matrix

W_{0} = A

and convergence of order two.

Another alternative is accelerated Newton’s method (

A N M

) that can be achieved via norm scaling parameter in (5) given by,

\{\begin{matrix} W_{0} & = A, \\ μ_{l} & = \sqrt{\frac{∥ W_{l}^{- 1} ∥}{∥ W_{l} ∥}}, \\ W_{l + 1} & = \frac{1}{2} (μ_{l} W_{l} + μ_{l}^{- 1} W_{l}^{- 1}) . \end{matrix}

(6)

In order to improve the convergence, acceleration and scaling of iteration (5), Kenney and Laub [3,4] suggested a family of matrix iterations by employing the Padé approximant to the function

H (x) = {(1 - x)}^{- 1 / 2},

(7)

where

x = 1 - z^{2}

and taking into the consideration

s i g n (z) = \frac{z}{{(1 - x)}^{1 / 2}} .

(8)

Let

(r, s)

-Padé approximant to

H (x)

with

r + s \geq 1

is denoted as

\frac{P_{r, s} (x)}{Q_{r, s} (x)}

, then

w_{l + 1} = \frac{w_{l} P_{r, s} (1 - w_{l}^{2})}{Q_{r, s} (1 - w_{l}^{2})} = ψ_{2 r + 1, 2 s},

(9)

converges to

\pm 1

r \geq s - 1

and it has convergence rate equal to

r + s + 1

. So, the

(r, s)

-Padé approximant iteration (

P A_{r, s}

) for computing

s i g n (A)

W_{l + 1} = W_{l} [P_{r, s} (1 - W_{l}^{2})] {[Q_{r, s} (1 - W_{l}^{2})]}^{- 1} .

(10)

These iterations were later investigated by Gomilko et al. [5]. Other popular iterative approaches for

s i g n (A)

, such as the Newton-Schulz method

(N S M)

[6]

W_{l + 1} = \frac{1}{2} W_{l} (3 I - W_{l}^{2}),

(11)

and Halley’s method

(H M)

[2]

W_{l + 1} = [I + 3 W_{l}^{2}] {[W_{l} (3 I + W_{l}^{2})]}^{- 1} .

(12)

are either members of the Padé family (10) or its reciprocal family.

A recent fifth order iterative method [7] (represented as Lotfi’s Method (

L M

)) for

s i g n (A)

can be investigated and is given by

W_{l + 1} = [W_{l} (18 I - 20 W_{l}^{2} - 30 W_{l}^{4})] {[5 I + 15 W_{l}^{2} - 45 W_{l}^{4} - 7 W_{l}^{6}]}^{- 1} .

(13)

For more recent literature on matrix sign function, one can go through references [8,9,10].

Motivated by constructing efficient iterative methods like the methods discussed above and to overcome on some of the drawbacks of the existing iterative methods, here we aim to propose a globally convergent iterative method for approximation of

s i g n (A)

. Matrix sign function has several application in scientific computing (see e.g., [11,12]). Another applications of matrix sign function include solving various matrix equations [13], solving stochastic differential equations [14] and many more.

In this article, we present a fifth order iterative scheme for evaluation of

s i g n (A)

. The outline for the rest of the article is defined here: In Section 2, the construction of new iterative scheme for evaluating matrix sign function has been concluded. In Section 3, basins of attraction guarantees the global convergence of developed method. Also, convergence analysis for the fifth rate of convergence is discussed along with stability analysis. Then the proposed method is extended to determine the generalized eigenvalues of a regular matrix pencil in Section 4. Computational complexity is examined in Section 5. Numerical examples are provided to evaluate the performance and to demonstrate the numerical precision of suggested method in Section 6. The final conclusion is presented in Section 7.

2. A Novel Iterative Method

In this section, we have developed an iterative method for evaluation of

s i g n (A)

. The primary goal is to develop the scheme by examining the factorization patterns and symmetry of existing schemes rather than using a non-linear equation solver [13]. We aim to design an iterative formula that does not conflict with the Padé family members (9) or its reciprocals.

If we assume the eigenvalues to converge to 1 or

- 1

, and we aim to get a higher convergence order five, then the iterative formula will satisfy the equation

\frac{w_{l + 1} - 1}{w_{l + 1} + 1} = \frac{{(w_{l} - 1)}^{5}}{{(w_{l} + 1)}^{5}} .

(14)

If we solve Equation (14) for

w_{l + 1}

, then we get

w_{l + 1} = \frac{w_{l} (5 + 10 w_{l}^{2} + w_{l}^{4})}{1 + 10 w_{l}^{2} + 5 w_{l}^{4}} .

(15)

Nevertheless, the formula mentioned above is not a novel one, as it is a member of Padé family (9). So, we make some modifications in (14) so that the formula becomes better in terms of converegnce as follows:

\frac{w_{l + 1} - 1}{w_{l + 1} + 1} = \frac{{(w_{l} - 1)}^{5} (- 4 + w_{l})}{{(w_{l} + 1)}^{5} (- 4 - w_{l})},

(16)

provided that

∥\frac{- 4 + w_{l}}{- 4 - w_{l}}∥ < 1

. Solving (16) for

w_{l + 1}

, we get

w_{l + 1} = \frac{w_{l} (21 + 50 w_{l}^{2} + 9 w_{l}^{4})}{4 + 45 w_{l}^{2} + 30 w_{l}^{4} + w_{l}^{6}},

(17)

which contains a polynomial of degree 5 in the numerator and a polynomial of degree 6 in the denominator. From (9), it is clear that Padé family has a unique member with a polynomial of degree

2 r + 1

in numerator and a polynomial of degree

2 s

in denominator for each

r, s \geq 0

. So

2 r + 1 = 5

and

2 s = 6

will result into

r = 2

and

s = 3

, i.e.,

P A_{2, 3}

and will have order

2 + 3 + 1 = 6

which is different from the order 5 that we are achieving. So, the above formula (17) is not a part of the family (9) or reciprocals of (9) and therefore our proposed method (

P M

) is as follows:

W_{l + 1} = [W_{l} (21 I + 50 W_{l}^{2} + 9 W_{l}^{4})] {[4 I + 45 W_{l}^{2} + 30 W_{l}^{4} + W_{l}^{6}]}^{- 1},

(18)

where

W_{0}

is the appropriate initial guess which we discuss in Section 3.

Remark 1.

Convergence can be slow if

∥ W_{l} ∥ > > 1

for some l. In that situation, it is advisable to introduce the scaling parameter

μ_{l}

[15] given as,

μ_{l} = \{\begin{matrix} \sqrt{\frac{∥ W_{l}^{- 1} ∥}{∥ W_{l} ∥}} (norm scaling), \\ \sqrt{\frac{ρ (W_{l}^{- 1})}{ρ (W_{l})}} (spectral scaling), \\ \sqrt{{|d e t (W_{l})|}^{- 1 / n}} (determinant scaling) . \end{matrix}

(19)

to speed up the convergence [16] of the scheme.

Hence, we also present the accelerated proposed method (

A P M

) which is accelerated version of (18) given by

\{\begin{matrix} W_{0} & = A, \\ μ_{l} & = Scaling parameter computed by (19), \\ W_{l + 1} & = [μ_{l} W_{l} (21 I + 50 μ_{l}^{2} W_{l}^{2} + 9 μ_{l}^{4} W_{l}^{4})] {[4 I + 45 μ_{l}^{2} Z_{l}^{2} + 30 μ_{l}^{4} W_{l}^{4} + μ_{l}^{6} W_{l}^{6}]}^{- 1}, \end{matrix}

(20)

where

{lim}_{l \to \infty} μ_{l} = 1

and

{lim}_{l \to \infty} W_{l} = S = s i g n (A)

3. Convergence and Stability Analysis

In this section, we discuss the convergence analysis to ensure that the sequence generated from (18) converges to

s i g n (A)

for a matrix

A \in C^{n \times n}

. We draw the basins of attraction to confirm the global convergence of the proposed scheme and compared with the existing schemes for solving

w^{2} - 1 = 0

[17].

We consider a region

R = [- 2, 2] \times [- 2, 2] \in C

. The region has grid size of

256 \times 256

and has two simple roots

- 1

and 1. The convergence was determined by the stopping criteria

∥ f (w_{l}) ∥ \leq 10^{- 2}

. The exact location of roots is shown by two white dots within the basins. The dark blue color (left side of basins) represent the convergence region corresponding to the root

w = - 1

and the light blue color (right side of basins) represent the convergence region corresponding to the root

w = 1

Figure 1 illustrates the basins for

N M

(5),

N S M

(11),

H M

(12),

P A_{r, s}

(10), and the proposed method

P M

(18). It is evident from Figure 1b,d that methods

N S M

and

P A_{3, 1}

exhibit local convergence, while the other methods exhibit the global convergence. Moreover, Figure 1e,f illustarte broader and lighter basins of attraction because of higher convergence order of

P A_{2, 2}

and

P M

as compared to the basins of

N M

and

H M

as shown in Figure 1a,c. One can notice the symmetry in convergence region of each method. The convergence regions are symmetric about imaginary axis.

After achieving global convergence of the proposed method, it is important to examine certain critical factors related to matrix iterations, including convergence rate and numerical stability.

3.1. Convergence Analysis

Theorem 1.

Suppose that the initial matrix

W_{0} = A \in C^{n \times n}

has no eigenvalue on the imaginary axis, then the sequence of matrices

{W_{l}}_{l = 0}^{\infty}

generated from (18) converges to

s i g n (A) = S

Proof.

Let

W = P J P^{- 1}

be Jordan decomposition of matrix

W \in C^{n \times n}

and consider

R (W) = P R (J) P^{- 1} .

(21)

where R is the rational function linked to (18). Consequently, this will lead to the mapping of an eigenvalue

λ

W_{l}

to an eigenvalue

R (λ)

W_{l + 1}

. Also, R will satisfy the properties given as follows:

$s i g n (R (w)) = s i g n (w)$ ∀ $w \in C$ .
$w_{l + 1} = R (w_{l})$ with $w_{0} = w$ converges to $s i g n (w)$ , where $R e (w) \neq 0$ .

To achieve this objective, consider

A = P Λ P^{- 1}

as the Jordan decomposition of A that can be modified as [2]

P^{- 1} A P = Λ = [\begin{matrix} C & 0 \\ 0 & N \end{matrix}],

(22)

where P is a non singular matrix while Jordan blocks C and N are associated to the eigenvalues in

C^{-}

and

C^{+}

respectively. Let the principal diagonals of blocks C and N have values

λ_{1}, λ_{2}, \dots, λ_{p}

and

λ_{p + 1}, λ_{p + 2}, \dots, λ_{n}

respectively. Then

s i g n (A) = P [\begin{matrix} - I_{p} & 0 \\ 0 & I_{n - p} \end{matrix}] P^{- 1} .

(23)

Therefore,

\begin{matrix} s i g n (Λ) & = s i g n (P^{- 1} A P) = P^{- 1} s i g n (A) P \\ = [\begin{matrix} s i g n (λ_{1}) \\ ⋱ \\ s i g n (λ_{p}) \\ s i g n (λ_{p + 1}) \\ ⋱ \\ s i g n (λ_{n}) \end{matrix}] \end{matrix}

(24)

Now, define a sequence

D_{l} = P^{- 1} W_{l} P

with

D_{0} = P^{- 1} A P

. Then from (18), we get

D_{l + 1} = [D_{l} (21 I + 50 D_{l}^{2} + 9 D_{l}^{4})] {[4 I + 45 D_{l}^{2} + 30 D_{l}^{4} + D_{l}^{6}]}^{- 1} .

(25)

Mathematical induction states that all

D_{l}

’s for

l = 1, 2, 3 \dots

will be diagonal if

D_{0}

is diagonal. When

D_{0}

is not diagonal, the proof may be considered in the similar methodology as stated at the end of the proof.

Now, (25) can be written for scalar case

f (w) = w^{2} - 1 = 0

as follows:

d_{l + 1}^{i} = \frac{d_{l}^{i} (21 + 50 {d_{l}^{i}}^{2} + 9 {d_{l}^{i}}^{4})}{4 + 45 {d_{l}^{i}}^{2} + 30 {d_{l}^{i}}^{4} + {d_{l}^{i}}^{6}},

(26)

where

d_{l}^{i} = {(D_{l})}_{i, i}

and

1 \leq i \leq n

. Furthermore,

s i g n (λ_{i}) = s_{i}

, with

s_{i} = \pm 1

and therefore from (26),

\frac{d_{l + 1}^{i} - s_{i}}{d_{l + 1}^{i} + s_{i}} = {(\frac{- s_{i} + d_{l}^{i}}{s_{i} + d_{l}^{i}})}^{5} (\frac{4 s_{i} - d_{l}^{i}}{4 s_{i} + d_{l}^{i}}) .

(27)

The factor

\frac{4 s_{i} - d_{l}^{i}}{4 s_{i} + d_{l}^{i}}

is bounded for

1 \leq i \leq n

and does not affect the convergence. Furthermore, by selecting the right initial matrix

W_{0} = A

and

|\frac{d_{0}^{i} - s_{i}}{d_{0}^{i} + s_{i}}| < 1

, we get

lim_{l \to \infty} |\frac{d_{l + 1}^{i} - s_{i}}{d_{l + 1}^{i} + s_{i}}| = 0,

(28)

implying that

{lim}_{l \to \infty} (d_{l}^{i}) = s_{i} = s i g n (λ_{i})

and hence

{lim}_{l \to \infty} (D_{l}) = s i g n (Λ)

In the case when

D_{0}

is not a diagonal matrix, it is essential that we investigate the relation between the eigenvalues of

W_{i}

’s which was briefly explained at the starting of the proof.

λ_{l + 1}^{i} = [λ_{l}^{i} (21 I + 50 {λ_{l}^{i}}^{2} + 9 {λ_{l}^{i}}^{4})] {[4 I + 45 {λ_{l}^{i}}^{2} + 30 {λ_{l}^{i}}^{4} + {λ_{l}^{i}}^{6}]}^{- 1} .

(29)

In a comparable approach, it is clear from (29) that

λ_{i}

’s will converge to

s_{i} = \pm 1

, i.e.,

lim_{l \to \infty} |\frac{λ_{l + 1}^{i} - s_{i}}{λ_{l + 1}^{i} + s_{i}}| = 0 .

(30)

As a result, it is easy to conclude that

lim_{l \to \infty} W_{l} = P (lim_{l \to \infty} D_{l}) P^{- 1} = P s i g n (Λ) P^{- 1} = s i g n (A),

(31)

which finishes the proof. □

Theorem 2.

Considering the assumptions of Theorem 1, the matrix sequence

{W_{l}}_{l = 0}^{\infty}

generated from (18) has convergence rate five to

s i g n (A) = S

Proof.

Since,

W_{l}

depends on matrix A∀

l \geq 0

and the matrix A commutes with S, so

W_{l}

also commutes with S∀

l \geq 0

. Taking into consideration the substitution

B_{l} = 4 I + 45 W_{l}^{2} + 30 W_{l}^{4} + W_{l}^{6},

(32)

we get,

\begin{matrix} W_{l + 1} - S & = [W_{l} (21 I + 50 W_{l}^{2} + 9 W_{l}^{4})] B_{l}^{- 1} - S \\ = [21 W_{l} + 50 W_{l}^{3} + 9 W_{l}^{5} - S B_{l}] B_{l}^{- 1} \\ = [21 W_{l} + 50 W_{l}^{3} + 9 W_{l}^{5} - 4 S - 45 S W_{l}^{2} - 30 S W_{l}^{4} - S W_{l}^{6}] B_{l}^{- 1} \\ = [- 4 S^{5} + 21 S^{4} W_{l} - 45 S^{3} W_{l}^{2} + 50 S^{2} W_{l}^{3} - 30 S W_{l}^{4} + 9 W_{l}^{5} - S W_{l}^{6}] B_{l}^{- 1} \\ = {(W_{l} - S)}^{5} (4 I - S W_{l}) B_{l}^{- 1} \end{matrix}

(33)

Apply any type of matrix norm to the both sides of (33),

∥W_{l + 1} - S∥ \leq (∥4 I - S W_{l}∥ ∥B_{l}^{- 1}∥) {∥W_{l} - S∥}^{5} .

(34)

This completes the proof. □

3.2. Stability Issues

Theorem 3.

Considering the assumptions of Theorem 1, the sequence of matrix iterations

{W_{l}}_{l = 0}^{\infty}

generated from (18) is stable asymptotically.

Proof.

The iterate obtained from (18) is a function of A because of

W_{0} = A

, and hence commutes with A. Let

Δ_{l}

be the perturbation generated at

l^{t h}

iterate in (18). Then

{\tilde{W}}_{l} = W_{l} + Δ_{l} .

(35)

Based on the results of the error analysis of order one, we may consider that

{(Δ_{l})}^{j} \approx 0

j \geq 2

. As long as

Δ_{l}

is small enough, this usual adjustment will work. Also suppose that

W_{l} \approx s i g n (A) = S

for larger l, then we get

\begin{matrix} {\tilde{W}}_{l + 1} & = [{\tilde{W}}_{l} (21 I + 50 {\tilde{W}}_{l}^{2} + 9 {\tilde{W}}_{l}^{4})] {[4 I + 45 {\tilde{W}}_{l}^{2} + 30 {\tilde{W}}_{l}^{4} + {\tilde{W}}_{l}^{6}]}^{- 1} \\ = [(W_{l} + Δ_{l}) (21 I + 50 {(W_{l} + Δ_{l})}^{2} + 9 {(W_{l} + Δ_{l})}^{4})] [4 I + 45 {(W_{l} + Δ_{l})}^{2} \\ + 30 {(W_{l} + Δ_{l})}^{4} + {(W_{l} + Δ_{l})}^{6}]^{- 1} \\ \approx [80 S + 148 Δ_{l} + 68 S Δ_{l} S] {[80 I + 108 S Δ_{l} + 108 Δ_{l} S]}^{- 1} \\ \approx [S + \frac{37}{20} Δ_{l} + \frac{17}{20} S Δ_{l} S] [I - \frac{27}{20} S Δ_{l} - \frac{27}{20} Δ_{l} S] \\ \approx S + \frac{1}{2} Δ_{l} - \frac{1}{2} S Δ_{l} S . \end{matrix}

(36)

Here, the identity [18] given by

{(B + C)}^{- 1} = B^{- 1} - B^{- 1} C B^{- 1} + O (C^{2}) .

(37)

has been applied, where B is any nonsingular matrix and C is any matrix. Also, considering

Δ_{l + 1} = {\tilde{W}}_{l + 1} - W_{l + 1}

, we get

Δ_{l + 1} \approx \frac{1}{2} Δ_{l} - \frac{1}{2} S Δ_{l} S .

(38)

At this point, it is possible to deduce that

Δ_{l + 1}

is bounded, i.e.,

∥Δ_{l + 1}∥ \leq \frac{1}{2} ∥Δ_{0} - S Δ_{0} S∥,

(39)

which completes the proof. □

4. Extension to Determine Generalized Eigenvalues of a Matrix Pencil

Let

A, B \in C^{n \times n}

. Then

λ \in C

is said to be a generalized eigenvalue of a matrix pencil

A - λ B

if there exists a non-zero vector

x \in C^{n}

such that

A x = λ B x .

(40)

Here

x

will be a generalized eigenvector corresponding to

λ

and (40) is referred as a generalized eigenvalue problem. This can be determined by finding the zeros of characteristic polynomial

d e t (A - λ B)

Complex problem solving in physics was the first scenario for the appearance of the characteristic polynomial and the eigenproblem. For a recent coverage of the topic one can see [19]. When more generality is needed, generalized eigenvalue problem is taking the place of the classical approach [20].

Problem (40) arises in several fields of numerical linear algebra. However, the dimensions of matrices A and B are typically quite enormous, hence further complicating the situation. There are a lot of different approaches that claim to be able to solve (40) including the necessary eigenvectors in a smaller subspace or series of matrix pencils [21,22,23] or using matrix decomposition [24,25,26]. Our interest is in employing the matrix sign function which involves evaluating certain matrix sign computations based on the given matrix pencil to compute the required decomposition and hence to solve the problem (40).

To solve the problem (40) for the case of regular matrix pencils using matrix sign function, we consider rank revealing

Q R

(R R Q R)

decomposition [27] after computing sign of two matrices depending on the given matrix pencil.

Theorem 4.

Let

A, B \in C^{n \times n}

such that

A - λ B

be a regular matrix pencil and let

U_{1}

and

U_{2}

are defined as

\begin{matrix} U_{1} = \frac{I - S_{1}}{2}, U_{2} = \frac{I - S_{2}}{2}, \end{matrix}

(41)

where

S_{1}

and

S_{2}

are sign functions given by

\begin{matrix} S_{1} & = s i g n ({(A - B)}^{- 1} (A + B)), \\ S_{2} & = s i g n ((A + B) {(A - B)}^{- 1}) . \end{matrix}

(42)

U_{1} Π_{1} = Q_{1} R_{1}

and

U_{2} Π_{2} = Q_{2} R_{2}

are

R R Q R

decompositions of

U_{1}

and

U_{2}

such that

Q_{2}^{*} A Q_{1}

and

Q_{2}^{*} B Q_{1}

are upper triangular, then

λ_{i} = \frac{α_{i i}}{β_{i i}},

(43)

where

λ_{i}

’s are eigenvalues of

A - λ B

α_{i i}

and

β_{i i}

are the diagonal entries of

Q_{2}^{*} A Q_{1}

and

Q_{2}^{*} B Q_{1}

respectively for

i = 1, 2, \dots, n

. Moreover, if B is non-invertible, then

β_{i i} = 0

for some i which implies that some eigenvalues of

A - λ B

will be infinity.

Let

τ

denotes tolerance value and

m a x

be the maximum iterations. Then, in accordance with the proposed method (18) and Theorem 4, we propose Algorithm 1 to find out the eigenvalues of a regular matrix pencil

A - λ B

Algorithm 1 Evaluating generalized eigenvalues of

A - λ B

Require:

A, B, m a x, τ

Ensure:

S_{1}, S_{2}, λ_{i}^{'} s

1:: Initialize $W_{0} = {(A - B)}^{- 1} (A + B)$ .
2:: while $l \leq m a x$ and $∥ W_{l}^{2} {- I ∥}_{*} > τ$ do
3:: Apply method (18).
4:: return $W_{\infty} = s i g n ({(A - B)}^{- 1} (A + B)) = S_{1}$ .
5:: end while
6:: Construct $R R Q R$ decomposition of $U_{1} = \frac{I - S_{1}}{2}$ as follows: $U_{1} Π_{1} = Q_{1} R_{1}$ .
7:: Initialize $W_{0} = (A + B) {(A - B)}^{- 1}$ .
8:: Repeat step 2–5 to return $W_{\infty} = s i g n ((A + B) {(A - B)}^{- 1}) = S_{2}$ .
9:: Construct $R R Q R$ decomposition of $U_{2} = \frac{I - S_{2}}{2}$ as follows: $U_{2} Π_{2} = Q_{2} R_{2}$ .
10:: Compute two matrices $Q_{2}^{*} A Q_{1} = [α_{i j}]$ and $Q_{2}^{*} B Q_{1} = [β_{i j}]$ , where $1 \leq i, j \leq n$ .
11:: if $α_{i j} = β_{i j} = 0$ ∀ $i > j$ then
12:: for i=1:n do
13:: $λ_{i} = \frac{α_{i i}}{β_{i i}}$ .
14:: end for
15:: end if

5. Computational Complexity

The computational efficiency index (

C E I

) [28] of an iterative method can be defined as follows:

C E I = p^{\frac{1}{C}} .

(44)

Here p denotes convergence order and C denotes computational cost per iteration.

The approximate cost of matrix multiplications in current programming packages is

2.373

[29]. For a matrix of size n, we therefore consider that the cost of one matrix-matrix multiplication is

n^{\frac{2373}{1000}}

and cost of one matrix inverse is

\frac{2 n^{3}}{3}

. Additionally, we do not take into account the cost per iteration of norm computation (in

A N M

(6)), addition and subtraction. Then,

C E I

’s for the methods

N M

(5),

A N M

(6),

H M

(12) and

P M

(18) to compute

s i g n (A)

are given in Table 1.

Figure 2 presents a comparison of the efficiencies of several approaches, demonstrating the greater efficiency index of the suggested method.

6. Numerical Aspects

In this section, we have examined a variety of examples to evaluate our proposed method by comparing it with various existing techniques. For this study, we have assumed the following stopping criteria:

∥ W_{l}^{2} {- I ∥}_{*} \leq τ,

(45)

where

τ

stands for the tolerance value and

{∥ . ∥}_{*}

will be suitable matrix norm. Also,

l_{\infty}

and

l_{2}

norms are considered for real and complex input matrices, respectively [14].

For the sake of comparison, we only consider the methods that exhibit global convergence and we do not include any methods that have local convergence. The methods that are being compared will be

N M

(5),

A N M

(6),

H M

(12),

P A_{2, 2}

((10) for

r = s = 2

L M

(13),

P M

(18) and

A P M

(20) with spectral scaling (19). All simulations are done in Mathematica using system specifications “Windows 11 Home Intel(R) Core(TM) i7-1255U CPU running at 1.70 GHz processor having 16.0 RAM (64-bit operating system)”.

Example 1.

Here, we perform the eigenvalue clustering of a random real

500 \times 500

matrix given as

SeedRandom[100]; A = RandomReal[{-5,5},{500,500}];

The outcomes are shown in Figure 3 using the stopping criteria (45) with

l_{\infty}

norm and

τ = 10^{- 4}

. We can see that eigenvalues of

W_{0}

are scattered at the initial stage as shown in Figure 3a, while proceeding for the next iterations

W_{i}

in Figure 3b–g for

i = 1, 2, \dots, 6

respectively using

P M

(18), the eigenvalues are converging towards 1 or

- 1

. Then in Figure 3h, all eigenvalues converge to 1 or

- 1

(based on the stopping criteria) implying that the sequence of matrices

{W_{l}}_{l = 0}^{\infty}

converges to

s i g n (A)

. The result of Theorem 1 is therefore confirmed.

Example 2.

In this example, convergence of many different methods for computing

s i g n (A)

has been compared, where A is a random

1000 \times 1000

complex square matrix given as follows:

SeedRandom[101];

A = RandomComplex[{-1 - 1.5 I, 1 + 1.5 I}, {1000, 1000}]

The comparison is shown in Figure 4 with

τ = 10^{- 4}

and

l_{2}

norm which illustrates the results with regard to the number of iterations and the absolute error and demonstrates the consistent behavior of both

P M

and

A P M

Example 3.

This example presents a comparison on convergence of many different methods for computing

s i g n (A)

, where A is a random

2000 \times 2000

real square matrix given as

SeedRandom[111];

A = RandomReal[{-25, 25}, {2000, 2000}];

The comparison is presented in Figure 5 with

τ = 10^{- 4}

and

l_{\infty}

norm which illustrates the results with respect to the number of iterations and the absolute error and demonstrates the consistent behavior of both

P M

and

A P M

for evaluating

s i g n (A)

Example 4.

Here, we examine the convergence of various methods with regard to the iteration number and the timings of CPU for evaluating the sign function of ten random matrices:

SeedRandom[121]; number = 10; tolerance = 10^(-5);

Table[A[j] = RandomComplex[{-3 - 2 I, 3 + 2 I}, {50j, 50j}];, {j, number}];

The comparative analysis based on iteration number and CPU timings for matrices of various sizes, are presented in Table 2 and Table 3 respectively using

τ = 10^{- 5}

and

l_{2}

-norm. The implementation of new methods

(P M)

and

(A P M)

result in an improvement in the average iteration number and the amount of time spent by CPU. The outcomes demonstrate that the suggested method

(P M)

is significantly better as compared to other fifth order methods like

P A_{2, 2}

and

L M

Now, in Examples 5 and 6, we aim to evaluate generalized eigenvalues of matrix pencils using matrix sign function for which we use Algorithm 1 for different methods by replacing (18) with existing methods in Step 3.

Example 5.

Let

A - λ B

be a matrix pencil, where matrices A and B are given by [30]

A = {[\begin{matrix} \frac{99}{100} & \frac{1}{100} & 0 & \dots & 0 \\ 0 & \frac{98}{100} & ⋱ & 0 \\ ⋮ & ⋱ & ⋱ & ⋱ & ⋮ \\ 0 & \dots & 0 & \frac{1}{100} & \frac{1}{100} \\ 0 & \dots & 0 & 0 & \frac{0}{100} \end{matrix}]}_{100 \times 100}, B = {[\begin{matrix} 0 \\ ⋱ \\ 0 \\ I_{20} \end{matrix}]}_{100 \times 100},

where

I_{20}

is the identity matrix of size 20.

The eigenvalues of given pencil

A - λ B

are

λ_{i} = \frac{i - 1}{100}

for

i = 1, 2, \dots 20

and remaining eigenvalues are infinite. We verify these eigenvalues using proposed method and Algorithm 1. We compute the matrices

Q_{2}^{*} A Q_{1} = {[α_{i j}]}_{100 \times 100}

and

Q_{2}^{*} B Q_{1} = {[β_{i j}]}_{100 \times 100}

and then calculate

λ_{i}^{'} s

from Algorithm 1. The outcomes are presented in Table 4 for

i = 1, 2, \dots 20

Moreover,

β_{i i} = 0

for

21 \leq i \leq 100

which implies that

λ_{i} = \infty

for

21 \leq i \leq 100

. Hence this example verifies the results given in Theorem 4.

Example 6 (Bounded finline dielectric waveguide problem).

The generalized eigenvalue problem (40) arises in the finite element analysis for finding the propagating modes of a rectangular waveguide. We consider three cases for the pairs of matrices A and B as

B F W 62 A, B F W 62 B

;

B F W 398 A, B F W 398 B

and

B F W 782 A, B F W 782 B

from the Matrix Market [31].

Table 5 shows the outcomes depending on the iteration numbers

l_{1}

and

l_{2}

for computing

S_{1} = s i g n ({(A - B)}^{- 1} (A + B))

and

S_{2} = s i g n ((A + B) {(A - B)}^{- 1})

with corresponding absolute errors

ϵ_{l_{1}} = {∥ W_{l_{1}}^{2} - I ∥}_{\infty}

and

ϵ_{l_{2}} = {∥ W_{l_{2}}^{2} - I ∥}_{\infty}

using stopping criterion (45) with

τ = 10^{- 10}

. Also, total computational time in seconds taken for computing

S_{1}

and

S_{2}

has been compared in Figure 6. The outcomes demonstrate that the suggested method is better with respect to errors and timings.

7. Conclusions

This paper has discussed an iterative approach for evaluation of sign of a matrix. Our proposed iterative method attains a convergence order of five, showcasing its efficiency in evaluation of the sign of complex matrices. Through comprehensive convergence analysis discussed globally and the demonstration of asymptotic stability, we have established the robustness and reliability of the proposed approach. Then an algorithm has been designed for determining the generalized eigenvalues of a regular matrix pencil utilizing matrix sign computation. Furthermore, extensive numerical experiments performed on matrices of different sizes have validated the effectiveness and superiority of our method over existing techniques.

Author Contributions

Conceptualization, M.K.; methodology, M.K.; software, P.S.; validation, M.K. and V.S.; formal analysis, L.J.; investigation, P.S.; resources, L.J.; data curation, L.J.; writing—original draft preparation, P.S.; writing—review and editing, L.J.; visualization, V.S.; supervision, M.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors would like to sincerely thank the editor and anonymous reviewers for their constructive and valuable suggestions that have helped in improving the quality of the manuscript.

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

Roberts, J.D. Linear model reduction and solution of the algebraic Riccati equation by use of the sign function. Int. J. Control 1980, 32, 677–687. [Google Scholar] [CrossRef]
Higham, N.J. Functions of Matrices: Theory and Computation; SIAM: Philadelphia, PA, USA, 2008. [Google Scholar] [CrossRef]
Kenney, C.; Laub, A.J. Rational iterative methods for the matrix sign function. SIAM J. Matrix Anal. Appl. 1991, 12, 273–291. [Google Scholar] [CrossRef]
Kenney, C.S.; Laub, A.J. The matrix sign function. IEEE Trans. Autom. Control 1995, 40, 1330–1348. [Google Scholar] [CrossRef]
Gomilko, O.; Greco, F.; Ziȩtak, K. A Padé family of iterations for the matrix sign function and related problems. Numer. Linear Algebra Appl. 2012, 19, 585–605. [Google Scholar] [CrossRef]
Schulz, G. Iterative berechung der reziproken matrix. J. Appl. Math. Mech. Angew. Math. Mech. 1933, 13, 57–59. [Google Scholar] [CrossRef]
Momenzadeh, M.; Lotfi, T. New iterative methods for finding matrix sign function: Derivation and application. Comput. Appl. Math. 2019, 38, 53. [Google Scholar] [CrossRef]
Ullah, M.Z.; Alaslani, S.M.; Mallawi, F.O.; Ahmad, F.; Shateyi, S.; Asma, M. A fast and efficient Newton-type iterative scheme to find the sign of a matrix. AIMS Math. 2023, 8, 19264–19274. [Google Scholar] [CrossRef]
Wang, X.; Cao, Y. A numerically stable high-order Chebyshev-Halley type multipoint iterative method for calculating matrix sign function. AIMS Math. 2023, 8, 12456–12471. [Google Scholar] [CrossRef]
Liu, T.; Ullah, M.Z.; Alshahrani, K.M.A.; Shateyi, S. From fractal behavior of iteration methods to an efficient solver for the sign of a matrix. Fractal Fract. 2023, 7, 32. [Google Scholar] [CrossRef]
Benner, P.; Quintana-Ortí, E.S. Solving stable generalized Lyapunov equations with the matrix sign function. Numer. Algorithms 1999, 20, 75–100. [Google Scholar] [CrossRef]
Howland, J.L. The sign matrix and the separation of matrix eigenvalues. Linear Algebra Its Appl. 1983, 49, 221–232. [Google Scholar] [CrossRef]
Soleymani, F.; Kumar, A. A fourth-order method for computing the sign function of a matrix with application in the Yang–Baxter-like matrix equation. Comput. Appl. Math. 2019, 38, 64. [Google Scholar] [CrossRef]
Soheili, A.R.; Toutounian, F.; Soleymani, F. A fast convergent numerical method for matrix sign function with application in SDEs. J. Comput. Appl. Math. 2015, 282, 167–178. [Google Scholar] [CrossRef]
Kenney, C.; Laub, A.J. On scaling Newton’s method for polar decomposition and the matrix sign function. SIAM J. Matrix Anal. Appl. 1992, 13, 688–706. [Google Scholar] [CrossRef]
Soleymani, F.; Stanimirović, P.S.; Stojanović, I. A novel iterative method for polar decomposition and matrix sign function. Discrete Dyn. Nat. Soc. 2015, 2015, 649423. [Google Scholar] [CrossRef]
Iannazzo, B. Numerical Solution of Certain Nonlinear Matrix Equations. Ph.D. Thesis, Dipartimento di Matematica, Università di Pisa, Pisa, Italy, 2007. [Google Scholar]
Henderson, H.V.; Searle, S.R. On deriving the inverse of a sum of matrices. SIAM Rev. 1981, 23, 53–60. [Google Scholar] [CrossRef]
Jäntschi, L. Eigenproblem Basics and Algorithms. Symmetry 2023, 15, 2046. [Google Scholar] [CrossRef]
Fischer, M.; Kostrzewa, B.; Ostmeyer, J.; Ottnad, K.; Ueding, M.; Urbach, C. On the generalised eigenvalue method and its relation to Prony and generalised pencil of function methods. Eur. Phys. J. A 2020, 56, 206. [Google Scholar] [CrossRef]
Ruhe, A. Rational Krylov sequence methods for eigenvalue computation. Linear Algebra Appl. 1984, 58, 391–405. [Google Scholar] [CrossRef]
Gallivan, K.; Grimme, E.; Van Dooren, P. Asymptotic waveform evaluation via a Lanczos method. Appl. Math. Lett. 1994, 7, 75–80. [Google Scholar] [CrossRef]
Gallivan, K.; Grimme, E.; Van Dooren, P. A rational Lanczos algorithm for model reduction. Numer. Algor. 1996, 12, 33–63. [Google Scholar] [CrossRef]
Moler, C.B.; Stewart, G.W. An algorithm for generalized matrix eigenvalue problems. SIAM J. Num. Anal. 1973, 10, 241–256. [Google Scholar] [CrossRef]
Bai, Z.; Demmel, J. Using the matrix sign function to compute invariant subspaces. SIAM J. Matrix Anal. Appl. 1998, 19, 205–225. [Google Scholar] [CrossRef]
Sun, X.; Quintana-Ortí, E. Spectral division methods for block generalized Schur decompositions. Math. Comput. 2004, 73, 1827–1847. [Google Scholar] [CrossRef]
Quintana-Ortí, G.; Sun, X.; Bischof, C.H. A BLAS-3 Version of the QR Factorization with Column Pivoting. SIAM J. Sci. Comp. 1998, 19, 1486–1494. [Google Scholar] [CrossRef]
Soleymani, F.; Khaksar Haghani, F.; Shateyi, S. Several numerical methods for computing unitary polar factor of a matrix. Adv. Differ. Equ. 2016, 2016, 4. [Google Scholar] [CrossRef]
Le Gall, F. Powers of tensors and fast matrix multiplication. In Proceedings of the 39th International Symposium on Symbolic and Algebraic Computation, Kobe, Japan, 23–25 July 2014; ACM: New York, NY, USA, 2014; pp. 296–303. [Google Scholar] [CrossRef]
Sakurai, T.; Sugiura, H. A projection method for generalized eigenvalue problems using numerical integration. J. Comput. Appl. Math. 2003, 159, 119–128. [Google Scholar] [CrossRef]
Matrix Market. Available online: https://math.nist.gov/MatrixMarket/ (accessed on 15 January 2024).

Figure 1. Attraction basins of various methods for equation

w^{2} - 1 = 0

Figure 1. Attraction basins of various methods for equation

w^{2} - 1 = 0

Figure 2. Comparative analysis of efficiency indices for different methods.

Figure 3. Eigenvalue clustering of

500 \times 500

matrix at every iterate in Example 1.

Figure 3. Eigenvalue clustering of

500 \times 500

matrix at every iterate in Example 1.

Figure 4. Convergence history of different methods for Example 2.

Figure 5. Convergence history of different methods for Example 3.

Figure 6. Comparison of CPU time for various methods to compute

S_{1}

and

S_{2}

in Example 6.

Figure 6. Comparison of CPU time for various methods to compute

S_{1}

and

S_{2}

in Example 6.

Table 1.

C E I

’s for various methods along with their convergence orders.

Table 1.

C E I

’s for various methods along with their convergence orders.

Method	Convergence Order	$CEI$
$N M$	2	$2^{\frac{1}{n^{2373 / 1000} + \frac{2 n^{3}}{3}}}$
$A N M$	$\sqrt{3} + 1$	${(\sqrt{3} + 1)}^{\frac{1}{n^{2373 / 1000} + \frac{2 n^{3}}{3}}}$
$H M$	3	$3^{\frac{1}{3 n^{2373 / 1000} + \frac{2 n^{3}}{3}}}$
$P M$	5	$5^{\frac{1}{5 n^{2373 / 1000} + \frac{2 n^{3}}{3}}}$

Table 2. Comparative analysis for the iteration number in Example 4.

Matrix	$NM$	$ANM$	$HM$	$LM$	${PA}_{2, 2}$	$PM$	$APM$
$A_{50 \times 50}$	12	11	8	6	6	5	4
$A_{100 \times 100}$	15	13	9	7	7	6	5
$A_{150 \times 150}$	15	15	10	8	7	6	5
$A_{200 \times 200}$	15	15	9	8	7	6	5
$A_{250 \times 250}$	15	17	10	8	7	6	5
$A_{300 \times 300}$	14	15	9	7	6	6	5
$A_{350 \times 350}$	15	16	10	8	7	6	5
$A_{400 \times 400}$	17	19	11	10	7	7	5
$A_{450 \times 450}$	14	15	9	7	6	6	5
$A_{500 \times 500}$	18	21	12	9	8	7	6
Mean	15	15.7	9.7	7.8	6.8	6.1	5

Table 3. Comparative analysis for the CPU-time (in seconds) in Example 4.

Matrix	$NM$	$ANM$	$HM$	$LM$	${PA}_{2, 2}$	$PM$	$APM$
$A_{50 \times 50}$	0.03	0.02	0.01	0.01	0.02	0.01	0.01
$A_{100 \times 100}$	0.07	0.15	0.05	0.05	0.09	0.05	0.10
$A_{150 \times 150}$	0.16	0.23	0.12	0.15	0.21	0.10	0.24
$A_{200 \times 200}$	0.29	0.40	0.19	0.26	0.34	0.18	0.40
$A_{250 \times 250}$	0.47	0.67	0.32	0.39	0.53	0.25	0.60
$A_{300 \times 300}$	0.58	0.87	0.44	0.55	0.63	0.44	1.02
$A_{350 \times 350}$	0.82	1.29	0.70	0.88	1.21	0.63	1.44
$A_{400 \times 400}$	1.20	2.04	1.46	1.78	1.62	0.99	1.70
$A_{450 \times 450}$	1.31	2.06	1.08	1.40	1.87	1.16	2.29
$A_{500 \times 500}$	2.06	3.58	1.91	2.32	3.07	1.60	3.32
Mean	0.70	1.13	0.63	0.78	0.96	0.54	1.11

Table 4. Eigenvalues (

λ_{i}

’s) of matrix pencil

A - λ B

in Example 5.

Table 4. Eigenvalues (

λ_{i}

’s) of matrix pencil

A - λ B

in Example 5.

i	1	2	3	4	5	6	7	8	9	10
$α_{i i}$	0.189762	0.18	−0.17	0.16	0.15	−0.14	0.13	−0.12	−0.11	0.1
$β_{i i}$	0.99875	1	−1	1	1	−1	1	−1	−1	1
$λ_{i}$	0.19	0.18	0.17	0.16	0.15	0.14	0.13	0.12	0.11	0.1
$i$	11	12	13	14	15	16	17	18	19	20
$α_{i i}$	−0.09	0.08	−0.07	0.05	−0.03	0.02	0.01	0	−0.06	0.04
$β_{i i}$	−1	1	−1	1	−1	1	1	−1	−1	1
$λ_{i}$	0.09	0.08	0.07	0.05	0.03	0.02	0.01	0	0.06	0.04

Table 5. Comparison of results obtained from various methods in Example 6.

Matrices	Method	$l_{1}$	$ϵ_{l_{1}}$	$l_{2}$	$ϵ_{l_{2}}$
BFW62A, BFW62B	$N M$	2	$5.79607 \times 10^{- 11}$	2	$5.79607 \times 10^{- 11}$
	$H M$	2	$4.44089 \times 10^{- 16}$	2	$4.44089 \times 10^{- 16}$
	$P A_{2, 2}$	1	$1.65201 \times 10^{- 13}$	1	$1.64757 \times 10^{- 13}$
	$L M$	1	$9.80327 \times 10^{- 13}$	1	$9.81437 \times 10^{- 13}$
	$P M$	1	$9.90319 \times 10^{- 14}$	1	$1.00364 \times 10^{- 13}$
BFW398A, BFW398B	$N M$	2	$6.18883 \times 10^{- 12}$	2	$6.18883 \times 10^{- 12}$
	$H M$	2	$4.44089 \times 10^{- 16}$	2	$4.44089 \times 10^{- 16}$
	$P A_{2, 2}$	1	$1.50990 \times 10^{- 14}$	1	$1.50990 \times 10^{- 14}$
	$L M$	1	$9.08162 \times 10^{- 14}$	1	$9.23706 \times 10^{- 14}$
	$P M$	1	$1.11022 \times 10^{- 14}$	1	$1.11022 \times 10^{- 14}$
BFW782A, BFW782B	$N M$	2	$5.82201 \times 10^{- 13}$	2	$5.82201 \times 10^{- 13}$
	$H M$	2	$4.44089 \times 10^{- 16}$	2	$4.44089 \times 10^{- 16}$
	$P A_{2, 2}$	1	$6.43929 \times 10^{- 15}$	1	$3.99680 \times 10^{- 15}$
	$L M$	1	$7.99361 \times 10^{- 15}$	1	$7.99361 \times 10^{- 15}$
	$P M$	1	$5.32907 \times 10^{- 16}$	1	$5.32907 \times 10^{- 16}$

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

Share and Cite

MDPI and ACS Style

Kansal, M.; Sharma, V.; Sharma, P.; Jäntschi, L. A Globally Convergent Iterative Method for Matrix Sign Function and Its Application for Determining the Eigenvalues of a Matrix Pencil. Symmetry 2024, 16, 481. https://doi.org/10.3390/sym16040481

AMA Style

Kansal M, Sharma V, Sharma P, Jäntschi L. A Globally Convergent Iterative Method for Matrix Sign Function and Its Application for Determining the Eigenvalues of a Matrix Pencil. Symmetry. 2024; 16(4):481. https://doi.org/10.3390/sym16040481

Chicago/Turabian Style

Kansal, Munish, Vanita Sharma, Pallvi Sharma, and Lorentz Jäntschi. 2024. "A Globally Convergent Iterative Method for Matrix Sign Function and Its Application for Determining the Eigenvalues of a Matrix Pencil" Symmetry 16, no. 4: 481. https://doi.org/10.3390/sym16040481

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Menu

A Globally Convergent Iterative Method for Matrix Sign Function and Its Application for Determining the Eigenvalues of a Matrix Pencil

Abstract

1. Introduction

2. A Novel Iterative Method

3. Convergence and Stability Analysis

3.1. Convergence Analysis

3.2. Stability Issues

4. Extension to Determine Generalized Eigenvalues of a Matrix Pencil

5. Computational Complexity

6. Numerical Aspects

7. Conclusions

Author Contributions

Funding

Data Availability Statement

Acknowledgments

Conflicts of Interest

References

Share and Cite

Article Metrics

Article Access Statistics

Further Information

Guidelines

MDPI Initiatives

Follow MDPI