From 39e0005e226ed6f04562e9d5d7548782fef81c20 Mon Sep 17 00:00:00 2001 From: John Glover Date: Tue, 21 Aug 2012 18:31:07 +0100 Subject: [sndobj] Update SndObj to use FFTW v3 (was using v2). Remove unused SndObj files. Whitespace clean up. --- src/sndobj/rfftw/generic.c | 102 --------------------------------------------- 1 file changed, 102 deletions(-) delete mode 100644 src/sndobj/rfftw/generic.c (limited to 'src/sndobj/rfftw/generic.c') diff --git a/src/sndobj/rfftw/generic.c b/src/sndobj/rfftw/generic.c deleted file mode 100644 index e0a2435..0000000 --- a/src/sndobj/rfftw/generic.c +++ /dev/null @@ -1,102 +0,0 @@ -/* - * Copyright (c) 1997-1999 Massachusetts Institute of Technology - * - * This program is free software; you can redistribute it and/or modify - * it under the terms of the GNU General Public License as published by - * the Free Software Foundation; either version 2 of the License, or - * (at your option) any later version. - * - * This program is distributed in the hope that it will be useful, - * but WITHOUT ANY WARRANTY; without even the implied warranty of - * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the - * GNU General Public License for more details. - * - * You should have received a copy of the GNU General Public License - * along with this program; if not, write to the Free Software - * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA - * - */ - -/* - * - * generic.c -- "generic" codelets. They work for all n (and they are - * slow) - */ -#include -#include - -void fftw_twiddle_generic(fftw_complex *A, const fftw_complex *W, - int m, int r, int n, int stride) -{ - int i, j, k; - const fftw_complex *jp; - fftw_complex *kp; - fftw_complex *tmp = (fftw_complex *) - fftw_malloc(r * sizeof(fftw_complex)); - - for (i = 0; i < m; ++i) { - for (k = 0, kp = tmp; k < r; ++k, kp++) { - fftw_real r0, i0, rt, it, rw, iw; - int l1 = i + m * k; - int l0; - - r0 = i0 = 0.0; - for (j = 0, jp = A + i * stride, l0 = 0; j < r; ++j, - jp += m * stride) { - rw = c_re(W[l0]); - iw = c_im(W[l0]); - rt = c_re(*jp); - it = c_im(*jp); - r0 += rt * rw - it * iw; - i0 += rt * iw + it * rw; - l0 += l1; - if (l0 >= n) - l0 -= n; - } - c_re(*kp) = r0; - c_im(*kp) = i0; - } - for (k = 0, kp = A + i * stride; k < r; ++k, kp += m * stride) - *kp = tmp[k]; - } - - fftw_free(tmp); -} - -void fftwi_twiddle_generic(fftw_complex *A, const fftw_complex *W, - int m, int r, int n, int stride) -{ - int i, j, k; - const fftw_complex *jp; - fftw_complex *kp; - fftw_complex *tmp = (fftw_complex *) - fftw_malloc(r * sizeof(fftw_complex)); - - for (i = 0; i < m; ++i) { - for (k = 0, kp = tmp; k < r; ++k, kp++) { - fftw_real r0, i0, rt, it, rw, iw; - int l1 = i + m * k; - int l0; - - r0 = i0 = 0.0; - for (j = 0, jp = A + i * stride, l0 = 0; j < r; ++j, - jp += m * stride) { - rw = c_re(W[l0]); - iw = c_im(W[l0]); - rt = c_re(*jp); - it = c_im(*jp); - r0 += rt * rw + it * iw; - i0 += it * rw - rt * iw; - l0 += l1; - if (l0 >= n) - l0 -= n; - } - c_re(*kp) = r0; - c_im(*kp) = i0; - } - for (k = 0, kp = A + i * stride; k < r; ++k, kp += m * stride) - *kp = tmp[k]; - } - - fftw_free(tmp); -} -- cgit v1.2.3