conatus.expert.numerical_abscissa¶

conatus.expert.numerical_abscissa(ss)[source]¶

ω(Ã) = max eig((Ã + Ãᴴ)/2) = d/dt ‖e^{tÃ}‖ at t = 0⁺, of this matrix.

On the full tau generator this is dominated by the unresolved top modes (the tau rows are not energy-consistent: measured ω grows with N). Use restrict(ss, cutoff=...) (declared truncation) or galerkin_numerical_abscissa (variational, converges) for the physical value.