1 files changed, 71 insertions, 0 deletions
diff --git a/ifftw~.c b/ifftw~.c
new file mode 100644
index 0000000..4818875
--- /dev/null
+++ b/ifftw~.c
@@ -0,0 +1,71 @@
+/* Copyright (c) 1997-1999 Miller Puckette and others.
+* For information on usage and redistribution, and for a DISCLAIMER OF ALL
+* WARRANTIES, see the file, "LICENSE.txt," in this distribution.  */
+
+#include "fftw~.h"
+
+#ifdef HAVE_LIBFFTW3F
+/* Support for fftw3 by Tim Blechmann                                       */
+
+/* ------------------------ ifft~ -------------------------------- */
+
+static t_class *sigifftw_class;
+
+typedef struct fftw
+{
+    t_object x_obj;
+    float x_f;
+
+    fftwf_plan plan;
+    fftwf_iodim dim;
+} t_sigifftw;
+
+static void *sigifftw_new(void)
+{
+    t_sigifftw *x = (t_sigifftw *)pd_new(sigifftw_class);
+    outlet_new(&x->x_obj, gensym("signal"));
+    outlet_new(&x->x_obj, gensym("signal"));
+    inlet_new(&x->x_obj, &x->x_obj.ob_pd, &s_signal, &s_signal);
+    x->x_f = 0;
+    return (x);
+}
+
+static void sigifftw_free(t_sigifftw * x)
+{
+    fftwf_destroy_plan(x->plan);
+}
+
+static t_int *sigifftw_perform(t_int *w)
+{
+    fftwf_execute(*(fftwf_plan *)w[1]);
+
+    return (w+2);
+}
+
+static void sigifftw_dsp(t_sigifftw *x, t_signal **sp)
+{
+    int n = sp[0]->s_n;
+    float *in1 = sp[0]->s_vec;
+    float *in2 = sp[1]->s_vec;
+    float *out1 = sp[2]->s_vec;
+    float *out2 = sp[3]->s_vec;
+    
+    x->dim.n=n;
+    x->dim.is=1;
+    x->dim.os=1;
+    x->plan = fftwf_plan_guru_split_dft(1, &(x->dim), 0, NULL, in2, 
+					in1, out2, out1, FFTW_ESTIMATE);
+    dsp_add(sigifftw_perform, 1, &x->plan);
+}
+
+EXTERNAL_SETUP ifftw_tilde_setup(void)
+{
+    sigifftw_class = class_new(gensym("ifftw~"), sigifftw_new, 
+			       (t_method) sigifftw_free, 
+			       sizeof(t_sigifftw), 0, 0);
+    CLASS_MAINSIGNALIN(sigifftw_class, t_sigifftw, x_f);
+    class_addmethod(sigifftw_class, (t_method)sigifftw_dsp, 
+		    gensym("dsp"), 0);
+}
+
+#endif