Jacobi-Davidson methods for generalized MHD eigenvalue problems

Abstract

A Jacobi-Davidson algorithm for computing selected eigenvalues and associated eigenvectors of the generalized eigenvalue problem Ax = Bx is presented. In this paper the emphasis is put on the case where one of the matrices, say the B-matrix, is Hermitian positive denite. The method is an inner-outer iterative scheme, in which the inner iteration process consists of solving linear systems to some accuracy. The factorization of either matrix is avoided. Numerical experiments are presented for problems arising in magnetohydrodynamics (MHD).