Function harmonics 

fn harmonics(lmax: i32) -> impl Iterator<Item = (i32, i32)>

Harmonic (p, m) pairs of degrees p <= lmax in table order; (p, m) sits at p * p + p + m.