Papers › Multi-output Gaussian process prediction of physical fields under linear equality constraints

Multi-output Gaussian process prediction of physical fields under linear equality constraints

26 Aug 2026arXiv:2608.25709added by Syntology

Mahamat Hamdan Nassouradine, Clément Gauchy, Pierre-Emmanuel Angeli, Sébastien da Veiga

Title, abstract, authors and date from arXiv's metadata (CC0); this paper is not in the Papers with Code archive (frozen 2025-07-28).

We address the simultaneous prediction of multiple high-dimensional physical fields governed by linear equality constraints, a setting that arises in many real-world applications in physics machine learning. Gaussian process (GP) regression is a widely used surrogate modeling approach due to its effectiveness in small-sample regimes and its ability to provide uncertainty quantification. However, applying GP models in this setting raises two major challenges: the high dimensionality of the discretized output fields and the enforcement of the physical constraint in predictions. For the latter, a common strategy consists in deducing one output from the others via the constraint relation. Through a benchmark, we show that this deductive approach is sensitive to the arbitrary choice of which output to deduce, affecting both predictive accuracy and uncertainty quantification. Consequently, there is a need for an approach that treats all fields symmetrically while strictly respecting the underlying physics. Motivated by these limitations, we propose a robust framework for jointly modeling constrained multi-field data. Our approach first leverages a specific PCA procedure for multi-field data, coined row-wise PCA, which has the interesting property of preserving the constraint in the latent space. Since standard PCA strategies for multi-field data do not preserve such constraints, we investigate theoretically the optimality of the row-wise choice. In a second step, we consider a linearly-constrained multi-output GP approach based on a specific kernel parametrization which is trained on the latent space of row-wise PCA. The proposed framework is validated on a population dynamics problem and on an industrial CFD application, which involves the prediction of Reynolds stress tensor components under the incompressibility constraint.

PaperPDFCode Syntology ran

In Syntology 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="2608.25709")

Code

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

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

nasrmht/paper_benchmark_repo found in paper text by Syntology report

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

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

10ran
1unverified

Licence: 11 of the 11 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 nasrmht/paper_benchmark_repo. “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_benchmark_config nasrmht/paper_benchmark_repo/benchmark_pca_gp/run_lotka_volterra.py found in paper text by Syntology ran no licence file found · pointer only · 279134da0d9d00e6 · report
build_tidy_df nasrmht/paper_benchmark_repo/constraint_benchmark/plot_ranking_percentages.py found in paper text by Syntology ran no licence file found · pointer only · 08caaf32ffd404e8 · report
denormalise_cross_covs nasrmht/paper_benchmark_repo/benchmark_pca_gp/models/base.py found in paper text by Syntology ran no licence file found · pointer only · 16e6b9ebacb9a5df · report
denormalise_means nasrmht/paper_benchmark_repo/benchmark_pca_gp/models/base.py found in paper text by Syntology ran no licence file found · pointer only · f9e3071d12ecc059 · report
generate_complex_data nasrmht/paper_benchmark_repo/constraint_benchmark/run_benchmark.py found in paper text by Syntology ran no licence file found · pointer only · 1cef51fea5f350b5 · report
load_and_aggregate_results nasrmht/paper_benchmark_repo/constraint_benchmark/plot_violinplot.py found in paper text by Syntology ran fingerprinted no licence file found · pointer only · 3977fc2078b3d93f · report
load_results nasrmht/paper_benchmark_repo/constraint_benchmark/plot_ranking_percentages.py found in paper text by Syntology ran fingerprinted no licence file found · pointer only · d7922f99dc136c87 · report
normalise_weights_per_mode nasrmht/paper_benchmark_repo/benchmark_pca_gp/models/base.py found in paper text by Syntology ran no licence file found · pointer only · 779de6fac9061a36 · report
scaled_branin nasrmht/paper_benchmark_repo/constraint_benchmark/run_benchmark.py found in paper text by Syntology ran fingerprinted no licence file found · pointer only · c12054bdf4376ade · report
scaled_ishigami nasrmht/paper_benchmark_repo/constraint_benchmark/run_benchmark.py found in paper text by Syntology ran fingerprinted no licence file found · pointer only · b175978eb8e073bd · report
compute_rank_percentages nasrmht/paper_benchmark_repo/constraint_benchmark/plot_ranking_percentages.py found in paper text by Syntology unverified no licence file found · pointer only · d93cbecbb17ff9aa · report

Results from the paper

The Papers with Code archive ends with its 2025-07-28 snapshot. This paper's arXiv identifier, 2608.25709, was issued in August 2026, after that date, so the archive has no leaderboard rows for it.

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