summaryrefslogtreecommitdiff
path: root/gr-digital/lib/pfb_clock_sync_ccf_impl.h
blob: 34a8bc981c429d3c0947727d77e3181cd0eba178 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
/* -*- c++ -*- */
/*
 * Copyright 2009,2010,2012 Free Software Foundation, Inc.
 *
 * This file is part of GNU Radio
 *
 * SPDX-License-Identifier: GPL-3.0-or-later
 *
 */

#ifndef INCLUDED_DIGITAL_PFB_CLOCK_SYNC_CCF_IMPL_H
#define INCLUDED_DIGITAL_PFB_CLOCK_SYNC_CCF_IMPL_H

#include <gnuradio/digital/pfb_clock_sync_ccf.h>

using namespace gr::filter;

namespace gr {
namespace digital {

class pfb_clock_sync_ccf_impl : public pfb_clock_sync_ccf
{
private:
    bool d_updated;
    double d_sps;
    float d_loop_bw;
    float d_damping;
    float d_alpha;
    float d_beta;

    const int d_nfilters;
    int d_taps_per_filter;
    std::vector<kernel::fir_filter_ccf> d_filters;
    std::vector<kernel::fir_filter_ccf> d_diff_filters;
    std::vector<std::vector<float>> d_taps;
    std::vector<std::vector<float>> d_dtaps;
    std::vector<float> d_updated_taps;

    float d_k;
    float d_rate;
    float d_rate_i;
    float d_rate_f;
    float d_max_dev;
    int d_filtnum;
    int d_osps;
    float d_error;
    int d_out_idx;

    uint64_t d_old_in, d_new_in, d_last_out;

    void create_diff_taps(const std::vector<float>& newtaps,
                          std::vector<float>& difftaps);

    void set_taps(const std::vector<float>& taps,
                  std::vector<std::vector<float>>& ourtaps,
                  std::vector<kernel::fir_filter_ccf>& ourfilter);

public:
    pfb_clock_sync_ccf_impl(double sps,
                            float loop_bw,
                            const std::vector<float>& taps,
                            unsigned int filter_size = 32,
                            float init_phase = 0,
                            float max_rate_deviation = 1.5,
                            int osps = 1);

    void setup_rpc();

    void update_gains();

    void forecast(int noutput_items, gr_vector_int& ninput_items_required) override;

    void update_taps(const std::vector<float>& taps);

    std::vector<std::vector<float>> taps() const;
    std::vector<std::vector<float>> diff_taps() const;
    std::vector<float> channel_taps(int channel) const;
    std::vector<float> diff_channel_taps(int channel) const;
    std::string taps_as_string() const;
    std::string diff_taps_as_string() const;

    void set_loop_bandwidth(float bw);
    void set_damping_factor(float df);
    void set_alpha(float alpha);
    void set_beta(float beta);
    void set_max_rate_deviation(float m) { d_max_dev = m; }

    float loop_bandwidth() const;
    float damping_factor() const;
    float alpha() const;
    float beta() const;
    float clock_rate() const;

    float error() const;
    float rate() const;
    float phase() const;

    /*******************************************************************
     *******************************************************************/

    bool check_topology(int ninputs, int noutputs) override;

    int general_work(int noutput_items,
                     gr_vector_int& ninput_items,
                     gr_vector_const_void_star& input_items,
                     gr_vector_void_star& output_items) override;
};

} /* namespace digital */
} /* namespace gr */

#endif /* INCLUDED_DIGITAL_PFB_CLOCK_SYNC_CCF_IMPL_H */