Skip to content
DraconiticResearch
Search Sources API
197,768 studies31 sourcesUpdated 1 h ago
PreprintarXiv licencearXiv

PDFP: GPU-Accelerated Feasibility Pump via PDHG

Abstract

The feasibility pump (FP) is among the most widely used primal heuristics for mixed-integer linear programming. It alternates between rounding an LP-feasible point and projecting the rounded point back onto the feasible region of the LP relaxation. Recent pump-based heuristics compute the projection with the primal--dual hybrid gradient method (PDHG), which is scalable on modern GPUs. They pose it as a linear program, which needs auxiliary variables and constraints for general integers. Once the projection runs on the GPU, the constraint-propagation rounding (CPR) of FP~2.0 takes a large share of each iteration. For the projection, we propose PDFP, a custom first-order method that solves the pump's projection problem directly, without auxiliary primal variables or constraints. We prove that it keeps the accelerated linear rate of restarted Halpern PDHG. For the rounding, we propose a new procedure that is more efficient than CPR and, under certain assumptions, still returns a CPR rounding. It also chooses rounded values by constraint violation, which can produce feasible roundings where nearest rounding cannot. Finally, our pump runs two independent trajectories concurrently on one GPU through CPU--GPU synchronization. This exploits the variance in trajectory length and lets the pump terminate sooner on average. Computational experiments on MIPLIB~2017 instances show early, promising results: these contributions speed up the discovery of the first incumbent solution. The code is publicly available at https://github.com/MIT-Lu-Lab/PDFP.

Related

  1. cuPDQPx: A GPU-Accelerated Restarted Halpern PDHG for Large-Scale Convex Quadratic Programming
    Haihao Lu, Zedong Peng, Jinwen Yang
    arXivPreprintarXiv licenceFull textOptimisation & controlarXiv
  2. H-PDLP: GPU-Accelerated Linear Programming via Generalized PDHG and Homogeneous Self-Dual Embedding
    Haihao Lu, Bo Tang
    arXivPreprintarXiv licenceFull textOptimisation & controlarXiv
  3. AT-SKM-Net: An Accelerated Trainable Sampling Kaczmarz-Motzkin Framework for Linear Hard-Constraint Feasibility on Dynamic Graphs
    Xiaochen Zhang, Haoyu Zhu, Yao Zhang +1
    arXivPreprintCC BY 4.0Full textOptimisation & controlarXiv
  4. Using a MIP Solver as a PDHG-Based MIP Heuristic
    Edward Rothberg
    arXivPreprintarXiv licenceFull textOptimisation & controlarXiv
  5. Inexact PDHG for Large-Scale Semidefinite Programming
    Karthik Prakhya, Zimeng Wang
    arXivPreprintarXiv licenceFull textOptimisation & controlarXiv
BibTeX
@misc{liu2026pdfp,
  title = {{PDFP: GPU-Accelerated Feasibility Pump via PDHG}},
  author = {Ruofeng Liu and Haihao Lu and Rahul Mazumder},
  howpublished = {arXiv},
  year = {2026},
  month = nov,
  doi = {10.48550/arxiv.2611.03629},
  eprint = {2611.03629},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2611.03629},
}
RIS
TY  - UNPB
TI  - PDFP: GPU-Accelerated Feasibility Pump via PDHG
AU  - Ruofeng Liu
AU  - Haihao Lu
AU  - Rahul Mazumder
PY  - 2026
DA  - 2026/11/02
PB  - arXiv
DO  - 10.48550/arxiv.2611.03629
UR  - https://arxiv.org/abs/2611.03629
AB  - The feasibility pump (FP) is among the most widely used primal heuristics for mixed-integer linear programming. It alternates between rounding an LP-feasible point and projecting the rounded point back onto the feasible region of the LP relaxation. Recent pump-based heuristics compute the projection with the primal--dual hybrid gradient method (PDHG), which is scalable on modern GPUs. They pose it as a linear program, which needs auxiliary variables and constraints for general integers. Once the projection runs on the GPU, the constraint-propagation rounding (CPR) of FP~2.0 takes a large share of each iteration. For the projection, we propose PDFP, a custom first-order method that solves the pump's projection problem directly, without auxiliary primal variables or constraints. We prove that it keeps the accelerated linear rate of restarted Halpern PDHG. For the rounding, we propose a new procedure that is more efficient than CPR and, under certain assumptions, still returns a CPR rounding. It also chooses rounded values by constraint violation, which can produce feasible roundings where nearest rounding cannot. Finally, our pump runs two independent trajectories concurrently on one GPU through CPU--GPU synchronization. This exploits the variance in trajectory length and lets the pump terminate sooner on average. Computational experiments on MIPLIB~2017 instances show early, promising results: these contributions speed up the discovery of the first incumbent solution. The code is publicly available at https://github.com/MIT-Lu-Lab/PDFP.
ER  -