Reproducible, Accurately Rounded and Efficient BLAS
Reproducible, Accurately Rounded and Efficient BLAS
复制标题
可重复、精确舍入且高效的 BLAS
DOI:
--
复制
发表时间:
2016
期刊:
影响因子:
--
通讯作者:
David Parello
中科院分区:
文献类型:
--
作者:
Chemseddine Chohra;P. Langlois;David Parello
Numerical reproducibility failures rise in parallel computation because floating-point summation is non-associative. Massively parallel and optimized executions dynamically modify the floating-point operation order. Hence, numerical results may change from one run to another. We propose to ensure reproducibility by extending as far as possible the IEEE-754 correct rounding property to larger operation sequences. We introduce our RARE-BLAS (Reproducible, Accurately Rounded and Efficient BLAS) that benefits from recent accurate and efficient summation algorithms. Solutions for level 1 (asum, dot and nrm2) and level 2 (gemv) routines are presented. Their performance is studied compared to the Intel MKL library and other existing reproducible algorithms. For both shared and distributed memory parallel systems, we exhibit an extra-cost of 2( imes ) in the worst case scenario, which is satisfying for a wide range of applications. For Intel Xeon Phi accelerator a larger extra-cost (4( imes ) to 6( imes )) is observed, which is still helpful at least for debugging and validation steps.