2022/01/07 by Joy Azzam, Daniel J. Henderson, Daniel Henderson +7
Computer Science · Engineering · Mathematics · #65F35 #68W2 #68W20 #90C53 #Advanced Image Processing Techniques #Advanced Numerical Analysis Techniques #Advanced Optimization Algorithms Research #FOS: Mathematics #Optimization and Control (math.OC) #math.OC #msc:65F35 #msc:68W2 #msc:68W20 #msc:90C53
paper · pdf · doi:10.48550/arxiv.2201.02608
arxiv created 2022/01/07 · openalex publication_date 2022/01/07 · arxiv updated 2022/01/10 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
This article explores how to effectively incorporate curvature information generated using SIMD-parallel forward-mode Algorithmic Differentiation (AD) into unconstrained Quasi-Newton (QN) minimization of a smooth objective function, f. Specifically, forward-mode AD can be used to generate block Hessian samples Y=∇2 f(x) S whenever the gradient is evaluated. Block QN algorithms then update approximate inverse Hessians, Hk ≈ ∇2 f(xk), with these Hessian samples. Whereas standard line-search based BFGS algorithms carefully filter and correct secant-based approximate curvature information to maintain positive definite approximations, our algorithms directly incorporate Hessian samples to update indefinite inverse Hessian approximations without filtering. The sampled directions supplement the standard QN two-dimensional trust-region sub-problem to generate a moderate dimensional subproblem which can exploit negative curvature. The resulting quadratically-constrained quadratic program is solved accurately with a generalized eigenvalue algorithm and the step advanced using standard trust region step acceptance and radius adjustments. The article aims to avoid serial bottlenecks, exploit accurate positive and negative curvature information, and conduct a preliminary evaluation of selection strategies for S.