From 97ff51156b492a1e52420a1f12cd75eccf2235ac Mon Sep 17 00:00:00 2001
From: Tom Rondeau <trondeau@vt.edu>
Date: Thu, 9 Dec 2010 16:53:06 -0500
Subject: volk: Removing erroneous SSE function that actually usese an SSE3
 intrin (mm_hadd_ps).

---
 .../volk/volk_16sc_magnitude_32f_aligned16.h       | 66 ----------------------
 1 file changed, 66 deletions(-)

(limited to 'volk/include')

diff --git a/volk/include/volk/volk_16sc_magnitude_32f_aligned16.h b/volk/include/volk/volk_16sc_magnitude_32f_aligned16.h
index 29c58ceb88..8ba5737e84 100644
--- a/volk/include/volk/volk_16sc_magnitude_32f_aligned16.h
+++ b/volk/include/volk/volk_16sc_magnitude_32f_aligned16.h
@@ -70,72 +70,6 @@ static inline void volk_16sc_magnitude_32f_aligned16_sse3(float* magnitudeVector
 }
 #endif /* LV_HAVE_SSE3 */
 
-#if LV_HAVE_SSE
-#include <xmmintrin.h>
-/*!
-  \brief Calculates the magnitude of the complexVector and stores the results in the magnitudeVector
-  \param complexVector The vector containing the complex input values
-  \param magnitudeVector The vector containing the real output values
-  \param scalar The data value to be divided against each input data value of the input complex vector
-  \param num_points The number of complex values in complexVector to be calculated and stored into cVector
-*/
-static inline void volk_16sc_magnitude_32f_aligned16_sse(float* magnitudeVector, const lv_16sc_t* complexVector, const float scalar, unsigned int num_points){
-  unsigned int number = 0;
-  const unsigned int quarterPoints = num_points / 4;
-
-  const int16_t* complexVectorPtr = (const int16_t*)complexVector;
-  float* magnitudeVectorPtr = magnitudeVector;
-
-  const float iScalar = 1.0 / scalar;
-  __m128 invScalar = _mm_set_ps1(iScalar);
-
-  __m128 cplxValue1, cplxValue2, result;
-
-  float inputFloatBuffer[8] __attribute__((aligned(128)));
-
-  for(;number < quarterPoints; number++){
-
-    inputFloatBuffer[0] = (float)(complexVectorPtr[0]);
-    inputFloatBuffer[1] = (float)(complexVectorPtr[1]);
-    inputFloatBuffer[2] = (float)(complexVectorPtr[2]);
-    inputFloatBuffer[3] = (float)(complexVectorPtr[3]);
-      
-    inputFloatBuffer[4] = (float)(complexVectorPtr[4]);
-    inputFloatBuffer[5] = (float)(complexVectorPtr[5]);
-    inputFloatBuffer[6] = (float)(complexVectorPtr[6]);
-    inputFloatBuffer[7] = (float)(complexVectorPtr[7]);
-
-    cplxValue1 = _mm_load_ps(&inputFloatBuffer[0]);
-    cplxValue2 = _mm_load_ps(&inputFloatBuffer[4]);
-
-    complexVectorPtr += 8;
-
-    cplxValue1 = _mm_mul_ps(cplxValue1, invScalar);
-    cplxValue2 = _mm_mul_ps(cplxValue2, invScalar);
-
-    cplxValue1 = _mm_mul_ps(cplxValue1, cplxValue1); // Square the values
-    cplxValue2 = _mm_mul_ps(cplxValue2, cplxValue2); // Square the Values
-
-    result = _mm_hadd_ps(cplxValue1, cplxValue2); // Add the I2 and Q2 values
-
-    result = _mm_sqrt_ps(result); // Square root the values
-
-    _mm_store_ps(magnitudeVectorPtr, result);
-      
-    magnitudeVectorPtr += 4;
-  }
-
-  number = quarterPoints * 4;
-  magnitudeVectorPtr = &magnitudeVector[number];
-  complexVectorPtr = (const int16_t*)&complexVector[number];
-  for(; number < num_points; number++){
-    float val1Real = (float)(*complexVectorPtr++) * iScalar;
-    float val1Imag = (float)(*complexVectorPtr++) * iScalar;
-    *magnitudeVectorPtr++ = sqrtf((val1Real * val1Real) + (val1Imag * val1Imag));
-  }
-}
-#endif /* LV_HAVE_SSE */
-
 #if LV_HAVE_GENERIC
 /*!
   \brief Calculates the magnitude of the complexVector and stores the results in the magnitudeVector
-- 
cgit v1.2.3