September 16, 2026 Dataset Open Access

Reproducibility package: Pressure-robustness by commuting interpolation operators for Stokes discretizations with continuous pressures

Philip Lederer

This repository contains the NGSolve/Netgen finite element code and raw convergence data used to produce the numerical examples (Section 8) of "Pressure-robustness by commuting interpolation operators for Stokes discretizations with continuous pressures" by P. L. Lederer and T. Vock. It includes two Python scripts (`convergence.py`, `convergence_3d.py`) that solve the Stokes equations on the unit square (resp. unit cube) with a manufactured smooth solution, using the MINI, stabilized equal-order $P^kP^k$, and Taylor--Hood finite elements, both with the standard method and with the pressure-robust method proposed in the paper (in which the right-hand side force is replaced by its Nédélec interpolant). It also includes the raw error/convergence data (`conv_data/`) generated by these scripts, from which all convergence tables in the paper's numerics section are directly computed. See README.md for software requirements, instructions to reproduce the data, and a file-to-table mapping.