From 6c61b72cc7be59138b59b1d64b64377fc5c54395 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?IOhannes=20m=20zm=C3=B6lnig?= Date: Wed, 20 Sep 2006 11:22:29 +0000 Subject: library providing FFT-objects using the (fast!) libFFTW svn path=/trunk/externals/fftw/; revision=5978 --- ifftw~.c | 71 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 71 insertions(+) create mode 100644 ifftw~.c (limited to 'ifftw~.c') 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 -- cgit v1.2.1