On Thu, Oct 25, 2012 at 10:12 AM,  <span dir="ltr"><<a href="mailto:l.buikstra@student.utwente.nl" target="_blank">l.buikstra@student.utwente.nl</a>></span> wrote:<br><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Dear all,<br>
<br>
I'm trying to get a feel for fourier transforms in Petsc, but I can find very little information about doing FTs through Petsc's interface. The examples all seem to just do a transform, then transform back to confirm the same function comes out again; that works fine over here, but I can't seem to figure out how the actual fourier transform itself is stored so it can be manipulated or used before transforming back. I wrote a short test program:<br>

<a href="http://pastebin.com/4saWgqLu" target="_blank">http://pastebin.com/4saWgqLu</a><br>
It fills a vector with a sine, then transforms it and displays the result, which looks like this:<br>
<a href="http://imageshack.us/a/img37/1135/whatph.png" target="_blank">http://imageshack.us/a/img37/1135/whatph.png</a><br>
(the original vector is 2000 points; I expect this is 4000 because it somehow represents a complex vector?)<br>
I would expect a smooth function with two peaks on the far left and right, but instead all it shows is a repetitively garbled mess in the first half of the output and mostly zeroes in the second. Am I just doing something wrong? If so, what, and if not, what exactly does this output represent? How does it relate to the output I'm expecting? Any help would be greatly appreciated.<br>
</blockquote><div><br></div><div>You want something this is a nice multiple of the fundamental of the box, so k pi i/N. Then you will get a spike.</div><div>You can find out which frequencies are where in the FFTW docs.</div>
<div><br></div><div>    Matt</div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Yours,<br>
<br>
Leon Buikstra</blockquote></div><br><br clear="all"><div><br></div>-- <br>What most experimenters take for granted before they begin their experiments is infinitely more interesting than any results to which their experiments lead.<br>
-- Norbert Wiener<br>