3 * Copyright 2004, 2009,2010 Free Software Foundation, Inc.
5 * This file is part of GNU Radio
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)
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.
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.
30 #include <gr_io_signature.h>
33 gr_make_@BASE_NAME@ (size_t vlen)
35 return gnuradio::get_initial_sptr (new @NAME@ (vlen));
38 @NAME@::@NAME@ (size_t vlen)
39 : gr_sync_block ("@BASE_NAME@",
40 gr_make_io_signature (1, -1, sizeof (@I_TYPE@)*vlen),
41 gr_make_io_signature (1, 1, sizeof (@O_TYPE@)*vlen)),
47 @NAME@::work (int noutput_items,
48 gr_vector_const_void_star &input_items,
49 gr_vector_void_star &output_items)
51 @O_TYPE@ *optr = (@O_TYPE@ *) output_items[0];
53 int ninputs = input_items.size ();
55 if (ninputs == 1){ // compute reciprocal
56 for (int i = 0; i < noutput_items*d_vlen; i++)
57 *optr++ = (@O_TYPE@) ((@O_TYPE@) 1 /
58 ((@I_TYPE@ *) input_items[0])[i]);
62 for (int i = 0; i < noutput_items*d_vlen; i++){
63 @I_TYPE@ acc = ((@I_TYPE@ *) input_items[0])[i];
64 for (int j = 1; j < ninputs; j++)
65 acc /= ((@I_TYPE@ *) input_items[j])[i];
67 *optr++ = (@O_TYPE@) acc;