199 lines
5 KiB
C++
199 lines
5 KiB
C++
/*****************************************************************************
|
|
|
|
Upsampler2xNeon.hpp
|
|
Author: Laurent de Soras, 2016
|
|
|
|
--- Legal stuff ---
|
|
|
|
This program is free software. It comes without any warranty, to
|
|
the extent permitted by applicable law. You can redistribute it
|
|
and/or modify it under the terms of the Do What The Fuck You Want
|
|
To Public License, Version 2, as published by Sam Hocevar. See
|
|
http://sam.zoy.org/wtfpl/COPYING for more details.
|
|
|
|
*Tab=3***********************************************************************/
|
|
|
|
|
|
|
|
#if ! defined (hiir_Upsampler2xNeon_CODEHEADER_INCLUDED)
|
|
#define hiir_Upsampler2xNeon_CODEHEADER_INCLUDED
|
|
|
|
|
|
|
|
/*\\\ INCLUDE FILES \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
|
|
|
#include "hiir/StageProcNeon.h"
|
|
|
|
#include <arm_neon.h>
|
|
|
|
#include <cassert>
|
|
|
|
|
|
|
|
namespace hiir
|
|
{
|
|
|
|
|
|
|
|
/*\\\ PUBLIC \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
|
|
|
|
|
|
|
/*
|
|
==============================================================================
|
|
Name: ctor
|
|
Throws: Nothing
|
|
==============================================================================
|
|
*/
|
|
|
|
template <int NC>
|
|
Upsampler2xNeon <NC>::Upsampler2xNeon ()
|
|
: _filter ()
|
|
{
|
|
for (int i = 0; i < NBR_STAGES + 1; ++i)
|
|
{
|
|
_filter [i]._coef4 = vdupq_n_f32 (0);
|
|
}
|
|
if ((NBR_COEFS & 1) != 0)
|
|
{
|
|
const int pos = (NBR_COEFS ^ 1) & (STAGE_WIDTH - 1);
|
|
_filter [NBR_STAGES]._coef [pos] = 1;
|
|
}
|
|
|
|
clear_buffers ();
|
|
}
|
|
|
|
|
|
|
|
/*
|
|
==============================================================================
|
|
Name: set_coefs
|
|
Description:
|
|
Sets filter coefficients. Generate them with the PolyphaseIir2Designer
|
|
class.
|
|
Call this function before doing any processing.
|
|
Input parameters:
|
|
- coef_arr: Array of coefficients. There should be as many coefficients as
|
|
mentioned in the class template parameter.
|
|
Throws: Nothing
|
|
==============================================================================
|
|
*/
|
|
|
|
template <int NC>
|
|
void Upsampler2xNeon <NC>::set_coefs (const double coef_arr [NBR_COEFS])
|
|
{
|
|
assert (coef_arr != 0);
|
|
|
|
for (int i = 0; i < NBR_COEFS; ++i)
|
|
{
|
|
const int stage = (i / STAGE_WIDTH) + 1;
|
|
const int pos = (i ^ 1) & (STAGE_WIDTH - 1);
|
|
_filter [stage]._coef [pos] = float (coef_arr [i]);
|
|
}
|
|
}
|
|
|
|
|
|
|
|
/*
|
|
==============================================================================
|
|
Name: process_sample
|
|
Description:
|
|
Upsamples (x2) the input sample, generating two output samples.
|
|
Input parameters:
|
|
- input: The input sample.
|
|
Output parameters:
|
|
- out_0: First output sample.
|
|
- out_1: Second output sample.
|
|
Throws: Nothing
|
|
==============================================================================
|
|
*/
|
|
|
|
template <int NC>
|
|
void Upsampler2xNeon <NC>::process_sample (float &out_0, float &out_1, float input)
|
|
{
|
|
const float32x2_t spl_in = vdup_n_f32 (input);
|
|
const float32x2_t spl_mid = vget_low_f32 (_filter [NBR_STAGES]._mem4);
|
|
float32x4_t y = vcombine_f32 (spl_in, spl_mid);
|
|
float32x4_t mem = _filter [0]._mem4;
|
|
|
|
StageProcNeon <NBR_STAGES>::process_sample_pos (&_filter [0], y, mem);
|
|
_filter [NBR_STAGES]._mem4 = y;
|
|
|
|
out_0 = vgetq_lane_f32 (y, 3);
|
|
out_1 = vgetq_lane_f32 (y, 2);
|
|
}
|
|
|
|
|
|
|
|
/*
|
|
==============================================================================
|
|
Name: process_block
|
|
Description:
|
|
Upsamples (x2) the input sample block.
|
|
Input and output blocks may overlap, see assert() for details.
|
|
Input parameters:
|
|
- in_ptr: Input array, containing nbr_spl samples.
|
|
- nbr_spl: Number of input samples to process, > 0
|
|
Output parameters:
|
|
- out_0_ptr: Output sample array, capacity: nbr_spl * 2 samples.
|
|
Throws: Nothing
|
|
==============================================================================
|
|
*/
|
|
|
|
template <int NC>
|
|
void Upsampler2xNeon <NC>::process_block (float out_ptr [], const float in_ptr [], long nbr_spl)
|
|
{
|
|
assert (out_ptr != 0);
|
|
assert (in_ptr != 0);
|
|
assert (out_ptr >= in_ptr + nbr_spl || in_ptr >= out_ptr + nbr_spl);
|
|
assert (nbr_spl > 0);
|
|
|
|
long pos = 0;
|
|
do
|
|
{
|
|
process_sample (out_ptr [pos * 2], out_ptr [pos * 2 + 1], in_ptr [pos]);
|
|
++ pos;
|
|
}
|
|
while (pos < nbr_spl);
|
|
}
|
|
|
|
|
|
|
|
/*
|
|
==============================================================================
|
|
Name: clear_buffers
|
|
Description:
|
|
Clears filter memory, as if it processed silence since an infinite amount
|
|
of time.
|
|
Throws: Nothing
|
|
==============================================================================
|
|
*/
|
|
|
|
template <int NC>
|
|
void Upsampler2xNeon <NC>::clear_buffers ()
|
|
{
|
|
for (int i = 0; i < NBR_STAGES + 1; ++i)
|
|
{
|
|
_filter [i]._mem4 = vdupq_n_f32 (0);
|
|
}
|
|
}
|
|
|
|
|
|
|
|
/*\\\ PROTECTED \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
|
|
|
|
|
|
|
/*\\\ PRIVATE \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|
|
|
|
|
|
|
|
} // namespace hiir
|
|
|
|
|
|
|
|
#endif // hiir_Upsampler2xNeon_CODEHEADER_INCLUDED
|
|
|
|
|
|
|
|
/*\\\ EOF \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\*/
|