Version 1.4.2
6th April 2022 User documentation
 Recent Changes

Code Download
 Single
 Double
 Single Complex
 Double Complex
HSL_EA19: Sparse symmetric or Hermitian: leftmost eigenpairs
HSL_EA19
uses a subspace iteration method to compute the leftmost eigenvalues and corresponding eigenvectors of a real symmetric (or Hermitian) operator \(\mathbf{A}\) acting in the \(n\)dimensional real (or complex) Euclidean space \(R^n\), or, more generally, of the problem \[\begin{aligned}
\mathbf{Ax} = \lambda \mathbf{Bx},\end{aligned}\] where \(\mathbf{B}\) a real symmetric (or Hermitian) positivedefinite operator. By applying HSL_EA19
to \(\mathbf{A}\), the user can compute the rightmost eigenvalues of \(\mathbf{A}\) and the corresponding eigenvectors. HSL_EA19
does not perform factorizations of \(\mathbf{A}\) or \(\mathbf{B}\) and thus is suitable for solving largescale problems for which a sparse direct solver for factorizing \(\mathbf{A}\) or \(\mathbf{B}\) is either not available or is too expensive.
The convergence may be accelerated by the provision of a symmetric positivedefinite operator \(\mathbf{T}\) that approximates the inverse of \((\mathbf{A}  \sigma \mathbf{B})\) for a value of \(\sigma\) that does not exceed the leftmost eigenvalue. Computation time may also be reduced by supplying vectors that are good approximations to some of the eigenvectors.