search

astro-informatics / purify

Next-generation radio interferometric imaging.
https://astro-informatics.github.io/purify
GNU General Public License v2.0
17 stars 13 forks source link

Replaced Eigen FFT operator with direct FFTW call. #113

Closed Luke-Pratley closed 6 years ago

Luke-Pratley commented 6 years ago

I have updated the code so that the functionality to not use Eigen for 2D FFTs exists. This is needed if we want to move to up to date versions of Eigen (stop using the separate fork).

One of the changes involves the FFTW plans being set using an enum class rather than a string. Also, the main FFTW operator in operators.h calls FFTW directly.

At the time of making this, the unit tests pass! Though, the FFT normalisaiton is something that had to be chosen that might be different.

Luke-Pratley commented 6 years ago

Okay, now there should be no dependency on the FFTOperator.h class/code. We should look at removing this and the old MeasurementOperator.h soon. Then we can update Eigen.