Linear operators for discretizations of differential equations and scientific machine learning (SciML)

ExaDG - High-Order Discontinuous Galerkin for the Exa-Scale

Modular and scalable computational imaging in Python with GPU/out-of-core computing.

PyTorch implementation of the Hessian-free optimizer

Stability, bifurcation, and sensitivity analysis for Nek5000

MFMG is an open-source library implementing matrix-free multigrid methods.

Second-Order Methods in JAX

Source and reproducibility code for the paper 'No 3D Matrices: A Unified Tensor-Product View of Matrix-Free Cartesian PDE Solvers' (Bay & Yearick)

One matrix-free spectral primitive for an operator's whole spectral world — densities, moments, free convolution, pseudospectra, the Brown measure, dequantization diagnostics. Spans 18+ fields of math & physics. pip install resona

Matrix-free implementation of accelerated projected gradient descent for large-scale non-negative least-squares problems.

C++, R and Python package for the Synthetic MR imaging program symr. For the R and python programs, see the R and python branches respectively.

Practical session on implementing probabilistic linear solvers at the Probabilistic Numerics Spring School 2024

Matrix-free 3D SIMP topology optimization with fused gather-GEMM-scatter CUDA kernels on NVIDIA RTX 4090. Companion code for arXiv:2604.18020.

Fully implicit density-based solver for compressible viscous flows on unstructured meshes. The project focuses on matrix-free approaches to implicit solver implementation

Jeremy L Thompson PhD dissertation

Iterative matrix-free eigensolver (Davidson-Liu algorithm) with ndarray

A library to build up lazily evaluated expressions of linear transforms for efficient scientific computing. To cite this Original Software Publication: https://www.sciencedirect.com/science/article/pii/S2352711022000218

GPU Library for Accelerated Solvers

Do we really need quantum computing? A runnable essay: one matrix-free dial — the effective rank Φ₁ — flags which quantum speedups survive a classical attack, and which were never quantum at all.

Quantum hamiltonian simulation with matrix-free linear operators