From a7dc893fc42d62dcbc84447a8be991f75bac26db Mon Sep 17 00:00:00 2001
From: Marcus Müller <marcus@hostalia.de>
Date: Thu, 9 Jun 2016 15:04:22 +0200
Subject: VOLKized the divide block. ~ 10x speedup
---
gr-blocks/lib/divide_cc_impl.cc | 66 +++++++++++++++++++++++++++++++++++++++++
1 file changed, 66 insertions(+)
create mode 100644 gr-blocks/lib/divide_cc_impl.cc
(limited to 'gr-blocks/lib/divide_cc_impl.cc')
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 */
--
cgit v1.2.3