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