Premium
Computing subdominant unstable modes of turbulent plasma with a parallel Jacobi–Davidson eigensolver
Author(s) -
Romero Eloy,
Roman Jose E.
Publication year - 2011
Publication title -
concurrency and computation: practice and experience
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.309
H-Index - 67
eISSN - 1532-0634
pISSN - 1532-0626
DOI - 10.1002/cpe.1740
Subject(s) - preconditioner , parallel computing , eigenvalues and eigenvectors , context (archaeology) , computer science , scalability , mathematics , factorization , computation , supercomputer , computational science , mathematical optimization , iterative method , algebra over a field , algorithm , physics , pure mathematics , paleontology , quantum mechanics , database , biology
SUMMARY In the numerical solution of large‐scale eigenvalue problems, Davidson‐type methods are an increasingly popular alternative to Krylov eigensolvers. The main motivation is to avoid the expensive factorizations that are often needed by Krylov solvers when the problem is generalized or interior eigenvalues are desired. In Davidson‐type methods, the factorization is replaced by iterative linear solvers that can be accelerated by a smart preconditioner. Jacobi–Davidson is one of the most effective variants. However, parallel implementations of this method are not widely available, particularly for non‐symmetric problems. We present a parallel implementation that has been included in SLEPc, the Scalable Library for Eigenvalue Problem Computations, and test it in the context of a highly scalable plasma turbulence simulation code. We analyze its parallel efficiency and compare it with a Krylov–Schur eigensolver. Copyright © 2011 John Wiley & Sons, Ltd.