Community
cancel
Showing results for 
Search instead for 
Did you mean: 
Highlighted
Beginner
17 Views

no consistent result with mkl functions

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.
0 Kudos
4 Replies
Highlighted
Employee
17 Views

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.
0 Kudos
Highlighted
Beginner
17 Views

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?
0 Kudos
Highlighted
Black Belt
17 Views

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.
0 Kudos
Highlighted
Beginner
17 Views

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?


0 Kudos