PU speedup (tril(K) vs tril(K)): %f\n',time_h/time_d);
%% Difference between results

fprintf('Difference between results:\n');
K2 = tril(K_f);
fprintf('\tCPU vs CPU (K vs tril(K)): %e\n',norm(K2(:)-K_s(:)));
K_d2 = gather(K_d);
fprintf('\tCPU vs GPU (tril(K) vs tril(K)): %e\n',norm(K_s(:)-K_d2(:)));
##### SOURCE END #####
--></body></html>