Xiao-Qing Jin

Teng-Teng Yao, Zheng-Jian Bai, Xiao-Qing Jin et al.

This paper is concerned with the least squares inverse eigenvalue problem of reconstructing a linear parameterized real symmetric matrix from the prescribed partial eigenvalues in the sense of least squares, which was originally proposed by Chen and Chu [SIAM J. Numer. Anal., 33 (1996), pp. 2417--2430]. We provide a Riemannian inexact Gausss-Newton method for solving the least squares inverse eigenvalue problem. The global and local convergence analysis of the proposed method is discussed. Also, a preconditioned conjugate gradient method with an efficient preconditioner is proposed for solving the Riemannian Gauss-Newton equation. Finally, some numerical tests, including an application in the inverse Sturm-Liouville problem, are reported to illustrate the efficiency of the proposed method.

Chun-Yueh Chiang, Matthew M. Lin, Xiao-Qing Jin

Inverse eigenvalue and singular value problems have been widely discussed for decades. The well-known result is the Weyl-Horn condition, which presents the relations between the eigenvalues and singular values of an arbitrary matrix. This result by Weyl-Horn then leads to an interesting inverse problem, i.e., how to construct a matrix with desired eigenvalues and singular values. In this work, we do that and more. We propose an eclectic mix of techniques from differential geometry and the inexact Newton method for solving inverse eigenvalue and singular value problems as well as additional desired characteristics such as nonnegative entries, prescribed diagonal entries, and even predetermined entries. We show theoretically that our method converges globally and quadratically, and we provide numerical examples to demonstrate the robustness and accuracy of our proposed method. {Having theoretical interest, we provide in the appendix a necessary and sufficient condition for the existence of a $2\times 2$ real matrix, or even a nonnegative matrix, with prescribed eigenvalues, singular values, and main diagonal entries.

Zhi Zhao, Zheng-Jian Bai, Xiao-Qing Jin

This paper is concerned with the nonnegative inverse eigenvalue problem of finding a nonnegative matrix such that its spectrum is the prescribed self-conjugate set of complex numbers. We first reformulate the nonnegative inverse eigenvalue problem as an under-determined constrained nonlinear matrix equation over several matrix manifolds. Then we propose a Riemannian inexact Newton-CG method for solving the nonlinear matrix equation. The global and quadratic convergence of the proposed method is established under some mild conditions. We also extend the proposed method to the case of prescribed entries. Finally, numerical experiments are reported to illustrate the efficiency of the proposed method.

Zhi Zhao, Xiao-Qing Jin, Matthew M. Lin

In this paper we want to propose practical numerical methods to solve a class of initial-boundary problem of space-time fractional advection-diffusion equations. To start with, an implicit method based on two-sided Grünwald formulae is proposed with a discussion of the stability and consistency. Then, the preconditioned generalized minimal residual (preconditioned GMRES) method and the preconditioned conjugate gradient normal residual ({preconditioned} CGNR) method, with an easily constructed preconditioner, are developed. Importantly, because the resulting systems are Topelitz-like, the fast Fourier transform can be applied to significantly reduce the computational cost. Numerical experiments are implemented to show the efficiency of our preconditioner, even with cases of variable coefficients.