SIMD vectorization for simultaneous solution of locally varying linear systems with multiple right-hand sides
Martin J. Kühn (Deutsches Zentrum für Luft- und Raumfahrt (DLR))
Johannes Holke (Deutsches Zentrum für Luft- und Raumfahrt (DLR))
Annette Lutz (Technische Universität Darmstadt)
Jonas Thies (TU Delft - Electrical Engineering, Mathematics and Computer Science)
Melven Röhrig-Zöllner (Deutsches Zentrum für Luft- und Raumfahrt (DLR))
Alexander Bleh (Deutsches Zentrum für Luft- und Raumfahrt (DLR))
Jan Backhaus (Deutsches Zentrum für Luft- und Raumfahrt (DLR))
Achim Basermann (Deutsches Zentrum für Luft- und Raumfahrt (DLR))
More Info
expand_more
Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.
Abstract
Developments in numerical simulation of flows and high-performance computing influence one another. More detailed simulation methods create a permanent need for more computational power, while new hardware developments often require changes to the software to exploit new hardware features. This dependency is very pronounced in the case of vector-units which are featured by all modern processors to increase their numerical throughput but require vectorization of the software to be used efficiently. We study the vectorization of a simulation method that exhibits an inherent level of vector-parallelism. This is of particular interest as SIMD operations will hopefully be available with std::simd in a future C++ standard. The simulation method considered here results in the simultaneous solution of multiple sparse linear systems of equations which only differ by their main diagonal and right-hand sides. Such structure arises in the simulation of unsteady flow in turbomachinery by means of a frequency domain approach called harmonic balance.
help_outline