AnsweredAssumed Answered

Non-deterministic results from dpotrf when re-running problem on the same machine

Question asked by rvstebbing on Sep 6, 2013
Latest reply on Sep 18, 2013 by rvstebbing

Dear List,


I am using ACML 5.3.1 for Windows (ifort64, 32-bit int build) with Microsoft Visual Studio 2012 (Update 3). I have a program to do the in-place Cholesky factorisation of a 126x126 matrix using dpotrf and I am linking against libacml_mp_dll.lib. When I run the program and dump the raw doubles (as hex) to file, I don't get consistent output between runs on the same machine. Specifically, there are differences in the lower-bits of some entries. However, when I run the factorisation multiple times within a single program instance the results are consistent.


Is this considered a bug?


The C++ to generate the dumps is here: ACML 5.3.1 dpotrf example - 8b959cb7

The Python code to run new instances until a difference is found is here: - dc2d3dca

Example output files are ACML_dpotrf_0.txt - 6f0435ad and ACML_dpotrf_1.txt - e8d61509


I appreciate any feedback you can give on this issue.