Working fft. Addresses issue #60.

JuliaLang · Jun 22, 2011 · 3697dae · 3697dae
1 parent 52df227
commit 3697dae
Showing 1 changed file with 22 additions and 0 deletions.
diff --git a/j/fft.j b/j/fft.j
@@ -0,0 +1,22 @@
+fftw = dlopen("libfftw3")
+
+FFTW_FORWARD = int32(-1)
+FFTW_ESTIMATE = uint32(64)
+
+function fft(X::Vector{Complex128})
+
+    Y = copy(X)
+
+    plan = ccall(dlsym(fftw, "fftw_plan_dft_1d"), 
+              Ptr{Void}, 
+              (Int32, Ptr{Complex128}, Ptr{Complex128}, Int32, Uint32, ),
+              length(X), X, Y, FFTW_FORWARD, FFTW_ESTIMATE)
+
+    ccall(dlsym(fftw, "fftw_execute"),
+          Void,
+          (Ptr{Void},),
+          plan)
+
+    return Y
+
+end