English

Efficiently preconditioned Inexact Newton methods for large symmetric eigenvalue problems

Numerical Analysis 2013-12-06 v1

Abstract

In this paper we propose an efficiently preconditioned Newton method for the computation of the leftmost eigenpairs of large and sparse symmetric positive definite matrices. A sequence of preconditioners based on the BFGS update formula is proposed, for the Preconditioned Conjugate Gradient solution of the linearized Newton system to solve Au=q(u)uA \mathbf{u} = q(\mathbf{u}) \mathbf{u}, q(u)q(\mathbf{u}) being the Rayleigh Quotient. We give theoretical evidence that the sequence of preconditioned Jacobians remains close to the identity matrix if the initial preconditioned Jacobian is so. Numerical results onto matrices arising from various realistic problems with size up to one million unknowns account for the efficiency of the proposed algorithm which reveals competitive with the Jacobi-Davidson method on all the test problems.

Keywords

Cite

@article{arxiv.1312.1553,
  title  = {Efficiently preconditioned Inexact Newton methods for large symmetric eigenvalue problems},
  author = {Luca Bergamaschi and Angeles Martinez},
  journal= {arXiv preprint arXiv:1312.1553},
  year   = {2013}
}

Comments

20 pages, 6 figures. Submitted to Optimization Methods and Software

R2 v1 2026-06-22T02:21:37.230Z