GPU and CPU computing results different

Question

peter 2011-7-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/11521-gpu-and-cpu-computing-results-different

Hi, I have tested FFT2 function for both workspace and GPU data. The difference in result is very small but still disturbs me. Does anyone have why that is and any suggested solutions?

Thanks,

im = phantom(512);
k = fft2(im);
kim = fft2(gpuArray(complex(im)));
OK = isequal(k,kim);

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ben Tordoff 2011-7-13

4
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/11521-gpu-and-cpu-computing-results-different#answer_15717

在 MATLAB Online 中打开

if you look at the magnitude of the difference compared to the magnitude of the values, it is pretty small:

   >> max(max(gather(abs(k-kim)))) / max(max(abs(k)))
   ans = 1.0144e-016

Infact, it is similar to the limit of accuracy for any calculation on the input:

   >> eps(max(max(abs(im))))
   ans = 2.2204e-016

(the two are not measuring the same thing, but hopefully you get the idea).

The CPU and GPU implementations of FFT/FFT2 are necessarily quite different in order to take best advantage of the hardware. The GPU version needs to be massively parallel. I believe that the difference you are seeing is well within what one should expect from any two different implementations of FFT/FFT2.

Whilst we can aspire to identical results from the CPU and GPU versions, sometimes it simply isn't possible without massive compromises on speed.