Merged eb/gcell -r8215:8243 into trunk. This adds gr-gcell, the GNU

[debian/gnuradio] / gcell / src / lib / wrapper / gcp_fft_1d_r2.cc

diff --git a/gcell/src/lib/wrapper/gcp_fft_1d_r2.cc b/gcell/src/lib/wrapper/gcp_fft_1d_r2.cc
index f92ee42c88aaaa0644ba99050e9847a62f34204c..07267e303a89ede8d9ba7a0bb7966ac63ee5e136 100644 (file)
--- a/gcell/src/lib/wrapper/gcp_fft_1d_r2.cc
+++ b/gcell/src/lib/wrapper/gcp_fft_1d_r2.cc
@@ -30,18 +30,19 @@ static void
 init_jd(gc_job_desc *jd,
        gc_proc_id_t proc_id,
        unsigned log2_fft_length,
-       bool forward,
+       bool shift,
        std::complex<float> *out,
        const std::complex<float> *in,
-       const std::complex<float> *W)
+       const std::complex<float> *twiddle,
+       const float *window)
 {
   jd->proc_id = proc_id;
   jd->input.nargs = 2;
   jd->output.nargs = 0;
-  jd->eaa.nargs = 3;
+  jd->eaa.nargs = 4;
 
   jd->input.arg[0].u32 = log2_fft_length;
-  jd->input.arg[1].u32 = forward;
+  jd->input.arg[1].u32 = shift;
   unsigned int fft_length = 1 << log2_fft_length;
 
   jd->eaa.arg[0].ea_addr = ptr_to_ea(out);
@@ -52,19 +53,28 @@ init_jd(gc_job_desc *jd,
   jd->eaa.arg[1].direction = GCJD_DMA_GET;
   jd->eaa.arg[1].get_size = sizeof(std::complex<float>) * fft_length;
 
-  jd->eaa.arg[2].ea_addr = ptr_to_ea(const_cast<std::complex<float>*>(W));
+  jd->eaa.arg[2].ea_addr = ptr_to_ea(const_cast<std::complex<float>*>(twiddle));
   jd->eaa.arg[2].direction = GCJD_DMA_GET;
   jd->eaa.arg[2].get_size = sizeof(std::complex<float>) * fft_length / 4;
-}
 
+  jd->eaa.arg[3].ea_addr = ptr_to_ea(const_cast<float*>(window));
+  jd->eaa.arg[3].direction = GCJD_DMA_GET;
+  if (window == 0)
+    jd->eaa.arg[3].get_size = 0;
+  else
+    jd->eaa.arg[3].get_size = sizeof(float) * fft_length;
+}
 
-gc_job_desc *
+  
+gc_job_desc_sptr
 gcp_fft_1d_r2_submit(gc_job_manager_sptr mgr,
                     unsigned int log2_fft_length,
                     bool forward,
+                    bool shift,
                     std::complex<float> *out,
                     const std::complex<float> *in,
-                    const std::complex<float> *W)
+                    const std::complex<float> *twiddle,
+                    const float *window)
 {
   unsigned int fft_length = 1 << log2_fft_length;
   if (fft_length > 4096)
@@ -74,29 +84,36 @@ gcp_fft_1d_r2_submit(gc_job_manager_sptr mgr,
     throw gc_bad_align("out");
   if ((intptr_t)in & 0xf)
     throw gc_bad_align("in");
-  if ((intptr_t)W & 0xf)
-    throw gc_bad_align("W");
+  if ((intptr_t)twiddle & 0xf)
+    throw gc_bad_align("twiddle");
+  if ((intptr_t)window & 0xf)
+    throw gc_bad_align("window");
+
+  std::string proc_name;
+  if (forward)
+    proc_name = "fwd_fft_1d_r2";
+  else
+    proc_name = "inv_fft_1d_r2";
 
-  gc_proc_id_t fft_id = mgr->lookup_proc("fft_1d_r2");
-  gc_job_desc *jd = mgr->alloc_job_desc();
-  init_jd(jd, fft_id, log2_fft_length, forward, out, in, W);
-  if (!mgr->submit_job(jd)){
+  gc_proc_id_t fft_id = mgr->lookup_proc(proc_name);
+  gc_job_desc_sptr jd = gc_job_manager::alloc_job_desc(mgr);
+  init_jd(jd.get(), fft_id, log2_fft_length, shift, out, in, twiddle, window);
+  if (!mgr->submit_job(jd.get())){
     gc_job_status_t s = jd->status;
-    mgr->free_job_desc(jd);
-    throw gc_bad_submit("fft_1d_r2", s);
+    throw gc_bad_submit(proc_name, s);
   }
   return jd;
 }
 
 void
-gcp_fft_1d_r2_twiddle(unsigned int log2_fft_length, std::complex<float> *W)
+gcp_fft_1d_r2_twiddle(unsigned int log2_fft_length, std::complex<float> *twiddle)
 {
   unsigned int n = 1 << log2_fft_length;
 
-  W[0].real() = 1.0;
-  W[0].imag() = 0.0;
+  twiddle[0].real() = 1.0;
+  twiddle[0].imag() = 0.0;
   for (unsigned i=1; i < n/4; i++){
-    W[i].real() =  cos(i * 2*M_PI/n);
-    W[n/4 - i].imag() = -W[i].real();
+    twiddle[i].real() =  cos(i * 2*M_PI/n);
+    twiddle[n/4 - i].imag() = -twiddle[i].real();
   }
 }