Adding SSE version of fft filters. Complex (ccc) version working.
[debian/gnuradio] / gnuradio-core / src / lib / filter / gr_fft_filter_ccc.cc
1 /* -*- c++ -*- */
2 /*
3  * Copyright 2005 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
18  * along with GNU Radio; see the file COPYING.  If not, write to
19  * the Free Software Foundation, Inc., 51 Franklin Street,
20  * Boston, MA 02110-1301, USA.
21  */
22
23 /*
24  * WARNING: This file is automatically generated by generate_gr_fft_filter_XXX.py
25  * Any changes made to this file will be overwritten.
26  */
27
28 #ifdef HAVE_CONFIG_H
29 #include "config.h"
30 #endif
31
32 #include <gr_fft_filter_ccc.h>
33 #include <gri_fft_filter_ccc_sse.h>
34 #include <gr_io_signature.h>
35 #include <gri_fft.h>
36 #include <math.h>
37 #include <assert.h>
38 #include <stdexcept>
39 #include <gr_firdes.h>
40
41 #include <cstdio>
42 #include <iostream>
43 #include <string.h>
44
45 gr_fft_filter_ccc_sptr gr_make_fft_filter_ccc (int decimation, const std::vector<gr_complex> &taps)
46 {
47   return gr_fft_filter_ccc_sptr (new gr_fft_filter_ccc (decimation, taps));
48 }
49
50
51 gr_fft_filter_ccc::gr_fft_filter_ccc (int decimation, const std::vector<gr_complex> &taps)
52   : gr_sync_decimator ("fft_filter_ccc",
53                        gr_make_io_signature (1, 1, sizeof (gr_complex)),
54                        gr_make_io_signature (1, 1, sizeof (gr_complex)),
55                        decimation),
56     d_updated(false)
57 {
58   set_history(1);
59   d_filter = new gri_fft_filter_ccc_sse(decimation, taps);
60   d_nsamples = d_filter->set_taps(taps);
61   set_output_multiple(d_nsamples);
62 }
63
64 gr_fft_filter_ccc::~gr_fft_filter_ccc ()
65 {
66   delete d_filter;
67 }
68
69 void
70 gr_fft_filter_ccc::set_taps (const std::vector<gr_complex> &taps)
71 {
72   d_new_taps = taps;
73   d_updated = true;
74 }
75
76 int
77 gr_fft_filter_ccc::work (int noutput_items,
78                          gr_vector_const_void_star &input_items,
79                          gr_vector_void_star &output_items)
80 {
81   const gr_complex *in = (const gr_complex *) input_items[0];
82   gr_complex *out = (gr_complex *) output_items[0];
83
84   if (d_updated){
85     d_nsamples = d_filter->set_taps(d_new_taps);
86     d_updated = false;
87     set_output_multiple(d_nsamples);
88     return 0;                           // output multiple may have changed
89   }
90
91   assert(noutput_items % d_nsamples == 0);
92
93   d_filter->filter(noutput_items, in, out);
94
95   //assert((out - (gr_complex *) output_items[0]) == noutput_items);
96
97   return noutput_items;
98 }