help-gsl
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

[Help-gsl] Does anybody know how to use FFT to compute numerical integra


From: Michael
Subject: [Help-gsl] Does anybody know how to use FFT to compute numerical integration?
Date: Wed, 27 Jun 2007 20:41:16 -0700

The FFT used for numerical integral has the problem of non-adaptive sample
points, which is inefficient.

However I have to use it in my 2D integration problem, for the following
reason:

My 2D integration is:

Integrate( F(v) * FourierTransform[g(t)] (v), v from -infinity to
+infinity ).

------------

g(t)'s evaluation is costly. So I plan to fix the parameter for g(t), and
compute the Fourier Transform of g(t) only once, and store in memroy. And
then in calibration loop, I only vary the parameters for F(v). And each time
for each different set of parameters of F(v), I compute the dot-product of
F(v) sample points and the FT[g(t)] sample points to obtain approximation to
the integral.

My question is: how to improve the accuracy of FFT-based integration? I know
it's inefficient, but is there any remedy at least?

Moreover, is there a better adaptive quadature based "smart" integration
method that can help me deal with the above situation efficiently? I am
thinking of doing a cache for the Fourier Transform of g(t), which is
FT[g(t)](v), since adaptive quadature based integration may sample different
point of FT[g(t)](v) each time... but perhaps the overhead introduced in the
cache may outweigh the smart adaptive integration itself...

Any suggestions? Thanks a lot!




reply via email to

[Prev in Thread] Current Thread [Next in Thread]