Nonlinear Optimization System Using Primal-Dual Barrier Methods

Resolve Bottlenecks,
Find Innovative Solutions
Generate Solutions

Solution Overview

Problem

Classical log-barrier methods for handling inequality constraints in optimization problems face challenges such as difficulty in finding a strictly feasible starting point and ill-conditioning of the Hessian matrix, leading to numerical difficulties and poor performance near the solution.

Innovation Solution

The proposed method employs a modified log-barrier function with a primal-dual framework, updating barrier parameters and Lagrange multipliers using complementarity slackness tests and a Newton-like method to ensure convergence, and applies quadratic extrapolation to manage infeasibility, thereby improving convergence properties and handling singularities.

Engineering Contradictions & Design Principles

VSEngineering Contradiction Analysis

1Reliability

If a modified log-barrier function is used to eliminate deficiencies of classical log-barrier methods, then convergence properties are improved, but numerical difficulties arise due to singularities at the boundary

Engineering Contradiction:
Improveconvergence propertiesVSAvoidnumerical difficulties from singularities
Core Design Contradiction:
ReliabilityVSObject-affected harmful factors

Solution Approach 1:

The patent introduces a smoothing parameter μ that acts as an intermediary to regularize the modified log-barrier function. This parameter prevents the function from becoming singular at the boundary by adding a regularization term, thereby maintaining numerical stability while preserving the improved convergence properties of the modified log-barrier approach

Inventive Principle:
Principle #24Intermediary (Mediator)

Solution Approach 2:

The patent dynamically adjusts the smoothing parameter μ during the optimization process. By changing this parameter from its initial value to approach zero progressively, the method maintains numerical stability in early iterations and achieves accurate convergence in later iterations, effectively resolving the singularity issue while maintaining reliability

Inventive Principle:
Principle #35Parameter changes

2Productivity

If a Hessian matrix of the log-barrier function is used in classical methods, then optimization can proceed, but the Hessian matrix becomes increasingly ill-conditioned as the barrier parameter converges to zero

Engineering Contradiction:
Improveoptimization progressVSAvoidHessian matrix conditioning
Core Design Contradiction:
ProductivityVSReliability

Solution Approach 1:

The patent transforms the Hessian matrix by applying a parameter change to its diagonal elements. Specifically, it adds a multiple of the identity matrix to the Hessian, where the multiple is determined by the smallest eigenvalue. This transformation maintains the positive definiteness and improves conditioning while preserving the optimization trajectory

Inventive Principle:
Principle #35Parameter changes

Solution Approach 2:

The patent anticipates the ill-conditioning problem by pre-computing the eigenvalue decomposition of the Hessian matrix before it becomes severely ill-conditioned. By doing the decomposition earlier when the matrix is still well-conditioned, the method avoids numerical instability while maintaining optimization progress

Inventive Principle:
Principle #11Beforehand cushioning (Prior cushioning)

Data Source

PatentUS11062219B1Nonlinear optimization system
Publication Date: 2021.07.13 SAS INSTITUTE INC
  • US11062219B1 patent drawing
  • US11062219B1 patent drawing
  • US11062219B1 patent drawing

AI summary

A computer solves a nonlinear optimization problem. An optimality check is performed for a current solution to an objective function that is a nonlinear equation with constraint functions on decision variables. When the performed optimality check indicates that the current solution is not an optimal solution, a barrier parameter value is updated, and a Lagrange multiplier value is updated for each constraint function based on a result of a complementarity slackness test. The current solution to the objective function is updated using a search direction vector determined by solving a primal-dual linear system that includes a dual variable for each constraint function and a step length value determined for each decision variable and for each dual variable. The operations are repeated until the optimality check indicates that the current solution is the optimal solution or a predefined number of iterations has been performed.