Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.
Announcements
Welcome to the Intel Community. If you get an answer you like, please mark it as an Accepted Solution to help others. Thank you!
6441 Discussions

no consistent result with mkl functions

Emily_Zhu
Beginner
111 Views
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
Todd_R_Intel
Employee
111 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.
Emily_Zhu
Beginner
111 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?
TimP
Black Belt
111 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.
Emily_Zhu
Beginner
111 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?


Reply