TransformCQ31.cpp 13 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499
  1. #include "TransformCQ31.h"
  2. #include <stdio.h>
  3. #include "Error.h"
  4. #include "Test.h"
  5. #define SNR_THRESHOLD 90
  6. void TransformCQ31::test_cfft_q31()
  7. {
  8. const q31_t *inp = input.ptr();
  9. q31_t *outfftp = outputfft.ptr();
  10. memcpy(outfftp,inp,sizeof(q31_t)*input.nbSamples());
  11. ASSERT_TRUE(status == ARM_MATH_SUCCESS);
  12. arm_cfft_q31(
  13. &(this->instCfftQ31),
  14. outfftp,
  15. this->ifft,
  16. 1);
  17. ASSERT_SNR(outputfft,ref,(float32_t)SNR_THRESHOLD);
  18. ASSERT_NEAR_EQ(outputfft,ref,(q31_t)32);
  19. ASSERT_EMPTY_TAIL(outputfft);
  20. }
  21. void TransformCQ31::test_cifft_q31()
  22. {
  23. const q31_t *inp = input.ptr();
  24. q31_t *outfftp = outputfft.ptr();
  25. q31_t *refp = ref.ptr();
  26. memcpy(outfftp,inp,sizeof(q31_t)*input.nbSamples());
  27. ASSERT_TRUE(status == ARM_MATH_SUCCESS);
  28. arm_cfft_q31(
  29. &(this->instCfftQ31),
  30. outfftp,
  31. this->ifft,
  32. 1);
  33. for(unsigned long i=0; i < outputfft.nbSamples();i++)
  34. {
  35. refp[i] >>= this->scaling;
  36. }
  37. ASSERT_SNR(outputfft,ref,(float32_t)SNR_THRESHOLD);
  38. ASSERT_EMPTY_TAIL(outputfft);
  39. }
  40. void TransformCQ31::setUp(Testing::testID_t id,std::vector<Testing::param_t>& paramsArgs,Client::PatternMgr *mgr)
  41. {
  42. (void)paramsArgs;
  43. switch(id)
  44. {
  45. case TransformCQ31::TEST_CFFT_Q31_1:
  46. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_16_Q31_ID,mgr);
  47. ref.reload( TransformCQ31::REF_CFFT_NOISY_16_Q31_ID,mgr);
  48. status=arm_cfft_init_q31(&instCfftQ31,16);
  49. this->ifft=0;
  50. break;
  51. case TransformCQ31::TEST_CIFFT_Q31_19:
  52. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_16_Q31_ID,mgr);
  53. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_16_Q31_ID,mgr);
  54. status=arm_cfft_init_q31(&instCfftQ31,16);
  55. this->ifft=1;
  56. this->scaling = 4;
  57. break;
  58. case TransformCQ31::TEST_CFFT_Q31_2:
  59. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_32_Q31_ID,mgr);
  60. ref.reload( TransformCQ31::REF_CFFT_NOISY_32_Q31_ID,mgr);
  61. status=arm_cfft_init_q31(&instCfftQ31,32);
  62. this->ifft=0;
  63. break;
  64. case TransformCQ31::TEST_CIFFT_Q31_20:
  65. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_32_Q31_ID,mgr);
  66. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_32_Q31_ID,mgr);
  67. status=arm_cfft_init_q31(&instCfftQ31,32);
  68. this->ifft=1;
  69. this->scaling = 5;
  70. break;
  71. case TransformCQ31::TEST_CFFT_Q31_3:
  72. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_64_Q31_ID,mgr);
  73. ref.reload( TransformCQ31::REF_CFFT_NOISY_64_Q31_ID,mgr);
  74. status=arm_cfft_init_q31(&instCfftQ31,64);
  75. this->ifft=0;
  76. break;
  77. case TransformCQ31::TEST_CIFFT_Q31_21:
  78. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_64_Q31_ID,mgr);
  79. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_64_Q31_ID,mgr);
  80. status=arm_cfft_init_q31(&instCfftQ31,64);
  81. this->ifft=1;
  82. this->scaling=6;
  83. break;
  84. case TransformCQ31::TEST_CFFT_Q31_4:
  85. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_128_Q31_ID,mgr);
  86. ref.reload( TransformCQ31::REF_CFFT_NOISY_128_Q31_ID,mgr);
  87. status=arm_cfft_init_q31(&instCfftQ31,128);
  88. this->ifft=0;
  89. break;
  90. case TransformCQ31::TEST_CIFFT_Q31_22:
  91. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_128_Q31_ID,mgr);
  92. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_128_Q31_ID,mgr);
  93. status=arm_cfft_init_q31(&instCfftQ31,128);
  94. this->ifft=1;
  95. this->scaling=7;
  96. break;
  97. case TransformCQ31::TEST_CFFT_Q31_5:
  98. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_256_Q31_ID,mgr);
  99. ref.reload( TransformCQ31::REF_CFFT_NOISY_256_Q31_ID,mgr);
  100. status=arm_cfft_init_q31(&instCfftQ31,256);
  101. this->ifft=0;
  102. break;
  103. case TransformCQ31::TEST_CIFFT_Q31_23:
  104. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_256_Q31_ID,mgr);
  105. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_256_Q31_ID,mgr);
  106. status=arm_cfft_init_q31(&instCfftQ31,256);
  107. this->ifft=1;
  108. this->scaling=8;
  109. break;
  110. case TransformCQ31::TEST_CFFT_Q31_6:
  111. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_512_Q31_ID,mgr);
  112. ref.reload( TransformCQ31::REF_CFFT_NOISY_512_Q31_ID,mgr);
  113. status=arm_cfft_init_q31(&instCfftQ31,512);
  114. this->ifft=0;
  115. break;
  116. case TransformCQ31::TEST_CIFFT_Q31_24:
  117. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_512_Q31_ID,mgr);
  118. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_512_Q31_ID,mgr);
  119. status=arm_cfft_init_q31(&instCfftQ31,512);
  120. this->ifft=1;
  121. this->scaling=9;
  122. break;
  123. case TransformCQ31::TEST_CFFT_Q31_7:
  124. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_1024_Q31_ID,mgr);
  125. ref.reload( TransformCQ31::REF_CFFT_NOISY_1024_Q31_ID,mgr);
  126. status=arm_cfft_init_q31(&instCfftQ31,1024);
  127. this->ifft=0;
  128. break;
  129. case TransformCQ31::TEST_CIFFT_Q31_25:
  130. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_1024_Q31_ID,mgr);
  131. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_1024_Q31_ID,mgr);
  132. status=arm_cfft_init_q31(&instCfftQ31,1024);
  133. this->ifft=1;
  134. this->scaling=10;
  135. break;
  136. case TransformCQ31::TEST_CFFT_Q31_8:
  137. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_2048_Q31_ID,mgr);
  138. ref.reload( TransformCQ31::REF_CFFT_NOISY_2048_Q31_ID,mgr);
  139. status=arm_cfft_init_q31(&instCfftQ31,2048);
  140. this->ifft=0;
  141. break;
  142. case TransformCQ31::TEST_CIFFT_Q31_26:
  143. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_2048_Q31_ID,mgr);
  144. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_2048_Q31_ID,mgr);
  145. status=arm_cfft_init_q31(&instCfftQ31,2048);
  146. this->ifft=1;
  147. this->scaling=11;
  148. break;
  149. case TransformCQ31::TEST_CFFT_Q31_9:
  150. input.reload(TransformCQ31::INPUTS_CFFT_NOISY_4096_Q31_ID,mgr);
  151. ref.reload( TransformCQ31::REF_CFFT_NOISY_4096_Q31_ID,mgr);
  152. status=arm_cfft_init_q31(&instCfftQ31,4096);
  153. this->ifft=0;
  154. break;
  155. case TransformCQ31::TEST_CIFFT_Q31_27:
  156. input.reload(TransformCQ31::INPUTS_CIFFT_NOISY_4096_Q31_ID,mgr);
  157. ref.reload( TransformCQ31::INPUTS_CFFT_NOISY_4096_Q31_ID,mgr);
  158. status=arm_cfft_init_q31(&instCfftQ31,4096);
  159. this->ifft=1;
  160. this->scaling=12;
  161. break;
  162. /* STEP FUNCTIONS */
  163. case TransformCQ31::TEST_CFFT_Q31_10:
  164. input.reload(TransformCQ31::INPUTS_CFFT_STEP_16_Q31_ID,mgr);
  165. ref.reload( TransformCQ31::REF_CFFT_STEP_16_Q31_ID,mgr);
  166. status=arm_cfft_init_q31(&instCfftQ31,16);
  167. this->ifft=0;
  168. break;
  169. case TransformCQ31::TEST_CIFFT_Q31_28:
  170. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_16_Q31_ID,mgr);
  171. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_16_Q31_ID,mgr);
  172. status=arm_cfft_init_q31(&instCfftQ31,16);
  173. this->ifft=1;
  174. this->scaling=4;
  175. break;
  176. case TransformCQ31::TEST_CFFT_Q31_11:
  177. input.reload(TransformCQ31::INPUTS_CFFT_STEP_32_Q31_ID,mgr);
  178. ref.reload( TransformCQ31::REF_CFFT_STEP_32_Q31_ID,mgr);
  179. status=arm_cfft_init_q31(&instCfftQ31,32);
  180. this->ifft=0;
  181. break;
  182. case TransformCQ31::TEST_CIFFT_Q31_29:
  183. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_32_Q31_ID,mgr);
  184. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_32_Q31_ID,mgr);
  185. status=arm_cfft_init_q31(&instCfftQ31,32);
  186. this->ifft=1;
  187. this->scaling=5;
  188. break;
  189. case TransformCQ31::TEST_CFFT_Q31_12:
  190. input.reload(TransformCQ31::INPUTS_CFFT_STEP_64_Q31_ID,mgr);
  191. ref.reload( TransformCQ31::REF_CFFT_STEP_64_Q31_ID,mgr);
  192. status=arm_cfft_init_q31(&instCfftQ31,64);
  193. this->ifft=0;
  194. break;
  195. case TransformCQ31::TEST_CIFFT_Q31_30:
  196. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_64_Q31_ID,mgr);
  197. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_64_Q31_ID,mgr);
  198. status=arm_cfft_init_q31(&instCfftQ31,64);
  199. this->ifft=1;
  200. this->scaling=6;
  201. break;
  202. case TransformCQ31::TEST_CFFT_Q31_13:
  203. input.reload(TransformCQ31::INPUTS_CFFT_STEP_128_Q31_ID,mgr);
  204. ref.reload( TransformCQ31::REF_CFFT_STEP_128_Q31_ID,mgr);
  205. status=arm_cfft_init_q31(&instCfftQ31,128);
  206. this->ifft=0;
  207. break;
  208. case TransformCQ31::TEST_CIFFT_Q31_31:
  209. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_128_Q31_ID,mgr);
  210. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_128_Q31_ID,mgr);
  211. status=arm_cfft_init_q31(&instCfftQ31,128);
  212. this->ifft=1;
  213. this->scaling=7;
  214. break;
  215. case TransformCQ31::TEST_CFFT_Q31_14:
  216. input.reload(TransformCQ31::INPUTS_CFFT_STEP_256_Q31_ID,mgr);
  217. ref.reload( TransformCQ31::REF_CFFT_STEP_256_Q31_ID,mgr);
  218. status=arm_cfft_init_q31(&instCfftQ31,256);
  219. this->ifft=0;
  220. break;
  221. case TransformCQ31::TEST_CIFFT_Q31_32:
  222. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_256_Q31_ID,mgr);
  223. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_256_Q31_ID,mgr);
  224. status=arm_cfft_init_q31(&instCfftQ31,256);
  225. this->ifft=1;
  226. this->scaling=8;
  227. break;
  228. case TransformCQ31::TEST_CFFT_Q31_15:
  229. input.reload(TransformCQ31::INPUTS_CFFT_STEP_512_Q31_ID,mgr);
  230. ref.reload( TransformCQ31::REF_CFFT_STEP_512_Q31_ID,mgr);
  231. status=arm_cfft_init_q31(&instCfftQ31,512);
  232. this->ifft=0;
  233. break;
  234. case TransformCQ31::TEST_CIFFT_Q31_33:
  235. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_512_Q31_ID,mgr);
  236. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_512_Q31_ID,mgr);
  237. status=arm_cfft_init_q31(&instCfftQ31,512);
  238. this->ifft=1;
  239. this->scaling=9;
  240. break;
  241. case TransformCQ31::TEST_CFFT_Q31_16:
  242. input.reload(TransformCQ31::INPUTS_CFFT_STEP_1024_Q31_ID,mgr);
  243. ref.reload( TransformCQ31::REF_CFFT_STEP_1024_Q31_ID,mgr);
  244. status=arm_cfft_init_q31(&instCfftQ31,1024);
  245. this->ifft=0;
  246. break;
  247. case TransformCQ31::TEST_CIFFT_Q31_34:
  248. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_1024_Q31_ID,mgr);
  249. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_1024_Q31_ID,mgr);
  250. status=arm_cfft_init_q31(&instCfftQ31,1024);
  251. this->ifft=1;
  252. this->scaling=10;
  253. break;
  254. case TransformCQ31::TEST_CFFT_Q31_17:
  255. input.reload(TransformCQ31::INPUTS_CFFT_STEP_2048_Q31_ID,mgr);
  256. ref.reload( TransformCQ31::REF_CFFT_STEP_2048_Q31_ID,mgr);
  257. status=arm_cfft_init_q31(&instCfftQ31,2048);
  258. this->ifft=0;
  259. break;
  260. case TransformCQ31::TEST_CIFFT_Q31_35:
  261. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_2048_Q31_ID,mgr);
  262. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_2048_Q31_ID,mgr);
  263. status=arm_cfft_init_q31(&instCfftQ31,2048);
  264. this->ifft=1;
  265. this->scaling=11;
  266. break;
  267. case TransformCQ31::TEST_CFFT_Q31_18:
  268. input.reload(TransformCQ31::INPUTS_CFFT_STEP_4096_Q31_ID,mgr);
  269. ref.reload( TransformCQ31::REF_CFFT_STEP_4096_Q31_ID,mgr);
  270. status=arm_cfft_init_q31(&instCfftQ31,4096);
  271. this->ifft=0;
  272. break;
  273. case TransformCQ31::TEST_CIFFT_Q31_36:
  274. input.reload(TransformCQ31::INPUTS_CIFFT_STEP_4096_Q31_ID,mgr);
  275. ref.reload( TransformCQ31::INPUTS_CFFT_STEP_4096_Q31_ID,mgr);
  276. status=arm_cfft_init_q31(&instCfftQ31,4096);
  277. this->ifft=1;
  278. this->scaling=12;
  279. break;
  280. }
  281. outputfft.create(ref.nbSamples(),TransformCQ31::OUTPUT_CFFT_Q31_ID,mgr);
  282. }
  283. void TransformCQ31::tearDown(Testing::testID_t id,Client::PatternMgr *mgr)
  284. {
  285. (void)id;
  286. outputfft.dump(mgr);
  287. }