Evaluating Accuracy and Efficiency of HPC Solvers for Sparse Linear Systems with Applications to PDEs

by   Antonella Galizia, et al.

Partial Differential Equations (PDEs) describe several problems relevant to many fields of applied sciences, and their discrete counterparts typically involve the solution of sparse linear systems. In this context, we focus on the analysis of the computational aspects related to the solution of large and sparse linear systems with HPC solvers, by considering the performances of direct and iterative solvers in terms of computational efficiency, scalability, and numerical accuracy. Our aim is to identify the main criteria to support application-domain specialists in the selection of the most suitable solvers, according to the application requirements and available resources. To this end, we discuss how the numerical solver is affected by the regular/irregular discretisation of the input domain, the discretisation of the input PDE with piecewise linear or polynomial basis functions, which generally result in a higher/lower sparsity of the coefficient matrix, and the choice of different initial conditions, which are associated with linear systems with multiple right-hand side terms. Finally, our analysis is independent of the characteristics of the underlying computational architectures, and provides a methodological approach that can be applied to different classes of PDEs or with approximation problems.


page 3

page 17


A Hybrid Iterative Numerical Transferable Solver (HINTS) for PDEs Based on Deep Operator Network and Relaxation Methods

Iterative solvers of linear systems are a key component for the numerica...

Learning Preconditioner for Conjugate Gradient PDE Solvers

Efficient numerical solvers for partial differential equations empower s...

Comparison of SeDuMi and SDPT3 Solvers for Stability of Continuous-time Linear System

SeDuMi and SDPT3 are two solvers for solving Semi-definite Programming (...

Modeling excitable cells with the EMI equations: spectral analysis and iterative solution strategy

In this work, we are interested in solving large linear systems stemming...

Chronos: A general purpose classical AMG solver for High Performance Computing

The numerical simulation of the physical systems has become in recent ye...

Natural Factor Based Solvers

We consider parametric families of partial differential equations–PDEs w...

Compressive Isogeometric Analysis

This work is motivated by the difficulty in assembling the Galerkin matr...

Please sign up or login with your details

Forgot password? Click here to reset