pg_utils.numerics.utils.eigen_drift

pg_utils.numerics.utils.eigen_drift(eig_base: ndarray, eig_comp: ndarray, waterlevel: float = 0.0, **opt_cluster)

Calculate eigenvalue drift ratio using Boyd’s method ([Boyd])

Parameters:

• eig_base (np.ndarray) – eigenvalue array used as a base

• eig_comp (np.ndarray) – eigenvalue array used for comparison

• waterlevel (float) – waterlevel for near-trivial eigenvalue to avoid division by zero; default to zero (assuming nontrivial eigenvalues).

• **opt_cluster – optional keyword arguments for clustering, see cluster_modes

[Boyd]

Boyd, Chebyshev and Fourier Spectral Methods.

Note

Be sure to pass in pre-sorted eigenvalues.