git.gag.com Git - debian/gnuradio/blob - gcell/src/lib/wrapper/spu/gcs_fft_1d_r2.c

   1 /* -*- c++ -*- */
   2 /*
   3  * Copyright 2008 Free Software Foundation, Inc.
   4  *
   5  * This file is part of GNU Radio
   6  *
   7  * GNU Radio is free software; you can redistribute it and/or modify
   8  * it under the terms of the GNU General Public License as published by
   9  * the Free Software Foundation; either version 3, or (at your option)
  10  * any later version.
  11  *
  12  * GNU Radio is distributed in the hope that it will be useful,
  13  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  14  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  15  * GNU General Public License for more details.
  16  *
  17  * You should have received a copy of the GNU General Public License along
  18  * with this program; if not, write to the Free Software Foundation, Inc.,
  19  * 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA.
  20  */
  21
  22 #include <gc_declare_proc.h>
  23 #include <libfft.h>
  24
  25 static void
  26 gcs_fft_1d_r2(const gc_job_direct_args_t *input,
  27               gc_job_direct_args_t *output __attribute__((unused)),
  28               const gc_job_ea_args_t *eaa)
  29 {
  30   vector float *out = (vector float *) eaa->arg[0].ls_addr;
  31   vector float *in = (vector float *) eaa->arg[1].ls_addr;
  32   vector float *W = (vector float *) eaa->arg[2].ls_addr;
  33   int log2_fft_length = input->arg[0].u32;
  34   int forward = input->arg[1].u32;      // non-zero if forward xform (FIXME use)
  35
  36   fft_1d_r2(out, in, W, log2_fft_length);
  37 }
  38
  39 GC_DECLARE_PROC(gcs_fft_1d_r2, "fft_1d_r2");