comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » Re: What are the errors in the FFT?
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: What are the errors in the FFT? [message #52441] Thu, 08 February 2007 10:24 Go to previous message
Kenneth Bowman is currently offline  Kenneth Bowman
Messages: 86
Registered: November 2006
Member
In article <1170953635.505681.59040@v33g2000cwv.googlegroups.com>,
monty@lanl.gov wrote:

> For a given function f(t) I am finding:
>
> FFT(FFT(f(t),-1),1) -f(t) varies between about 10^-7 to 1-^-8 for
> floating point
> and about 10^-14 to 10^-16 for double precision
>
> (I.e. the inverse transform of the transform deviates from the
> original function)
>
> Is this aproblem with the IDL implementation of the FFT, or is this a
> more fundamental issue with the algorithm itself?
>
> -Monty Wood

That is simply roundoff error. Unavoidable with floating-point
calculations, I'm afraid.

Ken Bowman
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Arrays of Structures
Next Topic: Re: READS as a speed improvement or simply style?

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Thu Oct 09 14:15:32 PDT 2025

Total time taken to generate the page: 1.03939 seconds