Abstract:This paper gives an approach to implementing a second-order primal-dual interior point method. It uses a Taylor polynomial of second order to approximate a primal-dual trajectory. The computations for the second derivative are combined with the computations for the centering direction. Computations in this approach do not require that primal and dual solutions be feasible. Expressions are given to compute all the higher-order derivatives of the trajectory of interest. The implementation ensures that a suitable potential function is reduced by a constant amount at each iteration. There are several salient features of this approach. An adaptive heuristic for estimating the centering parameter is given. The approach used to compute the step length is also adaptive. A new practical approach to compute the starting point is given. This approach treats primal and dual problems symmetrically. Computational results on a subset of problems available from netlib are given. On mutually tested problems the results show that the proposed method requires approximately 40 percent fewer iterations than the implementation proposed in Lustig, Marsten, and Shanno [Tech. Rep. TR J-89-11, Georgia Inst. of Technology, Atlanta, 1989]. It requires approximately 50 percent fewer iterations than the dual affine scaling method in Adler, Karmarkar, Resende, and Veiga [Math. Programming, 44 (1989), pp. 297336], and 35 percent fewer iterations than the second-order dual affine scaling method in the same paper. The new approach for estimating the centering parameter and finding the step length and the starting point have contributed to the reduction in the number of iterations. However, the contribution due to the use of second derivative is most significant. On the tested problems, on the average the implementation shown was found to be approximately two times faster than OBl (version 02/90) described in Lustig, Marsten, and Shanno and 2.5 times faster than MINOS 5.3 described in Murtagh and Saunders [Tech. Rep. SOL 83-20, Dept. of Operations Research, Stanford Univ., Stanford, CA, 1983].

A nonsmooth primal-dual method with interwoven PDE constraint solver

A nonsmooth primal-dual method with interwoven PDE constraint solver

An Efficient Interior-Point Algorithm with New Non-Monotone Line Search Filter Method for Nonlinear Constrained Programming

A Smooth Primal-Dual Optimization Framework for Nonsmooth Composite Convex Minimization

Inertial Accelerated Primal-Dual Methods for Linear Equality Constrained Convex Optimization Problems.

Solving Non-smooth Constrained Programs with Lower Complexity than O(1/epsilon): A Primal-Dual Homotopy Smoothing Approach

A Fast Single-Loop Primal-Dual Algorithm for Non-Convex Functional Constrained Optimization

Accelerated primal-dual methods with enlarged step sizes and operator learning for nonsmooth optimal control problems

On the Implementation of a Primal-Dual Interior Point Method

On the resolution of the non-smooth inverse Cauchy problem by the primal-dual method

On the Differentiability of the Primal-Dual Interior-Point Method

Solving Non-smooth Constrained Programs with Lower Complexity Than $\mathcal{o}(1/\varepsilon)$: A Primal-Dual Homotopy Smoothing Approach

Parametrization and convergence of a primal-dual block-coordinate approach to linearly-constrained nonsmooth optimization

A splitting primal-dual proximity algorithm for solving composite optimization problems

A distributed primal-dual interior-point method for loosely coupled problems using ADMM

A primal-dual interior-point relaxation method with global and rapidly local convergence for nonlinear programs

Fast Primal–dual Algorithm Via Dynamical System for a Linearly Constrained Convex Optimization Problem

A nested primal–dual iterated Tikhonov method for regularized convex optimization

A First-Order Primal-Dual Method for Nonconvex Constrained Optimization Based on the Augmented Lagrangian

A fast primal-dual algorithm via dynamical system with variable mass for linearly constrained convex optimization

A Primal-Dual Smoothing Framework for Max-Structured Non-Convex Optimization