[Numpy-discussion] Re: Weird (wrong?) real fft 2d behavior

Warren Focke focke at slac.stanford.edu
Thu Jan 26 18:36:10 CST 2006


On Wed, 25 Jan 2006, Andrew Jaffe wrote:

> Andrew Jaffe wrote:
> > If I start with what I thought was an appropriate (n, n/2+1) complex
> > matrix which should have a real inverse fft, and then take its real fft,
> > I don't get back the original matrix.
> >
> > Note the presence of very small but nonzero reals in the final matrix,

That's just roundoff.

> > and the fact that the 2d and 4th rows are duplicates. This seems to be a
> >  mistake somewhere.
> >
> > Or am I just misunderstanding/misremembering something about 2d real ffts?

It looks wrong to me, and I think I wrote those functions.  I get the same
results in Numeric.  I'll try to look into the problem.

> and I should point out that
> 	delta_rp = N.dft.real_fft2d(delta_kp)
> is 'allclose' to the original delta_r (which leads me to believe that I
> may indeed be misunderstanding something).

"Stable" does not neccessarily imply "correct".

w

>
> Andrew
>
>
> >
> > Andrew
> >
> >
> > In [41]: dims=(4,4)
> >
> > In [42]: delta_k = aniso.dft_realizn(dims, P)
> >          #### this just makes an appropriate matrix
> >
> > In [43]: delta_k
> > Out[43]:
> > array([[-0.    +0.j    ,  0.5904+0.0429j,  0.9063+0.j    ],
> >        [-0.221 +0.j    ,  0.4169+0.4602j, -0.6909+0.j    ],
> >        [ 0.9059+0.j    , -0.8037+0.2604j,  1.835 +0.j    ],
> >        [ 1.7436+0.j    , -0.2382+0.221j , -0.0607+0.j    ]])
> >        #### 16 real numbers
> >
> > In [44]: delta_r = N.dft.inverse_real_fft2d(delta_k)
> >
> > In [45]: delta_r
> > Out[45]:
> > array([[ 0.2718, -0.0956,  0.2805,  0.1505],
> >        [ 0.0297, -0.0533, -0.259 ,  0.0561],
> >        [ 0.1308, -0.2096,  0.2288, -0.3041],
> >        [ 0.0895,  0.1105, -0.3188, -0.1076]])
> >
> > In [46]: delta_kp = N.dft.real_fft2d(delta_r)
> >
> > In [47]: delta_kp
> > Out[47]:
> > array([[ 0.     +0.00e+00j,  0.5904 +4.2938e-02j,  0.9063 +0.0000e+00j],
> >        [ 0.7613 +5.5511e-17j,0.4169 +4.6020e-01j, -0.3758 +5.5511e-17j],
> >        [ 0.9059 +0.0000e+00j,-0.8037+2.6038e-01j, 1.835  +0.0000e+00j],
> >        [ 0.7613 -5.5511e-17j,-0.2382+2.2099e-01j,-0.3758-5.5511e-17j]])
> >
> > In [48]: N.__version__
> > Out[48]: '0.9.5.1982'
> > ------------------------
>
>
>
> -------------------------------------------------------
> This SF.net email is sponsored by: Splunk Inc. Do you grep through log files
> for problems?  Stop!  Download the new AJAX search engine that makes
> searching your log files as easy as surfing the  web.  DOWNLOAD SPLUNK!
> http://sel.as-us.falkag.net/sel?cmd=lnk&kid=103432&bid=230486&dat=121642
> _______________________________________________
> Numpy-discussion mailing list
> Numpy-discussion at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/numpy-discussion
>




More information about the Numpy-discussion mailing list