On Saturday, 27 June 2015 at 12:17:31 UTC, Timon Gehr wrote:


This computes a²·a̅ instead of a·a̅.

What is the source code for residualPowerOf2?

Also is there any performance issues? can I make this faster?


Probably you should use http://dlang.org/phobos/std_complex.html#.sqAbs instead. You then also don't need the final map to extract the real part.

You are absulately right instead
fftResult.zip(fftResult.map!(a => a * a.conj()))
I corrected it to
fftResult.zip(fftResult.map!(a => a.conj()))

Thanks a lot

Reply via email to