A High Performance QDWH-SVD Solver using Hardware Accelerators

Sukkari, Dalal E.; Ltaief, Hatem; Keyes, David E.

A High Performance QDWH-SVD Solver using Hardware Accelerators

Files

qdwh-svd-submitted.pdf (920.83 KB)

License

http://creativecommons.org/licenses/by/4.0/

Type

Technical Report

Authors

Sukkari, Dalal E.
Ltaief, Hatem

Keyes, David E.

KAUST Department

Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division

Online Publication Date

2016-08-15

Print Publication Date

2016-08-13

Date

2016-08-15

Abstract

This paper describes a new high performance implementation of the QR-based Dynamically Weighted Halley Singular Value Decomposition (QDWH-SVD) solver on multicore architecture enhanced with multiple GPUs. The standard QDWH-SVD algorithm was introduced by Nakatsukasa and Higham (SIAM SISC, 2013) and combines three successive computational stages: (1) the polar decomposition calculation of the original matrix using the QDWH algorithm, (2) the symmetric eigendecomposition of the resulting polar factor to obtain the singular values and the right singular vectors and (3) the matrix-matrix multiplication to get the associated left singular vectors. A comprehensive test suite highlights the numerical robustness of the QDWH-SVD solver. Although it performs up to two times more flops when computing all singular vectors compared to the standard SVD solver algorithm, our new high performance implementation on single GPU results in up to 3.8x improvements for asymptotic matrix sizes, compared to the equivalent routines from existing state-of-the-art open-source and commercial libraries. However, when only singular values are needed, QDWH-SVD is penalized by performing up to 14 times more flops. The singular value only implementation of QDWH-SVD on single GPU can still run up to 18% faster than the best existing equivalent routines. Integrating mixed precision techniques in the solver can additionally provide up to 40% improvement at the price of losing few digits of accuracy, compared to the full double precision floating point arithmetic. We further leverage the single GPU QDWH-SVD implementation by introducing the first multi-GPU SVD solver to study the scalability of the QDWH-SVD framework.

Citation

Sukkari, D., Ltaief, H., & Keyes, D. (2016). A High Performance QDWH-SVD Solver Using Hardware Accelerators. ACM Transactions on Mathematical Software, 43(1), 1–25. doi:10.1145/2894747

Acknowledgements

This work was supported by the Extreme Computing Research Center at KAUST. The authors would like to thank Ahmad Abdelfattah for his help to integrate KBLAS into QDWH-SVD and NVIDIA for the hardware donations.

Version History

You are currently viewing version 2 of the item.

Now showing 1 - 2 of 2

Version	Date	Summary
2*	2018-04-30 13:20:53	Published with DOI
1	2015-04-08 12:08:12

* Selected version