Force Deterministic Algorithm in FFTW?

Questions and postings pertaining to the development of ImageMagick, feature enhancements, and ImageMagick internals. ImageMagick source code and algorithms are discussed here. Usage questions which are too arcane for the normal user list should also be posted here.
Post Reply
blackhaz

Force Deterministic Algorithm in FFTW?

Post by blackhaz »

Is it possible to force Magick++ to use deterministic algorithm while doing Fourier transforms?

It looks like my Magick++ app produces slightly different images, sporadically. I suspect this is the fftw issue. I am doing deconvolution so most probably small errors multiply and I am able to notice the difference. I believe it has something to do with fftw running in non-deterministic mode:
http://www.fftw.org/faq/section3.html#nondeterministic

Thank you!
User avatar
fmw42
Posts: 25562
Joined: 2007-07-02T17:14:51-07:00
Authentication code: 1152
Location: Sunnyvale, California, USA

Re: Force Deterministic Algorithm in FFTW?

Post by fmw42 »

This was resolved as a non-issue, I am told. The FFT is in deterministic mode. The problem was really a bad OpenMP that failed sporatically on +ift on Dual Core Intel-based Mac. (So OpenMP needed to be disabled).
Post Reply