Move the aec_rdft* files to a more proper place beneath APM and make them thread-safe.

webrtc/modules/audio_processing/aec/aec_core_sse2.cc

 /*
  *  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.
  */
 
 /*
  * The core AEC algorithm, SSE2 version of speed-critical functions.
  */
 
 #include <emmintrin.h>
 #include <math.h>
 #include <string.h>  // memset
 
 extern "C" {
 #include "webrtc/common_audio/signal_processing/include/signal_processing_library.h"
 }
 #include "webrtc/modules/audio_processing/aec/aec_common.h"
 #include "webrtc/modules/audio_processing/aec/aec_core_optimized_methods.h"
-#include "webrtc/modules/audio_processing/aec/aec_rdft.h"
+#include "webrtc/modules/audio_processing/utility/ooura_fft.h"
 
 namespace webrtc {
 
 __inline static float MulRe(float aRe, float aIm, float bRe, float bIm) {
   return aRe * bRe - aIm * bIm;
 }
 
 __inline static float MulIm(float aRe, float aIm, float bRe, float bIm) {
   return aRe * bIm + aIm * bRe;
 }
@@ skipping 98 matching lines @@
       }
 
       // Stepsize factor
       ef[0][i] *= mu;
       ef[1][i] *= mu;
     }
   }
 }
 
 static void FilterAdaptationSSE2(
+    const OouraFft& ooura_fft,
     int num_partitions,
     int x_fft_buf_block_pos,
     float x_fft_buf[2][kExtendedNumPartitions * PART_LEN1],
     float e_fft[2][PART_LEN1],
     float h_fft_buf[2][kExtendedNumPartitions * PART_LEN1]) {
   float fft[PART_LEN2];
   int i, j;
   for (i = 0; i < num_partitions; i++) {
     int xPos = (i + x_fft_buf_block_pos) * (PART_LEN1);
     int pos = i * PART_LEN1;
@@ skipping 23 matching lines @@
       const __m128 h = _mm_unpackhi_ps(e, f);
       // Store
       _mm_storeu_ps(&fft[2 * j + 0], g);
       _mm_storeu_ps(&fft[2 * j + 4], h);
     }
     // ... and fixup the first imaginary entry.
     fft[1] =
         MulRe(x_fft_buf[0][xPos + PART_LEN], -x_fft_buf[1][xPos + PART_LEN],
               e_fft[0][PART_LEN], e_fft[1][PART_LEN]);
 
-    aec_rdft_inverse_128(fft);
+    ooura_fft.InverseFft(fft);
     memset(fft + PART_LEN, 0, sizeof(float) * PART_LEN);
 
     // fft scaling
     {
       float scale = 2.0f / PART_LEN2;
       const __m128 scale_ps = _mm_load_ps1(&scale);
       for (j = 0; j < PART_LEN; j += 4) {
         const __m128 fft_ps = _mm_loadu_ps(&fft[j]);
         const __m128 fft_scale = _mm_mul_ps(fft_ps, scale_ps);
         _mm_storeu_ps(&fft[j], fft_scale);
       }
     }
-    aec_rdft_forward_128(fft);
+    ooura_fft.Fft(fft);
 
     {
       float wt1 = h_fft_buf[1][pos];
       h_fft_buf[0][pos + PART_LEN] += fft[1];
       for (j = 0; j < PART_LEN; j += 4) {
         __m128 wtBuf_re = _mm_loadu_ps(&h_fft_buf[0][pos + j]);
         __m128 wtBuf_im = _mm_loadu_ps(&h_fft_buf[1][pos + j]);
         const __m128 fft0 = _mm_loadu_ps(&fft[2 * j + 0]);
         const __m128 fft4 = _mm_loadu_ps(&fft[2 * j + 4]);
         const __m128 fft_re =
@@ skipping 531 matching lines @@
   WebRtcAec_FilterAdaptation = FilterAdaptationSSE2;
   WebRtcAec_Overdrive = OverdriveSSE2;
   WebRtcAec_Suppress = SuppressSSE2;
   WebRtcAec_ComputeCoherence = ComputeCoherenceSSE2;
   WebRtcAec_UpdateCoherenceSpectra = UpdateCoherenceSpectraSSE2;
   WebRtcAec_StoreAsComplex = StoreAsComplexSSE2;
   WebRtcAec_PartitionDelay = PartitionDelaySSE2;
   WebRtcAec_WindowData = WindowDataSSE2;
 }
 }  // namespace webrtc