OLD | NEW |
1 /*********************************************************************** | 1 /*********************************************************************** |
2 Copyright (c) 2006-2011, Skype Limited. All rights reserved. | 2 Copyright (c) 2017 Google Inc. |
3 Redistribution and use in source and binary forms, with or without | 3 Redistribution and use in source and binary forms, with or without |
4 modification, are permitted provided that the following conditions | 4 modification, are permitted provided that the following conditions |
5 are met: | 5 are met: |
6 - Redistributions of source code must retain the above copyright notice, | 6 - Redistributions of source code must retain the above copyright notice, |
7 this list of conditions and the following disclaimer. | 7 this list of conditions and the following disclaimer. |
8 - Redistributions in binary form must reproduce the above copyright | 8 - Redistributions in binary form must reproduce the above copyright |
9 notice, this list of conditions and the following disclaimer in the | 9 notice, this list of conditions and the following disclaimer in the |
10 documentation and/or other materials provided with the distribution. | 10 documentation and/or other materials provided with the distribution. |
11 - Neither the name of Internet Society, IETF or IETF Trust, nor the | 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the |
12 names of specific contributors, may be used to endorse or promote | 12 names of specific contributors, may be used to endorse or promote |
13 products derived from this software without specific prior written | 13 products derived from this software without specific prior written |
14 permission. | 14 permission. |
15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" | 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" |
16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE | 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE |
17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE | 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE |
18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE | 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE |
19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR | 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR |
20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF | 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF |
21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS | 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS |
22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN | 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN |
23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) | 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) |
24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE | 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE |
25 POSSIBILITY OF SUCH DAMAGE. | 25 POSSIBILITY OF SUCH DAMAGE. |
26 ***********************************************************************/ | 26 ***********************************************************************/ |
27 | 27 |
28 #ifdef HAVE_CONFIG_H | 28 #ifndef SILK_WARPED_AUTOCORRELATION_FIX_ARM_H |
29 #include "config.h" | 29 # define SILK_WARPED_AUTOCORRELATION_FIX_ARM_H |
30 #endif | |
31 | 30 |
32 #include "main_FIX.h" | 31 # include "celt/arm/armcpu.h" |
33 | 32 |
34 #define QC 10 | 33 # if defined(FIXED_POINT) |
35 #define QS 14 | |
36 | 34 |
37 #if defined(MIPSr1_ASM) | 35 # if defined(OPUS_ARM_MAY_HAVE_NEON_INTR) |
38 #include "mips/warped_autocorrelation_FIX_mipsr1.h" | 36 void silk_warped_autocorrelation_FIX_neon( |
39 #endif | |
40 | |
41 | |
42 #ifndef OVERRIDE_silk_warped_autocorrelation_FIX | |
43 /* Autocorrelations for a warped frequency axis */ | |
44 void silk_warped_autocorrelation_FIX( | |
45 opus_int32 *corr, /* O
Result [order + 1]
*/ | 37 opus_int32 *corr, /* O
Result [order + 1]
*/ |
46 opus_int *scale, /* O
Scaling of the correlation vector
*/ | 38 opus_int *scale, /* O
Scaling of the correlation vector
*/ |
47 const opus_int16 *input, /* I
Input data to correlate
*/ | 39 const opus_int16 *input, /* I
Input data to correlate
*/ |
48 const opus_int warping_Q16, /* I
Warping coefficient
*/ | 40 const opus_int warping_Q16, /* I
Warping coefficient
*/ |
49 const opus_int length, /* I
Length of input
*/ | 41 const opus_int length, /* I
Length of input
*/ |
50 const opus_int order /* I
Correlation order (even)
*/ | 42 const opus_int order /* I
Correlation order (even)
*/ |
51 ) | 43 ); |
52 { | |
53 opus_int n, i, lsh; | |
54 opus_int32 tmp1_QS, tmp2_QS; | |
55 opus_int32 state_QS[ MAX_SHAPE_LPC_ORDER + 1 ] = { 0 }; | |
56 opus_int64 corr_QC[ MAX_SHAPE_LPC_ORDER + 1 ] = { 0 }; | |
57 | 44 |
58 /* Order must be even */ | 45 # if !defined(OPUS_HAVE_RTCD) && defined(OPUS_ARM_PRESUME_NEON) |
59 silk_assert( ( order & 1 ) == 0 ); | 46 # define OVERRIDE_silk_warped_autocorrelation_FIX (1) |
60 silk_assert( 2 * QS - QC >= 0 ); | 47 # define silk_warped_autocorrelation_FIX(corr, scale, input, warping_Q16, leng
th, order, arch) \ |
| 48 ((void)(arch), PRESUME_NEON(silk_warped_autocorrelation_FIX)(corr, scale, in
put, warping_Q16, length, order)) |
| 49 # endif |
| 50 # endif |
61 | 51 |
62 /* Loop over samples */ | 52 # if !defined(OVERRIDE_silk_warped_autocorrelation_FIX) |
63 for( n = 0; n < length; n++ ) { | 53 /*Is run-time CPU detection enabled on this platform?*/ |
64 tmp1_QS = silk_LSHIFT32( (opus_int32)input[ n ], QS ); | 54 # if defined(OPUS_HAVE_RTCD) && (defined(OPUS_ARM_MAY_HAVE_NEON_INTR) && !defi
ned(OPUS_ARM_PRESUME_NEON_INTR)) |
65 /* Loop over allpass sections */ | 55 extern void (*const SILK_WARPED_AUTOCORRELATION_FIX_IMPL[OPUS_ARCHMASK+1])(opus_
int32*, opus_int*, const opus_int16*, const opus_int, const opus_int, const opus
_int); |
66 for( i = 0; i < order; i += 2 ) { | 56 # define OVERRIDE_silk_warped_autocorrelation_FIX (1) |
67 /* Output of allpass section */ | 57 # define silk_warped_autocorrelation_FIX(corr, scale, input, warping_Q16, len
gth, order, arch) \ |
68 tmp2_QS = silk_SMLAWB( state_QS[ i ], state_QS[ i + 1 ] - tmp1_QS, w
arping_Q16 ); | 58 ((*SILK_WARPED_AUTOCORRELATION_FIX_IMPL[(arch)&OPUS_ARCHMASK])(corr, scale,
input, warping_Q16, length, order)) |
69 state_QS[ i ] = tmp1_QS; | 59 # elif defined(OPUS_ARM_PRESUME_NEON_INTR) |
70 corr_QC[ i ] += silk_RSHIFT64( silk_SMULL( tmp1_QS, state_QS[ 0 ] )
, 2 * QS - QC ); | 60 # define OVERRIDE_silk_warped_autocorrelation_FIX (1) |
71 /* Output of allpass section */ | 61 # define silk_warped_autocorrelation_FIX(corr, scale, input, warping_Q16, len
gth, order, arch) \ |
72 tmp1_QS = silk_SMLAWB( state_QS[ i + 1 ], state_QS[ i + 2 ] - tmp2_Q
S, warping_Q16 ); | 62 ((void)(arch), silk_warped_autocorrelation_FIX_neon(corr, scale, input, warp
ing_Q16, length, order)) |
73 state_QS[ i + 1 ] = tmp2_QS; | 63 # endif |
74 corr_QC[ i + 1 ] += silk_RSHIFT64( silk_SMULL( tmp2_QS, state_QS[ 0
] ), 2 * QS - QC ); | 64 # endif |
75 } | |
76 state_QS[ order ] = tmp1_QS; | |
77 corr_QC[ order ] += silk_RSHIFT64( silk_SMULL( tmp1_QS, state_QS[ 0 ] )
, 2 * QS - QC ); | |
78 } | |
79 | 65 |
80 lsh = silk_CLZ64( corr_QC[ 0 ] ) - 35; | 66 # endif /* end FIXED_POINT */ |
81 lsh = silk_LIMIT( lsh, -12 - QC, 30 - QC ); | 67 |
82 *scale = -( QC + lsh ); | 68 #endif /* end SILK_WARPED_AUTOCORRELATION_FIX_ARM_H */ |
83 silk_assert( *scale >= -30 && *scale <= 12 ); | |
84 if( lsh >= 0 ) { | |
85 for( i = 0; i < order + 1; i++ ) { | |
86 corr[ i ] = (opus_int32)silk_CHECK_FIT32( silk_LSHIFT64( corr_QC[ i
], lsh ) ); | |
87 } | |
88 } else { | |
89 for( i = 0; i < order + 1; i++ ) { | |
90 corr[ i ] = (opus_int32)silk_CHECK_FIT32( silk_RSHIFT64( corr_QC[ i
], -lsh ) ); | |
91 } | |
92 } | |
93 silk_assert( corr_QC[ 0 ] >= 0 ); /* If breaking, decrease QC*/ | |
94 } | |
95 #endif /* OVERRIDE_silk_warped_autocorrelation_FIX */ | |
OLD | NEW |