| 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 |