{"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/high-dimensional-gaussian-sampling-a-review","title":"High-dimensional Gaussian sampling: a review and a unifying approach based on a stochastic proximal point algorithm","arxiv_id":"2010.01510","date":"2020-10-04","proceeding":null,"authors":["Maxime Vono","Nicolas Dobigeon","Pierre Chainais"],"abstract":"Efficient sampling from a high-dimensional Gaussian distribution is an old but high-stake issue. Vanilla Cholesky samplers imply a computational cost and memory requirements which can rapidly become prohibitive in high dimension. To tackle these issues, multiple methods have been proposed from different communities ranging from iterative numerical linear algebra to Markov chain Monte Carlo (MCMC) approaches. Surprisingly, no complete review and comparison of these methods have been conducted. This paper aims at reviewing all these approaches by pointing out their differences, close relations, benefits and limitations. In addition to this state of the art, this paper proposes a unifying Gaussian simulation framework by deriving a stochastic counterpart of the celebrated proximal point algorithm in optimization. This framework offers a novel and unifying revisit of most of the existing MCMC approaches while extending them. Guidelines to choose the appropriate Gaussian simulation method for a given sampling problem in high dimension are proposed and illustrated with numerical examples.","url_abs":"https://arxiv.org/abs/2010.01510v2","url_pdf":"https://arxiv.org/pdf/2010.01510v2.pdf","source":{"archive":"pwc-archive (Hugging Face), CC BY-SA 4.0","snapshot":"2025-07-28","licence_url":"https://creativecommons.org/licenses/by-sa/4.0/legalcode","row_kind":"links_only","authors_date_abstract":"arXiv metadata, CC0 1.0 (https://info.arxiv.org/help/license), from the Kaggle arXiv metadata snapshot of 2026-09-12"},"code_links":[{"paper_slug":"high-dimensional-gaussian-sampling-a-review","repo_url":"https://github.com/mvono/PyGauss","is_official":1,"mentioned_in_paper":1,"mentioned_in_github":1,"framework":"none","reach":{"status":"ok","spdx":"MIT"}}],"tasks":[],"methods":[],"datasets_introduced":[],"methods_introduced":[],"results":[],"syntology":{"syntology_url":"https://syntology.ai/paper/2010.01510","atlas_url":"https://app.syntology.ai/?focus=2010.01510","mcp":{"get_harvested_code_for_paper":{"arxiv_id":"2010.01510"}},"developers":"https://syntology.ai/developers","read_at":"2026-09-25T09:33:49+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":[{"provenance":"external:paperswithcode_snapshot_2025-07-28","url":"https://github.com/mvono/PyGauss","reach":{"status":"ok","spdx":"MIT"}}],"summary":{"unverified":6},"by_repo_kind":{"official":{"samples":6,"ran":0,"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":0,"samples":[{"code_sha256_prefix":"568b3f2a3ea03bd8","entry":"CG","repo":"mvono/PyGauss","repo_kind":"official","path":"pygauss/utils.py","file_url":"https://github.com/mvono/PyGauss/blob/HEAD/pygauss/utils.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"MIT","inline_ok":true,"mcp_get_code":{"code_sha256":"568b3f2a3ea03bd8"}},{"code_sha256_prefix":"368718f126c7a20c","entry":"col_vector_norms","repo":"mvono/PyGauss","repo_kind":"official","path":"pygauss/utils.py","file_url":"https://github.com/mvono/PyGauss/blob/HEAD/pygauss/utils.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"MIT","inline_ok":true,"mcp_get_code":{"code_sha256":"368718f126c7a20c"}},{"code_sha256_prefix":"381a944879a3b4bd","entry":"diagonal_form","repo":"mvono/PyGauss","repo_kind":"official","path":"pygauss/utils.py","file_url":"https://github.com/mvono/PyGauss/blob/HEAD/pygauss/utils.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"MIT","inline_ok":true,"mcp_get_code":{"code_sha256":"381a944879a3b4bd"}},{"code_sha256_prefix":"2f9f15d43cbbaf03","entry":"sampler_band","repo":"mvono/PyGauss","repo_kind":"official","path":"pygauss/direct_sampling.py","file_url":"https://github.com/mvono/PyGauss/blob/HEAD/pygauss/direct_sampling.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"MIT","inline_ok":true,"mcp_get_code":{"code_sha256":"2f9f15d43cbbaf03"}},{"code_sha256_prefix":"4e0fcc70b3b86bf8","entry":"sampler_circulant","repo":"mvono/PyGauss","repo_kind":"official","path":"pygauss/direct_sampling.py","file_url":"https://github.com/mvono/PyGauss/blob/HEAD/pygauss/direct_sampling.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"MIT","inline_ok":true,"mcp_get_code":{"code_sha256":"4e0fcc70b3b86bf8"}},{"code_sha256_prefix":"bc624a1a1ecb2db7","entry":"sampler_factorization","repo":"mvono/PyGauss","repo_kind":"official","path":"pygauss/direct_sampling.py","file_url":"https://github.com/mvono/PyGauss/blob/HEAD/pygauss/direct_sampling.py","link_basis":"first_harvest_node","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"MIT","inline_ok":true,"mcp_get_code":{"code_sha256":"bc624a1a1ecb2db7"}}]},"arxiv_metadata":null,"syntology_extracted_results":null}