diff options
| author | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-05 11:12:15 -0300 |
|---|---|---|
| committer | Jaron Kent-Dobias <jaron@kent-dobias.com> | 2025-04-05 11:12:15 -0300 |
| commit | 96875ddafbf1ff821a4a1fd628c2e4758096700c (patch) | |
| tree | 16d00b3727c5273233be6af37c29361cdbf100a6 /fourier.cpp | |
| parent | 44ed3056e815910b7f61c62a9b71976a3c47f4b3 (diff) | |
| download | code-96875ddafbf1ff821a4a1fd628c2e4758096700c.tar.gz code-96875ddafbf1ff821a4a1fd628c2e4758096700c.tar.bz2 code-96875ddafbf1ff821a4a1fd628c2e4758096700c.zip | |
Save and read fftw wisdom to speed initialization
Diffstat (limited to 'fourier.cpp')
| -rw-r--r-- | fourier.cpp | 2 |
1 files changed, 2 insertions, 0 deletions
diff --git a/fourier.cpp b/fourier.cpp index 1bef4b9..0c22637 100644 --- a/fourier.cpp +++ b/fourier.cpp @@ -25,8 +25,10 @@ Real ddf(Real λ, unsigned p, unsigned s, Real q) { } FourierTransform::FourierTransform(unsigned n, Real Δω, Real Δτ, unsigned flags) : a(2 * n), â(n + 1), Δω(Δω), Δτ(Δτ) { + fftw_import_wisdom_from_filename("fftw.wisdom"); plan_r2c = fftw_plan_dft_r2c_1d(2 * n, a.data(), reinterpret_cast<fftw_complex*>(â.data()), flags); plan_c2r = fftw_plan_dft_c2r_1d(2 * n, reinterpret_cast<fftw_complex*>(â.data()), a.data(), flags); + fftw_export_wisdom_to_filename("fftw.wisdom"); } FourierTransform::~FourierTransform() { |