Line data    Source code

       1             : /*
       2             :  *  Copyright (c) 2013 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             : // Modified from the Chromium original:
      12             : // src/media/base/simd/sinc_resampler_sse.cc
      13             : 
      14             : #include "webrtc/common_audio/resampler/sinc_resampler.h"
      15             : 
      16             : #include <xmmintrin.h>
      17             : 
      18             : namespace webrtc {
      19             : 
      20           0 : float SincResampler::Convolve_SSE(const float* input_ptr, const float* k1,
      21             :                                   const float* k2,
      22             :                                   double kernel_interpolation_factor) {
      23             :   __m128 m_input;
      24           0 :   __m128 m_sums1 = _mm_setzero_ps();
      25           0 :   __m128 m_sums2 = _mm_setzero_ps();
      26             : 
      27             :   // Based on |input_ptr| alignment, we need to use loadu or load.  Unrolling
      28             :   // these loops hurt performance in local testing.
      29           0 :   if (reinterpret_cast<uintptr_t>(input_ptr) & 0x0F) {
      30           0 :     for (size_t i = 0; i < kKernelSize; i += 4) {
      31           0 :       m_input = _mm_loadu_ps(input_ptr + i);
      32           0 :       m_sums1 = _mm_add_ps(m_sums1, _mm_mul_ps(m_input, _mm_load_ps(k1 + i)));
      33           0 :       m_sums2 = _mm_add_ps(m_sums2, _mm_mul_ps(m_input, _mm_load_ps(k2 + i)));
      34             :     }
      35             :   } else {
      36           0 :     for (size_t i = 0; i < kKernelSize; i += 4) {
      37           0 :       m_input = _mm_load_ps(input_ptr + i);
      38           0 :       m_sums1 = _mm_add_ps(m_sums1, _mm_mul_ps(m_input, _mm_load_ps(k1 + i)));
      39           0 :       m_sums2 = _mm_add_ps(m_sums2, _mm_mul_ps(m_input, _mm_load_ps(k2 + i)));
      40             :     }
      41             :   }
      42             : 
      43             :   // Linearly interpolate the two "convolutions".
      44           0 :   m_sums1 = _mm_mul_ps(m_sums1, _mm_set_ps1(
      45           0 :       static_cast<float>(1.0 - kernel_interpolation_factor)));
      46           0 :   m_sums2 = _mm_mul_ps(m_sums2, _mm_set_ps1(
      47           0 :       static_cast<float>(kernel_interpolation_factor)));
      48           0 :   m_sums1 = _mm_add_ps(m_sums1, m_sums2);
      49             : 
      50             :   // Sum components together.
      51             :   float result;
      52           0 :   m_sums2 = _mm_add_ps(_mm_movehl_ps(m_sums1, m_sums1), m_sums1);
      53           0 :   _mm_store_ss(&result, _mm_add_ss(m_sums2, _mm_shuffle_ps(
      54             :       m_sums2, m_sums2, 1)));
      55             : 
      56           0 :   return result;
      57             : }
      58             : 
      59             : }  // namespace webrtc