Initial checkin for VOLK - Vector-Optimized Library of Kernels. This is a new SIMD library.

It currently stands by itself under the GNU Radio tree and can be used separately. We will integrate the build process into GNU Raio and start building off of its functionality over time.

1 files changed, 138 insertions, 0 deletions

diff --git a/volk/lib/qa_8sc_deinterleave_real_32f_aligned16.cc b/volk/lib/qa_8sc_deinterleave_real_32f_aligned16.cc
new file mode 100644
index 000000000..10e537cde
--- /dev/null
+++ b/volk/lib/qa_8sc_deinterleave_real_32f_aligned16.cc
@@ -0,0 +1,138 @@
+#include <volk/volk_runtime.h>
+#include <volk/volk.h>
+#include <qa_8sc_deinterleave_real_32f_aligned16.h>
+#include <volk/volk_8sc_deinterleave_real_32f_aligned16.h>
+#include <cstdlib>
+
+//test for sse
+
+#ifndef LV_HAVE_SSE4_1
+
+#ifndef LV_HAVE_SSE
+
+void qa_8sc_deinterleave_real_32f_aligned16::t1() {
+  printf("sse not available... no test performed\n");
+}
+
+#else
+
+void qa_8sc_deinterleave_real_32f_aligned16::t1() {
+  
+  volk_environment_init();
+  int ret;
+  clock_t start, end;
+  double total;
+  const int vlen = 3201;
+  const int ITERS = 100000;
+  std::complex<int8_t> input0[vlen] __attribute__ ((aligned (16)));
+  
+  float output_generic[vlen] __attribute__ ((aligned (16)));
+  float output_sse[vlen] __attribute__ ((aligned (16)));
+
+  int8_t* loadInput = (int8_t*)input0;
+  for(int i = 0; i < vlen*2; ++i) {   
+    loadInput[i] =((char)((((float) (rand() - (RAND_MAX/2))) / static_cast<float>((RAND_MAX/2))) * 128.0));
+  }
+  printf("8sc_deinterleave_real_32f_aligned\n");
+
+  start = clock();
+  for(int count = 0; count < ITERS; ++count) {
+    volk_8sc_deinterleave_real_32f_aligned16_manual(output_generic, input0, 32768.0, vlen, "generic");
+  }
+  end = clock();
+  total = (double)(end-start)/(double)CLOCKS_PER_SEC;
+  printf("generic_time: %f\n", total);
+  start = clock();
+  for(int count = 0; count < ITERS; ++count) {
+    volk_8sc_deinterleave_real_32f_aligned16_manual(output_sse, input0, 32768.0, vlen, "sse");
+  }
+  end = clock();
+  total = (double)(end-start)/(double)CLOCKS_PER_SEC;
+  printf("sse_time: %f\n", total);
+
+  for(int i = 0; i < 1; ++i) {
+    //printf("inputs: %d, %d\n", input0[i*2], input0[i*2 + 1]);
+    //printf("generic... %d, ssse3... %d\n", output0[i], output1[i]);
+  }
+  
+  for(int i = 0; i < vlen; ++i) {
+    //printf("%d...%d\n", output0[i], output01[i]);
+    CPPUNIT_ASSERT_DOUBLES_EQUAL(output_generic[i], output_sse[i], fabs(output_generic[i])*1e-4);
+  }
+}
+
+#endif /* LV_HAVE_SSE */
+
+#else
+
+void qa_8sc_deinterleave_real_32f_aligned16::t1() {
+  
+  
+  volk_runtime_init();
+
+  volk_environment_init();
+  int ret;
+  clock_t start, end;
+  double total;
+  const int vlen = 3201;
+  const int ITERS = 100000;
+  std::complex<int8_t> *input0;
+  
+  float* output_generic;
+  float* output_sse;
+  float* output_sse4_1;
+
+  ret = posix_memalign((void**)&input0, 16, 2*vlen * sizeof(int8_t));
+  ret = posix_memalign((void**)&output_generic, 16, vlen * sizeof(float));
+  ret = posix_memalign((void**)&output_sse, 16, vlen * sizeof(float));
+  ret = posix_memalign((void**)&output_sse4_1, 16, vlen * sizeof(float));
+
+  int8_t* loadInput = (int8_t*)input0;
+  for(int i = 0; i < vlen*2; ++i) {   
+    loadInput[i] =((char)(((float) (rand() - (RAND_MAX/2))) / static_cast<float>((RAND_MAX/2))) * 128.0);
+  }
+
+  printf("8sc_deinterleave_real_32f_aligned\n");
+
+  start = clock();
+  for(int count = 0; count < ITERS; ++count) {
+    volk_8sc_deinterleave_real_32f_aligned16_manual(output_generic, input0, 128.0, vlen, "generic");
+  }
+  end = clock();
+  total = (double)(end-start)/(double)CLOCKS_PER_SEC;
+  printf("generic_time: %f\n", total);
+  start = clock();
+  for(int count = 0; count < ITERS; ++count) {
+    volk_8sc_deinterleave_real_32f_aligned16_manual(output_sse, input0, 1288.0, vlen, "sse");
+  }
+  end = clock();
+  total = (double)(end-start)/(double)CLOCKS_PER_SEC;
+  printf("sse_time: %f\n", total);
+
+  start = clock();
+  for(int count = 0; count < ITERS; ++count) {
+    get_volk_runtime()->volk_8sc_deinterleave_real_32f_aligned16(output_sse4_1, input0, 128.0, vlen);
+  }
+  end = clock();
+  total = (double)(end-start)/(double)CLOCKS_PER_SEC;
+  printf("sse4_1_time: %f\n", total);
+
+
+  for(int i = 0; i < 1; ++i) {
+    //printf("inputs: %d, %d\n", input0[i*2], input0[i*2 + 1]);
+    //printf("generic... %d, ssse3... %d\n", output0[i], output1[i]);
+  }
+  
+  for(int i = 0; i < vlen; ++i) {
+    //printf("%d...%d\n", output0[i], output01[i]);
+    CPPUNIT_ASSERT_DOUBLES_EQUAL(output_generic[i], output_sse[i], fabs(output_generic[i])*1e-4);
+    CPPUNIT_ASSERT_DOUBLES_EQUAL(output_generic[i], output_sse4_1[i], fabs(output_generic[i])*1e-4);
+  }
+
+  free(input0);
+  free(output_generic);
+  free(output_sse);
+  free(output_sse4_1);
+}
+
+#endif /* LV_HAVE_SSE4_1 */