#include #include #include #include #include //test for sse #ifndef LV_HAVE_SSE void qa_32f_interleave_32fc_aligned16::t1() { printf("sse not available... no test performed\n"); } #else void qa_32f_interleave_32fc_aligned16::t1() { volk_environment_init(); clock_t start, end; double total; const int vlen = 3201; const int ITERS = 100000; float input0[vlen] __attribute__ ((aligned (16))); float input1[vlen] __attribute__ ((aligned (16))); std::complex output_generic[vlen] __attribute__ ((aligned (16))); std::complex output_sse[vlen] __attribute__ ((aligned (16))); for(int i = 0; i < vlen; ++i) { input0[i] = ((float) (rand() - (RAND_MAX/2))) / static_cast((RAND_MAX/2)); input1[i] = ((float) (rand() - (RAND_MAX/2))) / static_cast((RAND_MAX/2)); } printf("32f_interleave_32fc_aligned\n"); start = clock(); for(int count = 0; count < ITERS; ++count) { volk_32f_interleave_32fc_aligned16_manual(output_generic, input0, input1, 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_32f_interleave_32fc_aligned16_manual(output_sse, input0, input1, 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(std::real(output_generic[i]), std::real(output_sse[i]), fabs(std::real(output_generic[i]))*1e-4); CPPUNIT_ASSERT_DOUBLES_EQUAL(std::imag(output_generic[i]), std::imag(output_sse[i]), fabs(std::imag(output_generic[i]))*1e-4); } } #endif