{"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/pdlp-a-practical-first-order-method-for-large","title":"PDLP: A Practical First-Order Method for Large-Scale Linear Programming","arxiv_id":"2501.07018","date":"2025-01-13","proceeding":null,"authors":["David Applegate","Mateo Díaz","Oliver Hinder","Haihao Lu","Miles Lubin","Brendan O'Donoghue","Warren Schudy"],"abstract":"We present PDLP, a practical first-order method for linear programming (LP) designed to solve large-scale LP problems. PDLP is based on the primal-dual hybrid gradient (PDHG) method applied to the minimax formulation of LP. PDLP incorporates several enhancements to PDHG, including diagonal preconditioning, presolving, adaptive step sizes, adaptive restarting, and feasibility polishing. Our algorithm is implemented in C++, available in Google's open-source OR-Tools library, and supports multithreading. To evaluate our method, we introduce a new collection of eleven large-scale LP problems with sizes ranging from 125 million to 6.3 billion nonzeros. PDLP solves eight of these instances to optimality gaps of 1% (with primal and dual feasibility errors of less than $10^{-8}$) within six days on a single machine. We also compare PDLP with Gurobi barrier, primal simplex, and dual simplex implementations. These traditional methods are designed to solve linear programs to much tighter optimality gaps but struggle to solve these instances. Gurobi barrier solves only three instances, exceeding our 1TB RAM limit on the other eight. While primal and dual simplex are more memory-efficient than the barrier method, they are slower and solve only three instances within six days. Compared with the conference version of this work (in: Advances in Neural Information Processing Systems 34 (NeurIPS 2021)), the key new contributions are: (i) feasibility polishing, a technique that quickly finds solutions that are approximately optimal but almost exactly feasible (without which only two of the eleven problems can be solved); (ii) a multithreaded C++ implementation available in Google OR-Tools; and (iii) a new collection of large-scale LP problems. Note that the conference version should be referred to for comparisons with SCS and ablation studies, which we do not repeat in this paper.","url_abs":"https://arxiv.org/abs/2501.07018v1","url_pdf":"https://arxiv.org/pdf/2501.07018v1.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":"pdlp-a-practical-first-order-method-for-large","repo_url":"https://github.com/google-research/FirstOrderLp.jl","is_official":1,"mentioned_in_paper":1,"mentioned_in_github":0,"framework":"none","reach":{"status":"ok","spdx":"Apache-2.0"}}],"tasks":[],"methods":[],"datasets_introduced":[],"methods_introduced":[],"results":[],"syntology":{"atlas_url":"https://app.syntology.ai/?focus=2501.07018","mcp":{"get_harvested_code_for_paper":{"arxiv_id":"2501.07018"}},"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":[{"provenance":"external:paperswithcode_snapshot_2025-07-28","url":"https://github.com/google-research/FirstOrderLp.jl","reach":{"status":"ok","spdx":"Apache-2.0"}},{"provenance":"deterministic:regex_extraction","url":"https://github.com/ohinder/large-scale-LP-test-problems","reach":{"status":"ok","spdx":"MIT"}},{"provenance":"deterministic:regex_extraction","url":"https://github.com/google/or-tools","reach":{"status":"ok","spdx":"Apache-2.0"}}],"summary":{"ran_draft_wrong":1,"ran_fixture":1,"unverified":7},"by_repo_kind":{"found_in_text":{"samples":6,"ran":0,"repositories":1},"community":{"samples":3,"ran":2,"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":3,"samples":[{"code_sha256_prefix":"a6e67078a082f0df","entry":"decompress_mps","repo":"mmaaz-git/pdlp","repo_kind":"community","path":"benchmarks/mittelmann.py","file_url":"https://github.com/mmaaz-git/pdlp/blob/HEAD/benchmarks/mittelmann.py","link_basis":"first_harvest_node","language":"python","status":"ran_draft_wrong","verification_level":1,"contract_check":"OUTPUT_MISDECLARED","metamorphic_tier":"deterministic","behaviour_fingerprint":true,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"a6e67078a082f0df"}},{"code_sha256_prefix":"d4ec6c6df72cb4b0","entry":"solve","repo":"mmaaz-git/pdlp","repo_kind":"community","path":"pdlp.py","file_url":"https://github.com/mmaaz-git/pdlp/blob/HEAD/pdlp.py","link_basis":"first_harvest_node","language":"python","status":"ran_fixture","verification_level":1,"contract_check":"RAISES","metamorphic_tier":"none","behaviour_fingerprint":true,"licence":"NONE","inline_ok":false,"mcp_get_code":{"code_sha256":"d4ec6c6df72cb4b0"}},{"code_sha256_prefix":"fdadc7e6c0d5856d","entry":"deprecated","repo":"google/or-tools","repo_kind":"found_in_text","path":"ortools/sat/python/cp_model.py","file_url":"https://github.com/google/or-tools/blob/HEAD/ortools/sat/python/cp_model.py","link_basis":"harvester_set","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"Apache-2.0","inline_ok":true,"mcp_get_code":{"code_sha256":"fdadc7e6c0d5856d"}},{"code_sha256_prefix":"b064029eab8b1791","entry":"deprecated_method","repo":"google/or-tools","repo_kind":"found_in_text","path":"ortools/sat/python/cp_model.py","file_url":"https://github.com/google/or-tools/blob/HEAD/ortools/sat/python/cp_model.py","link_basis":"harvester_set","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"Apache-2.0","inline_ok":true,"mcp_get_code":{"code_sha256":"b064029eab8b1791"}},{"code_sha256_prefix":"03edfcca135ae5df","entry":"download_mps","repo":"mmaaz-git/pdlp","repo_kind":"community","path":"benchmarks/mittelmann.py","file_url":"https://github.com/mmaaz-git/pdlp/blob/HEAD/benchmarks/mittelmann.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":"03edfcca135ae5df"}},{"code_sha256_prefix":"a178873f45801076","entry":"is_a_number","repo":"google/or-tools","repo_kind":"found_in_text","path":"ortools/linear_solver/python/model_builder_numbers.py","file_url":"https://github.com/google/or-tools/blob/HEAD/ortools/linear_solver/python/model_builder_numbers.py","link_basis":"harvester_set","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"Apache-2.0","inline_ok":true,"mcp_get_code":{"code_sha256":"a178873f45801076"}},{"code_sha256_prefix":"91a1c170660ea3f6","entry":"is_integral","repo":"google/or-tools","repo_kind":"found_in_text","path":"ortools/linear_solver/python/model_builder_numbers.py","file_url":"https://github.com/google/or-tools/blob/HEAD/ortools/linear_solver/python/model_builder_numbers.py","link_basis":"harvester_set","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"Apache-2.0","inline_ok":true,"mcp_get_code":{"code_sha256":"91a1c170660ea3f6"}},{"code_sha256_prefix":"ae9dc9e14f1c761b","entry":"is_zero","repo":"google/or-tools","repo_kind":"found_in_text","path":"ortools/linear_solver/python/model_builder_numbers.py","file_url":"https://github.com/google/or-tools/blob/HEAD/ortools/linear_solver/python/model_builder_numbers.py","link_basis":"harvester_set","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"Apache-2.0","inline_ok":true,"mcp_get_code":{"code_sha256":"ae9dc9e14f1c761b"}},{"code_sha256_prefix":"3d245e14dadd311d","entry":"snake_case_to_camel_case","repo":"google/or-tools","repo_kind":"found_in_text","path":"ortools/sat/python/cp_model.py","file_url":"https://github.com/google/or-tools/blob/HEAD/ortools/sat/python/cp_model.py","link_basis":"harvester_set","language":"python","status":"unverified","verification_level":0,"contract_check":null,"metamorphic_tier":null,"behaviour_fingerprint":false,"licence":"Apache-2.0","inline_ok":true,"mcp_get_code":{"code_sha256":"3d245e14dadd311d"}}]},"arxiv_metadata":null,"syntology_extracted_results":null}