no consistent result with mkl functions

[Emily_Zhu]

Is anyone can solve the problem that mkl functions didn't get consistant result?

we use 10.3.3, It looks any function like DGEMM, DSYEV get different results

The different is big since for 1000*1000 matrix, the max difference for single element is

the matrix is 1.0e-14.

We also try 16 boundary alignment, but it doesn't work.

[Todd_R_Intel]

We've generally seen that if you are running your program on the same computer, with the same number of threads and your memory is aligned on 16-byte boundaries then you'll get consistent results from DGEMM from one run to the next.

We'd need to know a little more about your program and the system on which you are running.

[Emily_Zhu]

If I set each computer as mkl_sequential and use 16 boundary alignment.

but send whole computation to cluster, if the final result will be same?

[TimP]

Quoting Emily Zhu

If I set each computer as mkl_sequential and use 16 boundary alignment.

but send whole computation to cluster, if the final result will be same?

If each node has the same software installation and same CPU brand, it's reasonable to hope for identical results.

[Emily_Zhu]

Another question is that alought DGEMM is repeatable, if the result will shift with real result of matrix

multiplication. And how much it will shift?