1*a58d3d2aSXin Li /***********************************************************************
2*a58d3d2aSXin Li Copyright (c) 2006-2011, Skype Limited. All rights reserved.
3*a58d3d2aSXin Li Redistribution and use in source and binary forms, with or without
4*a58d3d2aSXin Li modification, are permitted provided that the following conditions
5*a58d3d2aSXin Li are met:
6*a58d3d2aSXin Li - Redistributions of source code must retain the above copyright notice,
7*a58d3d2aSXin Li this list of conditions and the following disclaimer.
8*a58d3d2aSXin Li - Redistributions in binary form must reproduce the above copyright
9*a58d3d2aSXin Li notice, this list of conditions and the following disclaimer in the
10*a58d3d2aSXin Li documentation and/or other materials provided with the distribution.
11*a58d3d2aSXin Li - Neither the name of Internet Society, IETF or IETF Trust, nor the
12*a58d3d2aSXin Li names of specific contributors, may be used to endorse or promote
13*a58d3d2aSXin Li products derived from this software without specific prior written
14*a58d3d2aSXin Li permission.
15*a58d3d2aSXin Li THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
16*a58d3d2aSXin Li AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
17*a58d3d2aSXin Li IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
18*a58d3d2aSXin Li ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
19*a58d3d2aSXin Li LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
20*a58d3d2aSXin Li CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
21*a58d3d2aSXin Li SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
22*a58d3d2aSXin Li INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
23*a58d3d2aSXin Li CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
24*a58d3d2aSXin Li ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
25*a58d3d2aSXin Li POSSIBILITY OF SUCH DAMAGE.
26*a58d3d2aSXin Li ***********************************************************************/
27*a58d3d2aSXin Li
28*a58d3d2aSXin Li #ifdef HAVE_CONFIG_H
29*a58d3d2aSXin Li #include "config.h"
30*a58d3d2aSXin Li #endif
31*a58d3d2aSXin Li
32*a58d3d2aSXin Li #include "main_FLP.h"
33*a58d3d2aSXin Li #include "tuning_parameters.h"
34*a58d3d2aSXin Li
silk_find_LTP_FLP(silk_float XX[MAX_NB_SUBFR * LTP_ORDER * LTP_ORDER],silk_float xX[MAX_NB_SUBFR * LTP_ORDER],const silk_float r_ptr[],const opus_int lag[MAX_NB_SUBFR],const opus_int subfr_length,const opus_int nb_subfr,int arch)35*a58d3d2aSXin Li void silk_find_LTP_FLP(
36*a58d3d2aSXin Li silk_float XX[ MAX_NB_SUBFR * LTP_ORDER * LTP_ORDER ], /* O Weight for LTP quantization */
37*a58d3d2aSXin Li silk_float xX[ MAX_NB_SUBFR * LTP_ORDER ], /* O Weight for LTP quantization */
38*a58d3d2aSXin Li const silk_float r_ptr[], /* I LPC residual */
39*a58d3d2aSXin Li const opus_int lag[ MAX_NB_SUBFR ], /* I LTP lags */
40*a58d3d2aSXin Li const opus_int subfr_length, /* I Subframe length */
41*a58d3d2aSXin Li const opus_int nb_subfr, /* I number of subframes */
42*a58d3d2aSXin Li int arch
43*a58d3d2aSXin Li )
44*a58d3d2aSXin Li {
45*a58d3d2aSXin Li opus_int k;
46*a58d3d2aSXin Li silk_float *xX_ptr, *XX_ptr;
47*a58d3d2aSXin Li const silk_float *lag_ptr;
48*a58d3d2aSXin Li silk_float xx, temp;
49*a58d3d2aSXin Li
50*a58d3d2aSXin Li xX_ptr = xX;
51*a58d3d2aSXin Li XX_ptr = XX;
52*a58d3d2aSXin Li for( k = 0; k < nb_subfr; k++ ) {
53*a58d3d2aSXin Li lag_ptr = r_ptr - ( lag[ k ] + LTP_ORDER / 2 );
54*a58d3d2aSXin Li silk_corrMatrix_FLP( lag_ptr, subfr_length, LTP_ORDER, XX_ptr, arch );
55*a58d3d2aSXin Li silk_corrVector_FLP( lag_ptr, r_ptr, subfr_length, LTP_ORDER, xX_ptr, arch );
56*a58d3d2aSXin Li xx = ( silk_float )silk_energy_FLP( r_ptr, subfr_length + LTP_ORDER );
57*a58d3d2aSXin Li temp = 1.0f / silk_max( xx, LTP_CORR_INV_MAX * 0.5f * ( XX_ptr[ 0 ] + XX_ptr[ 24 ] ) + 1.0f );
58*a58d3d2aSXin Li silk_scale_vector_FLP( XX_ptr, temp, LTP_ORDER * LTP_ORDER );
59*a58d3d2aSXin Li silk_scale_vector_FLP( xX_ptr, temp, LTP_ORDER );
60*a58d3d2aSXin Li
61*a58d3d2aSXin Li r_ptr += subfr_length;
62*a58d3d2aSXin Li XX_ptr += LTP_ORDER * LTP_ORDER;
63*a58d3d2aSXin Li xX_ptr += LTP_ORDER;
64*a58d3d2aSXin Li }
65*a58d3d2aSXin Li }
66