Parallel Arnoldi Method for the Construction of a Krylov Subspace Basis: AnApplication in Magnetohydrodynamics

摘要

In a recent article by Paige et al, a new method was proposed for computinginternal eigenvalues of symmetric matrices. In the present paper, we extend these ideas to non-hermitian eigenvalue problems and apply them to a practical example from the field of magnetohydrodynamics (MHD). The method is very suitable for an efficient parallel implementation. We give some results for the time-consuming kernels of the underlying orthogonalization process, the Arnoldi method, obtained for an MHD problem on a distributed memory multiprocessor.