| 123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473 |
- #include "TransformRF64.h"
- #include <stdio.h>
- #include "Error.h"
- #include "Test.h"
- #define SNR_THRESHOLD 250
- void TransformRF64::test_rfft_f64()
- {
- float64_t *inp = input.ptr();
- float64_t *tmp = inputchanged.ptr();
- float64_t *outp = outputfft.ptr();
- memcpy(tmp,inp,sizeof(float64_t)*input.nbSamples());
- arm_rfft_fast_f64(
- &this->instRfftF64,
- tmp,
- outp,
- this->ifft);
- ASSERT_SNR(outputfft,ref,(float64_t)SNR_THRESHOLD);
- ASSERT_EMPTY_TAIL(outputfft);
- }
- void TransformRF64::setUp(Testing::testID_t id,std::vector<Testing::param_t>& paramsArgs,Client::PatternMgr *mgr)
- {
- (void)paramsArgs;
- switch(id)
- {
- case TransformRF64::TEST_RFFT_F64_1:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_32_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_32_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,32);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_17:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_64_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_64_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_2:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_64_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_64_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_18:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_64_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_64_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_3:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_128_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_128_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_19:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_128_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_128_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_4:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_256_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_256_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_20:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_256_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_256_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_5:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_512_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_512_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_21:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_512_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_512_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_6:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_1024_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_1024_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_22:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_1024_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_1024_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_7:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_2048_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_2048_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_23:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_2048_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_2048_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_8:
- input.reload(TransformRF64::INPUTS_RFFT_NOISY_4096_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_NOISY_4096_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_24:
- input.reload(TransformRF64::INPUTS_RIFFT_NOISY_4096_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_NOISY_4096_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- /* STEP FUNCTIONS */
- case TransformRF64::TEST_RFFT_F64_9:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_32_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_32_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,32);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_25:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_64_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_64_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_10:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_64_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_64_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_26:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_64_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_64_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,64);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_11:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_128_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_128_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_27:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_128_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_128_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,128);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_12:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_256_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_256_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_28:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_256_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_256_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,256);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_13:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_512_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_512_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_29:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_512_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_512_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,512);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_14:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_1024_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_1024_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_30:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_1024_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_1024_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,1024);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_15:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_2048_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_2048_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_31:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_2048_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_2048_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,2048);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- case TransformRF64::TEST_RFFT_F64_16:
- input.reload(TransformRF64::INPUTS_RFFT_STEP_4096_F64_ID,mgr);
- ref.reload( TransformRF64::REF_RFFT_STEP_4096_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=0;
- break;
- case TransformRF64::TEST_RFFT_F64_32:
- input.reload(TransformRF64::INPUTS_RIFFT_STEP_4096_F64_ID,mgr);
- ref.reload( TransformRF64::INPUTS_RFFT_STEP_4096_F64_ID,mgr);
- arm_rfft_fast_init_f64(&this->instRfftF64 ,4096);
- inputchanged.create(input.nbSamples(),TransformRF64::TEMP_F64_ID,mgr);
- this->ifft=1;
- break;
- }
- outputfft.create(ref.nbSamples(),TransformRF64::OUTPUT_RFFT_F64_ID,mgr);
- }
- void TransformRF64::tearDown(Testing::testID_t id,Client::PatternMgr *mgr)
- {
- (void)id;
- outputfft.dump(mgr);
- }
|