src/third_party/webrtc/common_audio/resampler/sinc_resampler_sse.cc

   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 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   __m128 m_sums1 = _mm_setzero_ps();
  25   __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   if (reinterpret_cast<uintptr_t>(input_ptr) & 0x0F) {
  30     for (int i = 0; i < kKernelSize; i += 4) {
  31       m_input = _mm_loadu_ps(input_ptr + i);
  32       m_sums1 = _mm_add_ps(m_sums1, _mm_mul_ps(m_input, _mm_load_ps(k1 + i)));
  33       m_sums2 = _mm_add_ps(m_sums2, _mm_mul_ps(m_input, _mm_load_ps(k2 + i)));
  34     }
  35   } else {
  36     for (int i = 0; i < kKernelSize; i += 4) {
  37       m_input = _mm_load_ps(input_ptr + i);
  38       m_sums1 = _mm_add_ps(m_sums1, _mm_mul_ps(m_input, _mm_load_ps(k1 + i)));
  39       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   m_sums1 = _mm_mul_ps(m_sums1, _mm_set_ps1(
  45       static_cast<float>(1.0 - kernel_interpolation_factor)));
  46   m_sums2 = _mm_mul_ps(m_sums2, _mm_set_ps1(
  47       static_cast<float>(kernel_interpolation_factor)));
  48   m_sums1 = _mm_add_ps(m_sums1, m_sums2);
  49
  50   // Sum components together.
  51   float result;
  52   m_sums2 = _mm_add_ps(_mm_movehl_ps(m_sums1, m_sums1), m_sums1);
  53   _mm_store_ss(&result, _mm_add_ss(m_sums2, _mm_shuffle_ps(
  54       m_sums2, m_sums2, 1)));
  55
  56   return result;
  57 }
  58
  59 }  // namespace webrtc