pub fn decay_index( psi: &Array2<f64>, grid: &Grid2D, r: f64, z: f64, ) -> FusionResult<f64>
Calculate the decay index at a given position.
n = -(R/B_Z) ยท dB_Z/dR
Stable if 0 < n < 1.5.