Double definite integral (Fourier transform)

  • Thread starter secret2
  • Start date
I don't know if this question should be posted here, but I'll give it a shot anyways.

I am trying to find f(x,y), which can be obtain by doing the backward Fourier integral to F(\omega_x, \omega_y). I have 2 questions.

1. Is there any Fortran code that could evaluate the (numerical) Fourier integral?

2. Since the function f(x,y) is 2-dimensional, we have to do a double integral. Suppose that we evaluate first the x-integral. I have a polynimial in the denominator, but the roots of the polynomial will be functions of y. Then, how can I tell if the (simple) poles are in the upper half plane or not?


For the FFT, I generally recommend FFTW.

Alternatively, GSL includes FFT routines.

Another very good implementation (including a version in Fortran), which may be easier to use, is the one by Ooura.

Want to reply to this thread?

"Double definite integral (Fourier transform)" You must log in or register to reply here.

Physics Forums Values

We Value Quality
• Topics based on mainstream science
• Proper English grammar and spelling
We Value Civility
• Positive and compassionate attitudes
• Patience while debating
We Value Productivity
• Disciplined to remain on-topic
• Recognition of own weaknesses
• Solo and co-op problem solving