Update SILK code using the CELT range coder
[opus.git] / src_SigProc_FLP / SKP_Silk_SigProc_FLP.h
1 /***********************************************************************\r
2 Copyright (c) 2006-2010, Skype Limited. All rights reserved. \r
3 Redistribution and use in source and binary forms, with or without \r
4 modification, (subject to the limitations in the disclaimer below) \r
5 are permitted provided that the following conditions are met:\r
6 - Redistributions of source code must retain the above copyright notice,\r
7 this list of conditions and the following disclaimer.\r
8 - Redistributions in binary form must reproduce the above copyright \r
9 notice, this list of conditions and the following disclaimer in the \r
10 documentation and/or other materials provided with the distribution.\r
11 - Neither the name of Skype Limited, nor the names of specific \r
12 contributors, may be used to endorse or promote products derived from \r
13 this software without specific prior written permission.\r
14 NO EXPRESS OR IMPLIED LICENSES TO ANY PARTY'S PATENT RIGHTS ARE GRANTED \r
15 BY THIS LICENSE. THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND \r
16 CONTRIBUTORS ''AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING,\r
17 BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND \r
18 FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE \r
19 COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, \r
20 INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, BUT\r
21 NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; LOSS OF \r
22 USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON \r
23 ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT \r
24 (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE \r
25 OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.\r
26 ***********************************************************************/\r
27 \r
28 #ifndef _SKP_SILK_SIGPROC_FLP_H_\r
29 #define _SKP_SILK_SIGPROC_FLP_H_\r
30 \r
31 #include "SKP_Silk_SigProc_FIX.h"\r
32 #include <math.h>\r
33 \r
34 #ifdef  __cplusplus\r
35 extern "C"\r
36 {\r
37 #endif\r
38 \r
39 /********************************************************************/\r
40 /*                    SIGNAL PROCESSING FUNCTIONS                   */\r
41 /********************************************************************/\r
42 \r
43 /* first-order allpass filter */\r
44 void SKP_Silk_allpass_int_FLP(\r
45     const SKP_float *in,                /* I:   input signal [len]              */\r
46     SKP_float       *S,                 /* I/O: state [1]                       */\r
47     SKP_float       A,                  /* I:   coefficient     (0 <= A < 1)        */\r
48     SKP_float       *out,               /* O:   output signal [len]             */\r
49     const SKP_int32 len                 /* I:   number of samples               */\r
50 );\r
51 \r
52 /* downsample by a factor 2, coarser */\r
53 void SKP_Silk_decimate2_coarse_FLP(\r
54     const SKP_float     *in,                 /* I:  signal [2*len]            */\r
55     SKP_float       *S,                  /* I/O: state vector [2]         */\r
56     SKP_float       *out,                /* O:  decimated signal [len]    */\r
57     SKP_float       *scratch,            /* I:  scratch memory [3*len]    */\r
58     const SKP_int32 len                  /* I:  number of OUTPUT samples  */\r
59 );\r
60 \r
61 /* downsample by a factor 2, coarsest */\r
62 void SKP_Silk_decimate2_coarsest_FLP(\r
63     const SKP_float     *in,                 /* I:  signal [2*len]            */\r
64     SKP_float       *S,                  /* I/O: state vector [2]         */\r
65     SKP_float       *out,                /* O:  decimated signal [len]    */\r
66     SKP_float       *scratch,            /* I:  scratch memory [3*len]    */\r
67     const SKP_int32 len                  /* I:  number of OUTPUT samples  */\r
68 );\r
69 \r
70 /* Chirp (bw expand) LP AR filter */\r
71 void SKP_Silk_bwexpander_FLP( \r
72     SKP_float *ar,                     /* io   AR filter to be expanded (without leading 1)    */\r
73     const SKP_int d,                   /* i     length of ar                                       */\r
74     const SKP_float chirp              /* i     chirp factor (typically in range (0..1) )          */\r
75 );\r
76 \r
77 /* compute inverse of LPC prediction gain, and                                                  */\r
78 /* test if LPC coefficients are stable (all poles within unit circle)   */\r
79 /* this code is based on SKP_Silk_FLP_a2k()                                                             */\r
80 SKP_int SKP_Silk_LPC_inverse_pred_gain_FLP( /* O:   returns 1 if unstable, otherwise 0    */\r
81     SKP_float            *invGain,      /* O:   inverse prediction gain, energy domain    */\r
82     const SKP_float      *A,            /* I:   prediction coefficients [order]           */\r
83     SKP_int32            order          /* I:   prediction order                          */\r
84 );\r
85 \r
86 void SKP_Silk_schur_FLP(\r
87     SKP_float           refl_coef[],     /* O   reflection coefficients (length order)       */\r
88     const SKP_float     auto_corr[],     /* I   autotcorreation sequence (length order+1)    */\r
89     SKP_int             order            /* I   order                                        */\r
90 );\r
91 \r
92 void SKP_Silk_k2a_FLP(\r
93     SKP_float           *A,             /* O:   prediction coefficients [order]           */\r
94     const SKP_float     *rc,            /* I:   reflection coefficients [order]           */\r
95     SKP_int32           order           /* I:   prediction order                          */\r
96 );\r
97 \r
98 /* Solve the normal equations using the Levinson-Durbin recursion */\r
99 SKP_float SKP_Silk_levinsondurbin_FLP(  /* O    prediction error energy                                         */\r
100         SKP_float               A[],                            /* O    prediction coefficients [order]                         */\r
101         const SKP_float corr[],                         /* I    input auto-correlations [order + 1]                     */\r
102         const SKP_int   order                           /* I    prediction order                                                        */\r
103 );\r
104 \r
105 /* compute autocorrelation */\r
106 void SKP_Silk_autocorrelation_FLP( \r
107     SKP_float *results,                 /* o    result (length correlationCount)            */\r
108     const SKP_float *inputData,         /* i    input data to correlate                     */\r
109     SKP_int inputDataSize,              /* i    length of input                             */\r
110     SKP_int correlationCount            /* i    number of correlation taps to compute       */\r
111 );\r
112 \r
113 /* Pitch estimator */\r
114 #define SigProc_PE_MIN_COMPLEX        0\r
115 #define SigProc_PE_MID_COMPLEX        1\r
116 #define SigProc_PE_MAX_COMPLEX        2\r
117 \r
118 SKP_int SKP_Silk_pitch_analysis_core_FLP( /* O voicing estimate: 0 voiced, 1 unvoiced                       */\r
119     const SKP_float *signal,            /* I signal of length PE_FRAME_LENGTH_MS*Fs_kHz              */\r
120     SKP_int         *pitch_out,         /* O 4 pitch lag values                                             */\r
121     SKP_int         *lagIndex,          /* O lag Index                                                      */\r
122     SKP_int         *contourIndex,      /* O pitch contour Index                                            */\r
123     SKP_float       *LTPCorr,           /* I/O normalized correlation; input: value from previous frame     */\r
124     SKP_int         prevLag,            /* I last lag of previous frame; set to zero is unvoiced            */\r
125     const SKP_float search_thres1,      /* I first stage threshold for lag candidates 0 - 1                 */\r
126     const SKP_float search_thres2,      /* I final threshold for lag candidates 0 - 1                       */\r
127     const SKP_int   Fs_kHz,             /* I sample frequency (kHz)                                         */\r
128     const SKP_int   complexity,         /* I Complexity setting, 0-2, where 2 is highest                    */\r
129     const SKP_int   nb_subfr            /* I    number of 5 ms subframes                                    */\r
130 );\r
131 \r
132 #define PI               (3.1415926536f)\r
133 \r
134 void SKP_Silk_insertion_sort_decreasing_FLP(\r
135     SKP_float            *a,          /* I/O:  Unsorted / Sorted vector                */\r
136     SKP_int              *index,      /* O:    Index vector for the sorted elements    */\r
137     const SKP_int        L,           /* I:    Vector length                           */\r
138     const SKP_int        K            /* I:    Number of correctly sorted positions    */\r
139 );\r
140 \r
141 void SKP_Silk_insertion_sort_increasing_FLP(\r
142     SKP_float            *a,          /* I/O:  Unsorted / Sorted vector                */\r
143     SKP_int              *index,      /* O:    Index vector for the sorted elements    */\r
144     const SKP_int        L,           /* I:    Vector length                           */\r
145     const SKP_int        K            /* I:    Number of correctly sorted positions    */\r
146 );\r
147 \r
148 /* Laroia low complexity NLSF weights */\r
149 void SKP_Silk_NLSF_VQ_weights_laroia_FLP(       \r
150     SKP_float            *pXW,           /* 0: Pointer to input vector weights                [D x 1]    */\r
151     const SKP_float      *pX,            /* I: Pointer to input vector                        [D x 1]    */ \r
152     const SKP_int        D               /* I: Input vector dimension                                    */\r
153 );\r
154 \r
155 /* Compute reflection coefficients from input signal */\r
156 SKP_float SKP_Silk_burg_modified_FLP(       /* O    returns residual energy                                         */\r
157     SKP_float           A[],                /* O    prediction coefficients (length order)                          */\r
158     const SKP_float     x[],                /* I    input signal, length: nb_subfr*(D+L_sub)                        */\r
159     const SKP_int       subfr_length,       /* I    input signal subframe length (including D preceeding samples)   */\r
160     const SKP_int       nb_subfr,           /* I    number of subframes stacked in x                                */\r
161     const SKP_float     WhiteNoiseFrac,     /* I    fraction added to zero-lag autocorrelation                      */\r
162     const SKP_int       D                   /* I    order                                                           */\r
163 );\r
164 \r
165 /* multiply a vector by a constant */\r
166 void SKP_Silk_scale_vector_FLP( \r
167     SKP_float           *data1,\r
168     SKP_float           gain, \r
169     SKP_int             dataSize\r
170 );\r
171 \r
172 /* copy and multiply a vector by a constant */\r
173 void SKP_Silk_scale_copy_vector_FLP( \r
174     SKP_float           *data_out, \r
175     const SKP_float     *data_in, \r
176     SKP_float           gain, \r
177     SKP_int             dataSize\r
178 );\r
179 \r
180 /* inner product of two SKP_float arrays, with result as double */\r
181 double SKP_Silk_inner_product_FLP( \r
182     const SKP_float     *data1, \r
183     const SKP_float     *data2, \r
184     SKP_int             dataSize\r
185 );\r
186 \r
187 /* sum of squares of a SKP_float array, with result as double */\r
188 double SKP_Silk_energy_FLP( \r
189     const SKP_float     *data, \r
190     SKP_int             dataSize\r
191 );\r
192 \r
193 /********************************************************************/\r
194 /*                                MACROS                                */\r
195 /********************************************************************/\r
196 \r
197 #define SKP_min_float(a, b)                     (((a) < (b)) ? (a) :  (b)) \r
198 #define SKP_max_float(a, b)                     (((a) > (b)) ? (a) :  (b)) \r
199 #define SKP_abs_float(a)                        ((SKP_float)fabs(a))\r
200 \r
201 #define SKP_LIMIT_float( a, limit1, limit2)     ((limit1) > (limit2) ? ((a) > (limit1) ? (limit1) : ((a) < (limit2) ? (limit2) : (a))) \\r
202                                                                                                                              : ((a) > (limit2) ? (limit2) : ((a) < (limit1) ? (limit1) : (a))))\r
203 \r
204 /* sigmoid function */\r
205 SKP_INLINE SKP_float SKP_sigmoid(SKP_float x)\r
206 {\r
207     return (SKP_float)(1.0 / (1.0 + exp(-x)));\r
208 }\r
209 \r
210 /* floating-point to integer conversion (rounding) */\r
211 SKP_INLINE void SKP_float2short_array(\r
212     SKP_int16       *out, \r
213     const SKP_float *in, \r
214     SKP_int32       length\r
215\r
216 {\r
217     SKP_int32 k;\r
218     for (k = length-1; k >= 0; k--) {\r
219         double x = in[k];\r
220         out[k] = (SKP_int16)SKP_SAT16( ( x > 0 ) ? x + 0.5 : x - 0.5 );\r
221     }\r
222 }\r
223 \r
224 /* floating-point to integer conversion (rounding) */\r
225 SKP_INLINE SKP_int32 SKP_float2int(double x) \r
226 {\r
227     return (SKP_int32)( ( x > 0 ) ? x + 0.5 : x - 0.5 );\r
228 }\r
229 \r
230 /* integer to floating-point conversion */\r
231 SKP_INLINE void SKP_short2float_array(\r
232     SKP_float       *out, \r
233     const SKP_int16 *in, \r
234     SKP_int32       length\r
235\r
236 {\r
237     SKP_int32 k;\r
238     for (k = length-1; k >= 0; k--) {\r
239         out[k] = (SKP_float)in[k];\r
240     }\r
241 }\r
242 \r
243 #define SKP_round(x)            (SKP_float)((x)>=0 ? (SKP_int64)((x)+0.5) : (SKP_int64)((x)-0.5))\r
244 \r
245 #ifdef  __cplusplus\r
246 }\r
247 #endif\r
248 \r
249 #endif\r