efficient moving average filters, should speed up OFDM significantly
[debian/gnuradio] / gnuradio-core / src / lib / gengen / gr_moving_average_XX.cc.t
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
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 // @WARNING@
24
25 #ifdef HAVE_CONFIG_H
26 #include "config.h"
27 #endif
28
29 #include <@NAME@.h>
30 #include <gr_io_signature.h>
31
32 @SPTR_NAME@ 
33 gr_make_@BASE_NAME@ (int length, @O_TYPE@ scale, int max_iter)
34 {
35   return @SPTR_NAME@ (new @NAME@ (length, scale, max_iter));
36 }
37
38 @NAME@::@NAME@ (int length, @O_TYPE@ scale, int max_iter)
39   : gr_sync_block ("@BASE_NAME@",
40                    gr_make_io_signature (1, 1, sizeof (@I_TYPE@)),
41                    gr_make_io_signature (1, 1, sizeof (@O_TYPE@))),
42     d_length(length),
43     d_scale(scale),
44     d_max_iter(max_iter)
45 {
46   set_history(length);
47 }
48
49 @NAME@::~@NAME@ ()
50 {
51 }
52
53 int 
54 @NAME@::work (int noutput_items,
55               gr_vector_const_void_star &input_items,
56               gr_vector_void_star &output_items)
57 {
58   const @I_TYPE@ *in = (const @I_TYPE@ *) input_items[0];
59   @O_TYPE@ *out = (@O_TYPE@ *) output_items[0];
60
61   @I_TYPE@ sum = 0;
62   int num_iter = (noutput_items>d_max_iter) ? d_max_iter : noutput_items;
63   for (int i = 0; i < d_length-1 ; i++) {
64     sum += in[i];
65   }
66
67   for (int i = 0; i < num_iter; i++) {
68     sum += in[i+d_length-1];
69     out[i] = sum * d_scale;
70     sum -= in[i];
71   }
72
73   return num_iter;
74 }