mirror of
https://github.com/eried/portapack-mayhem.git
synced 2025-02-08 18:58:34 -05:00
Move IIR code into .cpp file.
A few hundred more text section bytes saved.
This commit is contained in:
parent
9fb22dfd1f
commit
316d5d433b
@ -148,6 +148,7 @@ CPPSRC = main.cpp \
|
|||||||
packet_builder.cpp \
|
packet_builder.cpp \
|
||||||
dsp_fft.cpp \
|
dsp_fft.cpp \
|
||||||
dsp_fir_taps.cpp \
|
dsp_fir_taps.cpp \
|
||||||
|
dsp_iir.cpp \
|
||||||
fxpt_atan2.cpp \
|
fxpt_atan2.cpp \
|
||||||
rssi.cpp \
|
rssi.cpp \
|
||||||
rssi_dma.cpp \
|
rssi_dma.cpp \
|
||||||
|
50
firmware/baseband/dsp_iir.cpp
Normal file
50
firmware/baseband/dsp_iir.cpp
Normal file
@ -0,0 +1,50 @@
|
|||||||
|
/*
|
||||||
|
* Copyright (C) 2015 Jared Boone, ShareBrained Technology, Inc.
|
||||||
|
*
|
||||||
|
* This file is part of PortaPack.
|
||||||
|
*
|
||||||
|
* This program 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 2, or (at your option)
|
||||||
|
* any later version.
|
||||||
|
*
|
||||||
|
* This program 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 this program; see the file COPYING. If not, write to
|
||||||
|
* the Free Software Foundation, Inc., 51 Franklin Street,
|
||||||
|
* Boston, MA 02110-1301, USA.
|
||||||
|
*/
|
||||||
|
|
||||||
|
#include "dsp_iir.hpp"
|
||||||
|
|
||||||
|
#include <hal.h>
|
||||||
|
|
||||||
|
void IIRBiquadFilter::execute(buffer_s16_t buffer_in, buffer_s16_t buffer_out) {
|
||||||
|
// TODO: Assert that buffer_out.count == buffer_in.count.
|
||||||
|
for(size_t i=0; i<buffer_out.count; i++) {
|
||||||
|
const int32_t output_sample = execute_sample(buffer_in.p[i]);
|
||||||
|
const int32_t output_sample_saturated = __SSAT(output_sample, 16);
|
||||||
|
buffer_out.p[i] = output_sample_saturated;
|
||||||
|
}
|
||||||
|
}
|
||||||
|
|
||||||
|
void IIRBiquadFilter::execute_in_place(buffer_s16_t buffer) {
|
||||||
|
execute(buffer, buffer);
|
||||||
|
}
|
||||||
|
|
||||||
|
float IIRBiquadFilter::execute_sample(const float in) {
|
||||||
|
x[0] = x[1];
|
||||||
|
x[1] = x[2];
|
||||||
|
x[2] = in;
|
||||||
|
|
||||||
|
y[0] = y[1];
|
||||||
|
y[1] = y[2];
|
||||||
|
y[2] = config.b[0] * x[2] + config.b[1] * x[1] + config.b[2] * x[0]
|
||||||
|
- config.a[1] * y[1] - config.a[2] * y[0];
|
||||||
|
|
||||||
|
return y[2];
|
||||||
|
}
|
@ -42,36 +42,15 @@ public:
|
|||||||
{
|
{
|
||||||
}
|
}
|
||||||
|
|
||||||
void execute(buffer_s16_t buffer_in, buffer_s16_t buffer_out) {
|
void execute(buffer_s16_t buffer_in, buffer_s16_t buffer_out);
|
||||||
// TODO: Assert that buffer_out.count == buffer_in.count.
|
void execute_in_place(buffer_s16_t buffer);
|
||||||
for(size_t i=0; i<buffer_out.count; i++) {
|
|
||||||
const int32_t output_sample = execute_sample(buffer_in.p[i]);
|
|
||||||
const int32_t output_sample_saturated = __SSAT(output_sample, 16);
|
|
||||||
buffer_out.p[i] = output_sample_saturated;
|
|
||||||
}
|
|
||||||
}
|
|
||||||
|
|
||||||
void execute_in_place(buffer_s16_t buffer) {
|
|
||||||
execute(buffer, buffer);
|
|
||||||
}
|
|
||||||
|
|
||||||
private:
|
private:
|
||||||
const iir_biquad_config_t config;
|
const iir_biquad_config_t config;
|
||||||
std::array<float, 3> x { { 0.0f, 0.0f, 0.0f } };
|
std::array<float, 3> x { { 0.0f, 0.0f, 0.0f } };
|
||||||
std::array<float, 3> y { { 0.0f, 0.0f, 0.0f } };
|
std::array<float, 3> y { { 0.0f, 0.0f, 0.0f } };
|
||||||
|
|
||||||
float execute_sample(const float in) {
|
float execute_sample(const float in);
|
||||||
x[0] = x[1];
|
|
||||||
x[1] = x[2];
|
|
||||||
x[2] = in;
|
|
||||||
|
|
||||||
y[0] = y[1];
|
|
||||||
y[1] = y[2];
|
|
||||||
y[2] = config.b[0] * x[2] + config.b[1] * x[1] + config.b[2] * x[0]
|
|
||||||
- config.a[1] * y[1] - config.a[2] * y[0];
|
|
||||||
|
|
||||||
return y[2];
|
|
||||||
}
|
|
||||||
};
|
};
|
||||||
|
|
||||||
#endif/*__DSP_IIR_H__*/
|
#endif/*__DSP_IIR_H__*/
|
||||||
|
Loading…
x
Reference in New Issue
Block a user