Naslov
A hybrid algorithm for solving shifted linear systems

Autori
Bujanović, Zvonimir

Vrsta, podvrsta i kategorija rada
Sažeci sa skupova, sažetak, znanstveni

Izvornik
6th Croatian Mathematical Congress / - , 2016, 1-1

Skup
6th Croatian Mathematical Congress

Mjesto i datum
Zagreb, Hrvatska, 14.06.2016. - 17.06.2016

Vrsta sudjelovanja
Predavanje

Vrsta recenzije
Nije recenziran

Ključne riječi
GPU ; Hessenberg matrix ; parallel solver ; shifted linear systems ; transfer function

Sažetak
In this talk, we propose a new algorithm for solving dense shifted linear systems with multiple right-hand sides and a large number of shifts. Such problems arise e.g. in control theory when computing the frequency response of a LTI system, and in many other applications as well. The new algorithm is designed for hybrid computer architectures that use classical multicore processors in combination with GPU accelerators with Nvidia CUDA technology. The algorithm consists of two phases: first, a hybrid CPU-GPU routine is used in order to transform the linear system to a so-called controller--Hessenberg form. This reduction is done only once, regardless of the number of shifts, and allows us to later solve the systems with far less computational effort. In the second phase, the transformed systems are solved by means of a pure GPU algorithm, simultaneously for a large batch of shifts. The solver combines custom made highly parallel kernels with the efficient cuBLAS routines. Such distribution of computational load shows significant performance benefits compared to classical CPU-bound algorithms, which we demonstrate by numerical experiments. This is joint work with Nela Bosner and Zlatko Drmač.

Izvorni jezik
Engleski

Znanstvena područja
Matematika

POVEZANOST RADA