diff options
| author | Marcus Müller <marcus@hostalia.de> | 2016-06-09 15:04:22 +0200 |
|---|---|---|
| committer | Marcus Müller <marcus@hostalia.de> | 2016-06-09 15:04:22 +0200 |
| commit | a7dc893fc42d62dcbc84447a8be991f75bac26db (patch) | |
| tree | 3b924c029b490f56c696db5281d928e30164b04a /gr-blocks/lib/divide_cc_impl.cc | |
| parent | 0f61e4e150359d4e68439085a5489f96df04f685 (diff) | |
VOLKized the divide block. ~ 10x speedup
Diffstat (limited to 'gr-blocks/lib/divide_cc_impl.cc')
| -rw-r--r-- | gr-blocks/lib/divide_cc_impl.cc | 66 |
1 files changed, 66 insertions, 0 deletions
diff --git a/gr-blocks/lib/divide_cc_impl.cc b/gr-blocks/lib/divide_cc_impl.cc new file mode 100644 index 0000000000..56ce09e04c --- /dev/null +++ b/gr-blocks/lib/divide_cc_impl.cc @@ -0,0 +1,66 @@ +/* -*- c++ -*- */ +/* + * Copyright 2004-2016 Free Software Foundation, Inc. + * + * This file is part of GNU Radio + * + * GNU Radio is free software; you can redistribute it and/or modify + * it under the terms of the GNU General Public License as published by + * the Free Software Foundation; either version 3, or (at your option) + * any later version. + * + * GNU Radio is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the + * GNU General Public License for more details. + * + * You should have received a copy of the GNU General Public License + * along with GNU Radio; see the file COPYING. If not, write to + * the Free Software Foundation, Inc., 51 Franklin Street, + * Boston, MA 02110-1301, USA. + */ + + +#ifdef HAVE_CONFIG_H +#include "config.h" +#endif + +#include <divide_cc_impl.h> +#include <volk/volk.h> +#include <gnuradio/io_signature.h> + +namespace gr { + namespace blocks { + + divide_cc::sptr divide_cc::make(size_t vlen) + { + return gnuradio::get_initial_sptr(new divide_cc_impl(vlen)); + } + + divide_cc_impl::divide_cc_impl(size_t vlen) + : sync_block ("divide_cc", + io_signature::make (2, -1, sizeof (gr_complex)*vlen), + io_signature::make (1, 1, sizeof (gr_complex)*vlen)), + d_vlen(vlen) + { + } + + int + divide_cc_impl::work(int noutput_items, + gr_vector_const_void_star &input_items, + gr_vector_void_star &output_items) + { + gr_complex *optr = (gr_complex *) output_items[0]; + size_t ninputs = input_items.size (); + gr_complex *numerator = (gr_complex *) input_items[0]; + for(size_t inp = 1; inp < ninputs; ++inp) + { + volk_32fc_x2_divide_32fc(optr, numerator, (gr_complex*) input_items[inp], noutput_items * d_vlen); + numerator = optr; + } + + return noutput_items; + } + + } /* namespace blocks */ +} /* namespace gr */ |