TransformRF64.cpp 13 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474
  1. #include "TransformRF64.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 250
  8. void TransformRF64::test_rfft_f64()
  9. {
  10. float64_t *inp = input.ptr();
  11. float64_t *tmp = inputchanged.ptr();
  12. float64_t *outp = outputfft.ptr();
  13. memcpy(tmp,inp,sizeof(float64_t)*input.nbSamples());
  14. arm_rfft_fast_f64(
  15. &this->instRfftF64,
  16. tmp,
  17. outp,
  18. this->ifft);
  19. ASSERT_SNR(outputfft,ref,(float64_t)SNR_THRESHOLD);
  20. ASSERT_EMPTY_TAIL(outputfft);
  21. }
  22. void TransformRF64::setUp(Testing::testID_t id,std::vector<Testing::param_t>& paramsArgs,Client::PatternMgr *mgr)
  23. {
  24. switch(id)
  25. {
  26. case TransformRF64::TEST_RFFT_F64_1:
  27. input.reload(TransformRF64::INPUTS_RFFT_NOISY_32_F64_ID,mgr);
  28. ref.reload( TransformRF64::REF_RFFT_NOISY_32_F64_ID,mgr);
  29. arm_rfft_fast_init_f64(&this->instRfftF64 ,32);
  30. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  31. this->ifft=0;
  32. break;
  33. case TransformRF64::TEST_RFFT_F64_17:
  34. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_64_F64_ID,mgr);
  35. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_64_F64_ID,mgr);
  36. arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
  37. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  38. this->ifft=1;
  39. break;
  40. case TransformRF64::TEST_RFFT_F64_2:
  41. input.reload(TransformRF64::INPUTS_RFFT_NOISY_64_F64_ID,mgr);
  42. ref.reload( TransformRF64::REF_RFFT_NOISY_64_F64_ID,mgr);
  43. arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
  44. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  45. this->ifft=0;
  46. break;
  47. case TransformRF64::TEST_RFFT_F64_18:
  48. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_64_F64_ID,mgr);
  49. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_64_F64_ID,mgr);
  50. arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
  51. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  52. this->ifft=1;
  53. break;
  54. case TransformRF64::TEST_RFFT_F64_3:
  55. input.reload(TransformRF64::INPUTS_RFFT_NOISY_128_F64_ID,mgr);
  56. ref.reload( TransformRF64::REF_RFFT_NOISY_128_F64_ID,mgr);
  57. arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
  58. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  59. this->ifft=0;
  60. break;
  61. case TransformRF64::TEST_RFFT_F64_19:
  62. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_128_F64_ID,mgr);
  63. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_128_F64_ID,mgr);
  64. arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
  65. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  66. this->ifft=1;
  67. break;
  68. case TransformRF64::TEST_RFFT_F64_4:
  69. input.reload(TransformRF64::INPUTS_RFFT_NOISY_256_F64_ID,mgr);
  70. ref.reload( TransformRF64::REF_RFFT_NOISY_256_F64_ID,mgr);
  71. arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
  72. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  73. this->ifft=0;
  74. break;
  75. case TransformRF64::TEST_RFFT_F64_20:
  76. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_256_F64_ID,mgr);
  77. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_256_F64_ID,mgr);
  78. arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
  79. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  80. this->ifft=1;
  81. break;
  82. case TransformRF64::TEST_RFFT_F64_5:
  83. input.reload(TransformRF64::INPUTS_RFFT_NOISY_512_F64_ID,mgr);
  84. ref.reload( TransformRF64::REF_RFFT_NOISY_512_F64_ID,mgr);
  85. arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
  86. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  87. this->ifft=0;
  88. break;
  89. case TransformRF64::TEST_RFFT_F64_21:
  90. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_512_F64_ID,mgr);
  91. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_512_F64_ID,mgr);
  92. arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
  93. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  94. this->ifft=1;
  95. break;
  96. case TransformRF64::TEST_RFFT_F64_6:
  97. input.reload(TransformRF64::INPUTS_RFFT_NOISY_1024_F64_ID,mgr);
  98. ref.reload( TransformRF64::REF_RFFT_NOISY_1024_F64_ID,mgr);
  99. arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
  100. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  101. this->ifft=0;
  102. break;
  103. case TransformRF64::TEST_RFFT_F64_22:
  104. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_1024_F64_ID,mgr);
  105. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_1024_F64_ID,mgr);
  106. arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
  107. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  108. this->ifft=1;
  109. break;
  110. case TransformRF64::TEST_RFFT_F64_7:
  111. input.reload(TransformRF64::INPUTS_RFFT_NOISY_2048_F64_ID,mgr);
  112. ref.reload( TransformRF64::REF_RFFT_NOISY_2048_F64_ID,mgr);
  113. arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
  114. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  115. this->ifft=0;
  116. break;
  117. case TransformRF64::TEST_RFFT_F64_23:
  118. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_2048_F64_ID,mgr);
  119. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_2048_F64_ID,mgr);
  120. arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
  121. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  122. this->ifft=1;
  123. break;
  124. case TransformRF64::TEST_RFFT_F64_8:
  125. input.reload(TransformRF64::INPUTS_RFFT_NOISY_4096_F64_ID,mgr);
  126. ref.reload( TransformRF64::REF_RFFT_NOISY_4096_F64_ID,mgr);
  127. arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
  128. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  129. this->ifft=0;
  130. break;
  131. case TransformRF64::TEST_RFFT_F64_24:
  132. input.reload(TransformRF64::INPUTS_RIFFT_NOISY_4096_F64_ID,mgr);
  133. ref.reload( TransformRF64::INPUTS_RFFT_NOISY_4096_F64_ID,mgr);
  134. arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
  135. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  136. this->ifft=1;
  137. break;
  138. /* STEP FUNCTIONS */
  139. case TransformRF64::TEST_RFFT_F64_9:
  140. input.reload(TransformRF64::INPUTS_RFFT_STEP_32_F64_ID,mgr);
  141. ref.reload( TransformRF64::REF_RFFT_STEP_32_F64_ID,mgr);
  142. arm_rfft_fast_init_f64(&this->instRfftF64 ,32);
  143. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  144. this->ifft=0;
  145. break;
  146. case TransformRF64::TEST_RFFT_F64_25:
  147. input.reload(TransformRF64::INPUTS_RIFFT_STEP_64_F64_ID,mgr);
  148. ref.reload( TransformRF64::INPUTS_RFFT_STEP_64_F64_ID,mgr);
  149. arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
  150. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  151. this->ifft=1;
  152. break;
  153. case TransformRF64::TEST_RFFT_F64_10:
  154. input.reload(TransformRF64::INPUTS_RFFT_STEP_64_F64_ID,mgr);
  155. ref.reload( TransformRF64::REF_RFFT_STEP_64_F64_ID,mgr);
  156. arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
  157. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  158. this->ifft=0;
  159. break;
  160. case TransformRF64::TEST_RFFT_F64_26:
  161. input.reload(TransformRF64::INPUTS_RIFFT_STEP_64_F64_ID,mgr);
  162. ref.reload( TransformRF64::INPUTS_RFFT_STEP_64_F64_ID,mgr);
  163. arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
  164. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  165. this->ifft=1;
  166. break;
  167. case TransformRF64::TEST_RFFT_F64_11:
  168. input.reload(TransformRF64::INPUTS_RFFT_STEP_128_F64_ID,mgr);
  169. ref.reload( TransformRF64::REF_RFFT_STEP_128_F64_ID,mgr);
  170. arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
  171. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  172. this->ifft=0;
  173. break;
  174. case TransformRF64::TEST_RFFT_F64_27:
  175. input.reload(TransformRF64::INPUTS_RIFFT_STEP_128_F64_ID,mgr);
  176. ref.reload( TransformRF64::INPUTS_RFFT_STEP_128_F64_ID,mgr);
  177. arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
  178. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  179. this->ifft=1;
  180. break;
  181. case TransformRF64::TEST_RFFT_F64_12:
  182. input.reload(TransformRF64::INPUTS_RFFT_STEP_256_F64_ID,mgr);
  183. ref.reload( TransformRF64::REF_RFFT_STEP_256_F64_ID,mgr);
  184. arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
  185. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  186. this->ifft=0;
  187. break;
  188. case TransformRF64::TEST_RFFT_F64_28:
  189. input.reload(TransformRF64::INPUTS_RIFFT_STEP_256_F64_ID,mgr);
  190. ref.reload( TransformRF64::INPUTS_RFFT_STEP_256_F64_ID,mgr);
  191. arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
  192. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  193. this->ifft=1;
  194. break;
  195. case TransformRF64::TEST_RFFT_F64_13:
  196. input.reload(TransformRF64::INPUTS_RFFT_STEP_512_F64_ID,mgr);
  197. ref.reload( TransformRF64::REF_RFFT_STEP_512_F64_ID,mgr);
  198. arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
  199. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  200. this->ifft=0;
  201. break;
  202. case TransformRF64::TEST_RFFT_F64_29:
  203. input.reload(TransformRF64::INPUTS_RIFFT_STEP_512_F64_ID,mgr);
  204. ref.reload( TransformRF64::INPUTS_RFFT_STEP_512_F64_ID,mgr);
  205. arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
  206. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  207. this->ifft=1;
  208. break;
  209. case TransformRF64::TEST_RFFT_F64_14:
  210. input.reload(TransformRF64::INPUTS_RFFT_STEP_1024_F64_ID,mgr);
  211. ref.reload( TransformRF64::REF_RFFT_STEP_1024_F64_ID,mgr);
  212. arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
  213. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  214. this->ifft=0;
  215. break;
  216. case TransformRF64::TEST_RFFT_F64_30:
  217. input.reload(TransformRF64::INPUTS_RIFFT_STEP_1024_F64_ID,mgr);
  218. ref.reload( TransformRF64::INPUTS_RFFT_STEP_1024_F64_ID,mgr);
  219. arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
  220. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  221. this->ifft=1;
  222. break;
  223. case TransformRF64::TEST_RFFT_F64_15:
  224. input.reload(TransformRF64::INPUTS_RFFT_STEP_2048_F64_ID,mgr);
  225. ref.reload( TransformRF64::REF_RFFT_STEP_2048_F64_ID,mgr);
  226. arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
  227. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  228. this->ifft=0;
  229. break;
  230. case TransformRF64::TEST_RFFT_F64_31:
  231. input.reload(TransformRF64::INPUTS_RIFFT_STEP_2048_F64_ID,mgr);
  232. ref.reload( TransformRF64::INPUTS_RFFT_STEP_2048_F64_ID,mgr);
  233. arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
  234. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  235. this->ifft=1;
  236. break;
  237. case TransformRF64::TEST_RFFT_F64_16:
  238. input.reload(TransformRF64::INPUTS_RFFT_STEP_4096_F64_ID,mgr);
  239. ref.reload( TransformRF64::REF_RFFT_STEP_4096_F64_ID,mgr);
  240. arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
  241. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  242. this->ifft=0;
  243. break;
  244. case TransformRF64::TEST_RFFT_F64_32:
  245. input.reload(TransformRF64::INPUTS_RIFFT_STEP_4096_F64_ID,mgr);
  246. ref.reload( TransformRF64::INPUTS_RFFT_STEP_4096_F64_ID,mgr);
  247. arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
  248. inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
  249. this->ifft=1;
  250. break;
  251. }
  252. outputfft.create(ref.nbSamples(),TransformRF64::OUTPUT_RFFT_F64_ID,mgr);
  253. }
  254. void TransformRF64::tearDown(Testing::testID_t id,Client::PatternMgr *mgr)
  255. {
  256. outputfft.dump(mgr);
  257. }