Newton's method in floating point arithmetic and iterative refinement of generalized eigenvalue problems

From MaRDI portal





In the symmetric generalized eigenvalue problem with indefinite \(A\) or \(B\) there exists a conflict between efficiency and backward stability. The paper advocates an iterative refinement of the computed eigenpairs by Newton's method. It is admitted that the linear solver is unstable and that the Jacobian itself is inaccurate. An extended precision is only assumed available in the computation of the residual. The key question which is analyzed concerns the sufficient precision of the residual. It is shown that and when an improvement of the forward and backward errors may be achieved in such a situation. Numerical illustrations confirm the theory using small matrices treated by the usual Cholesky-\(QR\) method.




Cited in
(25)


Describes a project that uses

Uses Software






This page was built for publication: Newton's method in floating point arithmetic and iterative refinement of generalized eigenvalue problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q2706305)