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 "pitch.h"
|
cannam@154
|
34
|
cannam@154
|
35 /* Copy and multiply a vector by a constant */
|
cannam@154
|
36 void silk_scale_copy_vector16(
|
cannam@154
|
37 opus_int16 *data_out,
|
cannam@154
|
38 const opus_int16 *data_in,
|
cannam@154
|
39 opus_int32 gain_Q16, /* I Gain in Q16 */
|
cannam@154
|
40 const opus_int dataSize /* I Length */
|
cannam@154
|
41 )
|
cannam@154
|
42 {
|
cannam@154
|
43 opus_int i;
|
cannam@154
|
44 opus_int32 tmp32;
|
cannam@154
|
45
|
cannam@154
|
46 for( i = 0; i < dataSize; i++ ) {
|
cannam@154
|
47 tmp32 = silk_SMULWB( gain_Q16, data_in[ i ] );
|
cannam@154
|
48 data_out[ i ] = (opus_int16)silk_CHECK_FIT16( tmp32 );
|
cannam@154
|
49 }
|
cannam@154
|
50 }
|
cannam@154
|
51
|
cannam@154
|
52 /* Multiply a vector by a constant */
|
cannam@154
|
53 void silk_scale_vector32_Q26_lshift_18(
|
cannam@154
|
54 opus_int32 *data1, /* I/O Q0/Q18 */
|
cannam@154
|
55 opus_int32 gain_Q26, /* I Q26 */
|
cannam@154
|
56 opus_int dataSize /* I length */
|
cannam@154
|
57 )
|
cannam@154
|
58 {
|
cannam@154
|
59 opus_int i;
|
cannam@154
|
60
|
cannam@154
|
61 for( i = 0; i < dataSize; i++ ) {
|
cannam@154
|
62 data1[ i ] = (opus_int32)silk_CHECK_FIT32( silk_RSHIFT64( silk_SMULL( data1[ i ], gain_Q26 ), 8 ) ); /* OUTPUT: Q18 */
|
cannam@154
|
63 }
|
cannam@154
|
64 }
|
cannam@154
|
65
|
cannam@154
|
66 /* sum = for(i=0;i<len;i++)inVec1[i]*inVec2[i]; --- inner product */
|
cannam@154
|
67 /* Note for ARM asm: */
|
cannam@154
|
68 /* * inVec1 and inVec2 should be at least 2 byte aligned. */
|
cannam@154
|
69 /* * len should be positive 16bit integer. */
|
cannam@154
|
70 /* * only when len>6, memory access can be reduced by half. */
|
cannam@154
|
71 opus_int32 silk_inner_prod_aligned(
|
cannam@154
|
72 const opus_int16 *const inVec1, /* I input vector 1 */
|
cannam@154
|
73 const opus_int16 *const inVec2, /* I input vector 2 */
|
cannam@154
|
74 const opus_int len, /* I vector lengths */
|
cannam@154
|
75 int arch /* I Run-time architecture */
|
cannam@154
|
76 )
|
cannam@154
|
77 {
|
cannam@154
|
78 #ifdef FIXED_POINT
|
cannam@154
|
79 return celt_inner_prod(inVec1, inVec2, len, arch);
|
cannam@154
|
80 #else
|
cannam@154
|
81 opus_int i;
|
cannam@154
|
82 opus_int32 sum = 0;
|
cannam@154
|
83 for( i = 0; i < len; i++ ) {
|
cannam@154
|
84 sum = silk_SMLABB( sum, inVec1[ i ], inVec2[ i ] );
|
cannam@154
|
85 }
|
cannam@154
|
86 return sum;
|
cannam@154
|
87 #endif
|
cannam@154
|
88 }
|
cannam@154
|
89
|
cannam@154
|
90 opus_int64 silk_inner_prod16_aligned_64_c(
|
cannam@154
|
91 const opus_int16 *inVec1, /* I input vector 1 */
|
cannam@154
|
92 const opus_int16 *inVec2, /* I input vector 2 */
|
cannam@154
|
93 const opus_int len /* I vector lengths */
|
cannam@154
|
94 )
|
cannam@154
|
95 {
|
cannam@154
|
96 opus_int i;
|
cannam@154
|
97 opus_int64 sum = 0;
|
cannam@154
|
98 for( i = 0; i < len; i++ ) {
|
cannam@154
|
99 sum = silk_SMLALBB( sum, inVec1[ i ], inVec2[ i ] );
|
cannam@154
|
100 }
|
cannam@154
|
101 return sum;
|
cannam@154
|
102 }
|