×

Sparse Hessian based semismooth Newton augmented Lagrangian algorithm for general \(\ell_1\) trend filtering. (English) Zbl 1507.90103

Summary: This paper investigates a semismooth Newton based augmented Lagrangian (Ssnal) algorithm for solving equivalent formulation of the general \(\ell_1\) trend filtering problem. The computational costs of a semismooth Newton (Ssn) algorithm for solving the subproblem in the Ssnal algorithm can be substantially reduced by exploiting the second order sparsity of Hessian matrix and some efficient techniques. The global convergence and the asymptotically superlinear local convergence of the Ssnal algorithm are given under mild conditions. Numerical comparisons between the Ssnal algorithm and other state-of-the-art algorithms on real and synthetic data sets validate that our algorithm has superior performance in both robustness and efficiency.

MSC:

90C06 Large-scale problems in mathematical programming
90C25 Convex programming
90C31 Sensitivity, stability, parametric optimization