Skip to content

GitLab

E
ESP32-audioI2S
Unverified Commit 33e09e72 authored by Wolle's avatar Wolle Committed by GitHub
Browse files
Options

Add files via upload

parent 5ab8fabe
Hide whitespace changes
Inline Side-by-side
Showing with 3943 additions and 241 deletions
src/opus_decoder/celt.cpp 0 → 100644
View file @ 33e09e72
This source diff could not be displayed because it is too large. You can view the blob instead.
src/opus_decoder/celt.h
View file @ 33e09e72
...@@ -37,11 +37,6 @@ ...@@ -37,11 +37,6 @@
#include "Arduino.h" #include "Arduino.h"
//#include "opus_decoder.h" //#include "opus_decoder.h"
#ifdef __cplusplus
extern "C" {
#endif
#define OPUS_OK 0 #define OPUS_OK 0
#define OPUS_BAD_ARG -1 #define OPUS_BAD_ARG -1
#define OPUS_BUFFER_TOO_SMALL -2 #define OPUS_BUFFER_TOO_SMALL -2
...@@ -51,78 +46,66 @@ extern "C" { ...@@ -51,78 +46,66 @@ extern "C" {
#define OPUS_INVALID_STATE -6 #define OPUS_INVALID_STATE -6
#define OPUS_ALLOC_FAIL -7 #define OPUS_ALLOC_FAIL -7
#define OPUS_RESET_STATE 4028 #define OPUS_RESET_STATE 4028
#define OPUS_GET_LOOKAHEAD_REQUEST 4027
#define OPUS_GET_SAMPLE_RATE_REQUEST 4029 #define OPUS_GET_SAMPLE_RATE_REQUEST 4029
#define OPUS_GET_FINAL_RANGE_REQUEST 4031
#define OPUS_GET_PITCH_REQUEST 4033
#define OPUS_SET_GAIN_REQUEST 4034
#define OPUS_GET_GAIN_REQUEST 4045 /* Should have been 4035 */
#define OPUS_SET_PHASE_INVERSION_DISABLED_REQUEST 4046
#define OPUS_GET_PHASE_INVERSION_DISABLED_REQUEST 4047
#define LEAK_BANDS 19 #define LEAK_BANDS 19
typedef struct { typedef struct {
int32_t valid; int32_t valid;
float tonality; float tonality;
float tonality_slope; float tonality_slope;
float noisiness; float noisiness;
float activity; float activity;
float music_prob; float music_prob;
float music_prob_min; float music_prob_min;
float music_prob_max; float music_prob_max;
int32_t bandwidth; int32_t bandwidth;
float activity_probability; float activity_probability;
float max_pitch_ratio; float max_pitch_ratio;
/* Store as Q6 char to save space. */ /* Store as Q6 char to save space. */
uint8_t leak_boost[LEAK_BANDS]; uint8_t leak_boost[LEAK_BANDS];
} AnalysisInfo; } AnalysisInfo;
/*OPT: ec_window must be at least 32 bits, but if you have fast arithmetic on a larger type, you can speed up the /*OPT: ec_window must be at least 32 bits, but if you have fast arithmetic on a larger type, you can speed up the
decoder by using it here.*/ decoder by using it here.*/
typedef uint32_t ec_window;
typedef struct ec_ctx ec_ctx;
typedef struct ec_ctx ec_enc;
typedef struct ec_ctx ec_dec;
typedef struct CELTMode CELTMode; typedef struct CELTMode CELTMode;
typedef struct CELTDecoder CELTDecoder; typedef struct CELTDecoder CELTDecoder;
struct ec_ctx { typedef struct _ec_ctx {
uint8_t *buf; /*Buffered input/output.*/ uint8_t *buf; /*Buffered input/output.*/
uint32_t storage; /*The size of the buffer.*/ uint32_t storage; /*The size of the buffer.*/
uint32_t end_offs; /*The offset at which the last byte containing raw bits was read/written.*/ uint32_t end_offs; /*The offset at which the last byte containing raw bits was read/written.*/
ec_window end_window; /*Bits that will be read from/written at the end.*/ uint32_t end_window; /*Bits that will be read from/written at the end.*/
int32_t nend_bits; /*Number of valid bits in end_window.*/ int32_t nend_bits; /*Number of valid bits in end_window.*/
int32_t nbits_total; int32_t nbits_total;
uint32_t offs; /*The offset at which the next range coder byte will be read/written.*/ uint32_t offs; /*The offset at which the next range coder byte will be read/written.*/
uint32_t rng; /*The number of values in the current range.*/ uint32_t rng; /*The number of values in the current range.*/
uint32_t val; uint32_t val;
uint32_t ext; uint32_t ext;
int32_t rem; /*A buffered input/output symbol, awaiting carry propagation.*/ int32_t rem; /*A buffered input/output symbol, awaiting carry propagation.*/
int32_t error; /*Nonzero if an error occurred.*/ int32_t error; /*Nonzero if an error occurred.*/
}; } ec_ctx_t;
struct band_ctx{ extern ec_ctx_t s_ec;
extern const uint8_t cache_bits50[392];
extern const int16_t cache_index50[105];
typedef struct _band_ctx{
int32_t encode; int32_t encode;
int32_t resynth; int32_t resynth;
const CELTMode *m;
int32_t i; int32_t i;
int32_t intensity; int32_t intensity;
int32_t spread; int32_t spread;
int32_t tf_change; int32_t tf_change;
ec_ctx *ec;
int32_t remaining_bits; int32_t remaining_bits;
const int32_t *bandE; const int32_t *bandE;
uint32_t seed; uint32_t seed;
int32_t arch;
int32_t theta_round; int32_t theta_round;
int32_t disable_inv; int32_t disable_inv;
int32_t avoid_split_noise; int32_t avoid_split_noise;
}; } band_ctx_t;
struct split_ctx{ struct split_ctx{
int32_t inv; int32_t inv;
...@@ -139,17 +122,14 @@ struct CELTDecoder { ...@@ -139,17 +122,14 @@ struct CELTDecoder {
int32_t channels; int32_t channels;
int32_t stream_channels; int32_t stream_channels;
int32_t downsample;
int32_t start, end; int32_t start, end;
int32_t signalling; int32_t signalling;
int32_t disable_inv; int32_t disable_inv;
int32_t arch;
uint32_t rng; uint32_t rng;
int32_t error; int32_t error;
int32_t last_pitch_index; int32_t last_pitch_index;
int32_t loss_count; int32_t loss_count;
int32_t skip_plc;
int32_t postfilter_period; int32_t postfilter_period;
int32_t postfilter_period_old; int32_t postfilter_period_old;
int16_t postfilter_gain; int16_t postfilter_gain;
...@@ -177,11 +157,6 @@ typedef struct { ...@@ -177,11 +157,6 @@ typedef struct {
int16_t i; int16_t i;
}kiss_twiddle_cpx; }kiss_twiddle_cpx;
typedef struct arch_fft_state{
int32_t is_supported;
void *priv;
} arch_fft_state;
#define MAXFACTORS 8 #define MAXFACTORS 8
typedef struct kiss_fft_state{ typedef struct kiss_fft_state{
...@@ -192,7 +167,6 @@ typedef struct kiss_fft_state{ ...@@ -192,7 +167,6 @@ typedef struct kiss_fft_state{
int16_t factors[2*MAXFACTORS]; int16_t factors[2*MAXFACTORS];
const int16_t *bitrev; const int16_t *bitrev;
const kiss_twiddle_cpx *twiddles; const kiss_twiddle_cpx *twiddles;
arch_fft_state *arch_fft;
} kiss_fft_state; } kiss_fft_state;
typedef struct { typedef struct {
...@@ -200,14 +174,9 @@ typedef struct { ...@@ -200,14 +174,9 @@ typedef struct {
int32_t maxshift; int32_t maxshift;
const kiss_fft_state *kfft[4]; const kiss_fft_state *kfft[4];
const int16_t * trig; const int16_t * trig;
} mdct_lookup; } mdct_lookup_t;
typedef struct {
int32_t size;
const int16_t *index;
const uint8_t *bits;
const uint8_t *caps;
} PulseCache;
/** Mode definition (opaque) /** Mode definition (opaque)
@brief Mode definition @brief Mode definition
...@@ -215,35 +184,24 @@ typedef struct { ...@@ -215,35 +184,24 @@ typedef struct {
struct CELTMode { struct CELTMode {
int32_t Fs; int32_t Fs;
int32_t overlap; int32_t overlap;
int32_t nbEBands; int32_t nbEBands;
int32_t effEBands; int32_t effEBands;
int16_t preemph[4]; int16_t preemph[4];
const int16_t *eBands; /**< Definition for each "pseudo-critical band" */
int32_t maxLM; int32_t maxLM;
int32_t nbShortMdcts; int32_t nbShortMdcts;
int32_t shortMdctSize; int32_t shortMdctSize;
int32_t nbAllocVectors; /**< Number of lines in the matrix below */ int32_t nbAllocVectors; /**< Number of lines in the matrix below */
const uint8_t *allocVectors; /**< Number of bits in each band for several rates */
const int16_t *logN;
const int16_t *window;
mdct_lookup mdct;
PulseCache cache;
}; };
extern const CELTMode m_CELTMode;
#define min(a,b) ((a)<(b)?(a):(b)) #define min(a,b) ((a)<(b)?(a):(b))
#define max(a,b) ((a)>(b)?(a):(b)) #define max(a,b) ((a)>(b)?(a):(b))
#define S_MUL(a,b) MULT16_32_Q15(b, a) inline int32_t S_MUL(int32_t a, int16_t b){return (int64_t)b * a >> 15;}
#define C_MUL(m,a,b) do{ (m).r = SUB32_ovflw(S_MUL((a).r,(b).r) , S_MUL((a).i,(b).i)); \ #define C_MUL(m,a,b) do{ (m).r = SUB32_ovflw(S_MUL((a).r,(b).r) , S_MUL((a).i,(b).i)); \
(m).i = ADD32_ovflw(S_MUL((a).r,(b).i) , S_MUL((a).i,(b).r)); }while(0) (m).i = ADD32_ovflw(S_MUL((a).r,(b).i) , S_MUL((a).i,(b).r)); }while(0)
#define C_MULC(m,a,b) do{ (m).r = ADD32_ovflw(S_MUL((a).r,(b).r) , S_MUL((a).i,(b).i)); \
(m).i = SUB32_ovflw(S_MUL((a).i,(b).r) , S_MUL((a).r,(b).i)); }while(0)
#define C_MULBYSCALAR( c, s ) do{ (c).r = S_MUL( (c).r , s ) ; (c).i = S_MUL( (c).i , s ) ; }while(0) #define C_MULBYSCALAR( c, s ) do{ (c).r = S_MUL( (c).r , s ) ; (c).i = S_MUL( (c).i , s ) ; }while(0)
#define DIVSCALAR(x,k) (x) = S_MUL( x, (32767-((k)>>1))/(k)+1 ) #define DIVSCALAR(x,k) (x) = S_MUL( x, (32767-((k)>>1))/(k)+1 )
...@@ -265,7 +223,7 @@ struct CELTMode { ...@@ -265,7 +223,7 @@ struct CELTMode {
#define VERY_LARGE16 ((int16_t)32767) #define VERY_LARGE16 ((int16_t)32767)
#define Q15_ONE ((int16_t)32767) #define Q15_ONE ((int16_t)32767)
#define EC_WINDOW_SIZE ((int32_t)sizeof(ec_window)*CHAR_BIT) #define EC_WINDOW_SIZE ((int32_t)sizeof(uint32_t)*CHAR_BIT)
#define EC_UINT_BITS (8) #define EC_UINT_BITS (8)
#define BITRES 3 #define BITRES 3
#define EC_MINI(_a,_b) ((_a)+(((_b)-(_a))&-((_b)<(_a)))) #define EC_MINI(_a,_b) ((_a)+(((_b)-(_a))&-((_b)<(_a))))
...@@ -284,7 +242,8 @@ inline int32_t MULT16_32_Q16(int64_t a, int64_t b){return (int32_t) (a * b) >> 1 ...@@ -284,7 +242,8 @@ inline int32_t MULT16_32_Q16(int64_t a, int64_t b){return (int32_t) (a * b) >> 1
/** 16x32 multiplication, followed by a 15-bit shift right. Results fits in 32 bits */ /** 16x32 multiplication, followed by a 15-bit shift right. Results fits in 32 bits */
#define MULT16_32_Q15(a,b) ((int32_t)((int64_t)((int16_t)(a))*(b) >> 15)) inline int32_t MULT16_32_Q15(int16_t a, int32_t b){return (int64_t)a * b >> 15;}
/** 32x32 multiplication, followed by a 31-bit shift right. Results fits in 32 bits */ /** 32x32 multiplication, followed by a 31-bit shift right. Results fits in 32 bits */
#define MULT32_32_Q31(a,b) ((int32_t)((int64_t)(a)*(int64_t)(b) >> 31)) #define MULT32_32_Q31(a,b) ((int32_t)((int64_t)(a)*(int64_t)(b) >> 31))
...@@ -309,7 +268,7 @@ inline int32_t MULT16_32_Q16(int64_t a, int64_t b){return (int32_t) (a * b) >> 1 ...@@ -309,7 +268,7 @@ inline int32_t MULT16_32_Q16(int64_t a, int64_t b){return (int32_t) (a * b) >> 1
#define SHL32(a,shift) ((int32_t)((uint32_t)(a)<<(shift))) #define SHL32(a,shift) ((int32_t)((uint32_t)(a)<<(shift)))
/** 32-bit arithmetic shift right with rounding-to-nearest instead of rounding down */ /** 32-bit arithmetic shift right with rounding-to-nearest instead of rounding down */
static inline int32_t PSHR(int32_t a, uint32_t shift){return (a + ((int32_t)1 << (shift >> 1))) >> shift;} inline int32_t PSHR(int32_t a, uint32_t shift){return (a + ((int32_t)1 << (shift >> 1))) >> shift;}
/** 32-bit arithmetic shift right where the argument can be negative */ /** 32-bit arithmetic shift right where the argument can be negative */
#define VSHR32(a, shift) (((shift)>0) ? SHR32(a, shift) : SHL32(a, -(shift))) #define VSHR32(a, shift) (((shift)>0) ? SHR32(a, shift) : SHL32(a, -(shift)))
...@@ -393,71 +352,67 @@ int32_t celt_rcp(int32_t x); ...@@ -393,71 +352,67 @@ int32_t celt_rcp(int32_t x);
#define LOG_MAX_PSEUDO 6 #define LOG_MAX_PSEUDO 6
#define ALLOC_NONE 1 #define ALLOC_NONE 1
extern const signed char tf_select_table[4][8];
extern const uint32_t SMALL_DIV_TABLE[129];
extern const uint8_t LOG2_FRAC_TABLE[24];
/* Prototypes and inlines*/ /* Prototypes and inlines*/
static inline int16_t SAT16(int32_t x) { inline int16_t SAT16(int32_t x) {
if(x > INT16_MAX) return INT16_MAX; if(x > INT16_MAX) return INT16_MAX;
if(x < INT16_MIN) return INT16_MIN; if(x < INT16_MIN) return INT16_MIN;
return (int16_t)x; return (int16_t)x;
} }
static inline int32_t celt_sudiv(int32_t n, int32_t d) { inline int32_t celt_sudiv(int32_t n, int32_t d) {
assert(d>0); return n/d; assert(d>0); return n/d;
} }
static inline int16_t sig2word16(int32_t x){ inline int16_t sig2word16(int32_t x){
x = PSHR(x, 12); x = PSHR(x, 12);
x = max(x, -32768); x = max(x, -32768);
x = min(x, 32767); x = min(x, 32767);
return (int16_t)(x); return (int16_t)(x);
} }
static inline int32_t ec_tell(ec_ctx *_this){ inline int32_t ec_tell(){
return _this->nbits_total-EC_ILOG(_this->rng); return s_ec.nbits_total-EC_ILOG(s_ec.rng);
} }
/* Atan approximation using a 4th order polynomial. Input is in Q15 format and normalized by pi/4. Output is in /* Atan approximation using a 4th order polynomial. Input is in Q15 format and normalized by pi/4. Output is in
Q15 format */ Q15 format */
static inline int16_t celt_atan01(int16_t x) { inline int16_t celt_atan01(int16_t x) {
return MULT16_16_P15( return MULT16_16_P15(
x, ADD32(32767, MULT16_16_P15(x, ADD32(-21, MULT16_16_P15(x, ADD32(-11943, MULT16_16_P15(4936, x))))))); x, ADD32(32767, MULT16_16_P15(x, ADD32(-21, MULT16_16_P15(x, ADD32(-11943, MULT16_16_P15(4936, x)))))));
} }
/* atan2() approximation valid for positive input values */ /* atan2() approximation valid for positive input values */
static inline int16_t celt_atan2p(int16_t y, int16_t x) { inline int16_t celt_atan2p(int16_t y, int16_t x) {
if (y < x) { if(y < x) {
int32_t arg; int32_t arg;
arg = celt_div(SHL32(EXTEND32(y), 15), x); arg = celt_div(SHL32(EXTEND32(y), 15), x);
if (arg >= 32767) arg = 32767; if(arg >= 32767) arg = 32767;
return SHR16(celt_atan01((int16_t)(arg)), 1); return SHR16(celt_atan01((int16_t)(arg)), 1);
} else { } else {
int32_t arg; int32_t arg;
arg = celt_div(SHL32(EXTEND32(x), 15), y); arg = celt_div(SHL32(EXTEND32(x), 15), y);
if (arg >= 32767) arg = 32767; if(arg >= 32767) arg = 32767;
return 25736 - SHR16(celt_atan01((int16_t)(arg)), 1); return 25736 - SHR16(celt_atan01((int16_t)(arg)), 1);
} }
} }
static inline int32_t celt_maxabs16(const int16_t *x, int32_t len) { inline int32_t celt_maxabs16(const int16_t *x, int32_t len) {
int32_t i; int32_t i;
int16_t maxval = 0; int16_t maxval = 0;
int16_t minval = 0; int16_t minval = 0;
for (i = 0; i < len; i++) { for(i = 0; i < len; i++) {
maxval = max(maxval, x[i]); maxval = max(maxval, x[i]);
minval = min(minval, x[i]); minval = min(minval, x[i]);
} }
return max(EXTEND32(maxval), -EXTEND32(minval)); return max(EXTEND32(maxval), -EXTEND32(minval));
} }
static inline int32_t celt_maxabs32(const int32_t *x, int32_t len) { inline int32_t celt_maxabs32(const int32_t *x, int32_t len) {
int32_t i; int32_t i;
int32_t maxval = 0; int32_t maxval = 0;
int32_t minval = 0; int32_t minval = 0;
for (i = 0; i < len; i++) { for(i = 0; i < len; i++) {
maxval = max(maxval, x[i]); maxval = max(maxval, x[i]);
minval = min(minval, x[i]); minval = min(minval, x[i]);
} }
...@@ -465,39 +420,35 @@ static inline int32_t celt_maxabs32(const int32_t *x, int32_t len) { ...@@ -465,39 +420,35 @@ static inline int32_t celt_maxabs32(const int32_t *x, int32_t len) {
} }
/** Integer log in base2. Undefined for zero and negative numbers */ /** Integer log in base2. Undefined for zero and negative numbers */
static inline int16_t celt_ilog2(int32_t x) { inline int16_t celt_ilog2(uint32_t x) {
assert(x > 0); assert(x > 0);
return EC_ILOG(x) - 1; return EC_ILOG(x) - 1;
} }
/** Integer log in base2. Defined for zero, but not for negative numbers */ /** Integer log in base2. Defined for zero, but not for negative numbers */
static inline int16_t celt_zlog2(int32_t x) { return x <= 0 ? 0 : celt_ilog2(x); } inline int16_t celt_zlog2(uint32_t x) { return x <= 0 ? 0 : celt_ilog2(x); }
/** Base-2 logarithm approximation (log2(x)). (Q14 input, Q10 output) */ /** Base-2 logarithm approximation (log2(x)). (Q14 input, Q10 output) */
static inline int16_t celt_log2(int32_t x) { inline int16_t celt_log2(int32_t x) {
int32_t i; int32_t i;
int16_t n, frac; int16_t n, frac, var1;
/* -0.41509302963303146, 0.9609890551383969, -0.31836011537636605, /* -0.41509302963303146, 0.9609890551383969, -0.31836011537636605, 0.15530808010959576, -0.08556153059057618 */
0.15530808010959576, -0.08556153059057618 */
static const int16_t C[5] = {-6801 + (1 << 3), 15746, -5217, 2545, -1401}; static const int16_t C[5] = {-6801 + (1 << 3), 15746, -5217, 2545, -1401};
if (x == 0) return -32767; if(x == 0) return -32767;
i = celt_ilog2(x); i = celt_ilog2(x);
n = VSHR32(x, i - 15) - 32768 - 16384; n = VSHR32(x, i - 15) - 32768 - 16384;
frac = ADD16( var1 = MULT16_16_Q15(n, ADD16(C[3], MULT16_16_Q15(n, C[4])));
C[0], frac = ADD16(C[0], MULT16_16_Q15(n, ADD16(C[1], MULT16_16_Q15(n, ADD16(C[2], var1)))));
MULT16_16_Q15(
n, ADD16(C[1], MULT16_16_Q15(n, ADD16(C[2], MULT16_16_Q15(n, ADD16(C[3], MULT16_16_Q15(n, C[4]))))))));
return SHL16(i - 13, 10) + SHR16(frac, 14 - 10); return SHL16(i - 13, 10) + SHR16(frac, 14 - 10);
} }
static inline int32_t celt_exp2_frac(int16_t x) { inline int32_t celt_exp2_frac(int16_t x) {
int16_t frac; int16_t frac = SHL16(x, 4);
frac = SHL16(x, 4); int16_t var1 = ADD16(14819, MULT16_16_Q15(10204, frac));
return ADD16(16383, return ADD16(16383, MULT16_16_Q15(frac, ADD16(22804, MULT16_16_Q15(frac, var1))));
MULT16_16_Q15(frac, ADD16(22804, MULT16_16_Q15(frac, ADD16(14819, MULT16_16_Q15(10204, frac))))));
} }
/** Base-2 exponential approximation (2^x). (Q10 input, Q16 output) */ /** Base-2 exponential approximation (2^x). (Q10 input, Q16 output) */
static inline int32_t celt_exp2(int16_t x) { inline int32_t celt_exp2(int16_t x) {
int32_t integer; int32_t integer;
int16_t frac; int16_t frac;
integer = SHR16(x, 10); integer = SHR16(x, 10);
...@@ -509,12 +460,12 @@ static inline int32_t celt_exp2(int16_t x) { ...@@ -509,12 +460,12 @@ static inline int32_t celt_exp2(int16_t x) {
return VSHR32(EXTEND32(frac), -integer - 2); return VSHR32(EXTEND32(frac), -integer - 2);
} }
static inline void dual_inner_prod(const int16_t *x, const int16_t *y01, const int16_t *y02, int32_t N, int32_t *xy1, inline void dual_inner_prod(const int16_t *x, const int16_t *y01, const int16_t *y02, int32_t N, int32_t *xy1,
int32_t *xy2) { int32_t *xy2) {
int32_t i; int32_t i;
int32_t xy01 = 0; int32_t xy01 = 0;
int32_t xy02 = 0; int32_t xy02 = 0;
for (i = 0; i < N; i++) { for(i = 0; i < N; i++) {
xy01 = MAC16_16(xy01, x[i], y01[i]); xy01 = MAC16_16(xy01, x[i], y01[i]);
xy02 = MAC16_16(xy02, x[i], y02[i]); xy02 = MAC16_16(xy02, x[i], y02[i]);
} }
...@@ -522,26 +473,24 @@ static inline void dual_inner_prod(const int16_t *x, const int16_t *y01, const i ...@@ -522,26 +473,24 @@ static inline void dual_inner_prod(const int16_t *x, const int16_t *y01, const i
*xy2 = xy02; *xy2 = xy02;
} }
/*We make sure a C version is always available for cases where the overhead of vectorization and passing around an inline uint32_t celt_inner_prod(const int16_t *x, const int16_t *y, int32_t N) {
arch flag aren't worth it.*/ int i;
static inline int32_t celt_inner_prod(const int16_t *x, const int16_t *y, int32_t N) { uint32_t xy = 0;
int32_t i; for (i = 0; i < N; i++) xy = (int32_t)x[i] * (int32_t)y[i] + xy;
int32_t xy = 0;
for (i = 0; i < N; i++) xy = MAC16_16(xy, x[i], y[i]);
return xy; return xy;
} }
static inline int32_t get_pulses(int32_t i){ inline int32_t get_pulses(int32_t i){
return i<8 ? i : (8 + (i&7)) << ((i>>3)-1); return i<8 ? i : (8 + (i&7)) << ((i>>3)-1);
} }
static inline int32_t bits2pulses(const CELTMode *m, int32_t band, int32_t LM, int32_t bits){ inline int32_t bits2pulses(int32_t band, int32_t LM, int32_t bits){
int32_t i; int32_t i;
int32_t lo, hi; int32_t lo, hi;
const uint8_t *cache; const uint8_t *cache;
LM++; LM++;
cache = m->cache.bits + m->cache.index[LM*m->nbEBands+band]; cache = cache_bits50 + cache_index50[LM * m_CELTMode.nbEBands + band];
lo = 0; lo = 0;
hi = cache[0]; hi = cache[0];
...@@ -561,125 +510,104 @@ static inline int32_t bits2pulses(const CELTMode *m, int32_t band, int32_t LM, i ...@@ -561,125 +510,104 @@ static inline int32_t bits2pulses(const CELTMode *m, int32_t band, int32_t LM, i
return hi; return hi;
} }
static inline int32_t pulses2bits(const CELTMode *m, int32_t band, int32_t LM, int32_t pulses){ inline int32_t pulses2bits(int32_t band, int32_t LM, int32_t pulses){
const uint8_t *cache; const uint8_t *cache;
LM++; LM++;
cache = m->cache.bits + m->cache.index[LM*m->nbEBands+band]; cache = cache_bits50 + cache_index50[LM * m_CELTMode.nbEBands + band];
return pulses == 0 ? 0 : cache[pulses]+1; return pulses == 0 ? 0 : cache[pulses]+1;
} }
_Pragma("GCC diagnostic push") void comb_filter_const(int32_t *y, int32_t *x, int32_t T, int32_t N, int16_t g10, int16_t g11, int16_t g12);
_Pragma("GCC diagnostic ignored \"-Wunused-function\"") void comb_filter(int32_t *y, int32_t *x, int32_t T0, int32_t T1, int32_t N, int16_t g0, int16_t g1, int32_t tapset0,
int32_t tapset1, int32_t overlap);
int32_t resampling_factor(int32_t rate); void init_caps(int32_t *cap, int32_t LM, int32_t C);
void comb_filter_const(int32_t *y, int32_t *x, int32_t T, int32_t N, int16_t g10, int16_t g11, int16_t g12);
void comb_filter(int32_t *y, int32_t *x, int32_t T0, int32_t T1, int32_t N, int16_t g0, int16_t g1, int32_t tapset0, int32_t tapset1,
const int16_t *window, int32_t overlap);
void init_caps(const CELTMode *m, int32_t *cap, int32_t LM, int32_t C);
uint32_t celt_lcg_rand(uint32_t seed); uint32_t celt_lcg_rand(uint32_t seed);
int16_t bitexact_cos(int16_t x); int16_t bitexact_cos(int16_t x);
int32_t bitexact_log2tan(int32_t isin, int32_t icos); int32_t bitexact_log2tan(int32_t isin, int32_t icos);
void denormalise_bands(const CELTMode *m, const int16_t * X, int32_t * freq, void denormalise_bands(const int16_t *X, int32_t *freq, const int16_t *bandLogE, int32_t start, int32_t end, int32_t M,
const int16_t *bandLogE, int32_t start, int32_t end, int32_t M, int32_t downsample, int32_t silence); int32_t silence);
void anti_collapse(const CELTMode *m, int16_t *X_, uint8_t *collapse_masks, int32_t LM, int32_t C, int32_t size, int32_t start, void anti_collapse(int16_t *X_, uint8_t *collapse_masks, int32_t LM, int32_t C, int32_t size, int32_t start,
int32_t end, const int16_t *logE, const int16_t *prev1logE, const int16_t *prev2logE, const int32_t *pulses, int32_t end, const int16_t *logE, const int16_t *prev1logE, const int16_t *prev2logE,
uint32_t seed, int32_t arch); const int32_t *pulses, uint32_t seed);
static void compute_channel_weights(int32_t Ex, int32_t Ey, int16_t w[2]); void compute_channel_weights(int32_t Ex, int32_t Ey, int16_t w[2]);
static void stereo_split(int16_t * X, int16_t * Y, int32_t N); void stereo_split(int16_t *X, int16_t *Y, int32_t N);
static void stereo_merge(int16_t * X, int16_t * Y, int16_t mid, int32_t N, int32_t arch); void stereo_merge(int16_t *X, int16_t *Y, int16_t mid, int32_t N);
static void deinterleave_hadamard(int16_t *X, int32_t N0, int32_t stride, int32_t hadamard); void deinterleave_hadamard(int16_t *X, int32_t N0, int32_t stride, int32_t hadamard);
static void interleave_hadamard(int16_t *X, int32_t N0, int32_t stride, int32_t hadamard); void interleave_hadamard(int16_t *X, int32_t N0, int32_t stride, int32_t hadamard);
void haar1(int16_t *X, int32_t N0, int32_t stride); void haar1(int16_t *X, int32_t N0, int32_t stride);
static int32_t compute_qn(int32_t N, int32_t b, int32_t offset, int32_t pulse_cap, int32_t stereo); int32_t compute_qn(int32_t N, int32_t b, int32_t offset, int32_t pulse_cap, int32_t stereo);
static void compute_theta(struct band_ctx *ctx, struct split_ctx *sctx, int16_t *X, int16_t *Y, int32_t N, int32_t *b, int32_t B, void compute_theta(struct split_ctx *sctx, int16_t *X, int16_t *Y, int32_t N, int32_t *b, int32_t B, int32_t __B0,
int32_t __B0, int32_t LM, int32_t stereo, int32_t *fill); int32_t LM, int32_t stereo, int32_t *fill);
static uint32_t quant_band_n1(struct band_ctx *ctx, int16_t *X, int16_t *Y, int32_t b, int16_t *lowband_out); uint32_t quant_band_n1(int16_t *X, int16_t *Y, int32_t b, int16_t *lowband_out);
static uint32_t quant_partition(struct band_ctx *ctx, int16_t *X, int32_t N, int32_t b, int32_t B, int16_t *lowband, int32_t LM, uint32_t quant_partition(int16_t *X, int32_t N, int32_t b, int32_t B, int16_t *lowband, int32_t LM, int16_t gain,
int16_t gain, int32_t fill); int32_t fill);
static uint32_t quant_band(struct band_ctx *ctx, int16_t *X, int32_t N, int32_t b, int32_t B, int16_t *lowband, int32_t LM, uint32_t quant_band(int16_t *X, int32_t N, int32_t b, int32_t B, int16_t *lowband, int32_t LM, int16_t *lowband_out,
int16_t *lowband_out, int16_t gain, int16_t *lowband_scratch, int32_t fill); int16_t gain, int16_t *lowband_scratch, int32_t fill);
static uint32_t quant_band_stereo(struct band_ctx *ctx, int16_t *X, int16_t *Y, int32_t N, int32_t b, int32_t B, int16_t *lowband, uint32_t quant_band_stereo(int16_t *X, int16_t *Y, int32_t N, int32_t b, int32_t B, int16_t *lowband, int32_t LM,
int32_t LM, int16_t *lowband_out, int16_t *lowband_scratch, int32_t fill); int16_t *lowband_out, int16_t *lowband_scratch, int32_t fill);
static void special_hybrid_folding(const CELTMode *m, int16_t *norm, int16_t *norm2, int32_t start, int32_t M, int32_t dual_stereo); void special_hybrid_folding(int16_t *norm, int16_t *norm2, int32_t start, int32_t M, int32_t dual_stereo);
void quant_all_bands(const CELTMode *m, int32_t start, int32_t end, int16_t *X_, int16_t *Y_, void quant_all_bands(int32_t start, int32_t end, int16_t *X_, int16_t *Y_, uint8_t *collapse_masks,
uint8_t *collapse_masks, const int32_t *bandE, int32_t *pulses, int32_t shortBlocks, int32_t spread, const int32_t *bandE, int32_t *pulses, int32_t shortBlocks, int32_t spread, int32_t dual_stereo,
int32_t dual_stereo, int32_t intensity, int32_t *tf_res, int32_t total_bits, int32_t balance, ec_ctx *ec, int32_t intensity, int32_t *tf_res, int32_t total_bits, int32_t balance, int32_t LM,
int32_t LM, int32_t codedBands, uint32_t *seed, int32_t complexity, int32_t arch, int32_t disable_inv); int32_t codedBands, uint32_t *seed, int32_t complexity, int32_t disable_inv);
int32_t opus_custom_decoder_get_size(const CELTMode *mode, int32_t channels); int32_t celt_decoder_get_size(int32_t channels);
int32_t celt_decoder_get_size(int32_t channels); int32_t celt_decoder_init(int32_t channels);
int32_t opus_custom_decoder_init(CELTDecoder *st, const CELTMode *mode, int32_t channels); void deemphasis_stereo_simple(int32_t *in[], int16_t *pcm, int32_t N, const int16_t coef0, int32_t *mem);
int32_t celt_decoder_init(CELTDecoder *st, int32_t sampling_rate, int32_t channels); void deemphasis(int32_t *in[], int16_t *pcm, int32_t N, int32_t C, const int16_t *coef, int32_t *mem);
static void deemphasis_stereo_simple(int32_t *in[], int16_t *pcm, int32_t N, const int16_t coef0, int32_t *mem); void celt_synthesis(int16_t *X, int32_t *out_syn[], int16_t *oldBandE, int32_t start, int32_t effEnd, int32_t C,
static void deemphasis(int32_t *in[], int16_t *pcm, int32_t N, int32_t C, int32_t downsample, const int16_t *coef, int32_t CC, int32_t isTransient, int32_t LM, int32_t silence);
int32_t *mem, int32_t accum); void tf_decode(int32_t start, int32_t end, int32_t isTransient, int32_t *tf_res, int32_t LM);
static void celt_synthesis(const CELTMode *mode, int16_t *X, int32_t *out_syn[], int16_t *oldBandE, int32_t start, int32_t celt_decode_with_ec(const uint8_t *inbuf, int32_t len, int16_t *outbuf, int32_t frame_size);
int32_t effEnd, int32_t C, int32_t CC, int32_t isTransient, int32_t LM, int32_t downsample, int32_t silence, int32_t arch); int32_t celt_decoder_ctl(int32_t request, ...);
static void tf_decode(int32_t start, int32_t end, int32_t isTransient, int32_t *tf_res, int32_t LM, ec_dec *dec); int32_t cwrsi(int32_t _n, int32_t _k, uint32_t _i, int32_t *_y);
static void celt_decode_lost(CELTDecoder * st, int32_t N, int32_t LM); int32_t decode_pulses(int32_t *_y, int32_t _n, int32_t _k);
int32_t celt_decode_with_ec(CELTDecoder *st, const uint8_t *data, int32_t len, int16_t * pcm, uint32_t ec_tell_frac();
int32_t frame_size, ec_dec *dec, int32_t accum); int32_t ec_read_byte();
int32_t celt_decoder_ctl(CELTDecoder * st, int32_t request, ...); int32_t ec_read_byte_from_end();
void celt_fir(const int16_t *x, const int16_t *num, int16_t *y, int32_t N, int32_t ord); void ec_dec_normalize();
void celt_iir(const int32_t *_x, const int16_t *den, int32_t *_y, int32_t N, int32_t ord, int16_t *mem, int32_t arch); void ec_dec_init(uint8_t *_buf, uint32_t _storage);
int32_t _celt_autocorr(const int16_t *x, int32_t *ac, const int16_t *window, int32_t overlap, int32_t lag, int32_t n, int32_t arch); uint32_t ec_decode(uint32_t _ft);
static int32_t cwrsi(int32_t _n, int32_t _k, uint32_t _i, int32_t *_y); uint32_t ec_decode_bin(uint32_t _bits);
int32_t decode_pulses(int32_t *_y, int32_t _n, int32_t _k, ec_dec *_dec); void ec_dec_update(uint32_t _fl, uint32_t _fh, uint32_t _ft);
uint32_t ec_tell_frac(ec_ctx *_this); int32_t ec_dec_bit_logp(uint32_t _logp);
static int32_t ec_read_byte(ec_dec *_this); int32_t ec_dec_icdf(const uint8_t *_icdf, uint32_t _ftb);
static int32_t ec_read_byte_from_end(ec_dec *_this); uint32_t ec_dec_uint(uint32_t _ft);
static void ec_dec_normalize(ec_dec *_this); uint32_t ec_dec_bits(uint32_t _bits);
void ec_dec_init(ec_dec *_this, uint8_t *_buf, uint32_t _storage); void kf_bfly2(kiss_fft_cpx *Fout, int32_t m, int32_t N);
uint32_t ec_decode(ec_dec *_this, uint32_t _ft); void kf_bfly4(kiss_fft_cpx *Fout, const size_t fstride, const kiss_fft_state *st, int32_t m, int32_t N, int32_t mm);
uint32_t ec_decode_bin(ec_dec *_this, uint32_t _bits); void kf_bfly3(kiss_fft_cpx *Fout, const size_t fstride, const kiss_fft_state *st, int32_t m, int32_t N, int32_t mm);
void ec_dec_update(ec_dec *_this, uint32_t _fl, uint32_t _fh, uint32_t _ft); void kf_bfly5(kiss_fft_cpx *Fout, const size_t fstride, const kiss_fft_state *st, int32_t m, int32_t N, int32_t mm);
int32_t ec_dec_bit_logp(ec_dec *_this, uint32_t _logp); void opus_fft_impl(const kiss_fft_state *st, kiss_fft_cpx *fout);
int32_t ec_dec_icdf(ec_dec *_this, const uint8_t *_icdf, uint32_t _ftb); uint32_t ec_laplace_get_freq1(uint32_t fs0, int32_t decay);
uint32_t ec_dec_uint(ec_dec *_this, uint32_t _ft); int32_t ec_laplace_decode(uint32_t fs, int32_t decay);
uint32_t ec_dec_bits(ec_dec *_this, uint32_t _bits);
static void kf_bfly2(kiss_fft_cpx *Fout, int32_t m, int32_t N);
static void kf_bfly4(kiss_fft_cpx *Fout, const size_t fstride, const kiss_fft_state *st, int32_t m, int32_t N, int32_t mm);
static void kf_bfly3(kiss_fft_cpx *Fout, const size_t fstride, const kiss_fft_state *st, int32_t m, int32_t N, int32_t mm);
static void kf_bfly5(kiss_fft_cpx *Fout, const size_t fstride, const kiss_fft_state *st, int32_t m, int32_t N, int32_t mm);
void opus_fft_impl(const kiss_fft_state *st, kiss_fft_cpx *fout);
static uint32_t ec_laplace_get_freq1(uint32_t fs0, int32_t decay);
int32_t ec_laplace_decode(ec_dec *dec, uint32_t fs, int32_t decay);
uint32_t isqrt32(uint32_t _val); uint32_t isqrt32(uint32_t _val);
int32_t frac_div32(int32_t a, int32_t b); int16_t celt_rsqrt_norm(int32_t x);
int16_t celt_rsqrt_norm(int32_t x); int32_t celt_sqrt(int32_t x);
int32_t celt_sqrt(int32_t x); int16_t celt_cos_norm(int32_t x);
int16_t celt_cos_norm(int32_t x); int32_t celt_rcp(int32_t x);
int32_t celt_rcp(int32_t x); void clt_mdct_backward(int32_t *in, int32_t *out, int32_t overlap, int32_t shift, int32_t stride);
void clt_mdct_backward(const mdct_lookup *l, int32_t *in, int32_t * out, const int16_t * window, int32_t overlap, void exp_rotation1(int16_t *X, int32_t len, int32_t stride, int16_t c, int16_t s);
int32_t shift, int32_t stride); void exp_rotation(int16_t *X, int32_t len, int32_t dir, int32_t stride, int32_t K, int32_t spread);
CELTMode *opus_custom_mode_create(int32_t Fs, int32_t frame_size, int32_t *error); void normalise_residual(int32_t *iy, int16_t *X, int32_t N, int32_t Ryy, int16_t gain);
static void exp_rotation1(int16_t *X, int32_t len, int32_t stride, int16_t c, int16_t s); uint32_t extract_collapse_mask(int32_t *iy, int32_t N, int32_t B);
void exp_rotation(int16_t *X, int32_t len, int32_t dir, int32_t stride, int32_t K, int32_t spread); uint32_t alg_unquant(int16_t *X, int32_t N, int32_t K, int32_t spread, int32_t B, int16_t gain);
static void normalise_residual(int32_t * iy, int16_t * X, int32_t N, int32_t Ryy, int16_t gain); void renormalise_vector(int16_t *X, int32_t N, int16_t gain);
static uint32_t extract_collapse_mask(int32_t *iy, int32_t N, int32_t B); int32_t interp_bits2pulses(int32_t start, int32_t end, int32_t skip_start, const int32_t *bits1, const int32_t *bits2,
uint32_t alg_unquant(int16_t *X, int32_t N, int32_t K, int32_t spread, int32_t B, ec_dec *dec, int16_t gain); const int32_t *thresh, const int32_t *cap, int32_t total, int32_t *_balance,
void renormalise_vector(int16_t *X, int32_t N, int16_t gain, int32_t arch); int32_t skip_rsv, int32_t *intensity, int32_t intensity_rsv, int32_t *dual_stereo,
int32_t dual_stereo_rsv, int32_t *bits, int32_t *ebits, int32_t *fine_priority, int32_t C,
int32_t celt_pitch_xcorr(const int16_t *_x, const int16_t *_y, int32_t *xcorr, int32_t len, int32_t max_pitch, int32_t arch); int32_t LM, int32_t prev, int32_t signalBandwidth);
int32_t clt_compute_allocation(int32_t start, int32_t end, const int32_t *offsets, const int32_t *cap,
static int32_t interp_bits2pulses(const CELTMode *m, int32_t start, int32_t end, int32_t skip_start, const int32_t *bits1, const int32_t *bits2, int32_t alloc_trim, int32_t *intensity, int32_t *dual_stereo, int32_t total,
const int32_t *thresh, const int32_t *cap, int32_t total, int32_t *_balance, int32_t skip_rsv, int32_t *balance, int32_t *pulses, int32_t *ebits, int32_t *fine_priority, int32_t C,
int32_t *intensity, int32_t intensity_rsv, int32_t *dual_stereo, int32_t dual_stereo_rsv, int32_t *bits, int32_t LM, int32_t prev, int32_t signalBandwidth);
int32_t *ebits, int32_t *fine_priority, int32_t C, int32_t LM, ec_ctx *ec, int32_t prev, int32_t signalBandwidth); void unquant_coarse_energy(int32_t start, int32_t end, int16_t *oldEBands, int32_t intra, int32_t C, int32_t LM);
int32_t clt_compute_allocation(const CELTMode *m, int32_t start, int32_t end, const int32_t *offsets, const int32_t *cap, int32_t alloc_trim, void unquant_fine_energy(int32_t start, int32_t end, int16_t *oldEBands, int32_t *fine_quant, int32_t C);
int32_t *intensity, int32_t *dual_stereo, int32_t total, int32_t *balance, int32_t *pulses, int32_t *ebits, void unquant_energy_finalise(int32_t start, int32_t end, int16_t *oldEBands, int32_t *fine_quant,
int32_t *fine_priority, int32_t C, int32_t LM, ec_ctx *ec, int32_t prev, int32_t signalBandwidth); int32_t *fine_priority, int32_t bits_left, int32_t C);
void unquant_coarse_energy(const CELTMode *m, int32_t start, int32_t end, int16_t *oldEBands, int32_t intra, ec_dec *dec, int32_t C, uint32_t celt_pvq_u_row(uint32_t row, uint32_t data);
int32_t LM);
void unquant_fine_energy(const CELTMode *m, int32_t start, int32_t end, int16_t *oldEBands, int32_t *fine_quant, ec_dec *dec, bool CELTDecoder_AllocateBuffers(void);
int32_t C); void CELTDecoder_FreeBuffers();
void unquant_energy_finalise(const CELTMode *m, int32_t start, int32_t end, int16_t *oldEBands, int32_t *fine_quant, void CELTDecoder_ClearBuffer(void);
int32_t *fine_priority, int32_t bits_left, ec_dec *dec, int32_t C);
static void xcorr_kernel(const int16_t *x, const int16_t *y, int32_t sum[4], int32_t len);
_Pragma("GCC diagnostic pop")
#ifdef __cplusplus
}
#endif
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment