Code for the paper: Incomplete double-cone factorizations of centrosymmetric matrices arising in spectralmethods
Authors: Chen Greif, Sarah Nataj, Manfred Trummer
Subjects: Numerical linear algebra, Spectral methods for PDEs
Pdf link: https://link.springer.com/article/10.1007/s11075-023-01612-y
Abstract: We develop structure-preserving incomplete LU type factorizations for preconditioning centrosymmetric matrices and use them to numerically solve centrosymmetric and nearly centrosymmetric linear systems arising from spectral methods for partial differential equations. Our algorithm builds in part on direct solution techniques previously developed for this type of linear systems, featuring double-cone factorizations. We illustrate our findings on discretizations of model problems involving the Poisson, diffusion, Helmholtz, and biharmonic equations in one, two, and three dimensions.