Journal of Computational Physics
ACM Transactions on Mathematical Software (TOMS)
Journal of Computational Physics
Towards a scalable fully-implicit fully-coupled resistive MHD formulation with stabilized FE methods
Journal of Computational Physics
Towards Adaptive Smoothed Aggregation ($\alpha$SA) for Nonsymmetric Problems
SIAM Journal on Scientific Computing
Stabilization and scalable block preconditioning for the Navier-Stokes equations
Journal of Computational Physics
A General Interpolation Strategy for Algebraic Multigrid Using Energy Minimization
SIAM Journal on Scientific Computing
On-the-Fly Adaptive Smoothed Aggregation Multigrid for Markov Chains
SIAM Journal on Scientific Computing
LSSC'09 Proceedings of the 7th international conference on Large-Scale Scientific Computing
Additive operator decomposition and optimization–based reconnection with applications
LSSC'09 Proceedings of the 7th international conference on Large-Scale Scientific Computing
Performance analysis of parallel Schwarz preconditioners in the LES of turbulent channel flows
Computers & Mathematics with Applications
A generalized Block FSAI preconditioner for nonsymmetric linear systems
Journal of Computational and Applied Mathematics
| Hi-index | 0.02 |
We propose a new variant of smoothed aggregation (SA) suitable for nonsymmetric linear systems. The new algorithm is based on two key generalizations of SA: restriction smoothing and local damping. Restriction smoothing refers to the smoothing of a tentative restriction operator via a damped Jacobi-like iteration. Restriction smoothing is analogous to prolongator smoothing in standard SA and in fact has the same form as the transpose of prolongator smoothing when the matrix is symmetric. Local damping refers to damping parameters used in the Jacobi-like iteration. In standard SA, a single damping parameter is computed via an eigenvalue computation. Here, local damping parameters are computed by considering the minimization of an energy-like quantity for each individual grid transfer basis function. Numerical results are given showing how this method performs on highly nonsymmetric systems.