A survey on direct solvers for Galerkin methods
MetadataShow full item record
In this paper we describe the history, performance, and design concepts of direct solvers for algebraic systems resulting from Galerkin discretizations of partial differential equations. Popular direct solver implementations of Gaussian elimination (also known as LU factorization) are introduced and briefly analyzed. We discuss three of the most relevant aspects influencing the performance of direct solvers on this kind of algebraic systems. First, the ordering of the degrees of freedom of the algebraic system has a significant impact on the solver performance, solution speed and memory requirements. The impact of unknowns ordering for elimination is exemplified and alternative ordering algorithms are described and compared. Second, the effect of round-off error on the simulation results is discussed. We detail this effect for uniform grids where the impact of round-off error on the solution is controlled by the condition number of the matrix in terms of the element size, but is independent of the polynomial order of approximation. Additionally, we discuss the link between unknown ordering and round-off error. Third, we describe the impact of the connectivity pattern (graph) of the basis functions on the performance of direct solvers. Variations in the connectivity structure of the resulting discrete system have severe impact on performance of the solver. That is, the resources needed to factorize the system strongly depend on its connectivity graph. Less connected graphs are cheaper to solve, that is, C0 finite element discretizations are cheaper to solve with direct solvers than Cp−1 discretizations.
Showing items related by title, author, creator and subject.
Garcia, D.; Pardo, D.; Dalcin, L.; Paszynski, M.; Collier, N.; Calo, Victor (2016)© 2016 Elsevier B.V.We propose the use of highly continuous finite element spaces interconnected with low continuity hyperplanes to maximize the performance of direct solvers. Starting from a highly continuous Isogeometric ...
Rodriguez Bernabeu, S.; Puzyrev, Vladimir; Hanzich, M.; Fernández, S. (2016)Frequency-domain seismic and electromagnetic modeling requires solving the linear systems resulting from the discretization of the corresponding time-harmonic equations. Geophysical inversion is typically performed using ...
Garcia, D.; Pardo, D.; Dalcin, L.; Calo, Victor (2018)© 2018 Elsevier B.V. Starting from a highly continuous Isogeometric Analysis (IGA) discretization, refined Isogeometric Analysis (rIGA) introduces C 0 hyperplanes that act as separators for the direct LU factorization ...