Papers › CoLA: Exploiting Compositional Structure for Automatic and Efficient Numerical Linear Algebra

CoLA: Exploiting Compositional Structure for Automatic and Efficient Numerical Linear Algebra

6 Sep 2023NeurIPS 2023 11arXiv:2309.03060archive 2025-07-28

Andres Potapczynski, Marc Finzi, Geoff Pleiss, Andrew Gordon Wilson

Many areas of machine learning and science involve large linear algebra problems, such as eigendecompositions, solving linear systems, computing matrix exponentials, and trace estimation. The matrices involved often have Kronecker, convolutional, block diagonal, sum, or product structure. In this paper, we propose a simple but general framework for large-scale linear algebra problems in machine learning, named CoLA (Compositional Linear Algebra). By combining a linear operator abstraction with compositional dispatch rules, CoLA automatically constructs memory and runtime efficient numerical algorithms. Moreover, CoLA provides memory efficient automatic differentiation, low precision computation, and GPU acceleration in both JAX and PyTorch, while also accommodating new objects, operations, and rules in downstream packages via multiple dispatch. CoLA can accelerate many algebraic operations, while making it easy to prototype matrix structures and algorithms, providing an appealing drop-in tool for virtually any computational effort that requires linear algebra. We showcase its efficacy across a broad range of applications, including partial differential equations, Gaussian processes, equivariant model construction, and unsupervised learning.

PaperPDFConference PDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2309.03060")

Code

Syntology Ran 17 of 17 code samples harvested from 1 repository linked to this paper; 0 have no recorded run. Of those that ran: 17 ran with no contract checked.

By repository: official repository: 17 samples from 1 repository, 17 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

wilson-labs/cola officialmentioned in papermentioned on GitHubjaxApache-2.0 report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

17 samples harvested; 17 ran; 0 honoured the contract we drafted; 0 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

17ran

Licence: 0 of the 17 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from wilson-labs/cola. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

build_tqdm wilson-labs/cola/cola/utils/jax_tqdm.py official repository ran Apache-2.0 (permissive) · 1222b814438ac04e · report
check_valid_dtype wilson-labs/cola/cola/backends/backends.py official repository ran Apache-2.0 (permissive) · 1505bc3ca64e8748 · report
combine wilson-labs/cola/cola/utils/custom_autodiff.py official repository ran Apache-2.0 (permissive) · e09ed6dad6d629b9 · report
eig wilson-labs/cola/cola/backends/np_fns.py official repository ran fingerprinted Apache-2.0 (permissive) · 194698227eb24ab9 · report
is_xnp_array wilson-labs/cola/cola/ops/operator_base.py official repository ran Apache-2.0 (permissive) · 79232c17bd3e80b9 · report
iterative_autograd wilson-labs/cola/cola/utils/custom_autodiff.py official repository ran Apache-2.0 (permissive) · 07094ab2794b77be · report
loop_tqdm wilson-labs/cola/cola/utils/jax_tqdm.py official repository ran Apache-2.0 (permissive) · d551c022a49dc59a · report
lstsq wilson-labs/cola/cola/backends/np_fns.py official repository ran fingerprinted Apache-2.0 (permissive) · 5dc26cca6fae567f · report
maybe_get_dtype wilson-labs/cola/cola/ops/operator_base.py official repository ran Apache-2.0 (permissive) · 12064a51bf1dd610 · report
product wilson-labs/cola/cola/ops/operators.py official repository ran fingerprinted Apache-2.0 (permissive) · 18ed40676168eff9 · report
scan_tqdm wilson-labs/cola/cola/utils/jax_tqdm.py official repository ran Apache-2.0 (permissive) · 06b12e67297b8604 · report
sparse_csr wilson-labs/cola/cola/backends/jax_fns.py official repository ran Apache-2.0 (permissive) · c56a560ef1f1c286 · report
sparse_csr wilson-labs/cola/cola/backends/torch_fns.py official repository ran Apache-2.0 (permissive) · 624ef6c3eeae7741 · report
svd wilson-labs/cola/cola/backends/jax_fns.py official repository ran Apache-2.0 (permissive) · a9583780ceda1e46 · report
to_np wilson-labs/cola/cola/backends/jax_fns.py official repository ran Apache-2.0 (permissive) · 0381ba678b73491b · report
while_loop wilson-labs/cola/cola/utils/torch_tqdm.py official repository ran Apache-2.0 (permissive) · 200e115745cfc25b · report
while_loop_winfo wilson-labs/cola/cola/utils/torch_tqdm.py official repository ran Apache-2.0 (permissive) · f71a8d09e632c999 · report

Tasks

Gaussian Processes

2 archive task tags without a task page not shown.

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Methods

COLA

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections