pub fn wilkinson_shift<N: Real>(tmm: N, tnn: N, tmn: N) -> N
Expand description

Computes the wilkinson shift, i.e., the 2x2 symmetric matrix eigenvalue to its tailing component tnn.

The inputs are interpreted as the 2x2 matrix: tmm tmn tmn tnn