English

Parallel matching-based AMG preconditioners for elliptic equations discretized by IgA

Numerical Analysis 2026-05-01 v2 Numerical Analysis

Abstract

Isogeometric analysis (IgA) offers enhanced approximation capabilities for the discretization of elliptic boundary-value problems, yet it results in large, sparse, and increasingly ill-conditioned linear systems due to higher interconnectivity among degrees of freedom. In particular, the discretization with tensor-product B-splines or NURBS of degree pp on a mesh with nn elements per parametric direction leads to symmetric positive-definite systems of the form Ku=FK\mathbf{u} = \mathbf{F}, where the matrix bandwidth and condition number scale unfavorably with both pp and spatial dimension dd. To address the computational challenges posed by such systems, especially in three-dimensional or high-order scenarios, Krylov subspace methods with specialized preconditioners become essential. This paper investigates the efficacy of algebraic multigrid (AMG) preconditioners tailored for IgA-based discretizations, with a focus on performance in modern high-performance computing (HPC) environments. Leveraging the Parallel Sparse Computation Toolkit (PSCToolkit), we explore distributed-memory and GPU-accelerated strategies for solving large-scale problems. The study assesses algorithmic efficiency and scalability across a range of benchmark tests. The results demonstrate that AMG preconditioners can achieve robust and scalable performance, confirming their potential as practical solvers for large IgA systems in engineering and scientific applications.

Keywords

Cite

@article{arxiv.2511.21268,
  title  = {Parallel matching-based AMG preconditioners for elliptic equations discretized by IgA},
  author = {Pasqua D'Ambra and Fabio Durastante and Salvatore Filippone},
  journal= {arXiv preprint arXiv:2511.21268},
  year   = {2026}
}
R2 v1 2026-07-01T07:55:58.599Z