annotate src/opus-1.3/silk/fixed/main_FIX.h @ 82:d0c2a83c1364

Add FFTW 3.3.8 source, and a Linux build
author Chris Cannam
date Tue, 19 Nov 2019 14:52:55 +0000
parents 7aeed7906520
children
rev   line source
Chris@69 1 /***********************************************************************
Chris@69 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
Chris@69 3 Redistribution and use in source and binary forms, with or without
Chris@69 4 modification, are permitted provided that the following conditions
Chris@69 5 are met:
Chris@69 6 - Redistributions of source code must retain the above copyright notice,
Chris@69 7 this list of conditions and the following disclaimer.
Chris@69 8 - Redistributions in binary form must reproduce the above copyright
Chris@69 9 notice, this list of conditions and the following disclaimer in the
Chris@69 10 documentation and/or other materials provided with the distribution.
Chris@69 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
Chris@69 12 names of specific contributors, may be used to endorse or promote
Chris@69 13 products derived from this software without specific prior written
Chris@69 14 permission.
Chris@69 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
Chris@69 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
Chris@69 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
Chris@69 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
Chris@69 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
Chris@69 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
Chris@69 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
Chris@69 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
Chris@69 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
Chris@69 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
Chris@69 25 POSSIBILITY OF SUCH DAMAGE.
Chris@69 26 ***********************************************************************/
Chris@69 27
Chris@69 28 #ifndef SILK_MAIN_FIX_H
Chris@69 29 #define SILK_MAIN_FIX_H
Chris@69 30
Chris@69 31 #include "SigProc_FIX.h"
Chris@69 32 #include "structs_FIX.h"
Chris@69 33 #include "control.h"
Chris@69 34 #include "main.h"
Chris@69 35 #include "PLC.h"
Chris@69 36 #include "debug.h"
Chris@69 37 #include "entenc.h"
Chris@69 38
Chris@69 39 #if ((defined(OPUS_ARM_ASM) && defined(FIXED_POINT)) \
Chris@69 40 || defined(OPUS_ARM_MAY_HAVE_NEON_INTR))
Chris@69 41 #include "fixed/arm/warped_autocorrelation_FIX_arm.h"
Chris@69 42 #endif
Chris@69 43
Chris@69 44 #ifndef FORCE_CPP_BUILD
Chris@69 45 #ifdef __cplusplus
Chris@69 46 extern "C"
Chris@69 47 {
Chris@69 48 #endif
Chris@69 49 #endif
Chris@69 50
Chris@69 51 #define silk_encoder_state_Fxx silk_encoder_state_FIX
Chris@69 52 #define silk_encode_do_VAD_Fxx silk_encode_do_VAD_FIX
Chris@69 53 #define silk_encode_frame_Fxx silk_encode_frame_FIX
Chris@69 54
Chris@69 55 #define QC 10
Chris@69 56 #define QS 13
Chris@69 57
Chris@69 58 /*********************/
Chris@69 59 /* Encoder Functions */
Chris@69 60 /*********************/
Chris@69 61
Chris@69 62 /* High-pass filter with cutoff frequency adaptation based on pitch lag statistics */
Chris@69 63 void silk_HP_variable_cutoff(
Chris@69 64 silk_encoder_state_Fxx state_Fxx[] /* I/O Encoder states */
Chris@69 65 );
Chris@69 66
Chris@69 67 /* Encoder main function */
Chris@69 68 void silk_encode_do_VAD_FIX(
Chris@69 69 silk_encoder_state_FIX *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 70 opus_int activity /* I Decision of Opus voice activity detector */
Chris@69 71 );
Chris@69 72
Chris@69 73 /* Encoder main function */
Chris@69 74 opus_int silk_encode_frame_FIX(
Chris@69 75 silk_encoder_state_FIX *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 76 opus_int32 *pnBytesOut, /* O Pointer to number of payload bytes; */
Chris@69 77 ec_enc *psRangeEnc, /* I/O compressor data structure */
Chris@69 78 opus_int condCoding, /* I The type of conditional coding to use */
Chris@69 79 opus_int maxBits, /* I If > 0: maximum number of output bits */
Chris@69 80 opus_int useCBR /* I Flag to force constant-bitrate operation */
Chris@69 81 );
Chris@69 82
Chris@69 83 /* Initializes the Silk encoder state */
Chris@69 84 opus_int silk_init_encoder(
Chris@69 85 silk_encoder_state_Fxx *psEnc, /* I/O Pointer to Silk FIX encoder state */
Chris@69 86 int arch /* I Run-time architecture */
Chris@69 87 );
Chris@69 88
Chris@69 89 /* Control the Silk encoder */
Chris@69 90 opus_int silk_control_encoder(
Chris@69 91 silk_encoder_state_Fxx *psEnc, /* I/O Pointer to Silk encoder state */
Chris@69 92 silk_EncControlStruct *encControl, /* I Control structure */
Chris@69 93 const opus_int allow_bw_switch, /* I Flag to allow switching audio bandwidth */
Chris@69 94 const opus_int channelNb, /* I Channel number */
Chris@69 95 const opus_int force_fs_kHz
Chris@69 96 );
Chris@69 97
Chris@69 98 /**************************/
Chris@69 99 /* Noise shaping analysis */
Chris@69 100 /**************************/
Chris@69 101 /* Compute noise shaping coefficients and initial gain values */
Chris@69 102 void silk_noise_shape_analysis_FIX(
Chris@69 103 silk_encoder_state_FIX *psEnc, /* I/O Encoder state FIX */
Chris@69 104 silk_encoder_control_FIX *psEncCtrl, /* I/O Encoder control FIX */
Chris@69 105 const opus_int16 *pitch_res, /* I LPC residual from pitch analysis */
Chris@69 106 const opus_int16 *x, /* I Input signal [ frame_length + la_shape ] */
Chris@69 107 int arch /* I Run-time architecture */
Chris@69 108 );
Chris@69 109
Chris@69 110 /* Autocorrelations for a warped frequency axis */
Chris@69 111 void silk_warped_autocorrelation_FIX_c(
Chris@69 112 opus_int32 *corr, /* O Result [order + 1] */
Chris@69 113 opus_int *scale, /* O Scaling of the correlation vector */
Chris@69 114 const opus_int16 *input, /* I Input data to correlate */
Chris@69 115 const opus_int warping_Q16, /* I Warping coefficient */
Chris@69 116 const opus_int length, /* I Length of input */
Chris@69 117 const opus_int order /* I Correlation order (even) */
Chris@69 118 );
Chris@69 119
Chris@69 120 #if !defined(OVERRIDE_silk_warped_autocorrelation_FIX)
Chris@69 121 #define silk_warped_autocorrelation_FIX(corr, scale, input, warping_Q16, length, order, arch) \
Chris@69 122 ((void)(arch), silk_warped_autocorrelation_FIX_c(corr, scale, input, warping_Q16, length, order))
Chris@69 123 #endif
Chris@69 124
Chris@69 125 /* Calculation of LTP state scaling */
Chris@69 126 void silk_LTP_scale_ctrl_FIX(
Chris@69 127 silk_encoder_state_FIX *psEnc, /* I/O encoder state */
Chris@69 128 silk_encoder_control_FIX *psEncCtrl, /* I/O encoder control */
Chris@69 129 opus_int condCoding /* I The type of conditional coding to use */
Chris@69 130 );
Chris@69 131
Chris@69 132 /**********************************************/
Chris@69 133 /* Prediction Analysis */
Chris@69 134 /**********************************************/
Chris@69 135 /* Find pitch lags */
Chris@69 136 void silk_find_pitch_lags_FIX(
Chris@69 137 silk_encoder_state_FIX *psEnc, /* I/O encoder state */
Chris@69 138 silk_encoder_control_FIX *psEncCtrl, /* I/O encoder control */
Chris@69 139 opus_int16 res[], /* O residual */
Chris@69 140 const opus_int16 x[], /* I Speech signal */
Chris@69 141 int arch /* I Run-time architecture */
Chris@69 142 );
Chris@69 143
Chris@69 144 /* Find LPC and LTP coefficients */
Chris@69 145 void silk_find_pred_coefs_FIX(
Chris@69 146 silk_encoder_state_FIX *psEnc, /* I/O encoder state */
Chris@69 147 silk_encoder_control_FIX *psEncCtrl, /* I/O encoder control */
Chris@69 148 const opus_int16 res_pitch[], /* I Residual from pitch analysis */
Chris@69 149 const opus_int16 x[], /* I Speech signal */
Chris@69 150 opus_int condCoding /* I The type of conditional coding to use */
Chris@69 151 );
Chris@69 152
Chris@69 153 /* LPC analysis */
Chris@69 154 void silk_find_LPC_FIX(
Chris@69 155 silk_encoder_state *psEncC, /* I/O Encoder state */
Chris@69 156 opus_int16 NLSF_Q15[], /* O NLSFs */
Chris@69 157 const opus_int16 x[], /* I Input signal */
Chris@69 158 const opus_int32 minInvGain_Q30 /* I Inverse of max prediction gain */
Chris@69 159 );
Chris@69 160
Chris@69 161 /* LTP analysis */
Chris@69 162 void silk_find_LTP_FIX(
Chris@69 163 opus_int32 XXLTP_Q17[ MAX_NB_SUBFR * LTP_ORDER * LTP_ORDER ], /* O Correlation matrix */
Chris@69 164 opus_int32 xXLTP_Q17[ MAX_NB_SUBFR * LTP_ORDER ], /* O Correlation vector */
Chris@69 165 const opus_int16 r_lpc[], /* I Residual signal after LPC */
Chris@69 166 const opus_int lag[ MAX_NB_SUBFR ], /* I LTP lags */
Chris@69 167 const opus_int subfr_length, /* I Subframe length */
Chris@69 168 const opus_int nb_subfr, /* I Number of subframes */
Chris@69 169 int arch /* I Run-time architecture */
Chris@69 170 );
Chris@69 171
Chris@69 172 void silk_LTP_analysis_filter_FIX(
Chris@69 173 opus_int16 *LTP_res, /* O LTP residual signal of length MAX_NB_SUBFR * ( pre_length + subfr_length ) */
Chris@69 174 const opus_int16 *x, /* I Pointer to input signal with at least max( pitchL ) preceding samples */
Chris@69 175 const opus_int16 LTPCoef_Q14[ LTP_ORDER * MAX_NB_SUBFR ],/* I LTP_ORDER LTP coefficients for each MAX_NB_SUBFR subframe */
Chris@69 176 const opus_int pitchL[ MAX_NB_SUBFR ], /* I Pitch lag, one for each subframe */
Chris@69 177 const opus_int32 invGains_Q16[ MAX_NB_SUBFR ], /* I Inverse quantization gains, one for each subframe */
Chris@69 178 const opus_int subfr_length, /* I Length of each subframe */
Chris@69 179 const opus_int nb_subfr, /* I Number of subframes */
Chris@69 180 const opus_int pre_length /* I Length of the preceding samples starting at &x[0] for each subframe */
Chris@69 181 );
Chris@69 182
Chris@69 183 /* Calculates residual energies of input subframes where all subframes have LPC_order */
Chris@69 184 /* of preceding samples */
Chris@69 185 void silk_residual_energy_FIX(
Chris@69 186 opus_int32 nrgs[ MAX_NB_SUBFR ], /* O Residual energy per subframe */
Chris@69 187 opus_int nrgsQ[ MAX_NB_SUBFR ], /* O Q value per subframe */
Chris@69 188 const opus_int16 x[], /* I Input signal */
Chris@69 189 opus_int16 a_Q12[ 2 ][ MAX_LPC_ORDER ], /* I AR coefs for each frame half */
Chris@69 190 const opus_int32 gains[ MAX_NB_SUBFR ], /* I Quantization gains */
Chris@69 191 const opus_int subfr_length, /* I Subframe length */
Chris@69 192 const opus_int nb_subfr, /* I Number of subframes */
Chris@69 193 const opus_int LPC_order, /* I LPC order */
Chris@69 194 int arch /* I Run-time architecture */
Chris@69 195 );
Chris@69 196
Chris@69 197 /* Residual energy: nrg = wxx - 2 * wXx * c + c' * wXX * c */
Chris@69 198 opus_int32 silk_residual_energy16_covar_FIX(
Chris@69 199 const opus_int16 *c, /* I Prediction vector */
Chris@69 200 const opus_int32 *wXX, /* I Correlation matrix */
Chris@69 201 const opus_int32 *wXx, /* I Correlation vector */
Chris@69 202 opus_int32 wxx, /* I Signal energy */
Chris@69 203 opus_int D, /* I Dimension */
Chris@69 204 opus_int cQ /* I Q value for c vector 0 - 15 */
Chris@69 205 );
Chris@69 206
Chris@69 207 /* Processing of gains */
Chris@69 208 void silk_process_gains_FIX(
Chris@69 209 silk_encoder_state_FIX *psEnc, /* I/O Encoder state */
Chris@69 210 silk_encoder_control_FIX *psEncCtrl, /* I/O Encoder control */
Chris@69 211 opus_int condCoding /* I The type of conditional coding to use */
Chris@69 212 );
Chris@69 213
Chris@69 214 /******************/
Chris@69 215 /* Linear Algebra */
Chris@69 216 /******************/
Chris@69 217 /* Calculates correlation matrix X'*X */
Chris@69 218 void silk_corrMatrix_FIX(
Chris@69 219 const opus_int16 *x, /* I x vector [L + order - 1] used to form data matrix X */
Chris@69 220 const opus_int L, /* I Length of vectors */
Chris@69 221 const opus_int order, /* I Max lag for correlation */
Chris@69 222 opus_int32 *XX, /* O Pointer to X'*X correlation matrix [ order x order ] */
Chris@69 223 opus_int32 *nrg, /* O Energy of x vector */
Chris@69 224 opus_int *rshifts, /* O Right shifts of correlations */
Chris@69 225 int arch /* I Run-time architecture */
Chris@69 226 );
Chris@69 227
Chris@69 228 /* Calculates correlation vector X'*t */
Chris@69 229 void silk_corrVector_FIX(
Chris@69 230 const opus_int16 *x, /* I x vector [L + order - 1] used to form data matrix X */
Chris@69 231 const opus_int16 *t, /* I Target vector [L] */
Chris@69 232 const opus_int L, /* I Length of vectors */
Chris@69 233 const opus_int order, /* I Max lag for correlation */
Chris@69 234 opus_int32 *Xt, /* O Pointer to X'*t correlation vector [order] */
Chris@69 235 const opus_int rshifts, /* I Right shifts of correlations */
Chris@69 236 int arch /* I Run-time architecture */
Chris@69 237 );
Chris@69 238
Chris@69 239 #ifndef FORCE_CPP_BUILD
Chris@69 240 #ifdef __cplusplus
Chris@69 241 }
Chris@69 242 #endif /* __cplusplus */
Chris@69 243 #endif /* FORCE_CPP_BUILD */
Chris@69 244 #endif /* SILK_MAIN_FIX_H */