TransformCQ31.cpp 13 KB

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