Intel® oneAPI Math Kernel Library
Ask questions and share information with other developers who use Intel® Math Kernel Library.

no consistent result with mkl functions

Emily_Zhu
Beginner
409 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
409 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
Emily_Zhu
Beginner
409 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
TimP
Honored Contributor III
409 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
Emily_Zhu
Beginner
409 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
Reply