| 74 | } |
| 75 | |
| 76 | void SpectrumStack::set_sample_num(uint64_t num) |
| 77 | { |
| 78 | _sample_num = num; |
| 79 | _xn.resize(_sample_num); |
| 80 | _xk.resize(_sample_num); |
| 81 | _power_spectrum.resize(_sample_num/2+1); |
| 82 | _fft_plan = fftw_plan_r2r_1d(_sample_num, _xn.data(), _xk.data(), |
| 83 | FFTW_R2HC, FFTW_ESTIMATE); |
| 84 | } |
| 85 | |
| 86 | int SpectrumStack::get_windows_index() |
| 87 | { |