pub fn kshell_factor(l_a: u8, l_b: u8) -> f64
Compute kshell factor for off-diagonal H0 scaling. GFN2 uses: kshell(k, l) = merge(2.0, (kdiag[l]+kdiag[k])/2, …)