dedalus.tools.sparse

Tools for working with sparse matrices.

Module Contents

STORE_LU
PERMC_SPEC
USE_UMFPACK
scipy_sparse_eigs(A, B, N, target, **kw)

Perform targeted eigenmode search using the scipy/ARPACK sparse solver for the reformulated generalized eigenvalue problem

A.x = λ B.x ==> (A - σB)B.x = (1/(λ-σ)) x

for eigenvalues λ near the target σ.

Parameters:
  • A, B (scipy sparse matrices) – Sparse matrices for generalized eigenvalue problem
  • N (int) – Number of eigenmodes to return
  • target (complex) – Target σ for eigenvalue search
  • Other keyword options passed to scipy.sparse.linalg.eigs.