Fast projection onto the intersection of simplex and singly linear constraint and its generalized Jacobian

Solving the distributional worst-case in the distributionally robust optimization problem is equivalent to finding the projection onto the intersection of simplex and singly linear inequality constraint, which is an important ingredient in the design of some first-order efficient algorithms. This paper focuses on efficient algorithms for computing the projection onto the intersection of simplex and singly linear inequality constraint. We apply an algorithm based on Lagrangian relaxation approach and secant method (LRSA) to find this projection. From the perspective of Lagrangian duality, we also develop an common algorithm based on semismooth Newton method (Ssn) for solving the dual of the projection problem. Numerical experiments demonstrate that LRSA outperforms Ssn algorithm and the state-of-the-art optimization solver called Gurobi. Moreover, we derive the explicit formulas for the generalized HS-Jacobian of the projection, which plays an important role in designing second-order nonsmooth Newton algorithms.