Resolvent4py: a parallel Python package for analysis, model reduction and control of large-scale linear systems

TL;DR

Resolvent4py is a Python package that enables efficient analysis, model reduction, and control of large-scale linear systems using parallel computing, suitable for various scientific fields.

Contribution

It introduces a user-friendly Python library that combines parallelism with advanced linear system analysis capabilities for large-scale problems.

Findings

Supports systems with millions or billions of degrees of freedom.

Enables efficient parallel computations in fluid, solid mechanics, and other fields.

Provides a Python-like interface with MPI-based parallelism.

Abstract

In this paper, we present resolvent4py, a parallel Python package for the analysis, model reduction and control of large-scale linear systems with millions or billions of degrees of freedom. This package provides the user with a friendly Python-like experience (akin to that of well-established libraries such as numpy and scipy), while enabling MPI-based parallelism through mpi4py, petsc4py and slepc4py. In turn, this allows for the development of streamlined and efficient Python code that can be used to solve several problems in fluid mechanics, solid mechanics, graph theory, molecular dynamics and several other fields.