1 /*
2 * Copyright (c) 2011 The WebRTC project authors. All Rights Reserved.
3 *
4 * Use of this source code is governed by a BSD-style license
5 * that can be found in the LICENSE file in the root of the source
6 * tree. An additional intellectual property rights grant can be found
7 * in the file PATENTS. All contributing project authors may
8 * be found in the AUTHORS file in the root of the source tree.
9 */
10
11
12 /*
13 * This file contains the function WebRtcSpl_DotProductWithScale().
14 * The description header can be found in signal_processing_library.h
15 *
16 */
17
18 #include "signal_processing_library.h"
19
WebRtcSpl_DotProductWithScale(WebRtc_Word16 * vector1,WebRtc_Word16 * vector2,int length,int scaling)20 WebRtc_Word32 WebRtcSpl_DotProductWithScale(WebRtc_Word16 *vector1, WebRtc_Word16 *vector2,
21 int length, int scaling)
22 {
23 WebRtc_Word32 sum;
24 int i;
25 #ifdef _ARM_OPT_
26 #pragma message("NOTE: _ARM_OPT_ optimizations are used")
27 WebRtc_Word16 len4 = (length >> 2) << 2;
28 #endif
29
30 sum = 0;
31
32 #ifndef _ARM_OPT_
33 for (i = 0; i < length; i++)
34 {
35 sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1++, *vector2++, scaling);
36 }
37 #else
38 if (scaling == 0)
39 {
40 for (i = 0; i < len4; i = i + 4)
41 {
42 sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
43 vector1++;
44 vector2++;
45 sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
46 vector1++;
47 vector2++;
48 sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
49 vector1++;
50 vector2++;
51 sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
52 vector1++;
53 vector2++;
54 }
55
56 for (i = len4; i < length; i++)
57 {
58 sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
59 vector1++;
60 vector2++;
61 }
62 }
63 else
64 {
65 for (i = 0; i < len4; i = i + 4)
66 {
67 sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
68 vector1++;
69 vector2++;
70 sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
71 vector1++;
72 vector2++;
73 sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
74 vector1++;
75 vector2++;
76 sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
77 vector1++;
78 vector2++;
79 }
80
81 for (i = len4; i < length; i++)
82 {
83 sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
84 vector1++;
85 vector2++;
86 }
87 }
88 #endif
89
90 return sum;
91 }
92