Function tri 

pub const fn tri(l: usize) -> usize

Triangular index base: q_l^m lives at tri(l) + m.

tri(L) + L = L(L+3)/2 < (L+1)^2 for every L, so the triangular tables always fit in the same N = (L+1)^2 allocation the flat output uses.