Abstract:Computational singular perturbation (CSP) is a method to analyze dynamical systems. It targets the decoupling of fast and slow dynamics using an alternate linear expansion of the right-hand side of the governing equations based on eigenanalysis of the associated Jacobian matrix. This representation facilitates diagnostic analysis, detection and control of stiffness, and the development of simplified models. We have implemented CSP in a C++ open-source library CSPlib 1 using the Kokkos 2 parallel programming model to address portability across diverse heterogeneous computing platforms, i.e., multi/many-core CPUs and GPUs. We describe the CSPlib implementation and present its computational performance across different computing platforms using several test problems. Specifically, we test the CSPlib performance for a constant pressure ignition reactor model on different architectures, including IBM Power 9, Intel Xeon Skylake, and NVIDIA V100 GPU. The size of the chemical kinetic mechanism is varied in these tests. As expected, the Jacobian matrix evaluation, the eigensolution of the Jacobian matrix, and matrix inversion are the most expensive computational tasks. When considering the higher throughput characteristic of GPUs, GPUs performs better for small matrices with higher occupancy rate. CPUs gain more advantages from the higher performance of well-tuned and optimized linear algebra libraries such as OpenBLAS. Program Program Title: CSPlib CPC Library link to program files: https://doi.org/10.17632/p9gb7z54sp.1 Developer's repository link: https://github.com/sandialabs/csplib Licensing provisions: BSD 2-clause Programming language: C++ Nature of problem: Dynamical systems can involve coupled processes with a wide range of time scales. The computational singular perturbation (CSP) method offers a reformulation of these systems which enables the use of dynamically-based diagnostic tools to better comprehend the dynamics by decoupling fast and slow processes. CSPlib is an open-source software library for analyzing general ordinary differential equation (ODE) and differential algebraic equation (DAE) systems, with specialized implementations for detailed chemical kinetic ODE/DAE systems. It relies on CSP for the analysis of these systems. CSPlib has been used in gas kinetic and heterogeneous catalytic kinetic models. Solution method: CSP analysis seeks a set of basis vectors to linearly decompose the right-hand side (RHS) of a dynamical system in a manner that decouples fast and slow processes. The CSP basis vectors are often well approximated with the right eigenvectors of the RHS Jacobian. And the left basis vectors are found by the inversion of the matrix, whose columns are the CSP basis vectors. Accordingly, the right and left CSP basis vectors are orthonormal. CSP defines mode amplitudes as the projections of the left basis vectors on the RHS; the time scales as the reciprocals of the RHS Jacobian eigenvalue magnitudes; and the CSP pointers, which are the element-wise multiplication of the transpose of the right CSP basis vectors with the left CSP basis vectors. For kinetic models that can be cast as the product of a generalized stoichiometric matrix and a rate of progress vector, CSP defines the participation index, which represents the contribution of a chemical reaction to each mode. Further, it defines the slow and fast importance indices, which describe the contribution of a chemical reaction to the slow and fast dynamics of a state variable, respectively. These indices are useful in diagnostic studies of dynamical systems and the construction of simplified models. Additional comments including restrictions and unusual features: CSPlib is a portable library that carries out many CSP analyses in parallel and can be used in modern high-performance platforms.

PyCPL: The ESO Common Pipeline Library in Python v1.0

binary_c-python: A Python-based stellar population synthesis tool and interface to binary_c

CSPlib: A performance portable parallel software toolkit for analyzing complex kinetic mechanisms

PyPDB: a Python API for the Protein Data Bank

rstoolbox - a Python library for large-scale analysis of computational protein design data and structural bioinformatics

PDBeCIF: an open-source mmCIF/CIF parsing and processing package

Post-processing CHARIS integral field spectrograph data with PyKLIP

PyMS: a Python toolkit for processing of gas chromatography-mass spectrometry (GC-MS) data. Application and comparative study of selected tools

PaPy: Parallel and Distributed Data-processing Pipelines in Python

The Public Photometry Pipelines for Exoplanets

prose: A Python framework for modular astronomical images processing

Parsl+CWL: Towards Combining the Python and CWL Ecosystems

PLCOjs, a FAIR GWAS web SDK for the NCI Prostate, Lung, Colorectal, and Ovarian Cancer Genetic Atlas Project

PHOTOMETRYPIPELINE: An Automated Pipeline for Calibrated Photometry

Pyleoclim: Paleoclimate Timeseries Analysis and Visualization with Python

opseestools: A Python library to streamline OpenSeesPy workflows

NMRPy: a novel NMR scripting system to implement artificial intelligence and advanced applications

pyEQUIB Python Package, an addendum to proEQUIB: IDL Library for Plasma Diagnostics and Abundance Analysis

Pyaneti: a fast and powerful software suite for multi-planet radial velocity and transit fitting

GalaPy: A highly optimised C++/Python spectral modelling tool for galaxies. I. Library presentation and photometric fitting

PyCPET - Computing Heterogeneous 3D Protein Electric Fields and Their Dynamics