This works if C is 2-dimensional, but gives an error otherwise (see e.g. this test, which fails since C is 3D and thus ij has 3 elements). Maybe this was meant to be ij[:-1] (drop last element, instead of invert the tuple)? Also to check carefully the need of using 'F' ordering in the previous line, since this basically switches the order of the 'i' and 'j' indices..
This works if
Cis 2-dimensional, but gives an error otherwise (see e.g. this test, which fails since C is 3D and thusijhas 3 elements). Maybe this was meant to beij[:-1](drop last element, instead of invert the tuple)? Also to check carefully the need of using 'F' ordering in the previous line, since this basically switches the order of the 'i' and 'j' indices..https://github.com/GO-Eratosthenes/dhdt/blob/8b8af330e3ec93212265d2be80fdee30d7d4fd5f/dhdt/processing/matching_tools.py#L76