Abstract
Matrix multiplication is a very important computation kernel both in its own right as a building block of many scientific applications and as a popular representative for other scientific applications. Cannon's algorithm which dates back to 1969 was the first efficient algorithm for parallel matrix multiplication providing theoretically optimal communication cost. However this algorithm requires a square number of processors. In the mid-1990s, the SUMMA algorithm was introduced. SUMMA overcomes the shortcomings of Cannon's algorithm as it can be used on a nonsquare number of processors as well. Since then the number of processors in HPC platforms has increased by two orders of magnitude making the contribution of communication in the overall execution time more significant. Therefore, the state of the art parallel matrix multiplication algorithms should be revisited to reduce the communication cost further. This paper introduces a new parallel matrix multiplication algorithm, Hierarchical SUMMA (HSUMMA), which is a redesign of SUMMA. Our algorithm reduces the communication cost of SUMMA by introducing a two-level virtual hierarchy into the two-dimensional arrangement of processors. Experiments on an IBM BlueGene/P demonstrate the reduction of communication cost up to 2.08 times on 2048 cores and up to 5.89 times on 16384 cores. © 2013 IEEE.
Original language | English (US) |
---|---|
Title of host publication | 2013 42nd International Conference on Parallel Processing |
Publisher | Institute of Electrical and Electronics Engineers (IEEE) |
Pages | 754-762 |
Number of pages | 9 |
ISBN (Print) | 9780769551173 |
DOIs | |
State | Published - Oct 2013 |
Externally published | Yes |
Bibliographical note
KAUST Repository Item: Exported on 2020-10-01Acknowledgements: The research in this paper was supported by IRCSET(IrishResearch Council for Science, Engineering and Technol-ogy) and IBM, grant numbers EPSG/2011/188 and EP-SPD/2011/207.Some of the experiments presented in this paper werecarried out using the Grid’5000 experimental testbed, beingdeveloped under the INRIA ALADDIN development actionwith support from CNRS, RENATER and several Universitiesas well as other funding bodies (see https://www.grid5000.fr)Another part of the experiments in this research were carriedout using the resources of the Supercomputing Laboratory atKing Abdullah University of Science&Technology (KAUST)in Thuwal, Saudi Arabia.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.