| 73 | |
| 74 | template <typename T_time,typename T_freq> |
| 75 | void fwd_inv(size_t nfft) |
| 76 | { |
| 77 | typedef typename NumTraits<T_freq>::Real Scalar; |
| 78 | vector<T_time> timebuf(nfft); |
| 79 | RandomFill(timebuf); |
| 80 | |
| 81 | vector<T_freq> freqbuf; |
| 82 | static FFT<Scalar> fft; |
| 83 | fft.fwd(freqbuf,timebuf); |
| 84 | |
| 85 | vector<T_time> timebuf2; |
| 86 | fft.inv(timebuf2,freqbuf); |
| 87 | |
| 88 | long double rmse = mag2(timebuf - timebuf2) / mag2(timebuf); |
| 89 | cout << "roundtrip rmse: " << rmse << endl; |
| 90 | } |
| 91 | |
| 92 | template <typename T_scalar> |
| 93 | void two_demos(int nfft) |
nothing calls this directly
no test coverage detected