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