<p>The primal-dual hybrid gradient (PDHG) algorithm is a prominent first-order primal-dual method designed to efficiently solve saddle point problems and associated convex optimization problems. Recently, the underlying mechanism of the PDHG algorithm has been elucidated through the high-resolution ordinary differential equation (ODE) and the implicit-Euler scheme as detailed in [<CitationRef CitationID="CR20">20</CitationRef>]. This insight has spurred the development of several accelerated variants of the PDHG algorithm, originally proposed by [<CitationRef CitationID="CR4">4</CitationRef>]. By employing discrete Lyapunov analysis, we establish that the PDHG algorithm with iteration-varying step sizes, converges at a rate near <InlineEquation ID="IEq1"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10957_2025_2835_Article_IEq1.gif" Format="GIF" Height="20" Rendition="HTML" Resolution="72" Type="Linedraw" Width="59" /> </InlineMediaObject> <EquationSource Format="TEX">\(O(1/k^2)\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>O</mi> <mo stretchy="false">(</mo> <mn>1</mn> <mo stretchy="false">/</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation>. Furthermore, for the specific setting where <InlineEquation ID="IEq2"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10957_2025_2835_Article_IEq2.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="83" /> </InlineMediaObject> <EquationSource Format="TEX">\(\tau _{k+1}\sigma _k = s^2\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi>τ</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <msub> <mi>σ</mi> <mi>k</mi> </msub> <mo>=</mo> <msup> <mi>s</mi> <mn>2</mn> </msup> </mrow> </math></EquationSource> </InlineEquation> and <InlineEquation ID="IEq3"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10957_2025_2835_Article_IEq3.gif" Format="GIF" Height="19" Rendition="HTML" Resolution="72" Type="Linedraw" Width="147" /> </InlineMediaObject> <EquationSource Format="TEX">\(\theta _k = \tau _{k+1}/\tau _k \in (0, 1)\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <msub> <mi>θ</mi> <mi>k</mi> </msub> <mo>=</mo> <msub> <mi>τ</mi> <mrow> <mi>k</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo stretchy="false">/</mo> <msub> <mi>τ</mi> <mi>k</mi> </msub> <mo>∈</mo> <mrow> <mo stretchy="false">(</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo stretchy="false">)</mo> </mrow> </mrow> </math></EquationSource> </InlineEquation> as proposed in [<CitationRef CitationID="CR4">4</CitationRef>], an even faster convergence rate of <InlineEquation ID="IEq4"> <InlineMediaObject> <ImageObject Color="BlackWhite" FileRef="10957_2025_2835_Article_IEq1.gif" Format="GIF" Height="20" Rendition="HTML" Resolution="72" Type="Linedraw" Width="59" /> </InlineMediaObject> <EquationSource Format="TEX">\(O(1/k^2)\)</EquationSource> <EquationSource Format="MATHML"><math> <mrow> <mi>O</mi> <mo stretchy="false">(</mo> <mn>1</mn> <mo stretchy="false">/</mo> <msup> <mi>k</mi> <mn>2</mn> </msup> <mo stretchy="false">)</mo> </mrow> </math></EquationSource> </InlineEquation> can be achieved. To substantiate these findings, we design a novel discrete Lyapunov function. This function is distinguished by its succinctness and straightforwardness, providing a clear and elegant proof of the enhanced convergence properties of the PDHG algorithm under the specified conditions. Finally, we utilize the discrete Lyapunov function to establish the optimal linear convergence rate when both the objective functions are strongly convex. The theoretical convergence results are validated by numerical experiments.</p>

错误:搜索内容不能为空,请输入英文关键词
错误:关键词超出字数限制,请精简
高级检索

A Lyapunov Analysis of Accelerated PDHG Algorithms

  • Xueying Zeng,
  • Bin Shi

摘要

The primal-dual hybrid gradient (PDHG) algorithm is a prominent first-order primal-dual method designed to efficiently solve saddle point problems and associated convex optimization problems. Recently, the underlying mechanism of the PDHG algorithm has been elucidated through the high-resolution ordinary differential equation (ODE) and the implicit-Euler scheme as detailed in [20]. This insight has spurred the development of several accelerated variants of the PDHG algorithm, originally proposed by [4]. By employing discrete Lyapunov analysis, we establish that the PDHG algorithm with iteration-varying step sizes, converges at a rate near \(O(1/k^2)\) O ( 1 / k 2 ) . Furthermore, for the specific setting where \(\tau _{k+1}\sigma _k = s^2\) τ k + 1 σ k = s 2 and \(\theta _k = \tau _{k+1}/\tau _k \in (0, 1)\) θ k = τ k + 1 / τ k ( 0 , 1 ) as proposed in [4], an even faster convergence rate of \(O(1/k^2)\) O ( 1 / k 2 ) can be achieved. To substantiate these findings, we design a novel discrete Lyapunov function. This function is distinguished by its succinctness and straightforwardness, providing a clear and elegant proof of the enhanced convergence properties of the PDHG algorithm under the specified conditions. Finally, we utilize the discrete Lyapunov function to establish the optimal linear convergence rate when both the objective functions are strongly convex. The theoretical convergence results are validated by numerical experiments.