{"about":{"site":"https://codewithpapers.app","non_affiliation":"Code with Papers and Syntology are not affiliated with, endorsed by, or sponsored by Papers with Code, Meta, or the pwc-archive mirror.","licence":"CC BY-SA 4.0","licence_url":"https://creativecommons.org/licenses/by-sa/4.0/legalcode","attribution":"https://codewithpapers.app/attribution","modified":"archive material modified by Syntology; see the attribution page"},"url":"/paper/arxiv-2608-25709","title":"Multi-output Gaussian process prediction of physical fields under linear equality constraints","arxiv_id":"2608.25709","date":"2026-08-26","proceeding":null,"authors":["Mahamat Hamdan Nassouradine","Clément Gauchy","Pierre-Emmanuel Angeli","Sébastien da Veiga"],"abstract":"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.","url_abs":"https://arxiv.org/abs/2608.25709","url_pdf":"https://arxiv.org/pdf/2608.25709","source":{"archive":null,"snapshot":"2025-07-28","note":"not in the Papers with Code archive (frozen at the snapshot)","row_kind":"graph","title_abstract_authors_date":"arXiv metadata, CC0 1.0 (https://info.arxiv.org/help/license)"},"code_links":[],"tasks":[],"methods":[],"datasets_introduced":[],"methods_introduced":[],"results":[],"syntology":{"atlas_url":null,"mcp":{"get_harvested_code_for_paper":{"arxiv_id":"2608.25709"}},"developers":"https://syntology.ai/developers","read_at":"2026-09-24T18:15:14+00:00","read_at_is":"when the build read Syntology's graph, not when any sample ran","claim":"Per-sample execution status on synthesized fixtures; not a correctness claim about the paper. Samples come from repositories linked to the paper, official or community; repo_kind says which.","repos":[{"mentioned_in_github":null,"is_official":null,"provenance":"deterministic:regex_extraction","mentioned_in_paper":null,"url":"https://github.com/nasrmht/paper_benchmark_repo","reach":null}],"summary":{"ran":10,"unverified":1},"by_repo_kind":{"found_in_text":{"samples":11,"ran":10,"repositories":1}},"repo_kind_vocabulary":{"official":"The archive marks this repository official for the paper","named_in_paper":"The archive records that the paper mentions this repository; it is not marked official","listed":"In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper","found_in_text":"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"},"n_pointer_only_for_licence":11,"samples":[{"code_sha256_prefix":"279134da0d9d00e6","entry":"build_benchmark_config","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"benchmark_pca_gp/run_lotka_volterra.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/benchmark_pca_gp/run_lotka_volterra.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"279134da0d9d00e6"}},{"code_sha256_prefix":"08caaf32ffd404e8","entry":"build_tidy_df","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/plot_ranking_percentages.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/plot_ranking_percentages.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"08caaf32ffd404e8"}},{"code_sha256_prefix":"16e6b9ebacb9a5df","entry":"denormalise_cross_covs","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"benchmark_pca_gp/models/base.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/benchmark_pca_gp/models/base.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"16e6b9ebacb9a5df"}},{"code_sha256_prefix":"f9e3071d12ecc059","entry":"denormalise_means","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"benchmark_pca_gp/models/base.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/benchmark_pca_gp/models/base.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"f9e3071d12ecc059"}},{"code_sha256_prefix":"1cef51fea5f350b5","entry":"generate_complex_data","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/run_benchmark.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/run_benchmark.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"1cef51fea5f350b5"}},{"code_sha256_prefix":"3977fc2078b3d93f","entry":"load_and_aggregate_results","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/plot_violinplot.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/plot_violinplot.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":true,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"3977fc2078b3d93f"}},{"code_sha256_prefix":"d7922f99dc136c87","entry":"load_results","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/plot_ranking_percentages.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/plot_ranking_percentages.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":true,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"d7922f99dc136c87"}},{"code_sha256_prefix":"779de6fac9061a36","entry":"normalise_weights_per_mode","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"benchmark_pca_gp/models/base.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/benchmark_pca_gp/models/base.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"779de6fac9061a36"}},{"code_sha256_prefix":"c12054bdf4376ade","entry":"scaled_branin","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/run_benchmark.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/run_benchmark.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":true,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"c12054bdf4376ade"}},{"code_sha256_prefix":"b175978eb8e073bd","entry":"scaled_ishigami","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/run_benchmark.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/run_benchmark.py","link_basis":"first_harvest_node","language":"python","status":"ran","verification_level":1,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":true,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"b175978eb8e073bd"}},{"code_sha256_prefix":"d93cbecbb17ff9aa","entry":"compute_rank_percentages","repo":"nasrmht/paper_benchmark_repo","repo_kind":"found_in_text","path":"constraint_benchmark/plot_ranking_percentages.py","file_url":"https://github.com/nasrmht/paper_benchmark_repo/blob/HEAD/constraint_benchmark/plot_ranking_percentages.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"d93cbecbb17ff9aa"}}]},"arxiv_metadata":{"licence":"arXiv metadata, CC0 1.0 (https://info.arxiv.org/help/license)","fields":["title","abstract","authors","date"],"primary_category":"cs.LG","source":"arxiv_2026.jsonl"},"syntology_extracted_results":{"kind":"leaderboard_placements","source":"Syntology's leaderboard-shaped extractor over the paper's own arXiv-HTML tables: a model pointed at a cell, the number was read from that cell and checked against the board's metric, dataset, split and scale, and an independent check accepted the entry; not reviewed by the paper's authors or the archive's editors","extractor_model":"global.anthropic.claude-sonnet-4-5-20250929-v1:0","verifier_model":null,"prompt_sha":"fa63d4bb9d755694","coverage":{"sentence":"Syntology has checked 6,264 of the 9,581 papers on this site that are newer than the archive; results from the others appear after they are checked.","papers_newer_than_archive":9581,"papers_checked":6264},"entries":[],"not_placed":{"boards":0,"rejected_by_independent_check":0,"refused_by_a_rule":0,"check_did_not_answer":0,"proposed_without_a_cell":0,"declined_by_site":0}}}