Is PZHEEVD deterministic? By that I mean, if I fix the input matrix and run PZHEEVD repeatedly on a fixed number of MPI tasks on the same architecture will I get the same result to machine precision every time. Let's assume everything is single-threaded, e.g. no threaded BLAS.
After some discussion, with a vendor (IBM) and the MPICH2 developers I am told that collectives are deterministic.
After taking a closer look at PZHEEVD, the only non-deterministic behavior I see is in the routine PZHETRD which uses BI_TreeComb which itself does not look deterministic from the use of MPI_ANY_SOURCE (define to the variable BANYNODE).
Is there a deterministic version of BI_TreeComb available?
Dr. Nichols A. Romero
Argonne Leadership Computing Facility