Mitigating GPU Moments Calculation Differences
By changing the algorithm that is used to calculate the eigenvalues and eigenvectors of the 3×3 matrix that gives the shower axis, we can improve the GPU-CPU differences, especially for Δα, Δφ and Δθ.
The new algorithm matches the non-iterative version used by Eigen, leveraging GPU intrinsics for better precision and performance.