Multithreading lapack numeric differences

Multithreading lapack numeric differences


I'm having different outputs with zgbsv LAPACK subroutine depending on the number of threads that are being used. My MKL dll version is This numeric differences are in the order of the 10th decimal number, but they exist. Is this normal?

Thanks in advance,


4 posts / 0 new
Last post
For more complete information about compiler optimizations, see our Optimization Notice.

yes, this is the expected behaviour with floating point operations. for your information, the latest version 11.0 introduces new functionality allowing users to have reproducibility in floating point calculations. Please look at this material regarding this topic:


Still, using mkl_cbwr_set(MKL_CBWR_COMPATIBLE) before LAPACK, it returns me different results with 3 or 4 threads for example. Am I using it fine?

MKL's CNR guarantees the bit to bit output results only if with the same #of threads.
Please participate this survey ( )
to influence how this feature will be extended in the future.

Leave a Comment

Please sign in to add a comment. Not a member? Join today