TransformRQ15.cpp 14 KB

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