TransformRQ15.cpp 14 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517
  1. #include "TransformRQ15.h"
  2. #include <stdio.h>
  3. #include "Error.h"
  4. #include "Test.h"
  5. #define SNR_THRESHOLD 40
  6. #define RIFFT_SNR_THRESHOLD 25
  7. void TransformRQ15::test_rfft_q15()
  8. {
  9. q15_t *inp = input.ptr();
  10. q15_t *tmp = inputchanged.ptr();
  11. q15_t *outp = outputfft.ptr();
  12. q15_t *overoutp = overheadoutputfft.ptr();
  13. memcpy(tmp,inp,sizeof(q15_t)*input.nbSamples());
  14. arm_rfft_q15(
  15. &this->instRfftQ15,
  16. tmp,
  17. overoutp);
  18. if (this->ifft)
  19. {
  20. for(unsigned long i = 0;i < overheadoutputfft.nbSamples(); i++)
  21. {
  22. overoutp[i] = overoutp[i] << this->scaling;
  23. }
  24. }
  25. memcpy(outp,overoutp,sizeof(q15_t)*outputfft.nbSamples());
  26. if (this->ifft)
  27. {
  28. ASSERT_SNR(outputfft,ref,(q15_t)RIFFT_SNR_THRESHOLD);
  29. }
  30. else
  31. {
  32. ASSERT_SNR(outputfft,ref,(q15_t)SNR_THRESHOLD);
  33. }
  34. ASSERT_EMPTY_TAIL(outputfft);
  35. }
  36. void TransformRQ15::setUp(Testing::testID_t id,std::vector<Testing::param_t>& paramsArgs,Client::PatternMgr *mgr)
  37. {
  38. (void)paramsArgs;
  39. switch(id)
  40. {
  41. case TransformRQ15::TEST_RFFT_Q15_1:
  42. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_32_Q15_ID,mgr);
  43. ref.reload( TransformRQ15::REF_RFFT_NOISY_32_Q15_ID,mgr);
  44. arm_rfft_init_q15(&this->instRfftQ15 ,32,0,1);
  45. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  46. this->ifft=0;
  47. break;
  48. case TransformRQ15::TEST_RFFT_Q15_17:
  49. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_32_Q15_ID,mgr);
  50. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_32_Q15_ID,mgr);
  51. arm_rfft_init_q15(&this->instRfftQ15 ,32,1,1);
  52. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  53. this->ifft=1;
  54. this->scaling=5;
  55. break;
  56. case TransformRQ15::TEST_RFFT_Q15_2:
  57. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_64_Q15_ID,mgr);
  58. ref.reload( TransformRQ15::REF_RFFT_NOISY_64_Q15_ID,mgr);
  59. arm_rfft_init_q15(&this->instRfftQ15 ,64,0,1);
  60. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  61. this->ifft=0;
  62. break;
  63. case TransformRQ15::TEST_RFFT_Q15_18:
  64. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_64_Q15_ID,mgr);
  65. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_64_Q15_ID,mgr);
  66. arm_rfft_init_q15(&this->instRfftQ15 ,64,1,1);
  67. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  68. this->ifft=1;
  69. this->scaling=6;
  70. break;
  71. case TransformRQ15::TEST_RFFT_Q15_3:
  72. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_128_Q15_ID,mgr);
  73. ref.reload( TransformRQ15::REF_RFFT_NOISY_128_Q15_ID,mgr);
  74. arm_rfft_init_q15(&this->instRfftQ15 ,128,0,1);
  75. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  76. this->ifft=0;
  77. break;
  78. case TransformRQ15::TEST_RFFT_Q15_19:
  79. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_128_Q15_ID,mgr);
  80. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_128_Q15_ID,mgr);
  81. arm_rfft_init_q15(&this->instRfftQ15 ,128,1,1);
  82. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  83. this->ifft=1;
  84. this->scaling=7;
  85. break;
  86. case TransformRQ15::TEST_RFFT_Q15_4:
  87. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_256_Q15_ID,mgr);
  88. ref.reload( TransformRQ15::REF_RFFT_NOISY_256_Q15_ID,mgr);
  89. arm_rfft_init_q15(&this->instRfftQ15 ,256,0,1);
  90. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  91. this->ifft=0;
  92. break;
  93. case TransformRQ15::TEST_RFFT_Q15_20:
  94. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_256_Q15_ID,mgr);
  95. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_256_Q15_ID,mgr);
  96. arm_rfft_init_q15(&this->instRfftQ15 ,256,1,1);
  97. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  98. this->ifft=1;
  99. this->scaling=8;
  100. break;
  101. case TransformRQ15::TEST_RFFT_Q15_5:
  102. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_512_Q15_ID,mgr);
  103. ref.reload( TransformRQ15::REF_RFFT_NOISY_512_Q15_ID,mgr);
  104. arm_rfft_init_q15(&this->instRfftQ15 ,512,0,1);
  105. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  106. this->ifft=0;
  107. break;
  108. case TransformRQ15::TEST_RFFT_Q15_21:
  109. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_512_Q15_ID,mgr);
  110. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_512_Q15_ID,mgr);
  111. arm_rfft_init_q15(&this->instRfftQ15 ,512,1,1);
  112. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  113. this->ifft=1;
  114. this->scaling=9;
  115. break;
  116. case TransformRQ15::TEST_RFFT_Q15_6:
  117. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_1024_Q15_ID,mgr);
  118. ref.reload( TransformRQ15::REF_RFFT_NOISY_1024_Q15_ID,mgr);
  119. arm_rfft_init_q15(&this->instRfftQ15 ,1024,0,1);
  120. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  121. this->ifft=0;
  122. break;
  123. case TransformRQ15::TEST_RFFT_Q15_22:
  124. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_1024_Q15_ID,mgr);
  125. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_1024_Q15_ID,mgr);
  126. arm_rfft_init_q15(&this->instRfftQ15 ,1024,1,1);
  127. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  128. this->ifft=1;
  129. this->scaling=10;
  130. break;
  131. case TransformRQ15::TEST_RFFT_Q15_7:
  132. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_2048_Q15_ID,mgr);
  133. ref.reload( TransformRQ15::REF_RFFT_NOISY_2048_Q15_ID,mgr);
  134. arm_rfft_init_q15(&this->instRfftQ15 ,2048,0,1);
  135. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  136. this->ifft=0;
  137. break;
  138. case TransformRQ15::TEST_RFFT_Q15_23:
  139. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_2048_Q15_ID,mgr);
  140. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_2048_Q15_ID,mgr);
  141. arm_rfft_init_q15(&this->instRfftQ15 ,2048,1,1);
  142. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  143. this->ifft=1;
  144. this->scaling=11;
  145. break;
  146. case TransformRQ15::TEST_RFFT_Q15_8:
  147. input.reload(TransformRQ15::INPUTS_RFFT_NOISY_4096_Q15_ID,mgr);
  148. ref.reload( TransformRQ15::REF_RFFT_NOISY_4096_Q15_ID,mgr);
  149. arm_rfft_init_q15(&this->instRfftQ15 ,4096,0,1);
  150. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  151. this->ifft=0;
  152. break;
  153. case TransformRQ15::TEST_RFFT_Q15_24:
  154. input.reload(TransformRQ15::INPUTS_RIFFT_NOISY_4096_Q15_ID,mgr);
  155. ref.reload( TransformRQ15::INPUTS_RFFT_NOISY_4096_Q15_ID,mgr);
  156. arm_rfft_init_q15(&this->instRfftQ15 ,4096,1,1);
  157. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  158. this->ifft=1;
  159. this->scaling=12;
  160. break;
  161. /* STEP FUNCTIONS */
  162. case TransformRQ15::TEST_RFFT_Q15_9:
  163. input.reload(TransformRQ15::INPUTS_RFFT_STEP_32_Q15_ID,mgr);
  164. ref.reload( TransformRQ15::REF_RFFT_STEP_32_Q15_ID,mgr);
  165. arm_rfft_init_q15(&this->instRfftQ15 ,32,0,1);
  166. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  167. this->ifft=0;
  168. break;
  169. case TransformRQ15::TEST_RFFT_Q15_25:
  170. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_32_Q15_ID,mgr);
  171. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_32_Q15_ID,mgr);
  172. arm_rfft_init_q15(&this->instRfftQ15 ,32,1,1);
  173. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  174. this->ifft=1;
  175. this->scaling=5;
  176. break;
  177. case TransformRQ15::TEST_RFFT_Q15_10:
  178. input.reload(TransformRQ15::INPUTS_RFFT_STEP_64_Q15_ID,mgr);
  179. ref.reload( TransformRQ15::REF_RFFT_STEP_64_Q15_ID,mgr);
  180. arm_rfft_init_q15(&this->instRfftQ15 ,64,0,1);
  181. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  182. this->ifft=0;
  183. break;
  184. case TransformRQ15::TEST_RFFT_Q15_26:
  185. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_64_Q15_ID,mgr);
  186. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_64_Q15_ID,mgr);
  187. arm_rfft_init_q15(&this->instRfftQ15 ,64,1,1);
  188. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  189. this->ifft=1;
  190. this->scaling=6;
  191. break;
  192. case TransformRQ15::TEST_RFFT_Q15_11:
  193. input.reload(TransformRQ15::INPUTS_RFFT_STEP_128_Q15_ID,mgr);
  194. ref.reload( TransformRQ15::REF_RFFT_STEP_128_Q15_ID,mgr);
  195. arm_rfft_init_q15(&this->instRfftQ15 ,128,0,1);
  196. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  197. this->ifft=0;
  198. break;
  199. case TransformRQ15::TEST_RFFT_Q15_27:
  200. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_128_Q15_ID,mgr);
  201. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_128_Q15_ID,mgr);
  202. arm_rfft_init_q15(&this->instRfftQ15 ,128,1,1);
  203. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  204. this->ifft=1;
  205. this->scaling=7;
  206. break;
  207. case TransformRQ15::TEST_RFFT_Q15_12:
  208. input.reload(TransformRQ15::INPUTS_RFFT_STEP_256_Q15_ID,mgr);
  209. ref.reload( TransformRQ15::REF_RFFT_STEP_256_Q15_ID,mgr);
  210. arm_rfft_init_q15(&this->instRfftQ15 ,256,0,1);
  211. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  212. this->ifft=0;
  213. break;
  214. case TransformRQ15::TEST_RFFT_Q15_28:
  215. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_256_Q15_ID,mgr);
  216. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_256_Q15_ID,mgr);
  217. arm_rfft_init_q15(&this->instRfftQ15 ,256,1,1);
  218. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  219. this->ifft=1;
  220. this->scaling=8;
  221. break;
  222. case TransformRQ15::TEST_RFFT_Q15_13:
  223. input.reload(TransformRQ15::INPUTS_RFFT_STEP_512_Q15_ID,mgr);
  224. ref.reload( TransformRQ15::REF_RFFT_STEP_512_Q15_ID,mgr);
  225. arm_rfft_init_q15(&this->instRfftQ15 ,512,0,1);
  226. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  227. this->ifft=0;
  228. break;
  229. case TransformRQ15::TEST_RFFT_Q15_29:
  230. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_512_Q15_ID,mgr);
  231. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_512_Q15_ID,mgr);
  232. arm_rfft_init_q15(&this->instRfftQ15 ,512,1,1);
  233. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  234. this->ifft=1;
  235. this->scaling=9;
  236. break;
  237. case TransformRQ15::TEST_RFFT_Q15_14:
  238. input.reload(TransformRQ15::INPUTS_RFFT_STEP_1024_Q15_ID,mgr);
  239. ref.reload( TransformRQ15::REF_RFFT_STEP_1024_Q15_ID,mgr);
  240. arm_rfft_init_q15(&this->instRfftQ15 ,1024,0,1);
  241. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  242. this->ifft=0;
  243. break;
  244. case TransformRQ15::TEST_RFFT_Q15_30:
  245. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_1024_Q15_ID,mgr);
  246. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_1024_Q15_ID,mgr);
  247. arm_rfft_init_q15(&this->instRfftQ15 ,1024,1,1);
  248. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  249. this->ifft=1;
  250. this->scaling=10;
  251. break;
  252. case TransformRQ15::TEST_RFFT_Q15_15:
  253. input.reload(TransformRQ15::INPUTS_RFFT_STEP_2048_Q15_ID,mgr);
  254. ref.reload( TransformRQ15::REF_RFFT_STEP_2048_Q15_ID,mgr);
  255. arm_rfft_init_q15(&this->instRfftQ15 ,2048,0,1);
  256. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  257. this->ifft=0;
  258. break;
  259. case TransformRQ15::TEST_RFFT_Q15_31:
  260. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_2048_Q15_ID,mgr);
  261. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_2048_Q15_ID,mgr);
  262. arm_rfft_init_q15(&this->instRfftQ15 ,2048,1,1);
  263. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  264. this->ifft=1;
  265. this->scaling=11;
  266. break;
  267. case TransformRQ15::TEST_RFFT_Q15_16:
  268. input.reload(TransformRQ15::INPUTS_RFFT_STEP_4096_Q15_ID,mgr);
  269. ref.reload( TransformRQ15::REF_RFFT_STEP_4096_Q15_ID,mgr);
  270. arm_rfft_init_q15(&this->instRfftQ15 ,4096,0,1);
  271. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  272. this->ifft=0;
  273. break;
  274. case TransformRQ15::TEST_RFFT_Q15_32:
  275. input.reload(TransformRQ15::INPUTS_RIFFT_STEP_4096_Q15_ID,mgr);
  276. ref.reload( TransformRQ15::INPUTS_RFFT_STEP_4096_Q15_ID,mgr);
  277. arm_rfft_init_q15(&this->instRfftQ15 ,4096,1,1);
  278. inputchanged.create(input.nbSamples(),TransformRQ15::TEMP_Q15_ID,mgr);
  279. this->ifft=1;
  280. this->scaling=12;
  281. break;
  282. }
  283. outputfft.create(ref.nbSamples(),TransformRQ15::OUTPUT_RFFT_Q15_ID,mgr);
  284. /*
  285. RFFT is writing more samples than it should.
  286. This is a temporary buffer allowing the test to pass.
  287. */
  288. overheadoutputfft.create(2*ref.nbSamples(),TransformRQ15::FULLOUTPUT_Q15_ID,mgr);
  289. }
  290. void TransformRQ15::tearDown(Testing::testID_t id,Client::PatternMgr *mgr)
  291. {
  292. (void)id;
  293. outputfft.dump(mgr);
  294. overheadoutputfft.dump(mgr);
  295. }