Function for_each_coupled 

fn for_each_coupled(
    destination: &Basis,
    source: &Basis,
    f: impl FnMut(usize, usize),
)

Call f(row, column) for each coupled entry of a destination and a source basis, column by column: equal axial wavenumbers and polarizations.