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