1 /*
2  *  Copyright (c) 2011 The WebRTC project authors. All Rights Reserved.
3  *
4  *  Use of this source code is governed by a BSD-style license
5  *  that can be found in the LICENSE file in the root of the source
6  *  tree. An additional intellectual property rights grant can be found
7  *  in the file PATENTS.  All contributing project authors may
8  *  be found in the AUTHORS file in the root of the source tree.
9  */
10 
11 /*
12  * Contains the core loop function for the lattice filter AR routine
13  * for iSAC codec.
14  *
15  */
16 
17 #include "common_audio/signal_processing/include/signal_processing_library.h"
18 #include "modules/audio_coding/codecs/isac/fix/source/settings.h"
19 
20 /* Filter ar_g_Q0[] and ar_f_Q0[] through an AR filter with coefficients
21  * cth_Q15[] and sth_Q15[].
22  */
WebRtcIsacfix_FilterArLoop(int16_t * ar_g_Q0,int16_t * ar_f_Q0,int16_t * cth_Q15,int16_t * sth_Q15,size_t order_coef)23 void WebRtcIsacfix_FilterArLoop(int16_t* ar_g_Q0,     // Input samples
24                                 int16_t* ar_f_Q0,     // Input samples
25                                 int16_t* cth_Q15,     // Filter coefficients
26                                 int16_t* sth_Q15,     // Filter coefficients
27                                 size_t order_coef) { // order of the filter
28   int n = 0;
29 
30   for (n = 0; n < HALF_SUBFRAMELEN - 1; n++) {
31     size_t k = 0;
32     int16_t tmpAR = 0;
33     int32_t tmp32 = 0;
34     int32_t tmp32_2 = 0;
35 
36     tmpAR = ar_f_Q0[n + 1];
37     for (k = order_coef; k > 0; k--) {
38       tmp32 = (cth_Q15[k - 1] * tmpAR - sth_Q15[k - 1] * ar_g_Q0[k - 1] +
39                16384) >> 15;
40       tmp32_2 = (sth_Q15[k - 1] * tmpAR + cth_Q15[k - 1] * ar_g_Q0[k - 1] +
41                  16384) >> 15;
42       tmpAR   = (int16_t)WebRtcSpl_SatW32ToW16(tmp32);
43       ar_g_Q0[k] = (int16_t)WebRtcSpl_SatW32ToW16(tmp32_2);
44     }
45     ar_f_Q0[n + 1] = tmpAR;
46     ar_g_Q0[0] = tmpAR;
47   }
48 }
49