Quadruple-Precision BLAS Using Bailey's Arithmetic with FMA Instruction: Its Performance and Applications.

When a floating-point arithmetic is executed on a processor unit, round-off and truncation errors occur every calculation. These errors cause a precision issue in a large simulation which requires a great number of calculations. Therefore, we have developed the quadruple-precision basic linear algebra subprograms (QPBLAS) based on Baileyu0027s double-double arithmetic. The multiplication operation of Baileyu0027s arithmetic is realized by 24 double-precision operations. When using an FMA (fused multiply-add) instruction, we can reduce the number of operations in about half. Therefore, we develop QPBLAS using the FMA instruction and evaluate its performance. The result shows that the QPBLAS using FMA is basically faster than QPBLAS. Moreover, when we replace QPBLAS of the quadruple-precision eigenvalue solver QPEigenK with QPBLAS using FMA, we can obtain about 1020% speedup.