summaryrefslogtreecommitdiffstats
path: root/src/CicEqualizer.cpp
blob: 52da55e0d3cd2b7c67c619b2393191ed2d458946 (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
/*
   Copyright (C) 2007, 2008, 2009, 2010, 2011 Her Majesty the Queen in
   Right of Canada (Communications Research Center Canada)
 */
/*
   This file is part of CRC-DADMOD.

   CRC-DADMOD 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 of the
   License, or (at your option) any later version.

   CRC-DADMOD 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 CRC-DADMOD.  If not, see <http://www.gnu.org/licenses/>.
 */

#include "CicEqualizer.h"
#include "PcDebug.h"

#include <stdio.h>
#include <stdexcept>


CicEqualizer::CicEqualizer(size_t nbCarriers, size_t spacing, int R) :
    ModCodec(ModFormat(nbCarriers * sizeof(complexf)),
            ModFormat(nbCarriers * sizeof(complexf))),
    myNbCarriers(nbCarriers),
    mySpacing(spacing)
{
    PDEBUG("CicEqualizer::CicEqualizer(%zu, %zu, %i) @ %p\n",
            nbCarriers, spacing, R, this);

    myFilter = new float[nbCarriers];
    const int M = 1;
    const int N = 4;
    const float pi = 4.0f * atanf(1.0f);
    for (size_t i = 0; i < nbCarriers; ++i) {
        int k = i < (nbCarriers + 1) / 2
            ? i + ((nbCarriers & 1) ^ 1)
            : i - (int)nbCarriers;
        float angle = pi * k / spacing;
        if (k == 0) {
            myFilter[i] = 1.0f;
	} else {
		myFilter[i] = sinf(angle / R) / sinf(angle * M);
		myFilter[i] = fabsf(myFilter[i]) * R * M;
		myFilter[i] = powf(myFilter[i], N);
	}
        PDEBUG("HCic[%zu -> %i] = %f (%f dB) -> angle: %f\n",
                i, k,myFilter[i], 20.0 * log10(myFilter[i]), angle);
    }
}


CicEqualizer::~CicEqualizer()
{
    PDEBUG("CicEqualizer::~CicEqualizer() @ %p\n", this);

    if (myFilter != NULL) {
        delete[] myFilter;
    }
}


int CicEqualizer::process(Buffer* const dataIn, Buffer* dataOut)
{
    PDEBUG("CicEqualizer::process(dataIn: %p, dataOut: %p)\n",
            dataIn, dataOut);

    dataOut->setLength(dataIn->getLength());

    const complexf* in = reinterpret_cast<const complexf*>(dataIn->getData());
    complexf* out = reinterpret_cast<complexf*>(dataOut->getData());
    size_t sizeIn = dataIn->getLength() / sizeof(complexf);
    size_t sizeOut = dataOut->getLength() / sizeof(complexf);

    if ((sizeIn % myNbCarriers) != 0) {
        PDEBUG("%zu != %zu\n", sizeIn, myNbCarriers);
        throw std::runtime_error(
                "CicEqualizer::process input size not valid!");
    }

    for (size_t i = 0; i < sizeOut; ) {
        for (size_t j = 0; j < myNbCarriers; ++j, ++i) {
            out[i] = in[i] * myFilter[j];
        }
    }

    return sizeOut;
}