conatus.expert.solve_sparse

conatus.expert.solve_sparse(L, M, target, n_modes=6, *, label='', tol=0.0)

n_modes eigenvalues nearest target by shift-invert Arnoldi on (L + target M)^{-1} M (k= is the deprecated former name of n_modes: k is a wavenumber everywhere else).

Works for singular, non-symmetric M (no M-inner product assumed): the standard problem -(L + s M)^{-1} M x = nu x has nu = 1 / (sigma - s).

Parameters:

• target (complex)

• n_modes (int)

• label (str)

• tol (float)

Return type:

EVPResult