Numerical solution to a stokes interface problem

A numerical algorithm is described for solving the generalized eigenvalue problem arising in the study of the spectrum of a preconditioned operator in the pressure equation derived from a Stokes interface problem. The algorithm is implemented for two finite element schemes. It is tested for a problem with an analytical solution and is applied to spectrum computations in the case of a piecewise constant viscosity. A large number of numerical experiments are analyzed, and recommendations are given for solving the Stokes interface problem in practice.