trunk/src/common_audio/signal_processing/dot_product_with_scale.c

/*
 *  Copyright (c) 2011 The WebRTC project authors. All Rights Reserved.
 *
 *  Use of this source code is governed by a BSD-style license
 *  that can be found in the LICENSE file in the root of the source
 *  tree. An additional intellectual property rights grant can be found
 *  in the file PATENTS.  All contributing project authors may
 *  be found in the AUTHORS file in the root of the source tree.
 */


/*
 * This file contains the function WebRtcSpl_DotProductWithScale().
 * The description header can be found in signal_processing_library.h
 *
 */

#include "signal_processing_library.h"

WebRtc_Word32 WebRtcSpl_DotProductWithScale(WebRtc_Word16 *vector1, WebRtc_Word16 *vector2,
                                            int length, int scaling)
{
    WebRtc_Word32 sum;
    int i;
#ifdef _ARM_OPT_
#pragma message("NOTE: _ARM_OPT_ optimizations are used")
    WebRtc_Word16 len4 = (length >> 2) << 2;
#endif

    sum = 0;

#ifndef _ARM_OPT_
    for (i = 0; i < length; i++)
    {
        sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1++, *vector2++, scaling);
    }
#else
    if (scaling == 0)
    {
        for (i = 0; i < len4; i = i + 4)
        {
            sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
            vector1++;
            vector2++;
            sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
            vector1++;
            vector2++;
            sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
            vector1++;
            vector2++;
            sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
            vector1++;
            vector2++;
        }

        for (i = len4; i < length; i++)
        {
            sum += WEBRTC_SPL_MUL_16_16(*vector1, *vector2);
            vector1++;
            vector2++;
        }
    }
    else
    {
        for (i = 0; i < len4; i = i + 4)
        {
            sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
            vector1++;
            vector2++;
            sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
            vector1++;
            vector2++;
            sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
            vector1++;
            vector2++;
            sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
            vector1++;
            vector2++;
        }

        for (i = len4; i < length; i++)
        {
            sum += WEBRTC_SPL_MUL_16_16_RSFT(*vector1, *vector2, scaling);
            vector1++;
            vector2++;
        }
    }
#endif

    return sum;
}