compute_frame() does the wrong thing when G.N == Nf

[pulquero]

trafficstars

Passing a 2D array to filter() is ambiguous when G.N == Nf. A 3D array should be used to be unambiguous, e.g.

s = np.identity(g.G.N) tig = g.filter(s[..., np.newaxis], **kwargs)