We develop a GPU accelerated variant of the blocked Bartels–Stewart algorithm for the generalized Lyapunov equation.
We analyze the influence of optimized data alignment on the GPU and discuss the mathematical problem behind it.
Reordering eigenvalues avoids alignment problems on the GPU and accelerates the computation.