Naslov
Parallel and heterogeneous m-Hessenberg-triangular-triangular reduction

Autori
Bosner, Nela ; Karlsson, Lars

Izvornik
SIAM journal on scientific computing (1064-8275) 39 (2017), 1; C29-C47

Vrsta, podvrsta i kategorija rada
Radovi u časopisima, članak, znanstveni

Ključne riječi
m-Hessenberg-triangular-triangular form, solving shifted linear systems, Givens rotations, heterogeneous CPU/GPU implementation

Sažetak
The m-Hessenberg-triangular-triangular (mHTT) reduction is a simultaneous orthogonal reduction of three matrices to condensed form. It has applications, for example, in solving shifted linear systems arising in various control theory problems. A new heterogeneous CPU/GPU implementation of the mHTT reduction is presented and evaluated against an existing CPU implementation. The algorithm offloads the compute-intensive matrix-matrix multiplications to the GPU and keeps the inner loop, which is memory intensive and has a complicated control flow, on the CPU. Experiments demonstrate that the heterogeneous implementation can be superior to the existing CPU implementation on a system with 2x8 CPU cores and one GPU. Future development should focus on improving the scalability of the CPU computations.

Izvorni jezik
Engleski

Znanstvena područja
Matematika, Računarstvo

POVEZANOST RADA