[Numpy-discussion] numpy.fft.irfftn fails apparently unexpectedly

[Warren Focke]

On Tue, 7 Feb 2012, Henry Gomersall wrote:

> On Tue, 2012-02-07 at 11:53 -0800, Warren Focke wrote:
>> You're not doing anything wrong.
>> irfftn takes complex input and returns real output.
>> The exception is a bug which is triggered because max(axes) >=
>> len(axes).
>
> Is this a bug I should register?

Yes.

It should work right if you replace
     s[axes[-1]] = (s[axes[-1]] - 1) * 2
with
     s[-1] = (a.shape[axes[-1]] - 1) * 2
but I'm not really in a position to test it right now.

Warren Focke