Source author record

Hong-kui Pang

Hong-kui Pang appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.NA

Catalog footprint

What is connected

2works

1topics

2close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2015arXiv

On the correction equation of the Jacobi-Davidson method

The Jacobi-Davidson method is one of the most popular approaches for iteratively computing a few eigenvalues and their associated eigenvectors of a large matrix. The key of this method is to expand the search subspace via solving the Jacobi-Davidson correction equation, whose coefficient matrix is singular. It is believed long by scholars that the Jacobi-Davidson correction equation is a consistent linear system. In this work, we point out that the correction equation may have a unique solution or have no solution at all, and we derive a computable necessary and sufficient condition for cheaply judging the existence and uniqueness of solution of the correction equation. Furthermore, we consider the difficulty of stagnation that bothers the Jacobi-Davidson method, and verify that if the Jacobi-Davidson method stagnates, then the corresponding Ritz value is a defective eigenvalue of the projection matrix. We provide a computable necessary and sufficient condition for expanding the search subspace successfully. The properties of the Jacobi-Davidson method with preconditioning and some alternative Jacobi-Davidson correction equations are also discussed.

preprint2014arXiv

Preconditioning the Restarted and Shifted Block FOM Algorithm for Matrix Exponential Computation

The approximation of $e^{tA}B$ where $A$ is a large sparse matrix and $B$ a rectangular matrix is the key ingredient in many scientific and engineering computations. A powerful tool to manage the matrix exponential function is to resort to a suitable rational approximation such as the Carath$\acute{\rm e}$odory-Fej$\acute{\rm e}$r approximation, whose core reduces to solve shifted linear systems with multiple right-hand sides. The restarted and shifted block FOM algorithm is a commonly used technique for this problem. However, determining good preconditioners for shifted systems that preserve the original structure is a difficult task. In this paper, we propose a new preconditioner for the restarted and shifted block FOM algorithm. The key is that the absolute values of the poles of the Carath$\acute{\rm e}$odory-Fej$\acute{\rm e}$r approximation are medium sized and can be much smaller than the norm of the matrix in question. The advantages of the proposed strategy are that we can precondition all the shifted linear systems simultaneously, and preserve the original structure of the shifted linear systems after restarting. Theoretical results are provided to show the rationality of our preconditioning strategy. Applications of the new approach to Toeplitz matrix exponential problem are also discussed. Numerical experiments illustrate the superiority of the new algorithm over many state-of-the-art algorithms for matrix exponential.