2022/02/18 by Denis Demidov, Demidov, Denis
Computer Science · Earth and Planetary Sciences · Engineering · #35-04 #65-04 #65Y05 #65Y10 #65Y15 #97N80 #Advanced Data Processing Techniques #Distributed and Parallel Computing Systems #FOS: Computer and information sciences #FOS: Mathematics #Geological Modeling and Analysis #Mathematical Software (cs.MS) #Numerical Analysis (math.NA)
paper · pdf · doi:10.48550/arxiv.2202.09056
openalex publication_date 2022/02/18 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
Efficient solution of 3D elasticity problems is an important part of many industrial and scientific applications. Smoothed aggregation algebraic multigrid using rigid body modes for the tentative prolongation operator construction is an efficient and robust choice for the solution of linear systems arising from the discretization of elasticity equations. The system matrices on every level of the multigrid hierarchy have block structure, so using block representation and block arithmetics should significantly improve the solver efficiency. However, the tentative prolongation operator construction may only be done using scalar representation. The paper proposes a couple of practical approaches for enabling the use of block arithmetics with smoothed aggregation algebraic multigrid based on the open-source AMGCL library. It is shown on the example of two real-world model problems that the suggested improvements may speed up the solution by 50% and reduce the memory requirements for the preconditioner by 30%. The implementation is straightforward and only requires a minimal amount of code.