I am trying to implement a phase correlation function using reikna (more info available under method section of this wikipedia article
). I was hoping that you might be able to give me a few tips as I'm getting started in order to avoid making silly decisions on program structure and usage of the reikna library.
For starters, one part of this task is to take the complex conjugate. I was thinking I'd use the split_complex transformation on the input, multiply the complex component by -1, and then recombine the previously split data using the combine_complex transformation. Is there a better way to do this?
I was thinking it may be most efficient and clean to do the whole operation (shown in below pictures), in one computation with several transformations connected, do you think this is a good strategy?
The phase correlation process also requires an FFT and IFFT, but these operations are both straightforward in the reikna library (although I'm unsure if these parts should also be connected to the previously discussed math).
One last thing: is there any way to efficiently compute the maximum using the reikna (or PyOpenCL) interface (I searched but wasn't able to find much of anything)?
Thank you for all your work here on the forum and on reikna itself!