umu.sePublications
Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
On the Tunability of a New Hessenberg Reduction Algorithm Using Parallel Cache Assignment
Umeå University, Faculty of Science and Technology, Department of Computing Science. (UMIT)
Umeå University, Faculty of Science and Technology, Department of Computing Science. (UMIT)
Umeå University, Faculty of Science and Technology, Department of Computing Science. (UMIT)
2017 (English)Conference paper, Published paper (Refereed)
Place, publisher, year, edition, pages
2017.
National Category
Computer Sciences Computational Mathematics
Identifiers
URN: urn:nbn:se:umu:diva-145342OAI: oai:DiVA.org:umu-145342DiVA, id: diva2:1186599
Conference
12th International Conference on Parallel Processing and Applied Mathematics (PPAM 2017)
Note

Accepted for publishing

Available from: 2018-02-28 Created: 2018-02-28 Last updated: 2018-06-09
In thesis
1. Identification and tuning of algorithmic parameters in parallel matrix computations: Hessenberg reduction and tensor storage format conversion
Open this publication in new window or tab >>Identification and tuning of algorithmic parameters in parallel matrix computations: Hessenberg reduction and tensor storage format conversion
2018 (English)Licentiate thesis, comprehensive summary (Other academic)
Abstract [en]

This thesis considers two problems in numerical linear algebra and high performance computing (HPC): (i) the parallelization of a new blocked Hessenberg reduction algorithm using Parallel Cache Assignment (PCA) and the tunability of its algorithm parameters, and (ii) storing and manipulating dense tensors on shared memory HPC systems.

The Hessenberg reduction appears in the Aggressive Early Deflation (AED) process for identifying converged eigenvalues in the distributed multishift QR algorithm (state-of-the-art algorithm for computing all eigenvalues for dense square matrices). Since the AED process becomes a parallel bottleneck it motivates a further study of AED components. We present a new Hessenberg reduction algorithm based on PCA which is NUMA-aware and targeting relatively small problem sizes on shared memory systems. The tunability of the algorithm parameters are investigated. A simple off-line tuning is presented and the performance of the new Hessenberg reduction algorithm is compared to its counterparts from LAPACK and ScaLAPACK. The new algorithm outperforms LAPACK in all tested cases and outperforms ScaLAPACK in problems smaller than order 1500, which are common problem sizes for AED in the context of the distributed multishift QR algorithm.

We also investigate automatic tuning of the algorithm parameters. The parameters span a huge search space and it is impractical to tune them using standard auto-tuning and optimization techniques. We present a modular auto-tuning framework which applies: search space decomposition, binning, and multi-stage search to enable searching the huge search space efficiently. The framework using these techniques exposes the underlying subproblems which allows using standard auto-tuning methods to tune them. In addition, the framework defines an abstract interface, which combined with its modular design, allows testing various tuning algorithms.

In the last part of the thesis, the focus is on the problem of storing and manipulating dense tensors. Developing open source tensor algorithms and applications is hard due to the lack of open source software for fundamental tensor operations. We present a software library dten, which includes tools for storing dense tensors in shared memory and converting a tensor storage format from one canonical form to another. The library provides two different ways to perform the conversion in parallel, in-place and out-of-place. The conversion involves moving blocks of contiguous data and are done to maximize the size of the blocks to move. In addition, the library supports tensor matricization for one or two tensors at the same time. The latter case is important in preparing tensors for contraction operations. The library is general purpose and highly flexible.

Place, publisher, year, edition, pages
Umeå: Umeå universitet, 2018. p. 15
Series
Report / UMINF, ISSN 0348-0542 ; UMINF 18.22
National Category
Computer Sciences
Identifiers
urn:nbn:se:umu:diva-145345 (URN)978-91-7601-843-9 (ISBN)
Supervisors
Available from: 2018-03-01 Created: 2018-02-28 Last updated: 2018-06-09Bibliographically approved

Open Access in DiVA

No full text in DiVA

Authority records BETA

Eljammaly, MahmoudKarlsson, LarsKågström, Bo

Search in DiVA

By author/editor
Eljammaly, MahmoudKarlsson, LarsKågström, Bo
By organisation
Department of Computing Science
Computer SciencesComputational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

urn-nbn

Altmetric score

urn-nbn
Total: 46 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf