To put it more bluntly, the quest for exact numerical reproducibility with floating point numbers is often very difficult, and in most cases not worth the effort. If small differences in input or differences that should be inconsequential (BLAS library, OS, architecture, etc) cause large differences in output, then either there’s a bug in one of the libraries, or your algorithm is unstable.
7 Likes