Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

Jonas_D_

Beginner

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

04-27-2017
12:07 PM

144 Views

MKL FFT 2D

Hello,

I am trying to compute the solution to the Laplace differential equation using a 2D FFT and MKL.

Are there sample codes to compute the solution to 2nd order differential equations using forward and inverse FFT?

The issue appears to be when I multiply the transformed values by the wavenumbers.

Thanks

Link Copied

2 Replies

Gennady_F_Intel

Moderator

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

04-27-2017
08:13 PM

144 Views

McCalpinJohn

Black Belt

- Mark as New
- Bookmark
- Subscribe
- Mute
- Subscribe to RSS Feed
- Permalink
- Email to a Friend
- Report Inappropriate Content

05-02-2017
11:03 AM

144 Views

I have implemented such solvers many times in the past, with many different FFT libraries, and it is quite easy to make mistakes....

I recommend starting with a very small 2D grid and right-hand-side composed of a single wave of unit amplitude in one direction. After the first set of transforms you should print out the transformed values and verify that you have a single non-negligible amplitude in Fourier space. The magnitude of this value will tell you whether the MKL routine scales the output on the forward transform. (Some libraries scale by 1/N on the inverse transform, but not on the forward transform. Many other variants exist.) The location of the value will tell you whether or not you understand how MKL is packing the transform of a real array into a half-sized complex array -- and there are many ways to do this, so read the documentation carefully and don't be shy about trying out lots of simple cases.

Repeat with the right-hand-side composed of a single wave of unit amplitude in the other direction.

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

For more complete information about compiler optimizations, see our Optimization Notice.