summaryrefslogtreecommitdiff
path: root/src/fftw3/dft/rader-omega.c
diff options
context:
space:
mode:
authorscuri <scuri>2009-08-20 12:35:06 +0000
committerscuri <scuri>2009-08-20 12:35:06 +0000
commit5d735255ddd3cb2f547abd3d03969af3fb7eb04e (patch)
tree8fb66510bc625bb1b08ccb41f1b83fb0f7cb8f19 /src/fftw3/dft/rader-omega.c
parent35733b87eed86e5228f12fa10c98a3d9d22a6073 (diff)
*** empty log message ***
Diffstat (limited to 'src/fftw3/dft/rader-omega.c')
-rw-r--r--src/fftw3/dft/rader-omega.c57
1 files changed, 0 insertions, 57 deletions
diff --git a/src/fftw3/dft/rader-omega.c b/src/fftw3/dft/rader-omega.c
deleted file mode 100644
index 7fd86b9..0000000
--- a/src/fftw3/dft/rader-omega.c
+++ /dev/null
@@ -1,57 +0,0 @@
-/*
- * Copyright (c) 2003 Matteo Frigo
- * Copyright (c) 2003 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
- *
- */
-
-#include "dft.h"
-
-static rader_tl *omegas = 0;
-
-R *X(dft_rader_mkomega)(plan *p_, int n, int ginv)
-{
- plan_dft *p = (plan_dft *) p_;
- R *omega;
- int i, gpower;
- trigreal scale;
-
- if ((omega = X(rader_tl_find)(n, n, ginv, omegas)))
- return omega;
-
- omega = (R *)MALLOC(sizeof(R) * (n - 1) * 2, TWIDDLES);
-
- scale = n - 1.0; /* normalization for convolution */
-
- for (i = 0, gpower = 1; i < n-1; ++i, gpower = MULMOD(gpower, ginv, n)) {
- omega[2*i] = X(cos2pi)(gpower, n) / scale;
- omega[2*i+1] = FFT_SIGN * X(sin2pi)(gpower, n) / scale;
- }
- A(gpower == 1);
-
- AWAKE(p_, 1);
- p->apply(p_, omega, omega + 1, omega, omega + 1);
- AWAKE(p_, 0);
-
- X(rader_tl_insert)(n, n, ginv, omega, &omegas);
- return omega;
-}
-
-void X(dft_rader_free_omega)(R **omega)
-{
- X(rader_tl_delete)(*omega, &omegas);
- *omega = 0;
-}