annotate src/opus-1.3/silk/resampler_down2_3.c @ 169:223a55898ab9 tip default

Add null config files
author Chris Cannam <cannam@all-day-breakfast.com>
date Mon, 02 Mar 2020 14:03:47 +0000
parents 4664ac0c1032
children
rev   line source
cannam@154 1 /***********************************************************************
cannam@154 2 Copyright (c) 2006-2011, Skype Limited. All rights reserved.
cannam@154 3 Redistribution and use in source and binary forms, with or without
cannam@154 4 modification, are permitted provided that the following conditions
cannam@154 5 are met:
cannam@154 6 - Redistributions of source code must retain the above copyright notice,
cannam@154 7 this list of conditions and the following disclaimer.
cannam@154 8 - Redistributions in binary form must reproduce the above copyright
cannam@154 9 notice, this list of conditions and the following disclaimer in the
cannam@154 10 documentation and/or other materials provided with the distribution.
cannam@154 11 - Neither the name of Internet Society, IETF or IETF Trust, nor the
cannam@154 12 names of specific contributors, may be used to endorse or promote
cannam@154 13 products derived from this software without specific prior written
cannam@154 14 permission.
cannam@154 15 THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
cannam@154 16 AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
cannam@154 17 IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE
cannam@154 18 ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE
cannam@154 19 LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR
cannam@154 20 CONSEQUENTIAL DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF
cannam@154 21 SUBSTITUTE GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
cannam@154 22 INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN
cannam@154 23 CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE)
cannam@154 24 ARISING IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
cannam@154 25 POSSIBILITY OF SUCH DAMAGE.
cannam@154 26 ***********************************************************************/
cannam@154 27
cannam@154 28 #ifdef HAVE_CONFIG_H
cannam@154 29 #include "config.h"
cannam@154 30 #endif
cannam@154 31
cannam@154 32 #include "SigProc_FIX.h"
cannam@154 33 #include "resampler_private.h"
cannam@154 34 #include "stack_alloc.h"
cannam@154 35
cannam@154 36 #define ORDER_FIR 4
cannam@154 37
cannam@154 38 /* Downsample by a factor 2/3, low quality */
cannam@154 39 void silk_resampler_down2_3(
cannam@154 40 opus_int32 *S, /* I/O State vector [ 6 ] */
cannam@154 41 opus_int16 *out, /* O Output signal [ floor(2*inLen/3) ] */
cannam@154 42 const opus_int16 *in, /* I Input signal [ inLen ] */
cannam@154 43 opus_int32 inLen /* I Number of input samples */
cannam@154 44 )
cannam@154 45 {
cannam@154 46 opus_int32 nSamplesIn, counter, res_Q6;
cannam@154 47 VARDECL( opus_int32, buf );
cannam@154 48 opus_int32 *buf_ptr;
cannam@154 49 SAVE_STACK;
cannam@154 50
cannam@154 51 ALLOC( buf, RESAMPLER_MAX_BATCH_SIZE_IN + ORDER_FIR, opus_int32 );
cannam@154 52
cannam@154 53 /* Copy buffered samples to start of buffer */
cannam@154 54 silk_memcpy( buf, S, ORDER_FIR * sizeof( opus_int32 ) );
cannam@154 55
cannam@154 56 /* Iterate over blocks of frameSizeIn input samples */
cannam@154 57 while( 1 ) {
cannam@154 58 nSamplesIn = silk_min( inLen, RESAMPLER_MAX_BATCH_SIZE_IN );
cannam@154 59
cannam@154 60 /* Second-order AR filter (output in Q8) */
cannam@154 61 silk_resampler_private_AR2( &S[ ORDER_FIR ], &buf[ ORDER_FIR ], in,
cannam@154 62 silk_Resampler_2_3_COEFS_LQ, nSamplesIn );
cannam@154 63
cannam@154 64 /* Interpolate filtered signal */
cannam@154 65 buf_ptr = buf;
cannam@154 66 counter = nSamplesIn;
cannam@154 67 while( counter > 2 ) {
cannam@154 68 /* Inner product */
cannam@154 69 res_Q6 = silk_SMULWB( buf_ptr[ 0 ], silk_Resampler_2_3_COEFS_LQ[ 2 ] );
cannam@154 70 res_Q6 = silk_SMLAWB( res_Q6, buf_ptr[ 1 ], silk_Resampler_2_3_COEFS_LQ[ 3 ] );
cannam@154 71 res_Q6 = silk_SMLAWB( res_Q6, buf_ptr[ 2 ], silk_Resampler_2_3_COEFS_LQ[ 5 ] );
cannam@154 72 res_Q6 = silk_SMLAWB( res_Q6, buf_ptr[ 3 ], silk_Resampler_2_3_COEFS_LQ[ 4 ] );
cannam@154 73
cannam@154 74 /* Scale down, saturate and store in output array */
cannam@154 75 *out++ = (opus_int16)silk_SAT16( silk_RSHIFT_ROUND( res_Q6, 6 ) );
cannam@154 76
cannam@154 77 res_Q6 = silk_SMULWB( buf_ptr[ 1 ], silk_Resampler_2_3_COEFS_LQ[ 4 ] );
cannam@154 78 res_Q6 = silk_SMLAWB( res_Q6, buf_ptr[ 2 ], silk_Resampler_2_3_COEFS_LQ[ 5 ] );
cannam@154 79 res_Q6 = silk_SMLAWB( res_Q6, buf_ptr[ 3 ], silk_Resampler_2_3_COEFS_LQ[ 3 ] );
cannam@154 80 res_Q6 = silk_SMLAWB( res_Q6, buf_ptr[ 4 ], silk_Resampler_2_3_COEFS_LQ[ 2 ] );
cannam@154 81
cannam@154 82 /* Scale down, saturate and store in output array */
cannam@154 83 *out++ = (opus_int16)silk_SAT16( silk_RSHIFT_ROUND( res_Q6, 6 ) );
cannam@154 84
cannam@154 85 buf_ptr += 3;
cannam@154 86 counter -= 3;
cannam@154 87 }
cannam@154 88
cannam@154 89 in += nSamplesIn;
cannam@154 90 inLen -= nSamplesIn;
cannam@154 91
cannam@154 92 if( inLen > 0 ) {
cannam@154 93 /* More iterations to do; copy last part of filtered signal to beginning of buffer */
cannam@154 94 silk_memcpy( buf, &buf[ nSamplesIn ], ORDER_FIR * sizeof( opus_int32 ) );
cannam@154 95 } else {
cannam@154 96 break;
cannam@154 97 }
cannam@154 98 }
cannam@154 99
cannam@154 100 /* Copy last part of filtered signal to the state for the next call */
cannam@154 101 silk_memcpy( S, &buf[ nSamplesIn ], ORDER_FIR * sizeof( opus_int32 ) );
cannam@154 102 RESTORE_STACK;
cannam@154 103 }