summaryrefslogtreecommitdiff
path: root/gr-digital/lib/cma_equalizer_cc_impl.cc
diff options
context:
space:
mode:
Diffstat (limited to 'gr-digital/lib/cma_equalizer_cc_impl.cc')
-rw-r--r--gr-digital/lib/cma_equalizer_cc_impl.cc111
1 files changed, 0 insertions, 111 deletions
diff --git a/gr-digital/lib/cma_equalizer_cc_impl.cc b/gr-digital/lib/cma_equalizer_cc_impl.cc
deleted file mode 100644
index d0dc1e0961..0000000000
--- a/gr-digital/lib/cma_equalizer_cc_impl.cc
+++ /dev/null
@@ -1,111 +0,0 @@
-/* -*- c++ -*- */
-/*
- * Copyright 2011,2012 Free Software Foundation, Inc.
- *
- * This file is part of GNU Radio
- *
- * SPDX-License-Identifier: GPL-3.0-or-later
- *
- */
-
-#ifdef HAVE_CONFIG_H
-#include "config.h"
-#endif
-
-#include "cma_equalizer_cc_impl.h"
-#include <gnuradio/io_signature.h>
-
-namespace gr {
-namespace digital {
-
-using namespace filter::kernel;
-
-cma_equalizer_cc::sptr
-cma_equalizer_cc::make(int num_taps, float modulus, float mu, int sps)
-{
- return gnuradio::make_block_sptr<cma_equalizer_cc_impl>(num_taps, modulus, mu, sps);
-}
-
-cma_equalizer_cc_impl::cma_equalizer_cc_impl(int num_taps,
- float modulus,
- float mu,
- int sps)
- : sync_decimator("cma_equalizer_cc",
- io_signature::make(1, 1, sizeof(gr_complex)),
- io_signature::make(1, 1, sizeof(gr_complex)),
- sps),
- fir_filter_ccc(std::vector<gr_complex>(num_taps, gr_complex(0, 0))),
- d_new_taps(num_taps, gr_complex(0, 0)),
- d_updated(false),
- d_error(gr_complex(0, 0))
-{
- set_modulus(modulus);
- set_gain(mu);
- if (num_taps > 0)
- d_new_taps[0] = 1.0;
- fir_filter_ccc::set_taps(d_new_taps);
-
- set_history(num_taps);
-}
-
-cma_equalizer_cc_impl::~cma_equalizer_cc_impl() {}
-
-void cma_equalizer_cc_impl::set_taps(const std::vector<gr_complex>& taps)
-{
- d_new_taps = taps;
- d_updated = true;
-}
-
-std::vector<gr_complex> cma_equalizer_cc_impl::taps() const { return d_taps; }
-
-gr_complex cma_equalizer_cc_impl::error(const gr_complex& out)
-{
- gr_complex error = out * (norm(out) - d_modulus);
- float re = gr::clip(error.real(), 1.0);
- float im = gr::clip(error.imag(), 1.0);
- return gr_complex(re, im);
-}
-
-void cma_equalizer_cc_impl::update_tap(gr_complex& tap, const gr_complex& in)
-{
- // Hn+1 = Hn - mu*conj(Xn)*zn*(|zn|^2 - 1)
- tap -= d_mu * conj(in) * d_error;
-}
-
-int cma_equalizer_cc_impl::work(int noutput_items,
- gr_vector_const_void_star& input_items,
- gr_vector_void_star& output_items)
-{
- const gr_complex* in = (const gr_complex*)input_items[0];
- gr_complex* out = (gr_complex*)output_items[0];
-
- if (d_updated) {
- d_taps = d_new_taps;
- set_history(d_taps.size());
- d_updated = false;
- return 0; // history requirements may have changed.
- }
-
- int j = 0;
- size_t k, l = d_taps.size();
- for (int i = 0; i < noutput_items; i++) {
- out[i] = filter(&in[j]);
-
- // Adjust taps
- d_error = error(out[i]);
- for (k = 0; k < l; k++) {
- // Update tap locally from error.
- update_tap(d_taps[k], in[j + k]);
-
- // Update aligned taps in filter object.
- fir_filter_ccc::update_tap(d_taps[k], k);
- }
-
- j += decimation();
- }
-
- return noutput_items;
-}
-
-} /* namespace digital */
-} /* namespace gr */