Resolvent4py

Resolvent4py is a petsc4py- and slepc4py-based toolbox to perform analysis, model reduction and control of high-dimensional linear dynamical systems. The goal of this project is to provide users with a friendly python-like experience, while also leveraging the high-performance and parallel-computing capabilities of the PETSc and SLEPc library. Current functionalities include:

  • Right and left eigendecomposition

  • Resolvent analysis (algebraic and using time-stepping techniques)

  • Harmonic resolvent analysis

  • Balanced truncation (time-invariant)

If you use resolvent4py in your work, please cite the following arXiv preprint (see here for the pdf)

@misc{padovan2025resolvent4pyparallelpythonpackage,
   title={Resolvent4py: a parallel Python package for analysis, model reduction and control of large-scale linear systems},
   author={Alberto Padovan and Vishal Anantharaman and Clarence W. Rowley and Blaine Vollmer and Tim Colonius and Daniel J. Bodony},
   year={2025},
   eprint={2506.20539},
   archivePrefix={arXiv},
   primaryClass={physics.comp-ph},
   url={https://arxiv.org/abs/2506.20539},
}

Getting Started

Examples