Function dispatched to the following (C++) functions:
[1] (sigma: int, k_idx: int) -> ndarray[complex, 2]
P^σ(k)_mν, returned as a matrix in (m ν)