← 研踪 PaperTrailarXiv:2505.00311 · v4
A Practical GPU-Enhanced Matrix-Free Primal-Dual Method for Large-Scale Conic Programs
首次提交 2025-05-01 13:13
版本更新 2026-10-29 07:45
雷达观测 2026-10-30 11:25
摘要
In this paper, we introduce a practical GPU-enhanced matrix-free first-order method for solving large-scale conic programming problems, which we refer to as PDCS, standing for the Primal-Dual Conic Programming Solver. Problems that it solves include linear programs, second-order cone programs, convex quadratic programs, and exponential cone programs. The method avoids matrix factorizations and uses sparse matrix-vector multiplication as its core operation, which is memory-efficient and well suited for GPU acceleration. It builds on a variant of the primal-dual hybrid gradient method (PDHG) called reflected restarted Halpern PDHG with averaging, which we prove converges on general conic programs, and further incorporates several practical enhancements. Its GPU implementation, cuPDCS, uses customized computational schemes that assign cones of different types and sizes to different levels of the GPU architecture. Numerical experiments demonstrate that cuPDCS is generally more efficient than state-of-the-art commercial solvers and other first-order methods on large-scale conic program applications, including Fisher market equilibrium problems, Lasso regression, and multi-period portfolio optimization. cuPDCS also exhibits better scalability, efficiency, and robustness compared to other first-order methods on the conic program benchmark dataset CBLIB. These advantages are more pronounced in large-scale settings.
作者备注：40 pages, 7 figures
观测到的版本
相邻记录可能跳过中间版本；首次发现 v2 及以上也会标为版本更新。
v4 · 2026-10-29 07:45 · 最新观测版本
A Practical GPU-Enhanced Matrix-Free Primal-Dual Method for Large-Scale Conic Programs
In this paper, we introduce a practical GPU-enhanced matrix-free first-order method for solving large-scale conic programming problems, which we refer to as PDCS, standing for the Primal-Dual Conic Programming Solver. Problems that it solves include linear programs, second-order cone programs, convex quadratic programs, and exponential cone programs. The method avoids matrix factorizations and uses sparse matrix-vector multiplication as its core operation, which is memory-efficient and well suited for GPU acceleration. It builds on a variant of the primal-dual hybrid gradient method (PDHG) called reflected restarted Halpern PDHG with averaging, which we prove converges on general conic programs, and further incorporates several practical enhancements. Its GPU implementation, cuPDCS, uses customized computational schemes that assign cones of different types and sizes to different levels of the GPU architecture. Numerical experiments demonstrate that cuPDCS is generally more efficient than state-of-the-art commercial solvers and other first-order methods on large-scale conic program applications, including Fisher market equilibrium problems, Lasso regression, and multi-period portfolio optimization. cuPDCS also exhibits better scalability, efficiency, and robustness compared to other first-order methods on the conic program benchmark dataset CBLIB. These advantages are more pronounced in large-scale settings.
查看此版本解读 →查看此版本原文 ↗